scholarly journals Local Truncation Error-Informed Code Verification

2020 ◽  
Vol 5 (4) ◽  
Author(s):  
Aaron M. Krueger ◽  
Vincent A. Mousseau ◽  
Yassin A. Hassan
Keyword(s):  
Numerical Method ◽  
Truncation Error ◽  
Fluid Dynamic ◽  
Local Truncation Error ◽  
Code Verification ◽  
Verification Method ◽  
First Order ◽  
Manufactured Solutions ◽  
Equation Analysis ◽  
Modified Equation

Abstract The method of manufactured solutions (MMS) has become increasingly popular in conducting code verification studies on predictive codes, such as nuclear power system codes and computational fluid dynamic codes. The reason for the popularity of this approach is that it can be used when an analytical solution is not available. Using MMS, code developers are able to verify that their code is free of coding errors that impact the observed order of accuracy. While MMS is still an excellent tool for code verification, it does not identify coding errors that are of the same order as the numerical method. This paper presents a method that combines MMS with modified equation analysis (MEA), which calculates the local truncation error (LTE) to identify coding error up to and including the order of the numerical method. This method is referred to as modified equation analysis methd of manufactured solutions (MEAMMS). MEAMMS is then applied to a custom-built code, which solves the shallow water equations, to test the performance of the code verification method. MEAMMS is able to detect all coding errors that impact the implementation of the numerical scheme. To show how MEAMMS is different than MMS, they are both applied to the same first-order numerical method test problem with a first-order coding error. When there are first-order coding errors, only MEAMMS is able to identify them. This shows that MEAMMS is able to identify a larger set of coding errors while still being able to identify the coding errors MMS is able to identify.

scholarly journals Rigorous Code Verification: An Additional Tool to Use With the Method of Manufactured Solutions

2019 ◽  
Author(s):  
Aaron M. Krueger ◽  
Vincent A. Mousseau ◽  
Yassin A. Hassan
Keyword(s):  
Exact Solution ◽  
Truncation Error ◽  
Roundoff Error ◽  
Local Truncation Error ◽  
Code Verification ◽  
Order Of Accuracy ◽  
Manufactured Solutions ◽  
Equation Analysis ◽  
Method Of Manufactured Solutions ◽  
Modified Equation

Abstract The Method of Manufactured Solutions (MMS) has proven to be useful for completing code verification studies. MMS allows the code developer to verify that the observed order-of-accuracy matches the theoretical order-of accuracy. Even though the solution to the partial differential equation is not intuitive, it provides an exact solution to a problem that most likely could not be solved analytically. The code developer can then use the exact solution as a debugging tool. While the order-of-accuracy test has been historically treated as the most rigorous of all code verification methods, it fails to indicate code “bugs” that are of the same order as the theoretical order-of-accuracy. The only way to test for these types of code bugs is to verify that the theoretical local truncation error for a particular grid matches the difference between the manufactured solution (MS) and the solution on that grid. The theoretical local truncation error can be computed by using the modified equation analysis (MEA) with the MS and its analytic derivatives, which we call modified equation analysis method of manufactured solutions (MEAMMS). In addition to describing the MEAMMS process, this study shows the results of completing a code verification study on a conservation of mass code. The code was able to compute the leading truncation error term as well as additional higher-order terms. When the code verification process was complete, not only did the observed order-of-accuracy match the theoretical order-of-accuracy for all numerical schemes implemented in the code, but it was also able to cancel the discretization error to within roundoff error for a 64-bit system.

Numerical Method for The Time Fractional Fokker-Planck Equation

2012 ◽  
Vol 4 (06) ◽  
pp. 848-863 ◽  
Author(s):  
Xue-Nian Cao ◽  
Jiang-Li Fu ◽  
Hu Huang
Keyword(s):  
Theoretical Analysis ◽  
Numerical Method ◽  
Numerical Experiments ◽  
Truncation Error ◽  
Planck Equation ◽  
Local Truncation Error ◽  
Order Of Accuracy ◽  
Fokker Planck Equation ◽  
The Stability ◽  
Fokker Planck

