Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)



Hydrocode subcycling stability

Author: D. L. Hicks
Journal: Math. Comp. 37 (1981), 69-78
MSC: Primary 65M10; Secondary 76-04, 76-08
MathSciNet review: 616360
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: The method of artificial viscosity was originally designed by von Neumann and Richtmyer for calculating the propagation of waves in materials that were hydrodynamic and rate-independent (e.g., ideal gas law). However, hydrocodes (such as WONDY) based on this method continue to expand their repertoire of material laws even unto material laws that are rate-dependent (e.g., Maxwell’s material law). Restrictions on the timestep required for stability with material laws that are rate-dependent can be considerably more severe than restrictions of the Courant-Friedrichs-Lewy (CFL) type that are imposed in these hydrocodes. These very small timesteps can make computations very expensive. An alternative is to go ahead and integrate the conservation laws with the usual CFL timestep while subcycling (integrating with a smaller timestep) the integration of the stress-rate equation. If the subcycling is done with a large enough number of subcycles (i.e., with a small enough subcycle timestep), then the calculation is stable. Specifically, the number of subcycles must be one greater than the ratio of the CFL timestep to the relaxation time of the material.

References [Enhancements On Off] (What's this?)

Similar Articles

Retrieve articles in Mathematics of Computation with MSC: 65M10, 76-04, 76-08

Retrieve articles in all journals with MSC: 65M10, 76-04, 76-08

Additional Information

Article copyright: © Copyright 1981 American Mathematical Society