Solution of variational problems of gasdynamics of supersonic nonequilibrium flows

The problem of the optimization of the supersonic portion of a nozzle for gas flow in the case of certain nonequilibrium processes was solved in [1, 2]. The authors examined the flow scheme in which the closing Mach line of the first family arrives at the initial rarefaction wave fan. At the same time, in [3] in the solution of the analogous problem for the case of gas flow with foreign particles it was shown that it is advisable to consider also a different scheme, in which the closing characteristic arrives at the axis of symmetry outside the initial rarefaction wave fan. In the following we present results of a study of such a scheme for gas flow with nonequilibrium processes taking place. The necessary conditions which define the optimum contour are obtained and, in particular, the conditions which define the coordinate x and the magnitude of the angle at the corner points.     

Calculation of a supersonic jet emerging from a hole with planar walls

The coordinates of the sonic line are derived by means of Frankl's solution [1], while the supersonic part of the jet is considered by the method of characteristics. The numerical solution has been used to calculate the family of rarefaction waves and a family of nozzles having a corner point and a curvilinear transition surface. These calculations show that, when there is a corner point, the shape of the sonic line has hardly any effect on the velocity distribution along the symmetry axis. It is also shown that a positive pressure gradient arises on the surface of the nozzle beyond the corner point if that point lies upstream from the limiting characteristic of the first family.     

Wave structure induced by fluid-dynamic limits in the Broadwell model

Consider the fluid-dynamic limit problem for the Broadwell system of the kinetic theory of gases, for Maxwellian Riemann initial data. The formal limit is the Riemann problem for a pair of conservation laws and is invariant under dilations of coordinates. The approach of self-similar fluid-dynamic limits consists in replacing the mean free path in the Broadwell model so that the resulting problem preserves the invariance under dilations. The limiting procedure was justified in [ST]. Here, we study the structure of the emerging solutions. We show that they consist of two wave fans separated by a constant state. Each wave fan is associated with one of the characteristic fields and is either a rarefaction wave or a shock wave. The shocks satisfy the Lax shock conditions and have the internal structure of a Broadwell shock profile.     

Some geometric acoustic problems of one-dimensional unsteady and two-dimensional steady flows

The behavior of discontinuities (weak shocks) of the parameters of a disturbed flow and their interaction with the discontinuities of the basic flow in the geometric acoustics approximation, when the variation of the intensity of such shocks along the characteristics or the bicharacteristics is described by ordinary differential equations, has been investigated by many authors. Thus, Keller [1] considered the case when the undisturbed flow is three-dimensional and steady, and the external inputs do not depend on the flow parameters. An analogous study was made by Bazer and Fleischman for the MGD isentropic flow of an ideal conducting medium [2], while Lugovtsov [3] studied the three-dimensional steady flow of a gas of finite conductivity for small magnetic Reynolds numbers and no electric field. Several studies (for example, [4]) have considered the behavior of discontinuities of the solutions from the general positions of the theory of hyperbolic systems of quasilinear equations. Finally, the interaction of weak shocks (or the equivalent continuous disturbances) with shock waves was studied in [5–11].In what follows we consider one-dimensional (with plane, cylindrical, and spherical waves) and quasi-one-dimensional unsteady flows, and also plane and axisymmetric steady flows. Two problems are investigated: the variation of the intensity of weak shocks in the presence of inputs which depend on the stream parameters, and the interaction of weak shocks with strong discontinuities which differ from contact (tangential) discontinuities.The thermodynamic properties of the gas are considered arbitrary. We note that the resulting formulas for the interaction coefficients of the weak and strong discontinuities are also valid for nonequilibrium flow.     

Decay of a plane shock wave in a two-parameter medium with arbitrary equation of state

