Entropy in self-similar shock profiles

In this paper, we study the structure of a gaseous shock, and in particular the distribution of entropy within, in both a thermodynamics and a statistical mechanics context. The problem of shock structure has a long and distinguished history that we review. We employ the Navier–Stokes equations to construct a self-similar version of Becker’s solution for a shock assuming a particular (physically plausible) Prandtl number; and that solution reproduces the well-known result of Morduchow & Libby that features a maximum of the equilibrium entropy inside the shock profile. We then construct an entropy profile, based on gas kinetic theory, that is smooth and monotonically increasing. The extension of equilibrium thermodynamics to irreversible processes is based in part on the assumption of local thermodynamic equilibrium. We show that this assumption is not valid except for the weakest shocks. We conclude by hypothesizing a thermodynamic nonequilibrium entropy and demonstrating that it closely estimates the gas kinetic nonequilibrium entropy within a shock.     

Viscous regularization of the full set of nonequilibrium‐diffusion Grey Radiation‐Hydrodynamic equations

A viscous regularization technique, based on the local entropy residual, was proposed by Delchini et al. (2015) to stabilize the nonequilibrium‐diffusion Grey Radiation‐Hydrodynamic equations using an artificial viscosity technique. This viscous regularization is modulated by the local entropy production and is consistent with the entropy minimum principle. However, Delchini et al. (2015) only based their work on the hyperbolic parts of the Grey Radiation‐Hydrodynamic equations and thus omitted the relaxation and diffusion terms present in the material energy and radiation energy equations. Here, we extend the theoretical grounds for the method and derive an entropy minimum principle for the full set of nonequilibrium‐diffusion Grey Radiation‐Hydrodynamic equations. This further strengthens the applicability of the entropy viscosity method as a stabilization technique for radiation‐hydrodynamic shock simulations. Radiative shock calculations using constant and temperature‐dependent opacities are compared against semi‐analytical reference solutions, and we present a procedure to perform spatial convergence studies of such simulations. Copyright © 2017 John Wiley & Sons, Ltd.     

A physical model of the structure and attenuation of shock waves in metals

In this paper, a physical model of the structure and attenuation of shock waves in metals is presented. In order to establish the constitutive equations of materials under high velocity deformation and to study the structure of transition zone of shock wave, two independent approaches are involved. Firstly, the specific internal energy is decomposed into the elastic compression energy and elastic deformation energy, and the later is represented by an expansion to third-order terms in elastic strain and entropy, including the coupling effect of heat and mechanical energy. Secondly, a plastic relaxation function describing the behaviour of plastic flow under high temperature and high pressure is suggested from the viewpoint of dislocation dynamics. In addition, a group of ordinary differential equations has been built to determine the thermo-mechanical state variables in the transition zone of a steady shock wave and the thickness of the high pressure shock wave, and an analytical solution of the equations can be found provided that the entropy change across the shock is assumed to be negligible and Hugoniot compression modulus is used instead of the isentropic compression modulus. A quite approximate method for solving the attenuation of shock wave front has been proposed for the flat-plate symmetric impact problem.     

Nonexistence of a shock layer in gas dynamics with a nonconvex equation of state

A classical result of Gilbarg states that a simple shock wave solution of Euler's equations is compressive if and only if a corresponding shock layer solution of the Navier-Stokes equations exists, assuming, among other things, that the equation of state is convex. An entropy condition appropriate for weeding out unphysical shocks in the nonconvex case has been introduced by T.-P. Liu. For shocks satisfying his entropy condition, Liu showed that purely viscous shock layers exist (with zero heat conduction). Dropping the convexity assumption, but retaining many other reasonable restrictions on the equation of state, we construct an example of a (large amplitude) shock which satisfies Liu's entropy condition but for which a shock layer does not exist if heat conduction dominates viscosity. We also give a simple restriction, weaker than convexity, which does guarantee that shocks which satisfy Liu's entropy condition always admit shock layers.     

Numerical simulations of high-speed chemically reacting flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacing flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.This work has been supported by the Aerospace Corporation through a Corporate Fellowship, by NASA Dryden Flight Research Center under Grant NCC 2-374, by ONR Grant N00014-86-K-0691, by NSF Grant DMS 88-11863, and by a DARPA Grant in the ACMP Program.     

Entropy generation analysis of natural convective radiative second grade nanofluid flow between parallel plates in a porous medium

