The tag has no wiki summary.

learn more… | top users | synonyms

0
votes
1answer
78 views

Dual time stepping for fluid dynamics

I'm attempting to implement the Weiss and Smith preconditioner in an existing finite volume code and I am struggling with the idea of dual time stepping. My inner time steps are predictor-corrector, ...
1
vote
0answers
31 views

Stationary phase approximation for an integral with infinity saddle points?

I need a hand with the numerical evaluation, in Mathematica, for this integral: $$f(t)=\int_{-\infty}^\infty Exp\{it(\omega_H-\omega_l-\omega_k) - \sum _{j\neq(l,k)} S_j [1-e^{-it\omega_j}]\}\, dt$$ ...
1
vote
1answer
60 views

RATTLE numerical integrator example

I want to understand how the RATTLE algorithm works. Can somebody give me an example (in pseudocode or using any programming language like python or matlab) of how would I implement a numerical ...
6
votes
2answers
144 views

Spectral Methods in time

I was reading up on Spectral Methods for PDEs. In all the descriptions I read, while the position component is approximated via a Fourier series or other methods, the time component is still ...
5
votes
1answer
121 views

Stable time step limits for Velocity-Verlet integration

I'm implementing a mass-spring solid mechanics solver and I'd like to use the Velocity-Verlet time integration scheme. However, I cannot find anything about the maximum stable time step -- either ...
3
votes
0answers
75 views

Stability of the first-order exponential integrator method

The question is about the first-order exponential integration method described in this article. Consider a system of ordinary differential equations $$y'(t) = -A\,y(t) + \mathcal{N}(t, y), \qquad ...
6
votes
1answer
232 views

Why am I getting so much error for my Runge Kutta Fehlberg solver?

My current project is a reprogramming of a protein folding model involving the solution of thousands of ODEs in C++. I've been making some stop and start progress as I'm writing the solver to run ...
4
votes
2answers
153 views

Numerical instability of spherical pendulum

Problem statement I am trying to simulate a spherical pendulum, with rod length $r$ south-polar angle $\theta$ and azimuthal angle $\phi$ initial values $(\theta_0,\phi_0)= (0,0)$ My particular ...
3
votes
0answers
92 views

Time-stepping for coupled nonlinear PDEs

What are good references for time-stepping of the coupled incompressible Navier-Stokes-heat equation (Boussinesq flow), $$ \begin{cases} \rho\left(\dot{\mathbf{u}} + \mathbf{u}\cdot\nabla ...
7
votes
1answer
119 views

What are the differences between Parareal, PITA, and PFASST?

The Parareal, PITA, and PFASST algorithms are all across-the-domain techniques for parallelizing the solution of time-dependent problems in time. What are the guiding principles behind these ...
7
votes
1answer
181 views

Navier-Stokes solver: How to adjust the time step based on non-linear terms?

My code solves the incompressible Navier-Stokes equation in a conducting fluid, together with the induction equation: $ \partial_t u + u \nabla u + 2\Omega \times u = -\nabla p + \nu \Delta u + ...
2
votes
0answers
54 views

Strong stability preserving RK scheme

For the ODE $$ \dot{x} = f(x) $$ we have the 2-stage, second order SSP RK scheme (Shu, Osher, Gottlieb) $$ x^{(0)} = x^n $$ $$ x^{(1)} = x^{(0)} + \Delta t f(x^{(0)}) $$ $$ x^{(2)} = \frac{1}{2} ...
4
votes
2answers
132 views

Improving the time integration of implicit discretized PDE with a non-linear source term

This might be a naive question, but when applying a implicit discretization to a PDE with a source term, should the source be averaged in time? For example if we take the diffusion equation with a ...
5
votes
1answer
81 views

Energy Conservation

I'm working on a time integration scheme for my research. As a result, I have come across an interesting phenomenon. Somehow, the total energy of the scheme oscillates. At any given time the total ...
12
votes
2answers
648 views

What is pseudo time-stepping?

While reading some literature on PDE solvers I came across the term pseudo time-stepping today. It seems to be a common term, however I failed to find a good definition or an introductionary article ...
12
votes
1answer
215 views

What is the correct way of integrating in astronomy simulations?

I'm creating a simple astronomy simulator that should use Newtonian physics to simulate movement of planets in a system (or any objects, for that matter). All the bodies are circles in an Euclidean ...
2
votes
2answers
119 views

What are good examples of problems which are stiff due to very long interval of integration?

There is a class of stiff initial value problems for ODEs that have small Lipschitz constants, slowly-changing solutions, but very long interval of integration. The only practical example of such a ...
6
votes
1answer
107 views

Linearized implicit time stepping

Consider the general FD implicit time stepping scheme $\frac{x_{t+1} - x_t}{\Delta t} = f(x_{t+1})$, where $x$ is the vector variable of interest and $f$ is some function, generally non-linear. ...
7
votes
1answer
137 views

Demonstrating that the time step size is small enough in a code with automatic step size selection

I recently inherited a large body of legacy code that solves a very stiff, transient problem. I would like to demonstrate that the spatial and temporal step sizes are small enough that the ...
10
votes
5answers
277 views

Examples of PDE computations using parallelism in both space and time

In the numerical solution of initial boundary value PDEs, it is very common to employ parallelism in space. It is much less common to employ some form of parallelism in the time discretization, and ...
6
votes
1answer
66 views

How do you change the desired accuracy of a TS object in PETSc?

I'm currently getting very long propagation times when attempting to use the Time Stepping propagators in Petsc 3.2, and in the interest of speeding things up, I'm curious how I can reduce the ...
8
votes
3answers
409 views

Can I use an explicit time stepping scheme to determine numerically whether an ODE is stiff?

I have an ODE: $u'=-1000u+sin(t)$ $u(0)=-\frac{1}{1000001}$ I know that this particular ODE is stiff, analytically. I also know that if we use an explicit (forward) time stepping method ...