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 ◽  
...  

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.

1987 ◽  
Vol 109 (3) ◽  
pp. 213-217 ◽  
Author(s):  
S. Abdallah ◽  
R. E. Henderson

Quasi three dimensional blade-to-blade solutions for stators and rotors of turbomachines are obtained using the Streamline Curvature Method (SLCM). The first-order velocity gradient equation of the SLCM, traditionally solved for the velocity field, is reformulated as a second-order elliptic differential equation and employed in tracing the streamtubes throughout the flow field. The equation of continuity is then used to calculate the velocity. The present method has the following advantages. First, it preserves the ellipticity of the flow field in the solution of the second-order velocity gradient equation. Second, it eliminates the need for curve fitting and strong smoothing under-relaxation in the classical SLCM. Third, the prediction of the stagnation streamlines is a straightforward matter which does not complicate the present procedure. Finally, body-fitted curvilinear coordinates (streamlines and orthogonals or quasi-orthogonals) are naturally generated in the method. Numerical solutions are obtained for inviscid incompressible flow in rotating and non-rotating passages and the results are compared with experimental data.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
R. Naz ◽  
I. Naeem ◽  
F. M. Mahomed

This paper analyzes the first integrals and exact solutions of mathematical models of epidemiology via the partial Lagrangian approach by replacing the three first-order nonlinear ordinary differential equations by an equivalent system containing one second-order equation and a first-order equation. The partial Lagrangian approach is then utilized for the second-order ODE to construct the first integrals of the underlying system. We investigate the SIR and HIV models. We obtain two first integrals for the SIR model with and without demographic growth. For the HIV model without demography, five first integrals are established and two first integrals are deduced for the HIV model with demography. Then we utilize the derived first integrals to construct exact solutions to the models under investigation. The dynamic properties of these models are studied too. Numerical solutions are derived for SIR models by finite difference method and are compared with exact solutions.


2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Don Liu ◽  
Qin Chen ◽  
Yifan Wang

A system of coupled nonlinear partial differential equations with convective and dispersive terms was modified from Boussinesq-type equations. Through a special formulation, a system of nonlinear partial differential equations was solved alternately and explicitly in time without linearizing the nonlinearity. Coupled compact schemes of sixth order accuracy in space were developed to obtain numerical solutions. Within couple compact schemes, variables and their first and second derivatives were solved altogether. The sixth order accuracy in space is achieved with a memory-saving arrangement of state variables so that the linear system is banded instead of blocked. This facilitates solving very large systems. The efficiency, simplicity, and accuracy make this coupled compact method viable as variational and weighted residual methods. Results were compared with exact solutions which were obtained via devised forcing terms. Error analyses were carried out, and the sixth order convergence in space and second order convergence in time were demonstrated. Long time integration was also studied to show stability and error convergence rates.


Mathematics ◽  
2018 ◽  
Vol 6 (12) ◽  
pp. 306 ◽  
Author(s):  
Juan Liang ◽  
Linke Hou ◽  
Xiaowu Li ◽  
Feng Pan ◽  
Taixia Cheng ◽  
...  

Orthogonal projection a point onto a parametric curve, three classic first order algorithms have been presented by Hartmann (1999), Hoschek, et al. (1993) and Hu, et al. (2000) (hereafter, H-H-H method). In this research, we give a proof of the approach’s first order convergence and its non-dependence on the initial value. For some special cases of divergence for the H-H-H method, we combine it with Newton’s second order method (hereafter, Newton’s method) to create the hybrid second order method for orthogonal projection onto parametric curve in an n-dimensional Euclidean space (hereafter, our method). Our method essentially utilizes hybrid iteration, so it converges faster than current methods with a second order convergence and remains independent from the initial value. We provide some numerical examples to confirm robustness and high efficiency of the method.


Sign in / Sign up

Export Citation Format

Share Document