Questions tagged [error-estimation]
For questions about determining the error caused by specific computational procedures, approximations, or numerical representations.
171
questions
2
votes
0
answers
100
views
Can I reduce my simulation error with a staggered grid, postprocessing and compatibility equation feedback?
What I did
Using the finite difference method, I solved with a certain amount of error the following system of hyperbolic partial differential equations in cylindrical coordinates (the problem is ...
0
votes
0
answers
34
views
Order of Error - Confusion: Clarifying Constraints on Constants and Determining Order of Error
I'm struggling to determine the order of error when considering the error value denoted by $\text{err}$ in relation to the variable $h$. Specifically, I aim to ascertain the value of $x$ in the ...
1
vote
2
answers
123
views
How to estimate the stage error for Runge kutta method
Consider an ordinary differential equation (ODE) in the form $u_t=g(t,u(t))$ and apply the explicit Runge-Kutta method, as defined by the following Butcher tableau:
$$
\mathrm{RK}(s,p):\begin{array}{c|...
0
votes
0
answers
39
views
What is the minimum error achievable using gaussian process emulation?
I am interested in using Gaussian processes as emulators for other computational models, and I would like to characterize the expected numerical precision of the emulator.
Specifically, how small can $...
4
votes
2
answers
438
views
Understanding the diffusion error of numerical schemes
In chapter 7 of Numerical Computation of Internal and External Flows (Second Edition) (https://www.sciencedirect.com/science/article/abs/pii/B9780750665940500497) the author describes the diffusion ...
0
votes
1
answer
60
views
Reverse engineering phase shift and numerical damping
I've been trying to validate the physics behind a particle system framework, but I'm having some difficulties.
A particle system is a set of lumped masses connected by spring-damper elements. Linear ...
4
votes
1
answer
154
views
The error propagation in calculating the inverse using a matrix decomposition
I have been trying to calculate the matrix inverse of some large matrix with entries ranging by orders of magnitude. I tried to use the matrix decomposition to simplify the computation, where a matrix
...
1
vote
1
answer
164
views
Storing Raw Simulation Data or Truncated Data?
I have a simulation that can generate quite a bit of data when it runs, for example $650\cdot 400 \cdot 400$ floating point numbers. Without compression, that's a few gigabytes worth if I want to save ...
0
votes
1
answer
84
views
Finite difference approximation error
I was reading Scientific Computin, An Introductory Survey, by Michael Heath. In the Example 1.11, he madr a Finite Difference Aproximation, with the usual approxination : $f’(x)\neq \frac{f(x+h)-f(x)}{...
1
vote
0
answers
119
views
Which dense matrices are hard to invert?
Suppose I'm solving $Ax=b$ for dense $m\times d$ matrix $A$. For which $A$ is this hard to do?
More concretely, is there any work on estimating the error after $k$ steps of iterative solver, $k\le d$, ...
1
vote
1
answer
79
views
Geometrically nonlinear finite element problem and mesh distortion
In my research on topology optimization of fluid-structure interaction problems (2D), I am using a geometrically nonlinear model to represent the structure. Body/surface forces are extracted from the ...
0
votes
0
answers
42
views
How to compare the bias between the two contour plots?
I wish to compare the value of the normalized concentration (c+), between two contour plots (one is experimental and another is a simulation). If I don't have a c+ value from the same point on both ...
0
votes
0
answers
57
views
Finding optimal values from multiple parameter estimation runs
I've performed a parameter estimation repeat (i.e. 1000 parallel runs with the same
initial values of parameters). I am trying to estimate ~20 parameters using measurements from experiments.
After ...
-1
votes
1
answer
66
views
Convergence of FEM on curved boundaries, and inhomogenous boundary data
In a smooth domain in $\mathbb{R}^{2}$ or $\mathbb{R}^{3}$ let's consider $-\Delta u = f$ with $u=g$ on a part of the boundary and $\partial_\nu u = w$ on another part of the boundary, which is far ...
0
votes
1
answer
90
views
finding discretization error in Burger equation
I was reading the paper given in the link http://www.unige.ch/~hairer/preprints/parareal.pdf and I have a problem in understanding in page 10 for the Burger equation on implementing Parareal method ...
7
votes
1
answer
332
views
How does non-dimensionalization improve the behavior of ODE solvers?
I have a set of coupled ODEs that I'm solving numerically. The independent variable is time and runs from values of $10^{15}$ to $10^{17}$ in units of seconds. The state variables in their usual ...
1
vote
1
answer
86
views
A priori estimates in finite elements for inhomogeneous heat equation
Consider the problem
$$\partial_t u-\Delta u = f\\
u(\Sigma_1)=f_D\\
\partial_\nu u (\Sigma_2)=f_N\\u(0)=u_0$$
where the sides of the space-time cylinder $\Sigma_i$ are disjoint (one of them could be ...
3
votes
0
answers
60
views
How to maximize effectiveness of 2 check bits overseeing 6 data bits?
Having a single parity bit gives a Hamming distance of 2, so only one single bit can be corrupted to reliably detect an error.
The next step I know of would be a Hamming code, which uses at least 3 ...
2
votes
1
answer
137
views
"A posteriori" estimates for finite difference methods
Suppose I have a PDE in a rectangular domain that I am solving numerically via a finite difference method. How do I answer the question, "How fine do I need to make the grid to so that my ...
2
votes
1
answer
190
views
Accuracy loss in single-precision Euclidean norm computation
I do hydrodynamics simulations with Fortran and recently I met with this issue:
I have a single-precision array b of length ...
-1
votes
1
answer
48
views
What is the ERRCON parameter in rkqs?
Ive take a course in computational physics and was asked to implement some numerical methods to solve ODES. I was reading up on the algorithms described in the textbook: NUMERICAL RECIPES IN FORTRAN.
...
0
votes
2
answers
92
views
calculating the Laplacian of the field variable in estimating the local residual error in the finite element method
to perform adaptive refinement in the finite element method according to the explicit residual method, the quantity
$$\eta_K^2=h_K^2\left\lVert r\right\rVert_{L_2(K)}^2+h_K\left\lVert R\right\rVert_{...
0
votes
1
answer
55
views
Numerical Error source when dealing with integer series
I am currently trying to compute the value of the first Fibonacci number recursively. the idea is as follow:
Compute $f_{n}$ and $f_{n-1}$ for $n = 2,...,100$,
Compute $f_k$ for $k = n−2, n−3, \dots, ...
8
votes
0
answers
104
views
How do we approximate the numerical error a numerical scheme (e.g Runge Kutta, Euler etc) makes without having access to an analytical solution?
So I recently encountered this question in my head while taking my Scientific Computing class, where the lecturer talked about computing numerical error of a scheme.
My guess would be that we take a ...
4
votes
1
answer
99
views
Roundoff errors in FEM computations - generalized eigenvalues
This is a continuation of my previous question. I am trying to effectively compute a bound for the roundoff errors in some FEM computation (2d polygons, triangular meshes). Below I will write some of ...
2
votes
2
answers
578
views
Different sources of error in Finite Element computations
Consider the problem $-\Delta u = f$ in $\Omega$, with $u=0$ on $\partial \Omega$. Suppose that $\Omega$ is a polygon and that we approximate the solutions to the previous problem using Lagrange ...
3
votes
1
answer
198
views
How to measure the error of Finite Element approximation in satisfying the PDE?
In Galerkin methods, we seldom can measure the accuracy of an approximation by tracking the value of the residual. For example, take the wave equation:
$ u_{tt} = u_{xx}$, $(x,t) \in (0,L) \times (0,T)...
1
vote
1
answer
308
views
Question about step in the proof of standard discrete trace inequality
I'm studying from Guermond lecture notes available at https://www.math.tamu.edu/~guermond/M661_FALL_2019/chap12.pdf (see Lemma 12.8( Discrete trace inequality).)
Consider the simple case $p=r$, i.e. ...
4
votes
0
answers
85
views
Global reconstruction defined elementwise in a-posteriori error estimator
This question is a follow-up of this previous one. In "Error Control for Discontinuous Galerkin Methods for First Order Hyperbolic Problems" by Georgoulis et al., an error estimator is ...
0
votes
1
answer
109
views
How is the integral of a projection over an element $T$ computed in practice? (deal.II related)
I'm studying an error estimator for the equation $\nabla\cdot(\beta u) + cu = f$ and it contains the following term $$||f - cU_h - \Pi(f-c U_h) ||_T$$
where :
$\Pi$ is the local orthogonal $L^2$ ...
2
votes
1
answer
5k
views
Error in Simpson's 3/8 rule is higher than that of Simpson's 1/3 rule
For a given function $f(x)$, I have tried to find its numerical integral using Simpson's 1/3 and Simpson's 3/8 rules.
I then compare the solution from the numerical quadratures to the analytical ...
3
votes
0
answers
72
views
Typo in a-priori error estimate in a Discontinuous Galerkin paper
I'm looking at this famous paper which is available in the link below:
Franco Brezzi, LD Marini, Endre Süli, Discontinuous Galerkin methods for first-order hyperbolic problems, Mathematical Models ...
1
vote
1
answer
116
views
Finding the source of numerical instability in a electrostatic problem solved by conformal mapping
I'm using conformal mapping to solve a 2D electrostatic problem (calculating the potential $u(x,y)$ in the plane). Let $C_1$ and $C_2$ be two circles at an electric potential $U_1$ and $U_2$, ...
3
votes
1
answer
156
views
Proof of R. Verfürth paper on adaptive mesh and bubble functions
I'm studying adaptive meshes, and my professor wrote the following property for a bubble function ( see this scicomp post for the definition I'm using)$b_T$ defined on a triangle $T$.
$$||b_T \phi ||_{...
1
vote
1
answer
111
views
Classical global estimate for $H^1$ error
I'm having lots of troubles in understanding the proof the estimation of the classical $H^1$ error using finite elements of degree $r$.
$$||u-u_h||_{H^1(\Omega)} \leq \frac{M}{\alpha} C h^r |u|_{H^{r+...
3
votes
1
answer
305
views
How to find the optimum finite difference method for derivatives?
Related to: What are the negatives of using higher order finite diference schemes?
Problem: I have some discrete data of a trajectory $x_t$ with errors $\delta x_t$ of a physical system sampled at ...
4
votes
1
answer
154
views
Computing the residual in a Dual Weighted Residual (DWR) method
I am in the process if computing the Dual-Weighted Residual (DWR) for a linear PDE with a linear functional but I am struggling with the residual part of the calculation.
For example suppose we want ...
1
vote
2
answers
234
views
Stability condition FCTS method
The FTCS method comes from the discretization of a diffusion PDE like this:
$$
a^{2} \frac{u_{i+1}^{k}-2 u_{i}^{k}+u_{i-1}^{k}}{\Delta x^{2}}=\frac{u_{i}^{k+1}-u_{i}^{k}}{\Delta t}
$$
If I have the ...
0
votes
1
answer
53
views
Round-off error step choice
In the Numerical Recipes in section 5.7.- Numerical derivatives the choice of the step size $h$ in the numerical derivative should lead to a difference between $x$ and $x+h$ representable by an exact ...
0
votes
1
answer
218
views
Taylor expansion round-off error
In the Numerical Recipes in section 5.7.- Numerical derivatives it's introduced de roundoff error of:
$$
f^{\prime}(x) \approx \frac{f(x+h)-f(x)}{h}
$$
as (with $h$ an "exact" number):
$$
\...
2
votes
1
answer
65
views
Initial condition precision
Is there a way to have an estimate of the error propagated on an ode numerical solution by the error of the initial conditions? I suppose this depend on the numerical method used and on the problem ...
1
vote
1
answer
299
views
FEM solution for Poisson is not exact at nodes
Let's consider the usual Poisson problem for FEM $$-u''(x)=1 \quad x \in [0,1]$$ with homogeneous Dirichlet boundary conditions. The solution is $u(x)=-\frac{1}{2}x(x-1)$
I know that the FEM solution (...
0
votes
1
answer
212
views
Perturbation problem using Runge-Kutta 4
I'm trying to evaluate the perturbations magnitude between 2 body orbiting a central one in three dimensions. In order to do this I need to have an estimate of the error, which I did using Richardson ...
2
votes
1
answer
229
views
Comparing numerical solutions with very different time grids
I've read an article (Long-term integrations and stability of planetary orbits in our Solar system) in which the authors solved the problem of the absence of an analytical solution for the solar ...
0
votes
1
answer
232
views
Error too large in leapfrog method for solving the wave equation of a vibrating string
I have been trying to figure out what I did wrong for the last two days. I do not know if I actually did something wrong or if the error is supposed to be this large in usual leapfrog problems. I ...
0
votes
1
answer
224
views
What are the advantages and disadvantages of using norm error control in the MATLAB ODE suit?
In MATLAB's ODE suit, there seem to be two basic methods of controlling the Local Truncation Error (LTE) of the ODE which the user can choose from, namely:
The absolute error control (default), ...
4
votes
2
answers
89
views
Dividing functions over a wide range-
I try to solve a system of coupled equations, where a very nasty division operation occurs. In fact, I need to compute a derivative of two exponential decaying functions. Let's illustrate this with ...
1
vote
1
answer
654
views
How to compute the $L^{2}$ error of the gradient in the Finite Element Method
Let $\Omega\subset \mathbb{R}^{2}$ and $\tau_{h} = \{\Omega_{k}\}_{k=1}^{N}$ be a triangulation of $\Omega$. The $L^2$ error for a FEM approximation $u_{h}$ is given by:
$ || u-u_{h} ||_{L^2} = \sqrt{ ...
3
votes
1
answer
107
views
Cauchy Lorentzian simulation on FFT with oscillation
Recently I do simulation on Lorentzian Function with FFT
Lorentzian Function is 2a/(x**2+a**2)
...
0
votes
0
answers
136
views
Error in Monte Carlo integration
I am looking for a concise description to help me understand the error for Monte Carlo Integration using Uniform Sampling and Importance Sampling
For Importance Sampling I have that the error is just ...