2,534 research outputs found
Hyperbolic/parabolic development for the GIM-STAR code
Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated
A parallel nearly implicit time-stepping scheme
Across-the-space parallelism still remains the most mature, convenient and natural way to parallelize large scale problems. One of the major problems here is that implicit time stepping is often difficult to parallelize due to the structure of the system. Approximate implicit schemes have been suggested to circumvent the problem. These schemes have attractive stability properties and they are also very well parallelizable.\ud
The purpose of this article is to give an overall assessment of the parallelism of the method
Unconditional Stability for Multistep ImEx Schemes: Theory
This paper presents a new class of high order linear ImEx multistep schemes
with large regions of unconditional stability. Unconditional stability is a
desirable property of a time stepping scheme, as it allows the choice of time
step solely based on accuracy considerations. Of particular interest are
problems for which both the implicit and explicit parts of the ImEx splitting
are stiff. Such splittings can arise, for example, in variable-coefficient
problems, or the incompressible Navier-Stokes equations. To characterize the
new ImEx schemes, an unconditional stability region is introduced, which plays
a role analogous to that of the stability region in conventional multistep
methods. Moreover, computable quantities (such as a numerical range) are
provided that guarantee an unconditionally stable scheme for a proposed
implicit-explicit matrix splitting. The new approach is illustrated with
several examples. Coefficients of the new schemes up to fifth order are
provided.Comment: 33 pages, 7 figure
An extension of A-stability to alternating direction implicit methods
An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented
Application of Operator Splitting Methods in Finance
Financial derivatives pricing aims to find the fair value of a financial
contract on an underlying asset. Here we consider option pricing in the partial
differential equations framework. The contemporary models lead to
one-dimensional or multidimensional parabolic problems of the
convection-diffusion type and generalizations thereof. An overview of various
operator splitting methods is presented for the efficient numerical solution of
these problems.
Splitting schemes of the Alternating Direction Implicit (ADI) type are
discussed for multidimensional problems, e.g. given by stochastic volatility
(SV) models. For jump models Implicit-Explicit (IMEX) methods are considered
which efficiently treat the nonlocal jump operator. For American options an
easy-to-implement operator splitting method is described for the resulting
linear complementarity problems.
Numerical experiments are presented to illustrate the actual stability and
convergence of the splitting schemes. Here European and American put options
are considered under four asset price models: the classical Black-Scholes
model, the Merton jump-diffusion model, the Heston SV model, and the Bates SV
model with jumps
- …