Homogenization of the Navier-Stokes equations in open sets perforated with tiny holes I. Abstract framework,a volume distribution of holes

This paper treats the homogenization of the Stokes or Navier-Stokes equations with a Dirichlet boundary condition in a domain containing many tiny solid obstacles, periodically distributed in each direction of the axes. (For example, in the three-dimensional case, the obstacles have a size of ³ and are located at the nodes of a regular mesh of size .) A suitable extension of the pressure is used to prove the convergence of the homogenization process to a Brinkman-type law (in which a linear zero-order term for the velocity is added to a Stokes or Navier-Stokes equation).     

Scheme of Investigation of the Spectrum of a Family of Perturbed Operators and Its Application to Spectral Problems in Thick Junctions

We develop a scheme for the investigation of the asymptotic behavior of eigenvalues and eigenvectors of a family of self-adjoint compact operators {A: > 0} that act in different spaces and lose their compactness in the limit case 0. We prove the Hausdorff convergence of the spectrum of the operator A to the spectrum of the limit operator A₀, obtain asymptotic estimates for this convergence both to points of the discrete spectrum and to points of the essential spectrum of the operator A₀, and prove asymptotic estimates for eigenvectors of A. This scheme is applied to the investigation of the asymptotic behavior of eigenvalues and eigenfunctions of the Neumann problem in a thick singularly degenerate junction that consists of two domains connected by an -periodic system of thin rods of fixed length.     

Homogenization of the Process of Phase Transitions in Multidimensional Heterogeneous Periodic Media

The homogenization of the Stefan multidimensional problem is carried out in the case where the medium is a composite consisting of two different substances with an –periodic structure. The averaged problem is deduced by asymptotic methods. It is shown that its solution is the limit of solutions of –problems.     

A solution of the E-ε model for nearly homogeneous turbulence with a mean shear

An analytical solution of the E- model for the downstream evolution of a stationary and nearly homogeneous turbulent shear flow is presented. In case that the turbulent time scale has adjusted itself to the time scale imposed by the shear, an asymptotic solution can be derived from the full solution, which shows that both E and increase downstream exponentially. By comparing this asymptotic solution with experimental data a value for the unknown constant c _l , in the -equation, is derived. Moreover, we find an expression for the downstream development of the variance of a scalar, which is also compared with experimental data. The analytical solution shows that a homogeneous shear flow with a uniform velocity gradient can only be obtained if the shear is sufficiently small. In the experiments this condition is not always satisfied. A discussion is given on how a nearly homogeneous shear flow can be obtained over a limited downstream interval by changing the initial conditions in E and , and a comparison is made with experimental data. Finally it is shown that better transverse homogeneity can be obtained by taking an exponential velocity profile instead of a linear profile.     

Homogenization of the Robin problem in a thick multilevel junction

We consider a mixed boundary-value problem for the Poisson equation in a plane two-level junction that is the union of a domain ₀ and a large number 2N of thin rods with variable thickness of order = (N ^–1). The thin rods are divided into two levels, depending on their length. In addition, the thin rods from each level are -periodically alternated. We investigate the asymptotic behavior of the solution as 0 under the Robin conditions on the boundaries of the thin rods. By using some special extension operators, a convergence theorem is proved.Published in Neliniini Kolyvannya, Vol. 7, No. 3, pp. 336–355, July–September, 2004.     

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.     

Spin-up in a saturated porous medium

It is shown that, on the Brinkman model, spin-up is confined to boundary layers whose thickness is of order k ^1/2, and the spin-up is established in a time of order k/, where k, , and denote permeability, density, porosity and dynamic viscosity, respectively.     

On the Collapse of Masonry Arches

In this paper a constitutive equation for masonry arches is defined and its main properties are proven; in this equation to each pair of generalized strains (, ), with the extensional strain and the curvature change of the centre line, is assigned the pair of generalized internal forces (N,M), where N is the normal force and M the bending moment. Subsequently, the collapse of masonry arches is characterized and the static and kinematic theorems proven. Finally, a method for determining the collapse load in the case of circular arches subjected to their own weight and a vertical point load applied at a point of the extrados is presented. The results obtained, of interest in some applications, are summarized in a series of graphs.     

