Asymptotic adaptive methods for multi-scale problems in fluid mechanics

This paper reports on the results of a three-year research effort aimed at investigating and exploiting the role of physically motivated asymptotic an alysis in the design of numerical methods for singular limit problems in fl uid mechanics. Such problems naturally arise, among others, in combustion, magneto-hydrodynamics, and geophysical fluid mechanics. Typically, they are characterized by multiple-space and/or -time scales and by the disturbing fact that standard computational techniques fail entirely, are unacceptably expensive, or both. The challenge here is to construct numerical methods w hich are robust, uniformly accurate, and efficient through different asympt otic regimes and over a wide range of relevant applications. Summaries of m ultiple-scales asymptotic analyses for low-Mach-number flows, magneto-hydro dynamics at small Mach and Alfven numbers, and of multiple-scales atmospher ic flows are provided. These reveal singular balances between selected term s in the respective governing equations within the considered flow regimes. These singularities give rise to problems of severe stiffness, stability, or to dynamic-range issues in straight-forward numerical discretizations. A formal mathematical framework for the multiple scales asymptotics is then summarized by use of the example of multiple length-scale single-time-scale asymptotics for low-Mach-number flows. The remainder of the paper focuses on the construction of numerical discretizations for the respective full go verning equation systems. These discretizations avoid the pitfalls of singu lar balances by exploiting the asymptotic results. Importantly, the asympto tics are not used here to derive simplified equation systems, which are the n solved numerically. Rather, numerical integration of the full equation se ts is aimed at and the asymptotics are used only to construct discretizatio ns that do not deteriorate as a singular limit is approached. One important ingredient of this strategy is the numerical identification of a singular limit regime given a set of discrete numerical state variables. This proble m is addressed in an exemplary fashion for multiple-length single-time-scal e low-Mach-number flows in one space dimension. The strategy allows a dynam ic determination of an instantaneous relevant Mach number, and it can thus be used to drive the appropriate adjustment of the numerical discretization s when the singular limit regime is approached.