Interior Differentiability of Weak Solutions to the Equations of Stationary Motion of a Class of Non-Newtonian Fluids

In this paper, we consider weak solutions to the equations of stationary motion of a class of non-Newtonian fluids the constitutive law of which includes the power law model as special case. We prove the existence of second order derivatives of weak solutions to these equations.     

On Regularity of Stationary Stokes and Navier-Stokes Equations near Boundary

We obtain local estimates of the steady-state Stokes system without pressure near boundary. As an application of the local estimates, we prove the partial regularity up to the boundary for the stationary Navier-Stokes equations in a smooth domain in five dimension.     

About the Regularized Navier–Stokes Equations

The first goal of this paper is to study the large time behavior of solutions to the Cauchy problem for the 3-dimensional incompressible Navier–Stokes system. The Marcinkiewicz space L^3, is used to prove some asymptotic stability results for solutions with infinite energy. Next, this approach is applied to the analysis of two classical regularized Navier–Stokes systems. The first one was introduced by J. Leray and consists in mollifying the nonlinearity. The second one was proposed by J.-L. Lions, who added the artificial hyper-viscosity (–)^{/ 2}, > 2 to the model. It is shown in the present paper that, in the whole space, solutions to those modified models converge as t toward solutions of the original Navier–Stokes system.     

Spatial Simple Waves on a Shear Flow

The paper studies simple waves of the shallowwater equations describing threedimensional wave motions of a rotational liquid in a freeboundary layer. Simple wave equations are derived for the general case. The existence of unsteady or steady simple waves adjacent continuously to a given steady shear flow along a characteristic surface is proved. Exact solutions of the equations describing steady simple waves were found. These solutions can be treated as extension of Prandtl–Mayer waves for sheared flows. For shearless flows, a general solution of the system of equations describing unsteady spatial simple waves was found.     

Investigation of the flow past wings of infinite span with a blunt leading edge in the vorticity interaction regime

An asymptotic analysis of the Navier-Stokes equations is carried out for the case of hypersonic flow past wings of infinite span with a blunt leading edge when 0, Re , and M . Analytic solutions are obtained for an inviscid shock layer and inviscid boundary layer. The results of a numerical solution of the problems of vorticity interaction at the blunt edge and on the lateral surface of the wing are presented. These solutions are compared with the solution of the equations of a thin viscous shock layer and on the basis of this comparison the boundaries of the asymptotic regions are estimated.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–127, November–December, 1987.     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

Transformation and some solutions of the equations of motion of an ideal gas

The equations of one-dimensional and plane steady adiabatic motion of an ideal gas are transformed to a new form in which the role of the independent variables are played by the stream function and the function introduced by Martin [1, 2], It is shown that the function retains a constant value on a strong shock wave (and on a strong shock for plane flows). For one-dimensional isentropic motions the resulting transformation permits new exact solutions to be obtained from the exact solutions of the equations of motion. It is shown also that the one-dimensional motions of an ideal gas with the equation of state p=f(t) and the one-dimensional adiabatic motions of a gas for which p=f() are equivalent (t is time, is the stream function). It is shown that if k=s=–1, m and n are arbitrary (m+n0) and =1, the general solution of the system of equations which is fundamental in the theory of one-dimensional adiabatic self-similar motions [3] is found in parametric form with the aid of quadratures. Plane adiabatic motions of an ideal gas having the property that the pressure depends only on a single geometric coordinate are studied.     

Asymptotic Modeling of Nonlinear Wave Processes in Shock‐Loaded Elastoplastic Materials

Nonlinear wave processes in shockloaded elastoplastic materials are modeled. A comparison of the results obtained with experimental data and numerical solutions of exact systems of dynamic equations shows that the model equations proposed qualitatively describe the stressdistribution evolution in both the elasticflow and plasticflow regions and can be used to solve one and twodimensional problems of pulsed deformation and fracture of elastoplastic media.     

Exact solutions and mathematical properties of boundary‐value problems for dynamic‐diffusion boundary layers

The paper studies boundaryvalue problems for dynamicdiffusion boundary layers occurring near a vertical wall at high Schmidt numbers and for dynamic boundary layers whose inner edge is adjacent to the dynamicdiffusion layers. Exact solutions for boundary layers at small and large times are derived. The wellposedness of the boundaryvalue problem for a steady dynamicdiffusion layer is studied.     

Numerical investigation of laminar separation in the case of supersonic flow of viscous gas past spiked bodies

The complete Navier-Stokes equations for a compressible viscous perfect heat conducting gas have been used in a numerical investigation of laminar separation in the case of supersymmetric axisymmetric flow past cylinders with a conical nose and a spike at the front of finite thickness. The flow structure has been studied in its dependence on the length of the spike and the half-angle of the conical tip. For the considered free-stream parameters (2 M 6, 100 Re 500) and spike lengths, which do not exceed the diameter of the cylinder, the existence of steady flow regimes has been established and it has been shown that the spike in front of the body reduces its total drag and the heat flux to its surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–131, March–April, 1984.     

Equilibrium Configurations for a Floating Drop

If a drop of fluid of density ₁ rests on the surface of a fluid of density ₂ below a fluid of density ₀, ₀ < ₁ < ₂, the surface of the drop is made up of a sessile drop and an inverted sessile drop which match an external capillary surface. Solutions of this problem are constructed by matching solutions of the axisymmetric capillary surface equation. For general values of the surface tensions at the common boundaries of the three fluids the surfaces need not be graphs and the profiles of these axisymmetric surfaces are parametrized by their tangent angles. The solutions are obtained by finding the value of the tangent angle for which the three surfaces match. In addition the asymptotic form of the solution is found for small drops.     

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

