Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms=1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock     

冲击波作用下Air/SF6斜界面不稳定性实验和数值模拟研究

开展了Mach数为1.23和1.41的冲击波作用下的Air/SF6斜界面不稳定性激波管实验,并利用王涛等人发展的可压缩多介质粘性流体和湍流大涡模拟程序MVFT(multi-viscous-fluid and turbulence),对该激波管实验进行了数值模拟,二者相比较一致性较好,包括界面图像、湍流混合区TMZ(turbulent mixing zone)宽度、气泡和尖钉位移,确认了该计算代码对界面不稳定性问题模拟的可靠性和有效性.数值模拟再现了冲击波作用下,Air/SF6斜界面的演化过程及流动中复杂波系结构的发展如冲击波的传播、折射和反射.结果还显示冲击波Mach数较大时,冲击波和界面相互作用时混合区获得的能量也较大,扰动界面发展的也更快.     

固体平板上超薄水膜去湿润过程的分子动力学模拟

运用分子动力学模拟研究固体表面超薄水膜的失稳和破裂过程．结果表明薄膜中小扰动将失稳,并在初始阶段线性增长．但固体和液体的相互作用对扰动的初期增长影响较小．最小厚度的下降导致薄膜发生破裂．此后破裂边缘以一定的动态接触角后退．与宏观理论预测一致,边缘半径随时间的变化与时间平方根成正比．若固液相互作用较强,将引起破裂时间延迟,动态接触角减小,固体表面附近的液体密度增加．     

Global axisymmetric classical solutions of full compressible magnetohydrodynamic equations with vacuum free boundary and large initial data

In this paper,the global existence of the classical solution to the vacuum free boundary problem of full compressible magnetohydrodynamic equations with large initial data and axial symmetry is studied.The solutions to the system(1.6)–(1.8) are in the class of radius-dependent solutions,i.e.,independent of the axial variable and the angular variable.In particular,the expanding rate of the moving boundary is obtained.The main difficulty of this problem lies in the strong coupling of the magnetic field,velocity,temperature and the degenerate density near the free boundary.We overcome the obstacle by establishing the lower bound of the temperature by using different Lagrangian coordinates,and deriving the uniform-in-time upper and lower bounds of the Lagrangian deformation variable r;by weighted estimates,and also the uniform-in-time weighted estimates of the higher-order derivatives of solutions by delicate analysis.     

The growth of fractal dimension of an interface evolution from the interaction of a shock wave with a rectangular block of SF6

The interface between air and a rectangular block of sulphur hexafluoride (SF6), impulsively accelerated by the passage of a planar shock wave, undergoes Richtmyer–Meshkov instability and the flow becomes turbulent. The evolution of the interface was previously simulated using a multi-component model based on a thermodynamically consistent and fully conservative formulation and results were validated against available experimental data (Bates et al. Richtmyer–Meshkov instability induced by the interaction of a shock wave with a rectangular block of SF6, Phys Fluids, 2007; 19:036101). In this study, the CFD results are analyzed using the fractal theory approach and the evolution of fractal dimension of the interface during the transition of the flow into fully developed turbulence is measured using the standard box-counting method. It is shown that as the Richtmyer–Meshkov instability on the interface develops and the flow becomes turbulent, the fractal dimension of the interface increases asymptotically toward a value close to 1.39, which agrees well to those measured for classical shear and fully developed turbulences.     

一维Chaplygin气体欧拉方程组中波的相互作用

研究一维Chaplygin气体欧拉方程组中波的相互作用.方程组的波包含接触间断和在密度变量以及内能变量上同时具有狄拉克函数的狄拉克激波.根据这些波的不同组合,问题被分成了7种情形.通过详细地构造每种情形的整体解,获得了各种波相互作用的完整结果.特别地,对于一类初值,两个接触间断相互作用后,产生了一个狄拉克激波;然而,对于另外一类初值,一个狄拉克激波与一个接触间断相互作用后,狄拉克激波消失.这些都是波相互作用中非常特别的现象.     

Existence in the Large for Pressure-Gradient System

