Difference between revisions of "Method of characteristics"
(Importing text file) |
Ulf Rehmann (talk | contribs) m (MR/ZBL numbers added) |
||
Line 63: | Line 63: | ||
Using the method of characteristics one can solve stationary multi-dimensional problems in a domain of hyperbolicity (for gas dynamics — of supersonic flow). One can also determine the position of secondary shock waves at places where the characteristics of a single family intersect or touch. By the method of characteristics one can tackle only problems with a small number of discontinuities, since under an accumulation of singularities the computations become tedious. A computation by the method of characteristics consists of a number of elementary problems: the calculation of an interior point, of a point on the shock wave, or on a body around which the flow is taking place, etc. | Using the method of characteristics one can solve stationary multi-dimensional problems in a domain of hyperbolicity (for gas dynamics — of supersonic flow). One can also determine the position of secondary shock waves at places where the characteristics of a single family intersect or touch. By the method of characteristics one can tackle only problems with a small number of discontinuities, since under an accumulation of singularities the computations become tedious. A computation by the method of characteristics consists of a number of elementary problems: the calculation of an interior point, of a point on the shock wave, or on a body around which the flow is taking place, etc. | ||
− | One can construct numerical schemes for the method of characteristics, which make it possible to carry out the calculations in | + | One can construct numerical schemes for the method of characteristics, which make it possible to carry out the calculations in "layers" , the grid-characteristic method (see [[#References|[2]]]). |
====References==== | ====References==== | ||
− | <table><TR><TD valign="top">[1]</TD> <TD valign="top"> , ''An attempt at calculating plane and axially-symmetric supersonic flows of a gas by the method of characteristics'' , Moscow | + | <table><TR><TD valign="top">[1]</TD> <TD valign="top"> , ''An attempt at calculating plane and axially-symmetric supersonic flows of a gas by the method of characteristics'' , Moscow (1961) (In Russian)</TD></TR><TR><TD valign="top">[2]</TD> <TD valign="top"> K.M. Magomedov, A.S. Kholodov, "The constructions of difference schemes for hyperbolic equations based on characteristic relations" ''USSR Math. Math. Comp.'' , '''9''' : 2 (1969) pp. 158–176 ''Zh. Vyshch. Mat. i Mat. Fiz.'' , '''9''' : 2 (1969) pp. 373–386</TD></TR></table> |
Line 74: | Line 74: | ||
====References==== | ====References==== | ||
− | <table><TR><TD valign="top">[a1]</TD> <TD valign="top"> | + | <table><TR><TD valign="top">[a1]</TD> <TD valign="top"> W.F. Ames, "Numerical methods for partial differential equations" , Acad. Press (1977) {{MR|0519123}} {{ZBL|0577.65077}} </TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top"> J. Massau, "Mémoire sur l'intégration graphique des équations aux derivées partielles" , F. Mayer-van Loo , Ghent (1899) {{MR|0052203}} {{ZBL|30.0312.01}} </TD></TR></table> |
Revision as of 16:59, 15 April 2012
A method for the numerical integration of equations of hyperbolic type. In the domain of hyperbolicity there is a linear combination of initial equations in which there occur only interior derivatives along characteristic surfaces. Then the equations to be solved simplify substantially. In the method of characteristics the solution is computed on a characteristic grid, which is constructed in the process of calculation, and so the domain of dependency of the solution can be determined exactly. For the method of characteristics the existence of a solution and its convergence have been proved. The widest application of the method of characteristics is in the solution of problems of the mechanics of continuous media (see [1]). For example, the equations in characteristic form
![]() | (1) |
![]() | (2) |
are a linear combination of the traditional equations of gas dynamics: continuity, momentum and energy. Here and below is the density,
is the velocity,
is the intrinsic energy of a unit mass,
is the pressure,
is the temperature,
is the spatial coordinate, and
is the time. The Cauchy problem asks for a solution in the domain
for given data on the line
. Entropy is the name for an integral
of the equation
![]() |
Then (1) has the form
![]() | (1prm) |
On the left-hand sides of (1prm) and (2) there stand the derivatives ,
,
, taken in the directions
![]() | (3) |
and
![]() | (4) |
the so-called characteristics. The system (1), (2) has three families of real characteristics. Along the characteristics (3) the relation
![]() |
holds, and along the characteristic (4) the relations
![]() | (5) |
Figure: m063580a
hold. Through the point (see Fig.) there passes a characteristic (4) on the side of increasing
,
![]() |
Through a point close to on the right, say
, there passes a characteristic (4) of another family
![]() |
where is the point of intersection of the characteristics. When one replaces the differential relations (5), which hold along characteristics, by differences, one obtains the algebraic system
![]() |
![]() |
from which and
can be determined. From
one draws the characteristic (3),
![]() |
up to the intersection with at
. The value of the entropy
at
is determined by means of interpolation between the points
and
(here
). From the equations
![]() |
![]() |
one can find the values of the intrinsic energy and the density
at
. When the data at the two points
and
are known, one can find the solution at
for large values of
. This procedure of computations is repeated for every pair of points. Then, by using new points
instead of the original
and
, one takes the next step in
. The calculation is made for the required values of
. However, since the equations of gas dynamics are non-linear, the calculation may come to a stop at a certain moment of time if the characteristics of one family touch each other or intersect.
The difference scheme thus described is of the first order of accuracy (an analogue of Euler's method of polygonal lines for solving ordinary differential equations). An increase in accuracy can be achieved by repetitions of the calculations, etc.
Using the method of characteristics one can solve stationary multi-dimensional problems in a domain of hyperbolicity (for gas dynamics — of supersonic flow). One can also determine the position of secondary shock waves at places where the characteristics of a single family intersect or touch. By the method of characteristics one can tackle only problems with a small number of discontinuities, since under an accumulation of singularities the computations become tedious. A computation by the method of characteristics consists of a number of elementary problems: the calculation of an interior point, of a point on the shock wave, or on a body around which the flow is taking place, etc.
One can construct numerical schemes for the method of characteristics, which make it possible to carry out the calculations in "layers" , the grid-characteristic method (see [2]).
References
[1] | , An attempt at calculating plane and axially-symmetric supersonic flows of a gas by the method of characteristics , Moscow (1961) (In Russian) |
[2] | K.M. Magomedov, A.S. Kholodov, "The constructions of difference schemes for hyperbolic equations based on characteristic relations" USSR Math. Math. Comp. , 9 : 2 (1969) pp. 158–176 Zh. Vyshch. Mat. i Mat. Fiz. , 9 : 2 (1969) pp. 373–386 |
Comments
The method of characteristics goes back to J. Massau (see [a2]).
References
[a1] | W.F. Ames, "Numerical methods for partial differential equations" , Acad. Press (1977) MR0519123 Zbl 0577.65077 |
[a2] | J. Massau, "Mémoire sur l'intégration graphique des équations aux derivées partielles" , F. Mayer-van Loo , Ghent (1899) MR0052203 Zbl 30.0312.01 |
Method of characteristics. Encyclopedia of Mathematics. URL: http://encyclopediaofmath.org/index.php?title=Method_of_characteristics&oldid=15839