Calculation of shock problems by using four different schemes

Conference · Sun Jan 01 04:00:00 EST 1984

OSTI ID:6872148

Lee, W H; Whalen, P P

Results are shown of the use of several different shock treatments in one- and two-dimensional Lagrangian code calculations of strong shock problems with known solutions. The shock treatments are (A) von Neumann-Richtmyer artificial viscosity; (B) Fixed length artificial viscosity; (C) Artificial energy diffusivity combined with artificial voscosity; and (D) Modified Godunov. The similarity test problems are plane and spherical implosions followed through a reflection. The problems are generated in the codes by imposing an inward directed velocity at one boundary of an initially, quiescent, gamma law gas with the other boundary fixed. Results are shown for calculations on uniform and non-uniform meshes. On non-uniform meshes, no method gives good results although method D is probably superior. Method B produces good appearing results with much shock smearing for the initial shock transit but the worst results after shock reflection. On uniform meshes, method C does the best job of handling the effects of shock initiation at boundaries while method A produces the worst results. Method D results in the smoothest flow field with less over-or-under-shooting, Gibbs phenomena. 14 references.

• Research Organization: Los Alamos National Lab., NM (USA)
• DOE Contract Number: W-7405-ENG-36
• OSTI ID: 6872148
• Report Number(s): LA-UR-84-344; CONF-840720-1; ON: DE84006032
• Country of Publication: United States
• Language: English

Related Subjects

658000* -- Mathematical Physics-- (-1987)
71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
COMPARATIVE EVALUATIONS
COMPUTER CALCULATIONS
FLUID FLOW
FUNCTIONS
IDEAL FLOW
LAGRANGIAN FUNCTION
MATERIALS
NUMERICAL SOLUTION
SHOCK WAVES
VISCOSITY
WAVE PROPAGATION