Questions tagged [convergence]

Questions related to whether the sequence of iterates generated by an iterative method has one or more limit points, and if those limit points have the correct properties.

1
vote
0answers
9 views

Abaqus, ANSYS, and FVM solver for thermal expansion problem converges to different values

Is it reasonable for a FEM and FVM code to converge to slightly different solutions for the same physical problem (identical BCs, geometry, properties, etc...), provided stability constraints are ...
2
votes
1answer
44 views

Blowup of error in Conjugate Gradient method with periodic Dirichlet Poisson matrix

My problem is that the L2-Norm of the residual for the periodic Poisson matrix $P$ is initially decreasing but starts to blow up after a certain number of iterations. The blowup happens earlier the ...
0
votes
0answers
26 views

Fit exponential convergence

I'm working with a numerical algorithm whose output $y$ asymptotically approaches a certain unknown value $a$. I expect an exponential convergence, i.e. the data $y$ given by my algorithm should be ...
2
votes
0answers
33 views

How to determine the order of convergence of the Euler-Maruyama method?

This question is originally posted in Quant.StackExchange but has been unanswered for some time so I ask in here. To make this simple let us consider the Geometric Brownian Motions (GBM). My ...
1
vote
1answer
47 views

Simplest way to precondition Uzawa iteration

I have a diffusion problem with an internal circular dirichlet constraint and a side condition which shall enforce a certain global volume integral. $\nabla(D \nabla u(x)) = 0$ outer boundary ...
7
votes
1answer
197 views

Lack of quadratic convergence in Newton's method

It is well-known that Newton's method can converge quadratically, if initial guess is close enough and if the arising linear systems are solved accurately. I am applying Newton's method to highly ill-...
0
votes
0answers
42 views

Correct order of convergence

