Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated

A study of numerical methods for hyperbolic conservation laws with stiff source terms

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained

A quick-retrieval high-speed digital framing camera

A new high-speed digital framing camera is described. The design is built around a rotating polygon mirror that provides a framing rate of 24 000 frames/s. The camera electronics digitizes an image into a 32×104 grid of pixels, where the second dimension of the grid can be varied and is determined by the 8 bit computer-aided measurement and control digitizer sampling rate. Available digitizer memory provides for 314 frames at this horizontal resolution. The advantages over other available high-speed framing cameras are (1) low cost of the system provided the digitizers are available, (2) rapid retrieval of a recorded event, and (3) the ease with which the system can be used. Sample results from an application in high-power arc photography are given to illustrate the system's spatial and temporal resolution

The Atmosphere Explorer and the shuttle glow

Recent analyses of the Atmosphere Explorer data are discussed in which it is demonstrated that the satellite glows have two components, one at high altitudes which is consistent with excitation in single collisions of atmospheric oxygen atoms with the vehicle surface and the other at low altitudes which is consistent with double collisions of nitrogen molecules. Contrary to an earlier suggestion, the low-altitude data are not consistent with collisions of oxygen molecules. The separation of the two components strengthens the conclusion that the high-altitude glow arises from vibrationally excited OH molecules produced by a formation mechanism that is different from that leading to the normal atmospheric OH airglow. The spectrum is consistent with association of oxygen and hydrogen atoms at sites on the surface into the vibrational levels of OH. The low-altitude glow is consistent with the green mechanism but there are difficulties with it. The shuttle glows are different and have the spectral appearance of emission from NO2. The characteristics of the shuttle glows and the satellite glows will be contrasted and a tentative resolution of the differences in the Atmosphere Explorer and shuttle glows will be offered

From the chiral magnetic wave to the charge dependence of elliptic flow

The quark-gluon plasma formed in heavy ion collisions contains charged chiral fermions evolving in an external magnetic field. At finite density of electric charge or baryon number (resulting either from nuclear stopping or from fluctuations), the triangle anomaly induces in the plasma the Chiral Magnetic Wave (CMW). The CMW first induces a separation of the right and left chiral charges along the magnetic field; the resulting dipolar axial charge density in turn induces the oppositely directed vector charge currents leading to an electric quadrupole moment of the quark-gluon plasma. Boosted by the strong collective flow, the electric quadrupole moment translates into the charge dependence of the elliptic flow coefficients, so that $v_2(\pi^+) < v_2(\pi^-)$ (at positive net charge). Using the latest quantitative simulations of the produced magnetic field and solving the CMW equation, we make further quantitative estimates of the produced $v_2$ splitting and its centrality dependence. We compare the results with the available experimental data.Comment: Contains 12 pages, 6 figures, written as a proceeding for the talk of Y. Burnier at the conference "P and CP-odd Effects in Hot and Dense Matter 2012" held in BN

Hypersonic blunt body computations including real gas effects

The recently developed second-order explicit and implicit total variation diminishing (TVD) shock-capturing methods of the Harten and Yee, Yee, and van Leer types in conjunction with a generalized Roe's approximate Riemann solver of Vinokur and the generalized flux-vector splittings of Vinokur and Montagne for two-dimensional hypersonic real gas flows are studied. A previous study on one-dimensional unsteady problems indicated that these schemes produce good shock-capturing capability and that the state equation does not have a large effect on the general behavior of these methods for a wide range of flow conditions for equilibrium air. The objective of this paper is to investigate the applicability and shock resolution of these schemes for two-dimensional steady-state hypersonic blunt body flows. The main contribution of this paper is to identify some of the elements and parameters which can affect the convergence rate for high Mach numbers or real gases but have negligible effect for low Mach number cases for steady-state inviscid blunt body flows