In this paper, the authors use Glimm scheme to study the global existence of BV solutions to Cauchy problem of the pressure-gradient system with large initial data. To this end, some important properties of the shock curves of the pressure-gradient system in the Riemann invariant coordinate system and verify that the shock curves satisfy Diperna’s conditions (see [Diperna, R. J., Existence in the large for quasilinear hyperbolic conservation laws, Arch. Ration. Mech. Anal., 52(3), 1973, 244–257]) are studied. Then they construct the approximate solution sequence through Glimm scheme. By establishing accurate local interaction estimates, they prove the boundedness of the approximate solution sequence and its total variation.     

Reflection and refraction of a plane thermoelastic wave at a solid–solid interface under perfect boundary condition,in presence of normal initial stress

Using the basic governing equations for isotropic and homogeneous generalized thermo elastic media under initial stress, the reflection and refraction of thermo elastic plane waves at the interface of two dissimilar thermo elastic solid half-spaces has been investigated. The amplitude ratios of various reflected and refracted waves are obtained for an ideal boundary for the incidence of SV-wave. The numerical computations are carried out for a particular model. The effect of initial stress on the amplitude ratios are shown graphically after numerical calculation.     

Continuous compression waves in the two-dimensional Riemann problem

The interaction between a plane shock wave in a plate and a wedge is considered within the framework of the nondissipative compressible fluid dynamic equations. The wedge is filled with a material that may differ from that of the plate. Based on the numerical solution of the original equations, self-similar solutions are obtained for several versions of the problem with an iron plate and a wedge filled with aluminum and for the interaction of a shock wave in air with a rigid wedge. The behavior of the solids at high pressures is approximately described by a two-term equation of state. In all the problems, a two-dimensional continuous compression wave develops as a wave reflected from the wedge or as a wave adjacent to the reflected shock. In contrast to a gradient catastrophe typical of one-dimensional continuous compression waves, the spatial gradient of a two-dimensional compression wave decreases over time due to the self-similarity of the solution. It is conjectured that a phenomenon opposite to the gradient catastrophe can occur in an actual flow with dissipative processes like viscosity and heat conduction. Specifically, an initial shock wave is transformed over time into a continuous compression wave of the same amplitude.     

Potential theory for regular and mach reflection of a shock at a wedge

If a plane shock hits a wedge, a self-similar pattern of reflected shocks travels outward as the shock moves forward in time. The nature of the pattern is explored for weak incident shocks (strength b) and small wedge angles 2θ_w through potential theory, a number of different scalings, some study of mixed equations and matching asymptotics for the different scalings. The self-similar equations are of mixed type. A linearization gives a linear mixed flow valid away from a sonic curve. Near the sonic curve a shock solution is constructed in another scaling except near the zone of interaction between the incident shock and the wall where a special scaling is used. The parameter β = c₁θ²_w(γ + 1)b ranges from 0 to ∞. Here γ is the polytropic constant and C₁ is the sound speed behind the incident shock. For β > 2 regular reflection (weak or strong) can occur and the whole pattern is reconstructed to lowest order in shock strength. For β < 1/2 Mach reflection occurs and the flow behind the reflection is subsonic and can be constructed in principle (with an open elliptic problem) and matched. The case β = 0 can be solved. For 1/2 < β < 2 or even larger β the flow behind a Mach reflection may be transonic and further investigation must be made to determine what happens. The basic pattern of reflection is an almost semi-circular shock issuing, for regular reflection, from the reflection point on the wedge and for Mach reflection, matched with a local interaction flow. Assuming their nature, choosing the least entropy generation, the weak regular reflection will occur for β sufficiently large (von Neumann paradox). An accumulation point of vorticity occurs on the wedge above the leading point. © 1994 John Wiley & Sons, Inc.     

On finite time BV blow-up for the p-system

Abstract

The paper studies the possible blowup of the total variation for entropy weak solutions of the p-system, modeling isentropic gas dynamics. It is assumed that the density remains uniformly positive, while the initial data can have arbitrarily large total variation (measured in terms of Riemann invariants). Two main results are proved. (I) If the total variation blows up in finite time, then the solution must contain an infinite number of large shocks in a neighborhood of some point in the t-x plane. (II) Piecewise smooth approximate solutions can be constructed whose total variation blows up in finite time. For these solutions the strength of waves emerging from each interaction is exact, while rarefaction waves satisfy the natural decay estimates stemming from the assumption of genuine nonlinearity.     

