29,561 research outputs found
Application of the method of lines for solutions of the Navier-Stokes equations using a nonuniform grid distribution
The feasibility of the method of lines for solutions of physical problems requiring nonuniform grid distributions is investigated. To attain this, it is also necessary to investigate the stiffness characteristics of the pertinent equations. For specific applications, the governing equations considered are those for viscous, incompressible, two dimensional and axisymmetric flows. These equations are transformed from the physical domain having a variable mesh to a computational domain with a uniform mesh. The two governing partial differential equations are the vorticity and stream function equations. The method of lines is used to solve the vorticity equation and the successive over relaxation technique is used to solve the stream function equation. The method is applied to three laminar flow problems: the flow in ducts, curved-wall diffusers, and a driven cavity. Results obtained for different flow conditions are in good agreement with available analytical and numerical solutions. The viability and validity of the method of lines are demonstrated by its application to Navier-Stokes equations in the physical domain having a variable mesh
Method of lines transpose: High order L-stable O(N) schemes for parabolic equations using successive convolution
We present a new solver for nonlinear parabolic problems that is L-stable and
achieves high order accuracy in space and time. The solver is built by first
constructing a single-dimensional heat equation solver that uses fast O(N)
convolution. This fundamental solver has arbitrary order of accuracy in space,
and is based on the use of the Green's function to invert a modified Helmholtz
equation. Higher orders of accuracy in time are then constructed through a
novel technique known as successive convolution (or resolvent expansions).
These resolvent expansions facilitate our proofs of stability and convergence,
and permit us to construct schemes that have provable stiff decay. The
multi-dimensional solver is built by repeated application of dimensionally
split independent fundamental solvers. Finally, we solve nonlinear parabolic
problems by using the integrating factor method, where we apply the basic
scheme to invert linear terms (that look like a heat equation), and make use of
Hermite-Birkhoff interpolants to integrate the remaining nonlinear terms. Our
solver is applied to several linear and nonlinear equations including heat,
Allen-Cahn, and the Fitzhugh-Nagumo system of equations in one and two
dimensions
On the fourth-order accurate compact ADI scheme for solving the unsteady Nonlinear Coupled Burgers' Equations
The two-dimensional unsteady coupled Burgers' equations with moderate to
severe gradients, are solved numerically using higher-order accurate finite
difference schemes; namely the fourth-order accurate compact ADI scheme, and
the fourth-order accurate Du Fort Frankel scheme. The question of numerical
stability and convergence are presented. Comparisons are made between the
present schemes in terms of accuracy and computational efficiency for solving
problems with severe internal and boundary gradients. The present study shows
that the fourth-order compact ADI scheme is stable and efficient
On Asymptotic Global Error Estimation and Control of Finite Difference Solutions for Semilinear Parabolic Equations
The aim of this paper is to extend the global error estimation and control
addressed in Lang and Verwer [SIAM J. Sci. Comput. 29, 2007] for initial value
problems to finite difference solutions of semilinear parabolic partial
differential equations. The approach presented there is combined with an
estimation of the PDE spatial truncation error by Richardson extrapolation to
estimate the overall error in the computed solution. Approximations of the
error transport equations for spatial and temporal global errors are derived by
using asymptotic estimates that neglect higher order error terms for
sufficiently small step sizes in space and time. Asymptotic control in a
discrete -norm is achieved through tolerance proportionality and uniform
or adaptive mesh refinement. Numerical examples are used to illustrate the
reliability of the estimation and control strategies
Recommended from our members
The numerical solution of stefan problems with front-tracking and smoothing methods
Recommended from our members
A numerical study of heat flow in an elliptic cylinder with interior derivative boundary conditions
- …