AbstractIn this paper, a new numerical algorithm for solving the time fractional Fokker-Planck equation is proposed. The analysis of local truncation error and the stability of this method are investigated. Theoretical analysis and numerical experiments show that the proposed method has higher order of accuracy for solving the time fractional Fokker-Planck equation.

scholarly journals Direct construction of optimized stellarator shapes. Part 3. Omnigenity near the magnetic axis

2019 ◽  
Vol 85 (6) ◽  
Author(s):  
Gabriel G. Plunk ◽  
Matt Landreman ◽  
Per Helander
Keyword(s):  
Numerical Method ◽  
Plasma Phys ◽  
Magnetic Axis ◽  
Direct Construction ◽  
Numerical Construction ◽  
First Order

The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren & Boozer (Phys. Fluids B, vol. 3 (10), 1991a, pp. 2805–2821). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman & Sengupta (J. Plasma Phys., vol. 84 (6), 2018, 905840616); Landreman et al. (J. Plasma Phys., vol. 85 (1), 2019, 905850103), enables their efficient numerical construction.

scholarly journals Verification of Convergence Rates of Numerical Solutions for Parabolic Equations

2019 ◽  
Vol 2019 ◽  
pp. 1-10
Author(s):  
Darae Jeong ◽  
Yibao Li ◽  
Chaeyoung Lee ◽  
Junxiang Yang ◽  
Yongho Choi ◽  
...  
Keyword(s):  
Parabolic Equations ◽  
Convergence Rates ◽  
Numerical Solutions ◽  
Second Order ◽  
Order Convergence ◽  
Numerical Convergence ◽  
Verification Method ◽  
First Order ◽  
Convergence Test ◽  
Second Order Convergence

In this paper, we propose a verification method for the convergence rates of the numerical solutions for parabolic equations. Specifically, we consider the numerical convergence rates of the heat equation, the Allen–Cahn equation, and the Cahn–Hilliard equation. Convergence test results show that if we refine the spatial and temporal steps at the same time, then we have the second-order convergence rate for the second-order scheme. However, in the case of the first-order in time and the second-order in space scheme, we may have the first-order or the second-order convergence rates depending on starting spatial and temporal step sizes. Therefore, for a rigorous numerical convergence test, we need to perform the spatial and the temporal convergence tests separately.

A Difference Method for Plane Problems in Magnetoelastodynamics

1972 ◽  
Vol 39 (3) ◽  
pp. 689-695 ◽  
Author(s):  
W. W. Recker
Keyword(s):  
Partial Differential Equations ◽  
Differential Equations ◽  
Numerical Method ◽  
Difference Method ◽  
Necessary Condition ◽  
Two Dimensional ◽  
Symmetric Hyperbolic System ◽  
Von Neumann ◽  
First Order ◽  
Independent Variables

The two-dimensional equations of magnetoelastodynamics are considered as a symmetric hyperbolic system of linear first-order partial-differential equations in three independent variables. The characteristic properties of the system are determined and a numerical method for obtaining the solution to mixed initial and boundary-value problems in plane magnetoelastodynamics is presented. Results on the von Neumann necessary condition are presented. Application of the method to a problem which has a known solution provides further numerical evidence of the convergence and stability of the method.

scholarly journals Weak solvability of the unconditionally stable difference scheme for the coupled sine-Gordon system

2020 ◽  
Vol 25 (6) ◽  
pp. 997-1014
Author(s):  
Ozgur Yildirim ◽  
Meltem Uzun
Keyword(s):  
Difference Scheme ◽  
Numerical Method ◽  
Existence And Uniqueness ◽  
Numerical Experiments ◽  
Finite Difference Schemes ◽  
Order Of Accuracy ◽  
Unconditionally Stable ◽  
First Order ◽  
The Difference ◽  
Weak Solvability

In this paper, we study the existence and uniqueness of weak solution for the system of finite difference schemes for coupled sine-Gordon equations. A novel first order of accuracy unconditionally stable difference scheme is considered. The variational method also known as the energy method is applied to prove unique weak solvability.We also present a new unified numerical method for the approximate solution of this problem by combining the difference scheme and the fixed point iteration. A test problem is considered, and results of numerical experiments are presented with error analysis to verify the accuracy of the proposed numerical method.

Sign in / Sign up

Export Citation Format

Share Document