The present article explores the entropy generation of radiating viscoelastic second grade nanofluid in a porous channel confined between two parallel plates. The boundaries of the plates are maintained at distinct temperatures and concentrations while the fluid is being sucked and injected periodically through upper and lower plates. The buoyancy forces, thermophoresis and Brownian motion are also considered due to the temperature and concentration differences across the channel. The system of governing partial differential equations has been transferred into a system of ordinary differential equations(ODEs) by appropriate similarity relations, and a shooting method with the fourth-order Runge-Kutta scheme is used for the solutions. The results are analyzed in detail for dimensionless velocity components. The temperature, concentration distributions, the entropy generation number, and the Bejan number corresponding to various fluid and geometric parameters are shown graphically. The skin friction, heat and mass transfer rates are presented in the form of tables. It is noticed that the temperature profile of the fluid is enhanced with the Brownian motion, whereas the concentration profile of the fluid is decreased with the thermophoresis parameter, and the entropy and Bejan numbers exhibit the opposite trend for the suction and injection ratio.     

A model for adiabatic supersonic turbulent boundary layers

An asymptotic analysis of the equations describing supersonic turbulent flow over an adiabatic wall is carried out for high Reynolds numbers, Re, and mainstream Mach numbers, M _e=O(1). A general expression for the adiabatic-wall temperature is derived. The asymptotic theory constrains the types of turbulence models that are suitable to represent the effects of viscous dissipation. A simple algebraic turbulence model is proposed and comparisons with measured total enthalpy profile data show good agreement, capturing the overshoot observed in total enthalpy near the boundarylayer edge.This work was supported by NASA Langley Research Center under Grant NAG-1-832 and the Air Force Office of Scientific Research under Grants AFOSR-91-0069 and F49620-93-0130; Dr. Ruban was supported by a grant from United Technologies Corporation.     

Shock wave in a gas-liquid medium

The propagation of weak shock waves and the conditions for their existence in a gas-liquid medium are studied in [1]. The article [2] is devoted to an examination of powerful shock waves in liquids containing gas bubbles. The possibility of the existence in such a medium of a shock wave having an oscillatory pressure profile at the front is demonstrated in [3] based on the general results of nonlinear wave dynamics. It is shown in [4, 5] that a shock wave in a gas-liquid mixture actually has a profile having an oscillating pressure. The drawback of [3–5] is the necessity of postulating the existence of the shock waves. This is connected with the absence of a direct calculation of the dissipative effects in the fundamental equations. The present article is devoted to the theoretical and experimental study of the structure of a shock wave in a gas-liquid medium. It is shown, within the framework of a homogeneous biphasic model, that the structure of the shock wave can be studied on the basis of the Burgers-Korteweg-de Vries equation. The results of piezoelectric measurements of the pressure profile along the shock wave front agree qualitatively with the theoretical representations of the structure of the shock wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 65–69, May–June, 1973.     

Turbulence Modeling in Noninertial Frames of Reference

The effect of an arbitrary change of frame on the structure of turbulence models is examined from a theoretical standpoint. It is proven, as a rigorous consequence of the Navier-Stokes equations, that turbulence models must be form invariant under arbitrary translational accelerations of the reference frame and should only be affected by rotations through the intrinsic mean vorticity. A direct application of this invariance property along with the Taylor-Proudman theorem, material frame-indifference in the limit of two-dimensional turbulence, and Rapid Distortion Theory is shown to yield powerful constraints on the allowable form of turbulence models. Most of the commonly used turbulence models are demonstrated to be in violation of these constraints and consequently are inconsistent with the Navier-Stokes equations in noninertial frames. Alternative models with improved noninertial properties are developed and some simple applications to rotating turbulent flows are considered.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-18107 while the author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665, U.S.A. Partial support was also provided by the Office of Naval Research under Contract No. N00014-85-K-0238.     

On the smallest scale for the incompressible Navier-Stokes equations

We prove that, for solutions to the two- and three-dimensional incompressible Navier-Stokes equations, the minimum scale is inversely proportional to the square root of the Reynolds number based on the kinematic viscosity and the maximum of the velocity gradients. The bounds on the velocity gradients can be obtained for two-dimensional flows, but have to be assumed in three dimensions. Numerical results in two dimensions are given which illustrate and substantiate the features of the proof. Implications of the minimum scale result, to the decay rate of the energy spectrum are discussed.Research was supported in part by the National Acronautics and Space Administration under NASA Contract No. NAS1-18107, while the second author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. Additional support was provided by the National Science Foundation under Grant DMS-8312264 and the Office of Naval Research under Contract N-00014-83-K-0422.     

Thermal and mechanical analysis of material response to non-steady ramp and steady shock wave loading