I have a sequence of points which was obtained from an iterative algorithm, and I computed the order of convergence $p$ of the method using the formula $$ p \approx \frac{\log({\rm err}(k+2))-\log({\...
3
votes
3answers
143 views

$H^1$-convergence rate of finite element method for Poisson equation, depending on element order

I wanted to verify my FEM-program by applying the method of manufactured solutions, while solving the Poisson equation in two dimensions using the continuous Galerkin method $$-\nabla^2u=f$$ with $$u=...
4
votes
2answers
95 views

Analytical convergent sequence and numerical divergent sequence

Is it possible to construct a sequence that converges in theory but when computed numerically with a computer program is diverging. I feel that today our computer programs doesn't allow such ...
-1
votes
1answer
59 views

Finite element convergence rate and possion's ratio

I am running simulations of a cantilever beam where it is fixed on one end and negative force applied to the other end. The first simulation is with 4-node linear quadrilateral elements and the other ...
3
votes
2answers
146 views

Asymptotic error of forward Euler

I'm trying to understand the asymptotic error behaviour of forward Euler (finite difference method), as timesteps are decreased (refined), so I feel trust in the method of manufactured solutions (MMS)....
1
vote
1answer
89 views

Analytical testcase for 2D/3D anisotropic Diffusion (Heat Kernel)

I want to verify and compare different Discretizations of the anisotropic diffusion equation in 2D / 3D. In order to both test the timestepping and the spatial discretisations I had a look at using ...
3
votes
0answers
102 views

Unstable convergence of a Poisson equation

What could be the reason that the solution of a Poisson equation is smooth when obtained by an iterative solver, only if the maximum residual is set to a high value (e.g. 0.1)? When the maximum ...
4
votes
0answers
57 views

Preconditioned residual converges, but true residual doesn't

I'm using Albany w/ Trilinos to solve an elasticity problem with thermal expansion mismatch. I'm using block GMRES with MueLu preconditioning. It works for problem size of several million dofs, but ...
4
votes
1answer
118 views

Convergence rate Jacobi/Gauss-Seidel with mesh resolution

In the book A Multigrid Tutorial - Briggs, Henson. McCormick in the beginning of Chapter 3, it is mentioned that ...because the convergence factor behaves as 1-$O(h^{2})$, the coarse grid ...
0
votes
1answer
101 views

Can I convert CUDA core to CPU core and use it as cpu core while running any program?

I was using Metatrader5 and have designed a strategy for trading using MQL5 programming language. While I was running a Strategy Optimization process, I saw the it will need 10,00= iterations or ...
0
votes
1answer
46 views

Conflicting definition of limit point

This question was raised at a different place without sufficient answers. Definition 1: We say that a vector $x \in R^n$ is a limit point of a sequence $\{x_k\}$ in $R^n$ if there exists a ...
1
vote
0answers
51 views

Decrease in slope during convergence analysis

I am using the method of manufactured solutions to perform the order of accuracy testing. I am using a cube for the testing. The cube is size 1m on all sides. I used 5 refinements: $dx = dy = dz = ...
2
votes
1answer
343 views

Numerical Sensitivity in Density of States of Tight-binding model

I'm working with the tight-binding model, and I'm trying to learn the basics of how to compute the Density of States (DOS) $N(E)$ numerically. The DOS is given by $$N(E) = \frac{1}{N}\sum_k \delta(...
3
votes
1answer
135 views

Numerical Lax-Wendroff scheme order of convergence on Burgers equation

I was suggested to move that question here. The question to be as follows. Statement of the problem Is it possible to achieve the second order of convergence (OOC) of Lax-Wendroff (LxW) scheme ...
1
vote
0answers
42 views

Convergence rate assesment in space/time

I am solving a hyperbolic PDE (e.g. the shallow water equations) which depends upon $x$ and $t$. Typically, the overall convergence rate is calculated by comparing the numerical error in different ...
0
votes
1answer
68 views

How to report non-monotonic runtimes in convergence plots

Let's say I have an algorithm that can be tuned by a parameter $h>0$ and is expected to converge as $h\to 0$. I want to study the computational complexity of this algorithm, i.e., how the ...
1
vote
0answers
56 views

BFGS convergence problem

I would like first to state that it is beyond my capability to identify whether this is a BFGS issue or a R package problem. I've been doing some mixed logit regression using the R package "mlogit". ...
1
vote
0answers
72 views

Successive iteration method for solving eigenvalue ploblem

I have a question concerning the branch of successive iteration methods (Newton, Runge-Kutta). I definitely know (or can read in Wikipedia) the implementation of these methods. But I was wondering ...
3
votes
1answer
355 views

Scaling/nondimensionalization for numerical optimization

I have a numerical optimization problem that I am trying to scale appropriately, in order to allow for the solver to achieve faster and more accurate results. I found a paper here that had a short ...
2
votes
2answers
504 views

Convergence problem for Poisson equation with periodic BC

I have written Poisson solvers using two different methods: A classic Jacobi scheme and one using the multigrid solver Hypre. I made up a couple of test cases ensuring the validity of those solvers. ...
1
vote
2answers
101 views

First-order ODE scheme implementation giving less than first-order convergence?

I am solving the initial value problem $$ \frac{d}{dt} (E C_g) = -\delta, \quad E(0) = E_0, $$ for $E$, where $E$ and $C_g$ are functions of $t$, $C_g$ is completely known, and $\delta$ is a function ...
1
vote
0answers
149 views

Accuracy of finite difference method for heat equation on a disk

To study an approximation for the heat equation $$\frac{\partial^2 u}{\partial r^2}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{1}{r^2}\frac{\partial^2 u}{\partial\theta^2}=f(r,\theta)$$ on the ...
4
votes
0answers
313 views

Convergence rate of Picard iterations

Given a first order ODE $y'(x)=f(x,y)$ with the initial condition $y(x_0)=x_0$ such that it satisfies Picard thoerem of existence and uniquness, one can compute the solution by Picard iterations : $$ ...
0
votes
1answer
275 views

Line search bracketing for proximal gradient. Is it good idea?

Maybe my question is obvious but i cannot find any good source which answers it I trying to learn about proximal gradient. One thing which is not clear for me is particular algorithm for line search. ...
1
vote
2answers
93 views

Relationship between global and local error?

In some cases I have seen that if the local error is: $Err = O (\Delta t^{p+1})$ where the global error is p. So if local error is 3, global will be 2. Does somebody know where it comes? For ...
3
votes
2answers
117 views

In numerical methods, eg, finite differencing approaches, does there exist convergent schemes that are not both consistent and stable?

In a book that our course is following this semester, the theorem given is only in one direction: if the scheme is both consistent and stable, then the scheme is convergent. However, since this ...
1
vote
1answer
175 views

Converge rate analysis: issue with time convergence

I have written a code which solves the incompressible formulation of the Navier-Stokes equations. It uses high-order methods both for time and spatial derivatives. I have been conducting convergence ...
13
votes
1answer
242 views

Non-monotonic convergence in fixed-point problem

Background I am solving a variant of the Ornstein-Zernike equation from liquid theory. Abstractly, the problem can be represented as solving the fixed point problem $A c(r)=c(r)$, where $A$ is an ...
9
votes
1answer
337 views

$L^2$-convergence of finite element method when right hand side is only in $H^{-1}$ (Poisson eqn)

I know that the piecewise linear finite element approximation $u_h$ of $$ \Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U $$ satisfies $$ \|u-u_h\|_{H^1_0(U)}\leq Ch\|f\|_{L^2(U)} $...
2
votes
1answer
227 views

Computing rate and order of convergence

this is a follow up question to Convergence rate vs convergence order I guess the whole confusion about how rate vs. order of convergence also came from the implementations I saw. For example, for ...
7
votes
1answer
651 views

Convergence rate vs convergence order

I'm a bit confused about the concepts of convergence rate and convergence order. Let me first give you the definitions we use: [sorry for the English, it's all self translated] Let $x^{*}$ be our ...
3
votes
0answers
95 views

Finite difference method for coupled PDEs: optimizing performance (time step, iterations per step)

I'm solving coupled PDEs using finite difference method: Incompressible Navier-Stokes and the divergence-free induction equation (Maxwell's equations) with non-uniform electrical conductivity. The ...
1
vote
0answers
135 views

Diverged HDG solution for 2D incompressible Navier-Stokes test case at SMALL time step. Why?

I wrote a hybridizable Discontinuous Galerkin code for transient Navier-Stokes flow, with thanks to Martin Kronbichler and Scott Miller's step-51 code in Deal.II. The main algorithm is from Nguyen ...
2
votes
0answers
383 views

Numeric integration over Dirac delta

I'm trying to solve the following integral numerically. $$H(y) = \int dx \, f(x) \, \delta(g(x,y)).$$ For this I chose a representation of the delta-function and employ convergence with respect to $\...
0
votes
1answer
73 views

2nd order accurate finite difference method variable material properties near boundary

I'm aware that a 2nd order accurate finite-difference method using variable properties for central differencing can be written in a finite-volume type way: $$ \nabla \bullet (k \nabla f) = \frac{\...
1
vote
2answers
243 views

Does the convergence of finite element have limit?

I have been trying to solve a nonlinear PDE related to structural mechanics (nonlinear Timoshenko beam to be precise). I am doing both h-refinement and p-refinement to reach the solution. The ...
2
votes
0answers
160 views

Does applying the Newton-Raphson iteration for matrix reciprocal refine a matrix inverse from LU/GE?

This is a follow-up to this answer. Suppose you have a possibly very ill-conditioned matrix $A$, and you compute its inverse with LU/GE to get $X_{\text{lu}}\approx A^{-1}$. The Newton-Raphson ...
1
vote
0answers
50 views

Decreasing - increasing - stabilising $l_{2}$ norm

Let $\bar{x}$ denote the analytical solution of a PDE. Let $x^{(k)}$ be the solution at the $k^{th}$ iteration. The initial guess for the solution is $x^{(0)} = 0$. Let $r_{0} = ||\bar{x}-x^{(0)}||_{2}...
2
votes
2answers
222 views

$O(h^2)$ convergence for Elliptic PDE

I am trying to solve an elliptic PDE in 2-D: $$-\nabla^{2} u = 20tanh(10x-5)(10-10tanh(10x-5)^2) = f$$ I know that the solution is $u = tanh(10x-5)$ but I am unable to get $O(h^2)$ solution with a ...
2
votes
0answers
101 views

Strange convergence behavior of WENO5 for Hamilton--Jacobi equations

I have the following question. I have a code function that computes right-biased and left-biased approximations of the derivative of a function using WENO5 for Hamilton--Jacobi equations as described ...
6
votes
1answer
382 views

Global convergence in trust region algorithm

I was reading about TR methods and there are some terms, which are confusing for me. It says, method is globaly convergent. What does it really mean? Converges to global minima, or converges for ...
3
votes
1answer
236 views

Measure the convergence rate of a discretization of a wave equation

I'm currently trying to approximate the following type of wave equation (in weak formulation): Let $\Omega \subset \mathbb{R}^d$ ($d=2$) be some polygonal domain. We seek a function $u \in L^2\left(0,...
2
votes
1answer
91 views

Interpreting convergence study results, fixed CFL

I am trying to determine the order of my numerical method for resolving a fluid-structure interaction problem using the immersed boundary method. I am using Crank-Nicolson to resolve the fluid ...
2
votes
0answers
60 views

Dynamic Successive Over/under Relaxation (SOR) with several variables

I am solving a partial differential algebraic equation (PDAE) system which has the following dependent variables: $f=f(X,T)$ and $g=g(T)$, along with a few others My current method for coupling is ...