Special curves, called shock polars, are frequently used to determine the state of the gas behind an oblique shock wave from known parameters of the oncoming flow. For a perfect gas, these curves have been constructed and investigated in detail [1]. However, for the solution of problems associated with gas flow at high velocities and high temperatures it is necessary to use models of gases with complicated equations of state. It is therefore of interest to study the properties of oblique shocks in such media. In the present paper, a study is made of the form of the shock polars for two-parameter media with arbitrary equation of state, these satisfying the conditions of Cemplen's theorem. Some properties of oblique shocks in such media that are new compared with a perfect gas are established. On the basis of the obtained results, the existence of triple configurations in steady supersonic flows obtained by the decay of plane shock waves is considered. It is shown that D'yakov-unstable discontinuities decompose into an oblique shock and a centered rarefaction wave, while spontaneously radiating discontinuities decompose into two shocks or into a shock and a rarefaction wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 147–153, November–December, 1982.     

Numerical Study of Three-Dimensional Stability Problems for Laminated and Ribbon-Reinforced Composites

Classes of basic problems (classical problems and new problems derived by partial separation of variables) of the theory of elasticity and three-dimensional linearized stability are formulated in Cartesian coordinates. For each class of problems, the base scheme and the base system of equations are constructed and substantiated. They allow us to write explicitly the global mesh equations corresponding to an arbitrary problem from the class under consideration     

The laminar boundary layer near a body corner point

A method is developed for calculating the characteristics of a laminar boundary layer near a body contour corner point, in the vicinity of which the outer supersonic stream passes through a rarefaction flow. In the study we use the asymptotic solution of the Navier-Stokes equations in the region with large longitudinal gradients of the flow functions for large values of the Reynolds number, the general form of which was used in [1].The pressure, heat flux, and friction distributions along the body surface are obtained. For small pressure differentials near the corner the solution of the corresponding equations for small disturbances is obtained in analytic form.The conventional method for studying viscous gas flow near body surfaces for large values of the Reynolds number is the use of the Prandtl boundary layer theory. Far from the body the asymptotic solution of the Navier-Stokes equations in the first approximation reduces to the solution of the Euler equations, while near the body it reduces to the solution of the Prandtl boundary layer equations. The characteristic feature of the boundary layer region is the small variation of the flow functions in the longitudinal direction in comparison with their variation in the transverse direction. However, in many cases this condition is violated.The necessity arises for constructing additional asymptotic expansions for the region in which the longitudinal and transverse variations of the flow functions are quantities of the same order. The general method for constructing asymptotic solutions for such flows with the use of the known method of outer and inner expansions is presented in [1].In the following we consider the flow in a laminar boundary layer for the case of a viscous supersonic gas stream in the vicinity of a body corner point. Behind the corner the flow separates from the body surface and flows around a stagnant zone, in which the pressure differs by a specified amount from the pressure in the undisturbed flow ahead of the point of separation. A pressure (rarefaction) disturbance propagates in the subsonic portion of the boundary layer upstream for a distance which in order of magnitude is equal to several boundary layer thicknesses. In the disturbed region of the boundary layer the longitudinal and transverse pressure and velocity disturbances are quantities of the same order. In this study we construct additional asymptotic expansions in the first approximation and calculate the distributions of the pressure, friction stress, and thermal flux along the body surface.     

PML absorbing boundary conditions for the linearized and nonlinear Euler equations in the case of oblique mean flow

For the case of uniform mean flow in an arbitrary direction, perfectly matched layer (PML) absorbing boundary conditions are presented for both the linearized and nonlinear Euler equations. Although linear perfectly matched side layers with an oblique mean flow have been studied in previous works, we propose in the present paper a construction of corner layer equations that are dynamically stable. Stability issues are investigated by examining the dispersion relations of linear waves supported by the corner layer equations. For increased efficiency, a pseudo mean flow is included in the derivation of the PML equations for the nonlinear case. Numerical examples are given to support the validity of the proposed equations. Specifically, the linear PML formulation is tested for the case of acoustic, vorticity, and entropy waves traveling with an oblique mean flow. The nonlinear formulation is tested with an isentropic vortex moving diagonally with a constant velocity. Copyright © 2008 John Wiley & Sons, Ltd.     

Normal interaction of shock waves in the presence of vibrational relaxation