Ramp wave experiments on the Sandia Z accelerator provide a new approach to study the rapid compression response of materials at pressures, temperatures and stress or strain rates not attainable in conventional shock experiments. Due to its shockless nature, the ramp wave experiment is often termed as an isentropic (or quasi-isentropic) compression experiment (ICE). However, in reality there is always some entropy produced when materials are subjected to large amplitude compression even under shockless loading. The entropy production mechanisms that cause deformation to deviate from the isentropic process can be attributed to mechanical and thermal dissipations. The former is due to inelasticity associated with various deformation mechanisms and the rate effect that is inherent in all the deformation processes and the latter is due to irreversible heat conduction. The main purpose of the current study is to gain insights into the effects of ramp and shock loading on the entropy production and thermomechanical responses of materials. Another purpose is to investigate the role of heat conduction in the material response to both the non-steady ramp wave and steady shock.Numerical simulations are used to address the aforementioned research objectives. The thermomechanical response associated with a steady shock wave is investigated first by solving a set of nonlinear ordinary differential equations. Using the steady wave solutions as the reference, the material responses under non-steady ramp waves are then studied with numerical wave propagation simulation. It is demonstrated that the material response to ramp and shock loading is essentially a manifestation of the interaction between the time scale associated with the loading and the intrinsic time scales associated with mechanical deformation and heat transfer. At lower loading rates as encountered in ramp loading, the loading path is closer to an isentrope and results in lower entropy production. The reasonable ramp rate to obtain a quasi-isentropic state depends on the intrinsic time scales of the dissipation mechanisms which are strongly material dependent. Thus shockless loading does not necessarily produce an isentropic response. Between two equilibrium states, heat conduction was shown to have significant effect on the temperature history but it contributes little to the overall temperature change if the specific heat remains constant. It also affects the history of entropy, but only the irreversible part of heat conduction contributes to the net entropy change. The various types of thermomechanical responses of materials would manifest themselves more significantly in terms of the thermal history than the mechanical history. Thus temperature measurement appears to be an important experimental tool in distinguishing the various mechanisms for the thermomechancial responses of the materials.     

Wall boundary conditions for high-order finite-difference schemes in computational aeroacoustics

High-order finite-difference schemes are less dispersive and dissipative but, at the same time, more isotropic than low-order schemes. They are well suited for solving computational acoustics problems. High-order finite-difference equations, however, support extraneous wave solutions which bear no resemblance to the exact solution of the original partial differential equations. These extraneous wave solutions, which invariably degrade the quality of the numerical solutions, are usually generated when solid-wall boundary conditions are imposed. A set of numerical boundary conditions simulating the presence of a solid wall for high-order finite-difference schemes using a minimum number of ghost values is proposed. The effectiveness of the numerical boundary conditions in producing quality solutions is analyzed and demonstrated by comparing the results of direct numerical simulations and exact solutions.This work was supported by the NASA Lewis Research Center Grant NAG 3-1267 and in part by the NASA Langley Research Center Grant NAG 1-1479 and the Florida State University through time granted on its Cray-YMP Supercomputer.     

Properties of shock adiabats in the neighborhood of jouguet points

Two well-known properties of shock adiabats in a gas [1] are proved for shock adiabats corresponding to discontinuous solutions of hyperbolic systems of equations expressing conservation laws. If the state on one side of a discontinuity is fixed, then at the point of extremum of the discontinuity velocity on the shock adiabat the velocity of the discontinuity is equal to one of the velocities of the characteristics on the other side of the discontinuity and vice versa. If for the systems there is defined an entropy flux or mass density of entropy, then at the points of extremum of the velocity there is an extremum of the entropy production at the discontinuity and the entropy mass density. If the system is a symmetric hyperbolic system [2, 3], then the extrema of the entropy production at the discontinuity correspond to extrema of the velocity. These properties may be helpful in the study of discontinuities in complex media, since the sections of a shock adiabat whose points can correspond to actually existing discontinuities are frequently bounded by points corresponding to discontinuities whose velocity is equal to the velocity of a characteristic on one of the sides of the discontinuity (see, for example, [1, 4, 5]).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 184–186, March–April, 1979.     

Thermodynamic fluctuations in canonical shock–turbulence interaction: effect of shock strength

In this work, we use numerical simulation and linear inviscid theory to study the thermodynamic field generated by the interaction of a shock wave with homogeneous isotropic turbulence. Fluctuations in density, pressure, temperature and entropy can play an important role in shock-induced mixing, combustion and energy transfer processes. Data from shock-captured direct numerical simulations (scDNS) are used to investigate the variation of thermodynamic fluctuations for varying shock strengths, and the results are compared with linear interaction analysis (LIA). The density, pressure and temperature variances attain large values at the shock, followed by, in general, a rapid decay in the downstream flow. The rapid variation behind the shock makes it difficult to compare numerical results with theoretical predictions. A threshold method based on instantaneous shock dilatation is used to overcome this problem, and it gives excellent match between scDNS and LIA. We find cases with non-monotonic variation with Mach number as well as local peaks in density fluctuations behind the shock. These are explained in terms of the contribution of the post-shock acoustic and entropy modes in the LIA solution and their cross-correlation. Budget of the transport equations reveals interesting insight into the physics governing the thermodynamic field behind the shock wave. It is found that the variances are primarily determined by the competing effects of dilatational and dissipation mechanisms. The dominant mechanisms are identified for a range of conditions, and their implication for developing predictive models is highlighted.     