Collision of two plane free streams with division of one of them into two flows

The collision of two streams of ideal incompressible fluid with the division of one stream into two flows in opposite directions along the boundary of the other is considered. The second stream is shifted from its initial direction due to the collision. The ratio of the total pressure heads of the separated and shifted streams is 1. The flow regions corresponding to each of these streams map conformally onto the first and fourth quadrants of the complex plane so that the streamline separating the streams corresponds to the real semiaxis in both cases. By using the dynamic and kinematic conditions for this streamline, nonlinear relations are found linking the boundary values of the logarithms of the complex velocity in each of the two flow regions. These analytic functions are computed approximately by linearization at small values of .Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 3–7, January–February, 1975.     

The influence of blending polystyrenes of narrow molecular weight distribution on melt creep flow and creep recovery in elongation

Creep and creep recovery experiments in elongation were performed with melts of anionically polymerized polystyrenes (PS) and with their blends at a temperature of 150 °C. For stresses ₀ < 10 000 N/m² the samples with narrow molecular weight distribution show linear viscoelastic behavior up to the maximum Hencky strain = 3.5, achievable in a newly developed elongational rheometer for polymer melts. The compliances,D (t), of the blends are linear-viscoelastic only up to a strain limit _L . For strains beyond _L the compliance of each blend depends on the stress ₀. For a series of binary blends, prepared from the same components of narrow MWD, the linear-viscoelastic limit _L seems to be independent of the mixing ratio and stress. _L seems to be a function only of the molecular weights of the original components, the blends investigated were made from.Paper presented at the Annual Conference of the German Society of Rheology at Ulm, March 7–10, 1983.     

Lower semicontinuity of the attractor for a singularly perturbed hyperbolic equation

For a smooth, bounded domain R, n 3, and a real, positive parameter, we consider the hyperbolic equationu _tt +u _t –u=–f(u) –g in with Dirichlet boundary conditions. Under certain conditions onf, this equation has a global attractorA inH ₀ ¹ () ×L ²(). For=0, the parabolic equation also has a global attractor which can be naturally embedded into a compact setA ₀ inH ₀ ¹ () ×L ²(). If all of the equilibrium points of the parabolic equation are hyperbolic, it is shown that the setsA are lower semicontinuous at=0. Moreover, we give an estimate of the symmetric distance betweenA ₀ andA .     

Resonance Capture in a Three Degree-of-Freedom Mechanical System

We study the phenomena of resonance capture in a three degree-of-freedom dynamical system modeling the dynamics of an unbalanced rotor, subject to a small constant torque, supported by orthogonal, linearly elastic supports, which is constrained to move in the plane. In the physical system the resonance exists between translational motions of the frame and the angular velocity of the unbalanced rotor. These equations, valid in the neighborhood of the resonance, possess a small parameter which is related to the imbalance. In the limit 0, the unperturbed system possesses a homoclinic orbit which separates bounded periodic motion corresponding to resonant solutions from unbounded motion which corresponds to solutions passing through the resonance. Using a generalized Melnikov integral, we characterize the splitting distance between the invariant manifolds which govern capture and escape from resonance for 0. It is shown that as certain slowly varying parameters evolve, the separation distance alternates sign, indicating that both capture into, and escape from resonance occur. We find that although a measurable set of initial conditions enter into a sustained resonance, as the system further evolves the orientation of the manifolds reverses and many of these captured solutions will subsequently escape.     

Square and Pulse Waves with Two Delays

This work is concerned with the effects on the dynamics of a differential difference equation with two delays as the delays become unbounded in a fixed direction. This leads to a singularly perturbed delay differential equation with singular parameter and delays (1, d). We study in detail d=2 for the case when =0 yields the Hénon map. In a neighborhood of a generic period doubling point for the Hénon map, we show that there can be either a stable square wave or an unstable pulse wave even though the period two point for the map is always stable.     

Transcendentally small transversality in the rapidly forced pendulum

