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.     

Stochastic-Perturbation Analysis of a One-Dimensional Dispersion-Reaction Equation: Effects of Spatially-Varying Reaction Rates

We carry out a stochastic-perturbation analysis of a one-dimensional convection–dispersion-reaction equation for reversible first-order reactions. The Damköhler number, Da, is distributed randomly from a distribution that has an exponentially decaying correlation function, controlled by a correlation length, . Zeroth- and first-order approximations of the dispersion coefficient, D are computed from moments of the residence-time distribution obtained by solving a one-dimensional network model, in which each unit of the network represents a Darcy-level transport unit, and the solution of the transfer function in zeroth- and first-order approximations of the transport equation. In the zeroth-order approximation, the dispersion coefficient is calculated using the convection–dispersion-reaction equation with constant parameters, that is, perturbation corrections to the local equation are ignored. This zeroth-order dispersion coefficient is a linear function of the variance of the Damköhler number, (Da)². A similar result was reported in a two-dimensional network simulation. The zeroth-order approximation does not give accurate predictions of mixing or spreading of a plume when Damköhler numbers, Da 1 and its variance, (Da)² > 0.25 Da². On the other hand, the first-order theory leads to a dispersion coefficient that is independent of the reaction parameters and to equations that do accurately predict mixing and spreading for Damköhler numbers and variances in the range (Da)²/Da0.3     

Nonstationary filtration in partially saturated porous media

From the mathematical formulation of a one-dimensional flow through a partially saturated porous medium, we arrive at a nonlinear free boundary problem, the boundary being between the saturated and the unsaturated regions in the medium. In particular we obtain an equation which is parabolic in the unsaturated part of the domain and elliptic in the saturated part.Existence, uniqueness, a maximum principle and regularity properties are proved for weak solutions of a Cauchy-Dirichlet problem in the cylinder {(x,t): 0x1, t0} and the nature, in particular the regularity, of the free boundary is discussed.Finally, it is shown that solutions of a large class of Cauchy-Dirichlet problems converge towards a stationary solution as t and estimates are given for the rate of convergence.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

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.     

Nonsimilar Solutions for Mixed Convection in Non‐Newtonian Fluids Along a Vertical Plate in a Porous Medium

A nonsimilar boundary layer analysis is presented for the problem of mixed convection in powerlaw type nonNewtonian fluids along a vertical plate with powerlaw wall temperature distribution. The mixed convection regime is divided into two regions, namely,the forced convection dominated regime and the free convection dominated regime. The two solutions are matched. Numerical results are presented for the details of the velocity and temperature fields. A discussion is provided for the effect of viscosity index on the surface heat transfer rate.     

Electrode layers in flows of an inviscid plasma with good conductivity

The structure of the electromagnetic electrode layers that are produced in flows across a magnetic field by a completely ionized and inviscid plasma with good conductivity and a high magnetic Reynolds number is examined in a linear approximation. Flow past a corrugated wall and flow in a plane channel of slowly varying cross section with segmented electrodes are taken as specific examples. The possibility is demonstrated of the formation of nondissipative electrode layers with thicknesses on the order of the Debye distance or electron Larmor radius and of dissipative layers with thicknesses on the order of the skin thickness, as calculated from the diffusion rate in a magnetic field [2].In plasma flow in a transverse magnetic field, near the walls, along with the gasdynamie boundary layers, which owe their formation to viscosity, thermal conductivity, etc. (because of the presence of electromagnetic fields, their structures may vary considerably from that of ordinary gasdynamic layers), proper electromagnetic boundary layers may also be produced. An example of such layers is the Debye layer in which the quasi-neutrality of the plasma is upset. No less important, in a number of cases, is the quasi-neutral electromagnetic boundary layer, in which there is an abrupt change in the frozen-in parameter k=B/p (B is the magnetic field and p is the density of the medium). This layer plays a special role when we must explicitly allow for the Hall effect and the related formation of a longitudinal electric field (in the direction of the veloeiryv of the medium). We will call this the magnetic layer. The magnetic boundary layer can be dissipative as well as noudissipative (see below). The dissipative magnetic layer has been examined in a number of papers: for an incompressible medium with a given motion law in [1], for a compressible medium with good conductivity in [2], and with poor conductivity in [3]. In the present paper, particular attention will be devoted to nondissipative magnetic boundary layers.     

Some problems of nonisothermal steady flow of a viscous fluid