The effect of nonequilibrium physicochemical processes on the flow resulting from the normal collision and reflection of shock waves is studied by the example of nonequilibrium excitation of molecular oscillations in nitrogen. It is shown that the thermal effect of vibrational relaxation is small and the problem can be linearized around a known solution [1]. A similar approach to the solution of the problem of flow around a wedge and certain one-dimensional non-steady-state problems was used earlier in [2–4]. The solution of these problems was constructed in an angular domain, bounded by the shock wave and a solid wall (or the contact surface) and was reduced to a well-known functional equation [6]. The solution of this problem, because of the presence of two angular domains divided by a tangential discontinuity, reduces to a functional equation of more general form than in [6]. The results are obtained in finite form. In the special case of shocks of equal intensity, the normal reflection parameters are obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 90–96, July–August, 1976.     

Macroscopic description of steady and unsteady rarefaction effects in boundary value problems of gas dynamics

Four basic flow configurations are employed to investigate steady and unsteady rarefaction effects in monatomic ideal gas flows. Internal and external flows in planar geometry, namely, viscous slip (Kramer’s problem), thermal creep, oscillatory Couette, and pulsating Poiseuille flows are considered. A characteristic feature of the selected problems is the formation of the Knudsen boundary layers, where non-Newtonian stress and non-Fourier heat conduction exist. The linearized Navier–Stokes–Fourier and regularized 13-moment equations are utilized to analytically represent the rarefaction effects in these boundary-value problems. It is shown that the regularized 13-moment system correctly estimates the structure of Knudsen layers, compared to the linearized Boltzmann equation data.     

A dynamic theory for magnetizable elastic solids with thermal and electrical conduction

A general dynamical theory of magnetizable, electrically and thermally conducting media is developed for soft ferromagnetic or paramagnetic materials in external electromagnetic fields. The general equations are linearized by assuming infinitesimal strains, linear constitutive equations and that all field variables may be divided into two parts: a "rigid body state" and a "perturbation state". The former is the same as the one in rigid body electrodynamics, and the latter which accounts for electromagnetic interaction with the deformable continuum is coupled with stress and strain through linearized field equations. The theory is developed for general anisotropy but specialized for materials with uniaxial, or higher, symmetry.     

Stability of nonplane-parallel three-dimensional jet flows

