Questions tagged [stability]

The study of the propagation of errors in a numerical algorithm.

Filter by
Sorted by
Tagged with
2
votes
1answer
79 views

Assessing numerical error in solving a least squares problem

I have a linear system of the type $$Ax = b$$ I want to minimise $|b - Ax|^2$. I know there are different approaches to directly solve the system (Normal equation + Cholesky, QR decomposition, SVD ...
2
votes
1answer
398 views

Von Neumann's stability analysis on non linear and coupled equations

I'd like to know if is possible to make a Von Neumann's stability analysis on an system of coupled equations, featuring quadratics: $$\begin{aligned} \frac{\partial u_1}{\partial t}&=D_1\Delta ...
2
votes
2answers
202 views

von Neumann stability analysis for spatial variable flux

Can we use the von Neumann stability analysis to investigate the stability of the discrete form of the following problem? $$u_t+\frac{\partial(x^2u)}{\partial x}=S(u,x)\ .$$ Please give some hint ...
2
votes
1answer
189 views

What does it mean for a finite difference scheme to be $L^1$-stable?

I am trying to answer a question about a finite difference scheme. I need to show that the method is stable in the L1-norm. I can't find a single definition of what that means, so does anyone have a ...
2
votes
1answer
91 views

Step size and stability of Euler forward method

I'm trying to calculate the maximum step size that provides stability for the following nonlinear IVP using the Euler forward method: $u'(t) = -200tu(t)^2,\qquad u_0 = 1, \qquad t\in [0,3]$, with ...
2
votes
1answer
45 views

Stability of a finite-difference scheme for the reaction-diffusion equation

I currently need to solve numerically the following reaction-diffusion equation: $$\partial_tu=\partial^2_xu+u-u^2$$ For this purpose, I use the following numerical scheme (Crank-Nicolson??): $$ \...
2
votes
2answers
110 views

How to support or contradict a hypothesis on unconditional stability using numerical optimization

The main motivation behind my next question is that I think I derived a higher order numerical scheme for linear advection equation that is unconditionally stable using Von Neumann stability analysis. ...
2
votes
0answers
224 views

Von Neumann analysis for coupled PDE's

I am interested in understanding how to perform stability analysis for coupled (to keep things do-able, lets say linear) PDE's In the case of a single PDE, i understand the logic behind the VN ...
2
votes
0answers
53 views

How to classify chaotic systems from a stability perspective

I am wondering what chaotic systems are from the perspective of numerical analysis. I am talking about 'deterministic chaos' such as for instance the 'logistic map' exhibits it. That is, the solution ...
2
votes
0answers
144 views

what do positive real parts of eigenvalues mean?

I am solving a 1D advection problem of the the form $$d{Q}/dt = [A]{Q}$$ where {Q} is the vector of unknowns and [A] is the matrix of coefficients of spatial discretisation. I have worked out the ...
2
votes
1answer
100 views

Adjusting Keplerian orbits for thrust with numerical stability

I'm writing a mod for a game that models orbital physics (Kerbal Space Program, or KSP). I'm attempting to model the effects of thrust on spacecraft in certain states where the game only models them ...
1
vote
1answer
186 views

Method to solve linear, first order ODE of generalized matrix matrix form

The equation and its meaning: Consider two sets $(A)_{l=0,...,m_a},$,$(B)_{l=0,...,m_b}$ of hermitian matrices and a set of positive semidefinite matrices $(C)_{l=0,...,m_c}$. Each matrix has the ...
1
vote
1answer
837 views

Improve numeric stability of subtraction in C++ [closed]

I'm writing a matrix library in c++. After some debugging I found that a simple double difference is not zero for two "equals" numbers. This is due how double are represented in a computer of course. ...
1
vote
1answer
447 views

More Smearing with decreasing timestep in advection problems

I find it kind of counter intuitive, that the result of an advection gets more smeared out at the borders when decreasing the timestep (which should make it more accurate). Let there be a equally ...
1
vote
2answers
378 views

significance of energy equation and its conservation

I am interested to know the significance of the energy conservation laws when modelling fluids (or other materials). Am I correct in saying that if energy is conserved then stability is achieved. In ...
1
vote
4answers
2k views

Why are upwind schemes stable in convection flow calculation?

It is well known that upwind schemes are stable when calculating convection flows with $|\text{Pe}|>2$, $\text{Pe}$ is the Peclet number. Why is that, and why is central difference unstable? Is ...
1
vote
1answer
131 views

Stability condition for explicit/implicit via non negative coefficients

To make stability proofs simpler, I can consider an explicit scheme written as $$V(n+1,i)=aV(n,i-1)+bV(n,i)+cV(n,i+1)$$ and one can show that if $a,b,c\ge 0$ and $a+b+c\le1$, then the explicit method ...
1
vote
1answer
241 views

Calculate inverse of dense matrix with entries of very different magnitude

I need to calculate the inverse of a dense matrix, with some elements taking values as high as 1e9 and some around 1e2. What would be the best method to do it? Note: I am more concerned about the ...
1
vote
1answer
280 views

energy norm for transport equation

I asked this question before but did not have any luck with an answer. It might be a student level question but I need to understand that with possibly some help. I am considering the hyperbolic ...
1
vote
1answer
487 views

Stability of forward euler method

I am trying to understand the stability of the forward Euler method. I read there's a model problem to see the stability. $$y'(t) = \lambda y(t) \qquad t \in (0, \infty)$$ $$y(0) = 1$$ then the book ...
1
vote
1answer
19 views

Implementation of a ratio with a well-defined limit

I'm implementing a function $f(x_1, x_2, \dots x_n)$ given by the ratio of two closed form equations. $$ f(x_1, x_2, \dots, x_n) = \frac{g(x_1, x_2, \dots, x_n)}{h(x_1, x_2, \dots, x_n)} $$ ...
1
vote
1answer
69 views

How do I add some floating point numbers, keeping numerical accuracy in mind?

I am solving a problem involving the line with the set of points $(x_3,y_3)$ that are equidistant to two given points $(x_1,y_1)$ and $(x_2,y_2)$. The equation for this line is $$(x_3 - x_1)^2 + (y_3 ...
1
vote
1answer
61 views

Locally evaluate nonlinear dynamic system's stability using eigenvalues

I'm working with Computational Neuroscience. I have a large Synaptic Matrix (x axis: presynaptic NeuronID, y axis: postsynaptic NeuronID) in a Modular network. This matrix is close to a random one and ...
1
vote
1answer
618 views

Von Neumann Stability Analysis

I came across the following task recently: Use the von-Neumann stability analysis to investigate the stability of the discrete form of $\frac{\partial c}{\partial x} = \frac{\partial^2 c}{\partial ...
1
vote
0answers
25 views

Finite difference/element method : time step and spatial resolution close to a finite singularity

I'm using the finite element method (FEM), but my question is quite a global question. It's related to this question but it is not the same. Let's assume we have this equation : $$\partial_t c - u\...
1
vote
0answers
58 views

Numerically Approximating the Jacobian and Comparing the Eigenvalues With Analytical Form

I am trying to study the stability of numerical discretization schemes using the Jacobian matrix of the residues with respect to the vector of conserved variables. For a simple diffusion equation ...
1
vote
0answers
45 views

Stability of SVD, Eigendecompositions, and pseudoinverse procedures in modern LAPACK routines

I have proposed an optimisation algorithm which I claim has improved upon the previous algorithm in a number of ways. One of these claims is that my proposed solution requires no explicit SVD and ...
1
vote
0answers
37 views

Stability of different quadrature rules in 1st-kind Volterra integral equation

I am dealing with a integral equation $$ f'(t) = -\int_0^t K(s) f(t-s)\quad t\in [0,t_\max] \tag{1}$$ in which $f(t)$ and $f'(t)$ are known, well-behaved functions of $t$ and $K(t)$ is the unknown. In ...
1
vote
0answers
83 views

Kahan Summation for Three-Term Recurrences

Kahan summation applies to summation problems, but not to three-term recurrence relations. However, a three-term recurrence shares many of the features of a summation-albeit with a rescaling step at ...
1
vote
0answers
140 views

Stability of the explicit MacCormack Scheme to solve the Navier Stokes equations with Wilcox's K-Omega Turbulence Model

I am solving turbulent pipe flow with an explicit MacCormack scheme and Wilcox k-omega model. The laminar version of the code had three distinct stability criteria which worked fine after ...
1
vote
0answers
81 views

Methods to approximate discretized derivatives in PDEs

When solving a general PDE such as $$ \frac {\partial ^2 E}{\partial t ^2} = \frac {\partial ^2 E}{\partial z ^2} - \frac {\partial E}{\partial z} $$ this equation can be solved by the method of ...
1
vote
0answers
40 views

Stability in Discretization of 1D Stationary Boltzmann equation

I want to discretize and numerically solve the following PDE: \begin{equation} v(k)\dfrac{\partial f}{\partial x} + E(x)\dfrac{\partial f}{\partial k} = S\{f\} \end{equation} using finite volume (box ...
1
vote
0answers
338 views

Computing solutions with singularities using MATLAB ODE45

I am new to solving numerically ODES and thus it is difficult for me to judge the reliability/trustworthiness of the results that I have produced for the following problem: I am dealing with a 2nd ...
1
vote
0answers
109 views

Reducing oscillations a 3D Alternating direction explicit scheme for the diffusion equation?

Hi I have made a 3D alternating direction explicit scheme for solving the diffusion equation, which will eventually replace a FTCS scheme in model of bubble dynamics in tissue. I have been testing it ...
1
vote
0answers
167 views

Tips on improving stability in numerical scheme for non-linear PDE

I am solving a non-linear second order system of PDEs in two variables. The equations are too complicated to write out here, but an essential feature is that there is a propagating wave which then ...
1
vote
0answers
41 views

Error analysis and the Model Problem [closed]

In numerical methods for ODE's, the model problem y' = cy where c is complex is regarded as sufficient in performing error analysis for different methods in ...
1
vote
0answers
91 views

Stability question (finite difference): dealing with corner nodes

Consider one initial boundary value problem for sphere. $$\frac{\partial u}{\partial t}=\operatorname{div}A\nabla u +f$$ Here is explicit numerical scheme (we consider that it is stable): $$\frac{u^{...
1
vote
0answers
35 views

Simulation of asymmetric structures (occupancy = 0.5) unstable

I am trying to simulate a metal-organic framework in LAMMPS using the UFF potential. It's working quite well for some structures where all molecules have an occupancy of 1. However, when I have a ...
0
votes
1answer
506 views

Avoid arithmetic overflow in matrix multiplication

I am solving the following matrix equation for $\mathbf{x}$: $$(J^{\mathbf{T}}J)\mathbf{x}=J^{\mathbf{T}}\mathbf{r}$$ $J$ is $m\times n$ matrix $\mathbf{x}$ is vector of size $n$ $\mathbf{r}$ is ...
0
votes
1answer
273 views

Runge-Kutta Stability Regions

Based on this link, in particular Figure 1, what is the exact meaning of the plot? To my understanding, it implies that for a given differential equation: $$ \frac {dy}{dt} = \lambda y $$ that the ...
0
votes
1answer
69 views

How do we find the condition?

Suppose that we are given a numerical scheme. In order to find the CFL condition , we set $U_j^n= \lambda ^ne^{ik x_j}$ and put it into the numerical scheme. I have shown that the given method is ...
0
votes
1answer
495 views

The condition for stability using the leapfrog method

I have the ODE below $$\frac{d}{dt}\pmatrix{x\\ y} = \pmatrix{0 &1\\-a &0}\pmatrix{x\\ y} \enspace .$$ The $m=1$ leapfrog method is defined as: $$y_{n+1} = y_{n-1} + 2f_nh \enspace .$$ For ...
0
votes
1answer
70 views

Numerical Stability of a Generalized Spatial Discretization Scheme

After reading the matrix stability chapter (10) of Hirsch [1], I decided to dive in the reference list of the chapter. One of the papers [2], which is cited as reference shows an very interesting ...
0
votes
1answer
75 views

Stability of PDEs

I am currently trying to solve some PDEs with FiPy. At page 56, the manual mentions (https://www.ctcms.nist.gov/fipy/download/fipy-3.0.pdf). The largest stable timestep that can be taken for this ...
0
votes
1answer
646 views

Finite Difference Method Limitations/Stability Criteria

Is it possible to solve an equation with only a single derivative such as: $$\frac{\partial U(x,t)}{\partial t} = A - BU(x,t)$$ with finite difference methods? I ask as I am trying to solve the ...
0
votes
0answers
17 views

Online Parameter Estimation

I have a system $\ y = θ^*u+d(t) $ where $\ y $ is the output , $\ u $ is the input which satisfies $\ u\epsilon L_{\infty} $, $\ θ^* $ is an unknown but constant parameter and $\ d(t) $ is an unknown ...
0
votes
0answers
34 views

Testing the SUPG method and other methods for hyperbolic equations

I am interesting in integrating the simple equation $$ \frac{\partial \phi}{\partial t} + \mathbf{u}\cdot\nabla \phi = 0 $$ with a Dirichlet boundary condition at the influx boundary ($\mathbf{u} \...
0
votes
0answers
57 views

Bound for Expectation of Singular Value

In my case, $X_{\boldsymbol{\delta}}\in\mathbb{R}^{d\times M}$ is a function of Rademacher variables $\boldsymbol{\delta}\in\{1,-1\}^M$ with $\delta_i$ independent uniform random variables taking ...
0
votes
0answers
59 views

The final Boundary Condition is Unknown, Is Backward Euler is still valid to be implemented?

I am working on conductive polymer modeling and supposed to do one-dimensional diffusion model in the thickness of the polymer, however, due to the small value of thickness in micro, when I use the ...
0
votes
0answers
342 views

Stability of nonlinear partial differential equation

I want to find an expression for the stability of the nonlinear Poisson equation. I know about von Neumann stability analysis which applies to linear equations as far as I know. Any suggestion how to ...