# Tagged Questions

The tag has no usage guidance.

28 views

### Time integration for elastodynamics

I'd like to solve the elastodynamics equation for a problem with fast scales from an impact (although no shocks). $\rho\ddot{\mathbf{u}} - \nabla \cdot \sigma(\mathbf{u}) = \mathbf{f}$ In structural ...
143 views

52 views

### How does Multi Body Dynamics software work for flexible joints?

I need to model a "fishing rod" in 2 dimensions by joining several "rigid sticks" by flexible/elastic joints. The joints act as plate/torsion springs with different spring constants. The "fishing rod" ...
129 views

### Computational time not proportional to integration interval in ODE-solver?

I am running octave and i have been trying ode45, ode54, ode23 etc to integrate the equation  $$Q''(t) = B\cos(Q)\sin(\omega t)$$ $$Q(t=0)=0.$$ When the time interval to be integrated increases, the ...
53 views

### Lie Expansion Integration V/S Runge-Kutta Schemes

It is known that one can perform numerical integration of particular order by considering the Lie Expansion (http://adsabs.harvard.edu/full/1984A%26A...132..203H). For example \frac{...
18 views

### Second and Higher Order Order Corrector in Spectral Deferred Correction

I am trying to work out a second order or higher order correction for the method of Spectral deferred Correction (SDC). Specifically using as a corrector a second order or third order multi-step. In ...
47 views

### What does a negative time stepping mean? (Adaptive time stepping)

Summary behind the problem: The following code aims at solving a static elasto-plastic problem. Like a 2D square mesh based on an elasto-plastic constitutive model like Von-Mises or Drucker-Prager ...
153 views

### Time step size for transient simulations of vortex shedding

I am simulating unsteady flow around a circular cylinder (using FLUENT). I am facing the following problem. Kindly please help me. OBJECTIVE: To find out the Reynolds number at which vortex shedding ...
108 views

### Fourth order IMEX Runge-Kutta method

I am looking for the Butcher tableau of a fourth order accurate Runge-Kutta method with IMEX splitting. I have been reading the ''classical'' paper on the subject by Ascher, Ruuth and Spiteri as well ...
115 views

### Why do planets move at the wrong speed in my solar system model?

Not 100% sure where this question belongs, since I'm not sure if the problem is code related or not. I've written a program to model the solar system. I am now testing its accuracy. I've checked ...
162 views

### Which numerical methods preserve time reversal symmetry?

If I have a physical system which contains a time reversal symmetry (for example a Hamiltonian $H(x,p)=p^2/2m + V(x)$ with $V(x)$ real) and I want to solve the differential equations which describe ...
112 views

### Is it worth switching to timesteppers provided by PETSc if I can't write down a Jacobian for my problem? Case study with “the amoeba” toy problem

I am considering using petsc4py instead of scipy.integrate.odeint (which is a wrapper for Fortran solvers) for a problem ...
77 views

### Why not solve molecular dynamics trajectories exactly?

The basis of Molecular Dynamics simulation is the integration of newton's equations of motion: $$\frac{F}{m} = \frac{d^2r}{dt^2}$$ I understand that there are various methods based on Taylor's ...
100 views

### Time discretization of wave equation

I am trying to model the seismic wave equation and have therefore been reading about discretization schemes and their stability. I recently came across an insightful paper on 'Galerkin FEM methods for ...
95 views

### BDF methods for implicit-explicit method

Are there BDF formulas like the ones given here but one that can be used with implicit-explicit discretization? The right hand side in those formulas is supposed to be implicitly discretized at the ...
100 views

### How to physically understand time dependent boundary conditions?

I am a beginner in Computational science and FEM. I came across some PDEs which implement time dependent boundary conditions. I am not able to visualize exactly a physical scenario of how that would ...
80 views

### PDEs appropriate for adaptive time stepping algorithms

I'm looking for some physical phenomena for which an adaptive time stepping algorithm would be ideal. A PDE or ODE that showed very large gradients in time at a small period of time and smoother ...
40 views

### comparison of stability of two non-linear methods

I have solved a numerical problem using two different sets of non-linear governing equations. I want to get an understanding of the stability of the methods relative to each other. To do so, I solving ...
151 views

146 views

1k views

### Numerical integration and filtering of acceleration experimental data

I have a vector containing acceleration measurements and the corresponding vector of times in which measurements are taken. To obtain velocity and displacement I used the cumtrapz() function already ...
190 views

### How does the L-stability or A-stability of a scheme relate to its ability to preserve a quadratic invariant?

I am working with the simple example of an oscillator: $$(1) \; \; \ddot{u} + u = 0, \; \; u(0) = u_0$$ I know that Forward Euler does not preserve an invariant of the above system: (2) \; \; \dot{...