The problem of the stability of nonplane-parallel flows is one of the most difficult and least studied problems in the theory of hydrodynamic stability [1]. In contrast to the Heisenberg approximation [1], the basic state whose stability is investigated depends on several variables, and the stability problem reduces to the solution of an eigenvalue problem for partial differential equations in which the coefficients depend on several variables [2–7]. In the case of a periodic dependence of these coefficients on the time [2] or the spatial coordinates [3, 4], the analog of Floquet theory for the partial differential equations is constructed. With rare exceptions, the case of a nonperiodic dependence has usually been considered under the assumption of weak nonplane-parallelism, i.e., a fairly small deviation from the plane-parallel case has been assumed and the corresponding asymptotic expansions in the linear [6] and nonlinear [7] stability analyses considered. The present paper considers the case of an arbitrary dependence of the velocity profile of the basic flow on two spatial variables. The deviation from the plane-parallel case is not assumed to be small, and the corresponding eigenvalue problem for the partial differential equations is solved by means of the direct methods of [5], which were introduced for the first time and justified in the theory of hydrodynamic stability by Petrov [8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 21–28, May–June, 1987.     

Corner stress singularities in an FGM thin plate

Describing the behaviors of stress singularities correctly is essential for obtaining accurate numerical solutions of complicated problems with stress singularities. This analysis derives asymptotic solutions for functionally graded material (FGM) thin plates with geometrically induced stress singularities. The classical thin plate theory is used to establish the equilibrium equations for FGM thin plates. It is assumed that the Young’s modulus varies along the thickness and Poisson’s ratio is constant. The eigenfunction expansion method is employed to the equilibrium equations in terms of displacement components for an asymptotic analysis in the vicinity of a sharp corner. The characteristic equations for determining the stress singularity order at the corner vertex and the corresponding corner functions are explicitly given for different combinations of boundary conditions along the radial edges forming the sharp corner. The non-homogeneous elasticity properties are present only in the characteristic equations corresponding to boundary conditions involving simple support. Finally, the effects of material non-homogeneity following a power law on the stress singularity orders are thoroughly examined by showing the minimum real values of the roots of the characteristic equations varying with the material properties and vertex angle.     

Free interaction of a three-dimensional boundary layer with an exterior irrotational flow

Asymptotic equations describing the unsteady free interaction of a three-dimensional boundary layer with an exterior flow are derived. The orders of the independent variables and perturbations of the flow parameters are chosen in such a way that the pressure gradient that occurs in the equation of the wall layer is due to the displacement of streamlines situated near the surface of the body. The Fourier method is used to construct a solution to the linearized problem. A class of perturbations satisfying homogeneous boundary conditions on the surface of the body is found.     

具有非线性homologous变形约束的桁架结构形态分析

研究了具有非线性homologous变形约束条件的桁架结构形态分析问题。在已有的线性homologous变形约束桁架形态分析的基础上,将结构的节点分成三类:homologous变形约束节点,形状可变节点和边界点。运用Moore-Penrose广义逆矩阵性质,将基础方程组解的存在条件表示为包含形状可变节点未知坐标的非线性方程组,为采用Newton-Raphson方法求解非线性方程组,对AA (A为任意矩阵,A 为A的Moore-Penrose广义逆矩阵)求偏导数,找到了满足保型要求的形态,给出的桁架算例说明了本文方法的有效性。     

Wave interactions and variation estimates for self-similar zero-viscosity limits in systems of conservation laws

We consider the problem of self-similar zero-viscosity limits for systems ofN conservation laws. First, we give general conditions so that the resulting boundary-value problem admits solutions. The obtained existence theory covers a large class of systems, in particular the class of symmetric hyperbolic systems. Second, we show that if the system is strictly hyperbolic and the Riemann data are sufficiently close, then the resulting family of solutions is of uniformly bounded variation and oscillation. Third, we construct solutions of the Riemann problem via self-similar zero-viscosity limits and study the structure of the emerging solution and the relation of self-similar zero-viscosity limits and shock profiles. The emerging solution consists ofN wave fans separated by constant states. Each wave fan is associated with one of the characteristic fields and consists of a rarefaction, a shock, or an alternating sequence of shocks and rarefactions so that each shock adjacent to a rarefaction on one side is a contact discontinuity on that side. At shocks, the solutions of the self-similar zero-viscosity problem have the internal structure of a traveling wave.     

Recursive-operator method in vibration problems for rod systems

Using linear differential equations with constant coefficients describing one-dimensional dynamical processes as an example, we show that the solutions of these equations and systems are related to the solution of the corresponding numerical recursion relations and one does not have to compute the roots of the corresponding characteristic equations. The arbitrary functions occurring in the general solution of the homogeneous equations are determined by the initial and boundary conditions or are chosen from various classes of analytic functions. The solutions of the inhomogeneous equations are constructed in the form of integro-differential series acting on the right-hand side of the equation, and the coefficients of the series are determined from the same recursion relations. The convergence of formal solutions as series of a more general recursive-operator construction was proved in [1]. In the special case where the solutions of the equation can be represented in separated variables, the power series can be effectively summed, i.e., expressed in terms of elementary functions, and coincide with the known solutions. In this case, to determine the natural vibration frequencies, one obtains algebraic rather than transcendental equations, which permits exactly determining the imaginary and complex roots of these equations without using the graphic method [2, pp. 448–449]. The correctness of the obtained formulas (differentiation formulas, explicit expressions for the series coefficients, etc.) can be verified directly by appropriate substitutions; therefore, we do not prove them here.     

Gas Flows with Several Thermal Nonequilibrium Modes

We study gas flows with any finite number of thermal nonequilibrium modes. The equations describing such flows are a hyperbolic system with several relaxation equations. An important feature is entropy increase dictated by physics for any irreversible process. Under physical assumptions we obtain properties of thermodynamic variables relevant to stability. By the energy method we prove global existence and uniqueness for the Cauchy problem when the initial data are small perturbations of constant equilibrium states. We give a precise formulation of the fundamental solution for the linearized system, and use it to obtain large time behavior of solutions to the nonlinear system. In particular, we show that the entropy increases but stays bounded. The resulting asymptotic picture of nonequilibrium flows is in a pointwise sense both in space and in time.