On the Quasi‐Steadiness Hypothesis as Applied to Gas Exhaustion from a Receiver

A semiempirical method of determining the stabilization time for a quasisteady mode of gas exhaustion from a receiver after sudden opening of the nozzle and the time evolution of the real flow rate at the stage of the transitional process are considered. The numerical solution of the equations of exhaustion gas dynamics in a twodimensional formulation and the results of model experiments demonstrated that the method can be used to estimate the conditions of applicability of the quasisteadiness hypothesis and to determine the discharge coefficient of the nozzle with controlled accuracy.     

Unsteady forced convection heat/mass transfer from a flat plate

Numerical methods are used to investigate the transient, forced convection heat/mass transfer from a finite flat plate to a steady stream of viscous, incompressible fluid. The temperature/concentration inside the plate is considered uniform. The heat/mass balance equations were solved in elliptic cylindrical coordinates by a finite difference implicit ADI method. These solutions span the parameter ranges 10 Re 400 and 0.1 Pr 10. The computations were focused on the influence of the product (aspect ratio) × (volume heat capacity ratio/Henry number) on the heat/mass transfer rate. The occurrence on the plates surface of heat/mass wake phenomena was also studied.     

Numerical Analysis of Branched Shapes of Arches in Bending

Nonlinear boundaryvalue problems of plane bending of elastic arches under a uniformly distributed load are solved by the shooting method. The problems are formulated for a system of six firstorder ordinary differential equations with a finiterotation field independent of displacements. Simply supported and clamped cases are considered. Branching solutions of the boundaryvalue problems are obtained. For a simply supported arch, a set of solutions describes symmetric and nonsymmetric shapes of bending, which correspond to positive, negative, and zero loads. For a clamped arch, the set of solutions consists of symmetric shapes that occur only for positive loads.     

Group Classification of Equations of the Form y″ = f(x,y)

The problem of classification of ordinary differential equations of the form y = f(x,y) by admissible local Lie groups of transformations is solved. Standard equations are listed on the basis of the equivalence concept. The classes of equations admitting a oneparameter group and obtained from the standard equations by invariant extension are described.     

Flow stability of a Grad-model liquid layer on an inclined plane

A study is presented of the flow of stability of a Grad-model liquid layer [1, 2] flowing over an inclined plane under the influence of the gravity force.It is assumed that at every point of the considered material continuum, along with the conventional velocity vector v, there is defined an angular velocity vector , the internal moment stresses are negligibly small, and in the general case the force stress tensor _kj is asymmetric. The model is characterized by the usual Newtonian viscosity , the Newtonian rolling viscosity _r, and the relaxation time = J/4 _r, where J is a scalar constant of the medium with dimensions of moment of inertia per unit mass, is the density. It is assumed that the medium is incompressible, the coefficients , _r, J are constant [2].The exact solution of the equations of motion, corresponding to flow of a layer with a plane surface, coincides with the solution of the Navier-Stokes equations in the case of flow of a layer of Newtonian fluid. The equations for three-dimensional periodic disturbances differ considerably from the corresponding equations for the problem of the flow stability of a layer of a Newtonian medium. It is shown that the Squire theorem is valid for parallel flows of a Grad liquid.The flow stability of the layer with respect to long-wave disturbances is studied using the method of sequential approximations suggested in [3, 4].     

Numerical analysis of axisymmetric buckling of conical shells

Nonlinear boundaryvalue problems of axisymmetric buckling of conical shells under a uniformly distributed normal pressure are solved by the shooting method. The problems are formulated for a system of six firstorder ordinary differential equations with independent rotation and displacement fields. Simply supported and clamped cases are considered. Branching solutions of the boundaryvalue problems are studied for different pressures and geometrical parameters of the shells. The nonmonotonic and discontinuous curves of equilibrium states obtained show that collapse, i.e., snapthrough instability is possible. For a simply supported shell, multivalued solutions are obtained for both external and internal pressure. For a clamped thinwalled shell, theoretical results are compared with experimental data.     

Direct simulation of fluid particle motions

Continuum models of two-phase flows of solids and liquids use constitutive assumptions to close the equations. A more fundamental approach is a molecular dynamic simulation of flowing big particles based on reliable macroscopic equations for both solid and liquid. We developed a package that simulates the unsteady two-dimensional solid-liquid two-phase flows using the Navier-Stokes equations for the liquid and Newton's equations of motion for the solid particles. The Navier-Stokes equations are solved using a finite-element formulation and Newton's equations of motion are solved using an explicit-implicit scheme. We show that the simplest fully explicit scheme to update the particle motion using Newton's equations is unstable. To correct this instability we propose and implement and Explicit-Implicit Scheme in which, at each time step, the positions of the particles are updated explicitly, the computational domain is remeshed, the solution at the previous time is mapped onto the new mesh, and finally the nonlinear Navier-Stokes equation and the implicitly discretized Newton's equations for particle velocities are solved on the new mesh iteratively. The numerical simulation reveals the effect of vortex shedding on the motion of the cylinders and reproduces the drafting, kissing, and tumbling scenario which is the dominant rearrangement mechanism in two-phase flow of solids and liquids in beds of spheres which are constrained to move in only two dimensions.This work was supported by the National Science Foundation, the Department of Energy, and the Army Research Office, Mathematics.