The rapidly forced pendulum equation with forcing sin((t/), where =<₀^p,p = 5, for ₀, sufficiently small, is considered. We prove that stable and unstable manifolds split and that the splitting distanced(t) in the ( ,t) plane satisfiesd(t) = sin(t/) sech(/2) +O( ₀ exp(–/2)) (2.3a) and the angle of transversal intersection,, in thet = 0 section satisfies 2 tan/2 = 2S _s = (/2) sech(/2) +O(( ₀ /) exp(–/2)) (2.3b) It follows that the Melnikov term correctly predicts the exponentially small splitting and angle of transversality. Our method improves a previous result of Holmes, Marsden, and Scheuerle. Our proof is elementary and self-contained, includes a stable manifold theorem, and emphasizes the phase space geometry.     

Quasi-Stability of the Primary Flow in a Cone and Plate Viscometer

We investigate the flow between a shallow rotating cone and a stationary plate. This cone and plate device is used in rheometry, haemostasis as well as in food industry to study the properties of the flow w.r.t. shear stress. Physical experiments and formal computations show that close to the apex the flow is approximately azimuthal and the shear-stress is constant within the device, the quality of the approximation being controlled essentially by the single parameter Re ², where Re is the Reynolds number and the thinness of the cone-plate gap. We establish this fact by means of rigorous energy estimates and numerical simulations. Surprisingly enough, this approximation is valid though the primary flow is not itself a solution of the Navier-Stokes equations, and it does not even fulfill the correct boundary conditions, which are in this particular case discontinuous along a line, thus not allowing for a usual Leray solution. To overcome this difficulty we construct a suitable corrector.     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

The dynamics of resonant capture

Resonant capture describes the behavior of a weakly coupled multi-degree-of-freedom system when two or more of its uncoupled frequencies become locked in resonance. Flow on the region of phase space near the resonance (the resonance manifold) involves a region bounded by a separatrix in the uncoupled (=0) system. Capture corresponds to motions which appear to cross into the interior of the separated region for >0. We offer two approximate methods for estimating which initial conditions lead to capture: an energy method and a perturbation method based on invariant manifold theory. These methods are applied to a model problem involving the spinup of an unbalanced rotor attached to an elastic support.     

The eddy viscosity of turbulent equilibrium layers

Summary Similarity laws for the mean flow and scaling laws for the turbulent motion are used in an attempt to obtain a general expression for the eddy viscosity of equilibrium layers. It is found that =0.09 w ² /w*, in which w ² is a Reynolds stress representative for the region of overlap between the law of the wall and the velocity-defect law, while w* is the logarithmic slope of the mean velocity profile in that region. The distinction between w and w* is related to the strong inhomogeneity of the mean rate of strain in the inner layer. The results of the theory agree with experimental evidence obtained from transpired equilibrium layers.     

Use of turbulent energy balance equation in the theory of turbulent flows along walls

The turbulent flow of an incompressible fluid is considered in a plane channel, a circular tube, and the boundary layer on a flat plate. The system of equations describing the motion of the fluid consists of the Reynolds equations and the mean kinetic energy balance equation for turbulent fluctuations. On the basis of an analysis of experimental data, hypotheses are formulated with respect to the eddy kinematic viscosity and lengthl entering into the expression for specific dissipation of turbulent energy into heat. It is assumed that in the central (outer) region of the flow in a channel, andl are constants, and expressions are taken for them which are used for a free boundary layer; near the walll varies linearly and almost linearly. Results of calculations of the turbulent energy distribution, the mean velocity, and the drag coefficient are in good agreement with the existing experimental data. The values of two empirical coefficients, which enter into the system of equations as the result of the hypotheses, are close to those obtained for a free boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 25–33, May–June, 1973.     

On stability of the Hill equation with damping

We consider the Hill equation with damping describing the parametric oscillations of a torsional pendulum excited by varying the moment of inertia of the rotating body. Using the method of a small parameter, we analytically calculate a fundamental system of solutions of this equation in the form of power series in the excitation amplitude with accuracy O(²) and verify conditions for its stability. In the first-order approximation in , we prove that the resonance domain exists only if the excitation frequency is sufficiently close to the double natural frequency of the pendulum; the corresponding equation of the stability boundary is obtained.Published in Neliniini Kolyvannya, Vol. 7, No. 2, pp. 169–179, April–June, 2004.