The paper presents solutions to the problems of plane Couette flow, axial flow in an annulus between two infinite cylinders, and flow between two rotating cylinders. Taking into account energy dissipation and the temperature dependence of viscosity, as given by Reynolds's relation =₀ exp (–T) (₀, =const). Two types of boundary conditions are considered: a) the two surfaces are held at constant (but in general not equal) temperatures; b) one surface is held at a constant temperature, the other surface is insulated.Nonisothermal steady flow in simple conduits with dissipation of energy and temperature-dependent viscosity has been studied by several authors [1–11]. In most of these papers [1–6] viscosity was assumed to be a hyperbolic function of temperature, viz. =_m 1/1+²(T–T_m.Under this assumption the energy equation is linear in temperature and can he easily integrated. Couette flow with an exponential viscosity-temperature relation. =₀ e ^–T (₀, =const), (0.1) was studied in [7, 8]. Couette flow with a general (T) relation was studied in (9).Forced flow in a plane conduit and in a circular tube with a general (T) relation was studied in [10]. In particular, it has been shown in [10] that in the case of sufficiently strong dependence of viscosity on temperature there can exist a critical value of the pressure gradient, such that a steady flow is possible only for pressure gradients below this critical value.In a previous work [11] the authors studied Polseuille flow in a circular tube with an exponential (T) relation. This thermohydrodynamic problem was reduced to the problem of a thermal explosion in a cylindrical domain, which led to the existence of a critical regime. The critical conditions for the hydrodynamic thermal explosion and the temperature and velocity profiles were calculated.In this paper we treat the problems of Couette flow, pressureless axial flow in an annulus, and flow between two rotating cylinders taking into account dissipation and the variation of viscosity with temperature according to Reynolds's law (0.1). The treatment of the Couette flow problem differs from that given in [8] in that the constants of integration are found by elementary methods, whereas in [8] this step involved considerable difficulties. The solution to the two other problems is then based on the Couette problem.     

Asymptotic solution to the problem of hypersonic flow over bodies in the neighborhood of the separation point of a thin shock layer

A study is made of the problem of hypersonic flow of an inviscid perfect gas over a convex body with continuously varying curvature. The solution is sought in the framework of the asymptotic theory of a strongly compressed gas [1–4] in the limit M when the specific heat ratio tends to 1. Under these assumptions, the disturbed flow is situated in a thin shock layer between the body and the shock wave. At the point where the pressure found by the Newton-Buseman formula vanishes there is separation of the flow and formation of a free layer next to the shock wave [1–4]. The singularity of the asymptotic expansions with respect to the parameter ₁ = ( –1)/( + 1) associated with separation of the strongly compressed layer has been investigated previously by various methods [3–9]. Local solutions to the problem valid in the neighborhood of the singularity have been obtained for some simple bodies [3–7]. Other solutions [7, 9] eliminate the singularity but do not give the transition solution entirely. In the present paper, an asymptotic solution describing the transition from the attached to the free layer is constructed for a fairly large class of flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 99–105, January–February, 1982.     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

Almost analytic solutions and their tests of the horizontal diffusion equation for the movement of water in unsaturated soil

Itisalwaysdifficulttofindthesolutionsoftheequationforthemovementofwaterinunsaturatedsoi1.Theprimar}'reasonisthatthehydraulicconductivityK(T)orthediffusivityofsoiIwaterD(o)isfunctionofwaterpotential(W)orwatercontent'(o)'Atpresent,thegeneralwaystofindthesol…     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

Some problems concerning nonstationary currents in dense plasma systems confined by walls

Nonstationary currents are examined in a dense magnetized plasma with 1, in which energy release and heat loss by thermal conduction and radiation are possible. Solutions are found in two limiting cases: ¦f¦ ¦ div (T)¦ and ¦f¦ ¦ div(T)¦ (f is the radiation intensity, is the coefficient of heat conduction, and T is the temperature). In the first case a solution was obtained of some problems of the cooling and heating of a plasma illustrated in part by the evolution in time of the temperature profile in the boundary layer. In the second case an isomorphic solution was found for an arbitrary dependence of the coefficient of heat conduction on the temperature, pressure, and magnetic field.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 3–8, January–February, 1972.The author is grateful to G. I. Budker for formulating the problem.     

Various reciprocal theorems and variational principles in the theories of nonlocal micropolar linear elastic mediums

In the first part of our paper, we have extended the concepts of the classical convolution and the convolution scalar product given by I. Hlaváck and presented the concepts of the convolution vector and the convolution vector scalar product, which enable us to extend the initial value as well as the initial-boundary value problems for the equation with the operator coefficients to those for the system of equations with the operator coefficients.In the second part of this paper, based on the concepts of the convolution vector and the convolution vector scalar product, two fundamental types of reciprocal theorems of the non-local micropolar linear elastodynamics for inhomogeneous and anisotropic solids are derived.In the third part of this paper, based on the concepts and results in the first and second parts as well as the Lagrange multiplies method which is presented by W. Z. Chien, four main types of variational principles are given for the nonlocal micropolar linear elastodynamics for inhomogeneous and anisotropic solids. These are the counterparts of the variational principles of Hu-Washizu type, Hellinger-Reissner type and Gurtin type in classical elasticity as well as Hlaváck type and Iesan type in local micropolar and nonlocal elasticity. Finally, we have proved the equivalence of the last two main variational principles which are given in this paper.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

ON THE BOUNDEDNESS AND THE STABILITY RESULTS FOR THE SOLUTION OF CERTAIN FOURTH ORDER DIFFERENTIAL EQUATIONS VIA THE INTRINSIC METHOD

ＯＮＴＨＥＢＯＵＮＤＥＤＮＥＳＳＡＮＤＴＨＥＳＴＡＢＩＬＩＴＹＲＥＳＵＬＴＳＦＯＲＴＨＥＳＯＬＵＴＩＯＮＯＦＣＥＲＴＡＩＮＦＯＵＲＴＨＯＲＤＥＲＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＶＩＡＴＨＥＩＮＴＲＩＮＳＩＣＭＥＴＨＯＤＣｅｍｉｌＴＵＮＣ;Ａ...     

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.     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.