Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [stiffness]

Referring to ordinary differential equations that require an extremely small timestep to guarantee numerical stability.

2
votes
2answers
106 views

How does a stiff equation solver work?

I am trying to understand how stiff differential equations are solved. For instance the equation, $$\frac{\partial y}{\partial t} = \alpha\frac{\partial ^2 y}{\partial z^2}$$ can be solved using ...
3
votes
1answer
99 views

Derivation of backward differentiation formulas(BDF)

I have been reading upon numerical techniques that are used to solve stiff ordinary differential equations. From the description given here, I could follow the steps till equation (5). I am finding ...
2
votes
1answer
182 views

Step-size selection for an Trapezoidal Method ODE solver (ode23t)

I was reading the documentation of the MatLab ODE solver ode23t, and I've seen that the trapezoidal rule is used. Moreover, the error is estimated by ...
0
votes
1answer
69 views

Is the time step size of a Rosenbrock method for stiff systems iteratively calculated?

I have an ODE system of the general form y' = k(y)(x) + q(z)(x) x' = a(z)(x) + b(x)(x) where k,q,a and b are also dependent on the states x and y. The ...
4
votes
2answers
154 views

Failing integration with the radau5-implementation in DotNumerics, over a discontinuity

Summary We are trying to solve some models of emission of substances to the air. In these models, emission stops at a certain point. We are using the DotNumerics implementation of radau5 We had ...
1
vote
2answers
368 views

Checking if a given differential equation is stiff

i have to decide if the following differential equation is stiff: $$y''(t)=-201y'-200y^2 + 2, \quad t\in[0,20].$$ Sadly, I don't have any solutions. So, what I did was implementing explicit and ...
1
vote
2answers
778 views

stop condition in scipy.integrate.ode for stiff system

I'm using Python scipy.integrate.ode, and I want to stop my integration at a certain condition. So I use the integrator "dopri5" and use the method "set_solout" to specify a function for my stop ...
2
votes
2answers
2k views

Understanding the Courant–Friedrichs–Lewy condition

I understand these equations in particular can be solved easily without use of computational methods. Although right now I am concerned with trying to solve these equations using numerical integration ...
0
votes
1answer
274 views

Implementing odespy for system of PDEs

After trying to use RK4 to solve the below system of equations, it appears the output had large and fast oscillations even with an adaptive time step I incorporated using the Cash-Karp method. I am ...
0
votes
1answer
199 views

Applying Runge-Kutta to nonlinear system of PDEs

I am applying a 4th order Runge-Kutta code, using the method of lines, to solve the following: \begin{equation} \frac {\partial y_1}{\partial t} = y_2 y_3 - C_1 y_1 \end{equation} \begin{equation} \...
0
votes
1answer
436 views

How to properly implement Backwards Euler on a system of bodies

I have a system of two bodies, $b^i$ and $b^j$, each in position $\vec{p}^i = (x^i, y^i)$ and $\vec{p}^j = (x^j, y^j)$. The two are connected via spring and I'd like to know, given their states at ...
1
vote
0answers
140 views

Numerically solving a system of stiff nonlinear PDEs

I am attempting to numerically solve the following: \begin{align} \frac {\partial y_1}{\partial t} &= i(y_2y_3 - y_2^*y_3^*) - y_1 \tag{1}\\ \frac {\partial y_2}{\partial t} &= y_1^*y_3 - y_2 ...
5
votes
0answers
146 views

Will penalty-augmented stiffness matrix cause numerical issues in eigenvalue analysis?