GLOBAL SOLUTIONS TO THE NAVIER-STOKES EQUATIONS FOR COMPRESSIBLE HEAT-CONDUCTING FLOW WITH SYMMETRY AND FREE BOUNDARY

ABSTRACT

Global solutions of the multidimensional Navier-Stokes equations for compressible heat-conducting flow are constructed, with spherically symmetric initial data of large oscillation between a static solid core and a free boundary connected to a surrounding vacuum state. The free boundary connects the compressible heat-conducting fluids to the vacuum state with free normal stress and zero normal heat flux. The fluids are initially assumed to fill with a finite volume and zero density at the free boundary, and with bounded positive density and temperature between the solid core and the initial position of the free boundary. One of the main features of this problem is the singularity of solutions near the free boundary. Our approach is to combine an effective difference scheme to construct approximate solutions with the energy methods and the pointwise estimate techniques to deal with the singularity of solutions near the free boundary and to obtain the bounded estimates of the solutions and the free boundary as time evolves. The convergence of the difference scheme is established. It is also proved that no vacuum develops between the solid core and the free boundary, and the free boundary expands with finite speed.     

Initial Algebraic Growth of Small Angular Dependent Disturbances in Pipe Poiseuille Flow

The evolution of small, angular dependent velocity disturbances in laminar pipe flow is studied. In particular, streamwise independent perturbations are considered. To fully describe the flow field, two equations are required, one for the radial and the other for the streamwise velocity perturbation. Whereas the former is homogeneous, the latter has the radial velocity component as a forcing term. First, the normal modes of the system are determined and analytical solutions for eigenfunctions, damping rates, and phase velocities are calculated. As the azimuthal wave number (n) increases, the damping rate increases and the phase velocities decrease. Particularly interesting are results showing that the phase velocities associated with the streamwise eigenfunctions are independent of the radial mode index when n = 1, and when n = 5 the same is obtained for phase velocities associated with the eigenfunctions of the radial component. Then, the initial value problem is treated and the time development of the disturbances is determined. The radial and the azimuthal velocity components always decay but, owing to the forcing, the streamwise component shows an initial algebraic growth, followed by a decay. The kinetic energy density is used to characterize the induced streamwise disturbance. Its dependence on the Reynolds number, the radial mode, and the azimuthal wave number is investigated. With a normalized initial disturbance, n = 1 gives the largest amplification, followed by n = 2 etc. However, for small times, higher values of n are associated with the largest energy density. As n increases, the distribution of the streamwise velocity perturbation becomes more concentrated to the region near the pipe wall.     

A difference scheme with anti-diffusion terms to improve the computation of contact discontinuities and the study of weak discontinuities in compressible flow

The investigation of Mach reflection formed after the impingement of a weak plane shock wave on a wedge with shock Mach number M_s near 1, is still an open problem[12]. It's difficult for shock tube experiments with interferometer to detect contact discontinuities if it is too weak; also difficult to catch with due accuracy the transition condition between Mach reflection and regular reflection. The interest to this phenomenon is continuing, especially for weak shocks, because there was systematic discrepancy between simplified three shock theory of von Neumann [8] and shock tube results [15] which was named by G. Birkhoff as “von Neumann Paradox on three shock theory” [18].In 1972, K.O.Friedrichs called for more computational efforts on this problem. Recently it is known that for weak impinging shocks it's still difficult to get contact discontinuities and curved Mach stem with satisfactory accuracy. Recent numerical computation sometimes even fails to show reflected shock wave[6]. These explain why von Neumann paradox of the three shock theory in case of weak discontinuities is still a problem of interesting [9,12,14]. In this paper, on one hand, we investigate the numerical methods for Euler's equation for compressible inviscid flow, aiming at improving the computation of contact discontinuities, on the other hand, a methodology is suggested to correctly plot flow data from the massive information in storage. On this basis, all the reflected shock wave , contact discontinuities and the curved Mach stem are determined. We get Mach reflection under the condition when over-simplified shock theory predicts no such configuration[5].     