Lagrangian description of transport equations for shock waves in three-dimensional elastic solids

A set of transport equations for the growth or decay of theamplitudes of shock waves along an arbitrary propagation directionin three-dimensional nonlinear elastic solids is derived using theLagrangian coordinates.The transport equations obtained showthat the time derivative of the amplitude of a shock wave alongany propagation ray depends on (i) an unknown quantity immediatelybehind the shock wave,(ii) the two principal curvatures of theshock surface,(iii) the gradient taken on the shock surface ofthe normal shock wave speed and (iv) the inhomogeneous term.whichis related to the motion ahead of the shock surface.vanisheswhen the motion ahead of the shock surface is uniform.Severalchoices of the propagation vector are given for which the tran-sport equations can be simplified.Some universal relations,which relate the time derivatives of various jump quantities toeach other but which do not depend on the constitutive equationsof the material,are also presented.     

Numerical simulation of shock propagation in a polydisperse bubbly liquid

The effect of distributed bubble nuclei sizes on shock propagation in a bubbly liquid is numerically investigated. An ensemble-averaged technique is employed to derive the statistically averaged conservation laws for polydisperse bubbly flows. A finite-volume method is developed to solve the continuum bubbly flow equations coupled to a single-bubble-dynamic equation that incorporates the effects of heat transfer, liquid viscosity and compressibility. The one-dimensional shock computations reveal that the distribution of equilibrium bubble sizes leads to an apparent damping of the averaged shock dynamics due to phase cancellations in oscillations of the different-sized bubbles. If the distribution is sufficiently broad, the phase cancellation effect can dominate over the single-bubble-dynamic dissipation and the averaged shock profile is smoothed out.     

Globally convergent homotopy algorithms for nonlinear systems of equations

Probability-one homotopy methods are a class of algorithms for solving nonlinear systems of equations that are accurate, robust, and converge from an arbitrary starting point almost surely. These new globally convergent homotopy techniques have been successfully applied to solve Brouwer fixed point problems, polynomial systems of equations, constrained and unconstrained optimization problems, discretizations of nonlinear two-point boundary value problems based on shooting, finite differences, collocation, and finite elements, and finite difference, collocation, and Galerkin approximations to nonlinear partial differential equations. This paper introduces, in a tutorial fashion, the theory of globally convergent homotopy algorithms, deseribes some computer algorithms and mathematical software, and presents several nontrivial engineering applications.This work was supported in part by DOE Grant DE-FG05-88ER25068, NASA Grant NAG-1-1079, and AFOSR Grant 89-0497.     

三角翼的双襟翼控涡作用的数值模拟研究

对装有“前端襟翼”和“前缘襟翼”的７４°后掠三角翼的不可压缩流场作了数值模拟,以研究襟翼的旋涡控制作用．数值模拟是用拟压缩性方法求解一般曲线坐标系下的三维不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,时间离散用向后Ｅｕｌｅｒ差分,空间无粘项的离散用二阶迎风ＴＶＤ格式,所得的离散方程用对角化形式的近似隐式因子分解格式求解．湍流模型用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模式．计算了三种平面形状的机翼在迎角范围为１０°～５０°的绕流和气动特性．计算和实验的比较表明,襟翼向下偏转可以推迟旋涡破裂,且对提高机翼的减阻能力、升阻比和改善失速前后的气动特性有明显效果,双襟翼具有更佳的控涡效果．     

Entropy analysis in electrical magnetohydrodynamic(MHD) flow of nanofluid with effects of thermal radiation,viscous dissipation,and chemical reaction

The unsteady mixed convection flow of electrical conducting nanofluid and heat transfer due to a permeable linear stretching sheet with the combined effects of an electric field, magnetic field, thermal radiation, viscous dissipation, and chemical reaction have been investigated. A similarity transformation is used to transform the constitutive equations into a system of nonlinear ordinary differential equations.The resultant system of equations is then solved numerically using implicit finite difference method.The velocity, temperature, concentration, entropy generation, and Bejan number are obtained with the dependence of different emerging parameters examined. It is noticed that the velocity is more sensible with high values of electric field and diminished with a magnetic field. The radiative heat transfer and viscous dissipation enhance the heat conduction in the system. Moreover, the impact of mixed convection parameter and Buoyancy ratio parameter on Bejan number profile has reverse effects. A chemical reaction reduced the nanoparticle concentration for higher values.