In the finite element method, we often construct the constraints of the system by adding penalty-function terms ( which often are many many magnitudes, up to $10^6$ order bigger than the largest ...
3
votes
1answer
177 views

Numerical methods for coupled stiff PDEs

I'm dealing with a set of nonlinear coupled PDEs that have the form: \begin{align} \frac {\partial y_1}{\partial t} &= y_2y_3 - y_1 \tag{1}\\ \frac {\partial y_2}{\partial t} &= y_1y_3 - y_2 \...
2
votes
3answers
108 views

Solving ODE with multiple equilibriums

Consider an ODE of the form: $$ u'(t)=-\frac{1}{\varepsilon}u(u-\frac{1}{2})(u-1) $$ with the initial value $$ u(0)=u_0. $$ Here $\varepsilon>0$ is a constant. It is easy to verify that $u\equiv0$ ...
1
vote
0answers
52 views

Stiff ODEs coupled with PDEs (computational efficiency)

I am simulating in COMSOL a system of 3 coupled PDEs (parabolic & elliptic) along with 10 stiff ODEs. In order to have the system working, I am downsizing the time step size too much to achieve ...
2
votes
1answer
127 views

Applicable solvers for nonlinear coupled PDEs

I've been trying to find an applicable PDE solver for cases such as this: Although when dealing with stiff equations in the complex domain, applying existing packages has been problematic. I've ...
1
vote
0answers
106 views

Stiffness emerges as number of ODEs increases

I want to solve a system of ordinary differential equations with Matlab. I need this to solve a mechanical engineering related problem. If $n$ is the number of degrees of freedom of my mechanical ...
0
votes
1answer
248 views

FEM, Direct Stiffness Method with a nonlinear displacement constraint in one node

i have a question about a FE problem im working on. I made a finite element model of an linear elastic block of material (double striped block) attached with a rigid connection to the environment (...
0
votes
1answer
197 views

Convection diffusion reaction equation (stiffness, solver)

I am trying to solve the CDR-Equation in 2D: $$\frac{\partial c(x,y)}{\partial t} + \nabla \cdot ( -d\nabla c(x,y) + \vec{v}(x,y) c(x,y))+ a c(x,y)=0\,,$$ with Boundary Conditions (length of square ...
-3
votes
1answer
199 views

How can I solve stiff equations by method of line (MOL)?

I want to solve 7 coupled equations.I use method of line(MOL) and discrete the equation in Length and radius and convert them to a system of ODEs in time.and use ode15s to solve them in MATLAB. But an ...
2
votes
0answers
49 views

Resolving a stiff hyperbolic problem with Neumann boundary conditions

I am trying to numerically resolve the equation for an Euler-Bernoulli beam that is inextensible, unshearable, and subject to planar deformations: $$\rho I(s) \frac{\partial^2 \theta}{\partial t^2}(s,...
9
votes
2answers
5k views

In FEM, why is the stiffness matrix positive definite?

In FEM classes, it's usually taken for granted that the stiffness matrix is positive definite, but I just can't understand why. Could anyone give some explanation? For instance, we can consider the ...
3
votes
0answers
175 views

Spectral Collocation (or Weighted Residual) Methods to solve Stiff ODEs?

I have a system of ODEs which is (at least moderately) stiff. Consider the class of spectral collocation methods https://en.wikipedia.org/wiki/Spectral_method or the related class of weighted ...
3
votes
1answer
1k views

Does the global stiffness matrix size depend on the number of joints or the number of elements?

When assembling all the stiffness matrices for each element together, is the final matrix size equal to the number of joints or elements?
4
votes
1answer
531 views

What is the case of trade-off in different Runge Kutta methods

There are so many Runge Kutta methods, including Dormand-Prince 45 Cash-Karp 54 Fehlberge 78 Is there any comparison between them? What is each approach sacrificing? What is the general trade-off ...
-1
votes
2answers
100 views

RK4 for stiff IVP

I need a solver of stiff Inital-Value Problems (IVP) in python exploiting RK4 preferably explicit. I have been searching for past few days but could not find it. Following are my queries: Does the ...
1
vote
0answers
39 views

Resources for viscous behavior in simple FEM

I am working on a simple explicit-integration lumped-mass elastic FEM code which implements CST+DKT triangles (plate+shell) and constant-strain tetrahedra (http://woodem.eu/doc/theory/membrane-element....
2
votes
2answers
291 views

Computing element stiffness matrices with variable coefficients

I am trying to implement a simple FEM approach, using p1 triangular elements, for solving the diffusion equation with variable nodal diffusivities and I was wondering how to incorporate the variable ...
1
vote
1answer
84 views

Stiff Equations - What to plot as a qualitative or quantitative measure of stiffness

On a recommendation from Mathematica.SE, I am posting this on Computational Science.SE: I am trying to quantify stiffness of an ODE by relating it to the fine-ness with which NDSolve treats it's ...
3
votes
2answers
357 views

Stiffness ratio evolution for odes

I had a general query related to calculation of stiffness ratio evolution for a set of coupled odes over a certain time interval. My question is, while calculating the eigen values from the jacobian ...
7
votes
2answers
156 views

How do you numerically solve a multivariable ODE system with different time steps per state variable?

If you have a large multivariable ODE system, and certain processes occur at a much smaller time scale, how can you implement a solver that uses smaller time steps for state variables involved in fast ...
6
votes
2answers
400 views

ODE: How to measure stiffness if the Jacobian has zero eigenvalues?

Say you have a system of ODE's where the Jacobian has one zero eigenvalue; what does that tell you about the stiffness of the system? This case doesn't seem to be discussed in the cases I have been ...
2
votes
2answers
143 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 ...
10
votes
3answers
944 views

Explicit Euler method too slow for reaction-diffusion problem

I am solving Turing's reaction-diffusion system with following C++ code. It is too slow: for 128x128 pixel texture, acceptable number of iterations is 200 – which results in 2.5 seconds of delay. I ...
10
votes
3answers
816 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 (...
17
votes
4answers
4k views

The definition of stiff ODE system

Consider an IVP for ODE system $y'=f(x,y)$, $y(x_0)=y_0$. Most commonly this problem is considered stiff when Jacobi matrix $\frac{\partial f}{\partial y}(x_0,y_0)$ has both eigenvalues with very ...