The Stability of Large‐Amplitude Shallow Interfacial Non‐Boussinesq Flows

The system of equations describing the shallow‐water limit dynamics of the interface between two layers of immiscible fluids of different densities is formulated. The flow is bounded by horizontal top and bottom walls. The resulting equations are of mixed type: hyperbolic when the shear is weak and the behavior of the system is internal‐wave like, and elliptic for strong shear. This ellipticity, or ill‐posedness is shown to be a manifestation of large‐scale shear instability. This paper gives sharp nonlinear stability conditions for this nonlinear system of equations. For initial data that are initially hyperbolic, two different types of evolution may occur: the system may remain hyperbolic up to internal wave breaking, or it may become elliptic prior to wave breaking. Using simple waves that give a priori bounds on the solutions, we are able to characterize the condition preventing the second behavior, thus providing a long‐time well‐posedness, or nonlinear stability result. Our formulation also provides a systematic way to pass to the Boussinesq limit, whereby the density differences affect buoyancy but not momentum, and to recover the result that shear instability cannot occur from hyperbolic initial data in that case.     

Existence and continuous dependence of large solutions for the magnetohydrodynamic equations

Global solutions of the nonlinear magnetohydrodynamic (MHD) equations with general large initial data are investigated. First the existence and uniqueness of global solutions are established with large initial data in H ¹. It is shown that neither shock waves nor vacuum and concentration are developed in a finite time, although there is a complex interaction between the hydrodynamic and magnetodynamic effects. Then the continuous dependence of solutions upon the initial data is proved. The equivalence between the well-posedness problems of the system in Euler and Lagrangian coordinates is also showed.     

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach numberMc = 0.8 andRe = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of A -vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.     

Construction of Solutions and L^1-error Estimates of Viscous Methods for Scalar Conservation Laws with Boundary

This paper is concerned with an initial boundary value problem for strictly convex conservation laws whose weak entropy solution is in the piecewise smooth solution class consisting of finitely many discontinuities. By the structure of the weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux Nedelec, we give a construction method to the weak entropy solution of the initial boundary value problem. Compared with the initial value problem, the weak entropy solution of the initial boundary value problem includes the following new interaction type： an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary. According to the structure and some global estimates of the weak entropy solution, we derive the global L^1-error estimate for viscous methods to this initial boundary value problem by using the matching travelling wave solutions method. If the inviscid solution includes the interaction that an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary, or the inviscid solution includes some shock wave which is tangent to the boundary, then the error of the viscosity solution to the inviscid solution is bounded by O（ε^1/2） in L^1-norm; otherwise, as in the initial value problem, the L^1-error bound is O（ε｜ In ε｜）.     

Global existence of shock front solution to axially symmetric piston problem in compressible flow

In this paper we study the global existence of BV solution to two-dimensional piston problem in fluid dynamics. Different from previous results on related problems we remove the restriction on the strength of the leading shock and require the velocity of the piston is rather fast or the density is quite small instead. The main tool in our proof is Glimm Scheme with some improvement. To define the Glimm functional we derive more precise estimates for the interaction of elementary waves, particularly in the region near the leading shock.     

Existence and nonexistence of TV bounds for scalar conservation laws with discontinuous flux

For the scalar conservation laws with discontinuous flux, an infinite family of (A, B)‐interface entropies are introduced and each one of them is shown to form an L¹‐contraction semigroup (see [2]). One of the main unsettled questions concerning conservation law with discontinuous flux is boundedness of total variation of the solution. Away from the interface, boundedness of total variation of the solution has been proved in a recent paper [6]. In this paper, we discuss this particular issue in detail and produce a counterexample to show that the solution, in general, has unbounded total variation near the interface. In fact, this example illustrates that smallness of the BV norm of the initial data is immaterial. We hereby settle the question of determining for which of the aforementioned (A, B) pairs the solution will have bounded total variation in the case of strictly convex fluxes. © 2010 Wiley Periodicals, Inc.