Incompressible Fluids in Thin Domains with Navier Friction Boundary Conditions (I)

This study is motivated by problems arising in oceanic dynamics. Our focus is the Navier–Stokes equations in a three-dimensional domain Ω_ɛ, whose thickness is of order O(ɛ) as ɛ → 0, having non-trivial topography. The velocity field is subject to the Navier friction boundary conditions on the bottom and top boundaries of Ω_ɛ, and to the periodicity condition on its sides. Assume that the friction coefficients are of order O(ɛ^3/4) as ɛ → 0. It is shown that if the initial data, respectively, the body force, belongs to a large set of H¹(Ω_ɛ), respectively, L²(Ω_ɛ), then the strong solution of the Navier–Stokes equations exists for all time. Our proofs rely on the study of the dependence of the Stokes operator on ɛ, and the non-linear estimate in which the contributions of the boundary integrals are non-trivial.     

Influence of the material hardening and compressibility on the solution of elastoplastic deformation problems for a space with a cylindrical cavity

The present paper deals with plane deformation problems (ɛ _z = 0) concerned with elastoplastic deformation of a space with a cylindrical cavity in the case where the load is given either at infinity or on the cavity surface. It is assumed that the material obeys the relations of the theory of flow with isotropic hardening and the von Mises plasticity condition. The effects of the elastic compressibility (Poisson’s ratio) and the coefficient of linear hardening on the stress-strain state are studied. The influence of the linear hardening is shown to be small, while that of the elastic compressibility is shown to be quite significant.     

Reversible Bifurcation of Homoclinic Solutions in Presence of an Essential Spectrum

We consider bifurcations of a class of infinite dimensional reversible dynamical systems which possess a family of symmetric equilibria near the origin. We also assume that the linearized operator at the origin L_ɛ has an essential spectrum filling the entire real line, in addition to the simple eigenvalue at 0. Moreover, for parameter values ɛ < 0 there is a pair of imaginary eigenvalues which meet in 0 for ɛ = 0, and which disappear for ɛ > 0. The above situation occurs for example when one looks for travelling waves in a system of superposed perfect fluid layers, one being infinitely deep. We give quite general assumptions which apply in such physical examples, under which one obtains a family of bifurcating solutions homoclinic to every equilibrium near the origin. These homoclinics are symmetric and decay algebraically at infinity, being approximated at main order by the Benjamin–Ono homoclinic. For the water wave example, this corresponds to a family of solitary waves, such that at infinity the upper layer slides with a uniform velocity, over the bottom layer (at rest).     

New Explicit Expression of Barnett-Lothe Tensors for Anisotropic Linear Elastic Materials

The three Barnett-Lothe tensors H, L, S appear often in the Stroth formalism of two-dimensional deformations of anisotropic elastic materials [1–3]. They also appear in certain three-dimensional problems [4, 5]. The algebraic representation of H, L, S requires computation of the eigenvalues p^v(v=1,2,3) and the normalized eigenvectors (a, b). The integral representation of H, L, S circumvents the need for computing p ^v(v=1,2,3) and (a, b), but it is not simple to integrate the integrals except for special materials. Ting and Lee [6] have recently obtained an explicit expression of H for general anisotropic materials. We present here the remaining tensors L, S using the algebraic representation. They key to our success is the obtaining of the normalization factor for (a, b) in a simple form. The derivation of L and S then makes use of (a, b) but the final result does not require computation of (a, b), which makes the result attractive to numerical computation. Even though the tensor H given in [6] is in terms of the elastic stiffnesses C^μ v while the tensors L, S presented here are in terms of the reduced elastic compliances s′ _μv , the structure of L, S is similar to that of H. Following the derivation of H, we also present alternate expressions of L, S that remain valid for the degenerate cases p ₁ p ₂ and p₁=p₂ = p ₃. One may want to compute H, L, S using either C _μv or s′ _μv v, but not both. We show how an expression in C^μ v can be converted to an expression in s′ _μv v, and vice versa. As an application of the conversion, we present explicit expressions of the extic equation for p in C^μ v and s′ _μv v. This revised version was published online in August 2006 with corrections to the Cover Date.     

Elliptical inhomogeneity with polynomial eigenstrains embedded in orthotropic materials

A closed-form solution for elastic field of an elliptical inhomogeneity with polynomial eigenstrains in orthotropic media having complex roots is presented. The distribution of eigenstrains is assumed to be in the form of quadratic functions in Cartesian coordinates of the points of the inhomogeneity. Elastic energy of inhomogeneity–matrix system is expressed in terms of 18 real unknown coefficients that are analytically evaluated by means of the principle of minimum potential energy and the corresponding elastic field in the inhomogeneity is obtained. Results indicate that quadratic terms in the eigenstrains induce zeroth-order elastic strain components, which reflect the coupling effect of the zeroth- and second-order terms in the polynomial expressions on the elastic field. In contrast, the first-order terms in the eigenstrains only produce corresponding elastic fields in the form of the first-order terms. Numerical examples are given to demonstrate the normal and shear stresses at the interface between the inhomogeneity and the matrix. Furthermore, the solution reduces to known results for the special cases.     

Homogenization of a Boundary-Value Problem with Varying type of Boundary Conditions in a Thick Two-Level 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 thickness of order ɛ = (N ⁻¹). Depending on their lengths, the thin rods are divided into two levels. In addition, the rods from each level are ɛ-periodically alternated. Inhomogeneous Neumann boundary conditions are given on the vertical sides of the thin rods of the first level, and homogeneous Dirichlet boundary conditions are given on the vertical sides of the rods of the second level. We investigate the asymptotic behavior of a solution of this problem as ɛ → 0 and prove a convergence theorem and the convergence of the energy integral. __________ Translated from Neliniini Kolyvannya, Vol. 8, No. 2, pp. 241–257, April–June, 2005.     

Minimizing Sequences Selected via Singular Perturbations, and their Pattern Formation

We study the Cahn-Hilliard energy E _ɛ(u) over the unit square under the constraint of a constant mass m with (ɛ > 0) and without ɛ= 0) interfacial energy. Minimizers of E ₀(u) have no preferred pattern and we select patterns via sequences of conditionally critical points of E _ɛ(u) converging to minimizers as ɛ tends to zero. Those critical points are not minimizers if the singular limit has no minimal interface. We obtain them by a global bifurcation analysis of the Euler-Lagrange equations for E _ɛ(u) where the mass m is the bifurcation parameter. We make use of the symmetry of the unit square, and the elliptic maximum principle, in turn, implies that the location of maxima and minima is fixed for all solutions on global branches. This property is used to guarantee the existence of a singular limit and to verify the Weierstrass-Erdmann corner condition which proves its minimizing property. Accepted January 21, 2000?Published online November 24, 2000     

K-ε-T model of dense liquid-solid two-phase turbulent flow and its application to the pipe flow

To predict the characteristics of dense liquid-solid two-phase flow, K-ε-T model is established, in which the turbulent flow of fluid phase is described with fluid turbulent kinetic energy K_f and its dissipation rate ε_f, and the particles random motion is described with particle turbulent energy K_p and its dissipation rate ε_p and pseudothermal temperature T_p. The governing equations are also derived. With K-ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.     

A Sharp-Interface Limit for a Two-Well Problem in Geometrically Linear Elasticity

In the theory of solid-solid phase transitions the deformation of an elastic body is determined via a functional containing a nonconvex energy density and a singular perturbation. We study Frame indifference, within a linearized setting, requires that W depends only on the symmetric part of ∇u. The potential W is nonnegative and vanishes on two wells, i.e., for d = 2, on two lines in the space of matrices. We determine, for d = 2, the Gamma limit I₀ = Γ− lim _ɛ→0 I_ɛ. The limit I₀[u] is finite only for deformations u that fulfill W(∇u)=0 almost everywhere and have sharp interfaces where ∇u has jumps. For these u, I₀[u] equals the line integral over the interfaces of a surface energy density.     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999     

Identification of Heterogeneous Constitutive Parameters in a Welded Specimen: Uniform Stress and Virtual Fields Methods for Material Property Estimation

Local strain data obtained throughout the entire weld region encompassing both the weld nugget and heat affected zones (HAZs) are processed using two methodologies, uniform stress and virtual fields, to estimate specific heterogeneous material properties throughout the weld zone. Results indicate that (a) the heterogeneous stress–strain behavior obtained by using a relatively simple virtual fields model offers a theoretically sound approach for modeling stress–strain behavior in heterogeneous materials, (b) the local stress–strain results obtained using both a uniform stress assumption and a simplified uniaxial virtual fields model are in good agreement for strains ɛ _xx < 0.025, (c) the weld nugget region has a higher hardening coefficient, higher initial yield stress and a higher hardening exponent, consistent with the fact that the steel weld is overmatched and (d) for ɛ _xx > 0.025, strain localization occurs in the HAZ region of the specimen, resulting in necking and structural effects that complicate the extraction of local stress strain behavior using either of the relatively simple models.

Application of the grating objective speckle method to composite materials

The grating objective speckle method has been applied in whole-field strain analysis of carbon-fiber-woven polyimide-composite materials. A sequence of rosette fringe patterns was obtained from one single specklegram, indicating displacement components along four different directions with a 45-deg interval. The spatial frequencies, which represent the sensitivities of the fringe intervals, were 2400 ℓ/mm forU _x ,U _y and 1697 ℓ/mm forU ₄₅ andU ₁₃₅. Normal strain components, ɛ _x , ɛ _y , ɛ₄₅, and ɛ₁₃₅, were trasformed to quantitative topographic representations by coupling to a digital computer system. Shear strain can be determined by the rosette equations.     

Crack linkup: An experimental analysis

TheT _ɛ ^* integral was used to assess stable crack growth and crack linkup in 0.8 mm thick 2024-T3 aluminum tension specimens with multiple site damage (MSD) under monotonic and cyclic loads. TheT _ɛ ^* values were obtained directly from the recorded moiré fringes on the fracture specimens with and without MSD. TheT _ɛ ^* resistance curves of these fracture specimens of different geometries were in excellent agreement with each other. The results suggest thatT _ɛ ^* is a material parameter which can be used to characterize crack growth and linkup in the absence of large overloading.T _ɛ ^* based crack growth and net-section-yield based crack linkup criteria for MSD specimens are proposed. The crack tip opening angle (CTOA) criterion can also be used to correlate crack growth larger than 2 mm.     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

An Efficient Derivation of the Aronsson Equation

For 1<p<∞, the equation which characterizes minima of the functional u↦∫ _U |Du| ^p ,dx subject to fixed values of u on ∂U is −Δ _p u=0. Here −Δ _p is the well-known ``p-Laplacian'. When p=∞ the corresponding functional is u↦|| |Du|<sup>2</sup>|| <sub>L∞(U)</sub> . A new feature arises in that minima are no longer unique unless U is allowed to vary, leading to the idea of ``absolute minimizers'. Aronsson showed that then the appropriate equation is −Δ_∞ u=0, that is, u is ``infinity harmonic' as explained below. Jensen showed that infinity harmonic functions, understood in the viscosity sense, are precisely the absolute minimizers. Here we advance results of Barron, Jensen and Wang concerning more general functionals u↦||f(x,u,Du)|| _L∞(U) by giving a simplified derivation of the corresponding necessary condition under weaker hypotheses. (Accepted September 6, 2002) Published online April 14, 2003 Communicated by S. Muller     

An experimental study of dynamic tensile properties of glass fiber at low-temperatures

Tensile impact experiments of EC8.0−24×7 glass fiber bundles at different low temperaturesT(14°C, −40°C and −10°C) and strain rates ɛ were carried out, and complete stress-strain curves were obtained. Within the range of the experiment temperatures and strain rates, it is found that the initial modulusE, the ultimate strength σ_max and the unstable strain ɛ _b of the glass fiber bundles all increase with ɛ at an identicalT. At an identical ɛ, with the decrease ofT, E and σ_max increase; but ɛ _b increases when 10°C>T>−40°C and decreases when −40°C>T>−100°C. The strain-rate- and temperature-dependent bimodal Weibull statistical constitutive theory was adopted for the statistical analysis of the experimental results, and the Weibull parameters of single fiber were obtained. The results show that the bimodal Weibull distribution function is suitable to represent the strength distribution of the glass fiber at low temperature and different strain rates. The differences in the mechanical properties between EC8.0−24×7 and EC5.5−12 ×14 glass fiber bundles were also discussed. Project supported by the National Natural Science Foundation of China (No. 19772058).     

Effects of Large Scales Motion in Unstable Stratified Shear Flows

Homogeneous turbulence under unstable uniform stratification (N ² < 0) and vertical shear is investigated by using the linear theory (or the so-called rapid distortion theory, RDT) for an initial isotropic turbulence over a range −∞ ≤ R _i =N ²/S ² ≤ 0. The initial potential energy is zero and P _r =1 (i.e. the molecular Prandtl number).One-dimensional (streamwise) k ₁−spectra, especially Θ₃₃(k ₁) (i.e., that of the vertical kinetic energy, are investigated. In agreement with previous experiments, it is found that the unstable stratification affects the turbulence quantities at all scales. A significant increase of the vertical kinetic energy is observed at low wavenumbers k ₁ (i.e. large scales) due to an increase of the stratification . The effect of the shear (S) is appreciable only at high wavenumbers k ₁ (i.e. small scales).Based on the importance of the spectral components with k ₁ = 0, the asymptotic forms of Θ _ij (k ₁ = 0) or equivalently the so-called “two-dimensional” energy components (2DEC) are analyzed in detail. The asymptotic form for the ratio of 2DEC is compared to the long-time limit of the ratio of real energies. In the unstable stratified shearless case (S=0,N ² ≠ 0) the variances and the covariances of the velocity and the density are derived analytically in terms of the Weber functions, while when S ≠ 0 and N ² ≠ 0 they are obtained numerically (−100 ≤ R _i <0 and . The results are discussed in connection to previous experimental results in unstable stratified open channel flows cooled from above by Komori et al. Phy Fluids 25, 1539–1546 (1982).It is shown that the Richardson number dependence of the long-time limit of the ratios of real energies is well described by this “simple” model (i.e. the dependence of the long-time limit of 2DEC on R _i ). For example, the ratio of the potential energy to the kinetic energy (q ²/2), approaches −R _i /(1−R _i ), the ratio of turbulent energy production by buoyancy forces to production by shearing forces (i.e. the flux Richardson number, R _if ), approaches R _i . Also, the Richardson number dependence of the principal angle (β) of the Reynolds stress tensor and the angle (β_ρ) of the scalar flux vector is fairly predicted by this model .On the other hand, it is found that the above ratios are insensitive to viscosity, while the ratios ɛ /q ² and , depend on the viscosity and they evolve asymptotically like t ⁻¹. The turbulent Froude number, F _rt =(L _Oz /L _E )^2/3, where L _Oz and L _E are the Ozmidov length scale and the Ellison length scale, respectively, evolves asymptotically like t ^−1/3.     

On the Navier�CStokes equations with rotating effect and prescribed outflow velocity

We consider the equations of Navier–Stokes modeling viscous fluid flow past a moving or rotating obstacle in \mathbb R^d{\mathbb R^d} subject to a prescribed velocity condition at infinity. In contrast to previously known results, where the prescribed velocity vector is assumed to be parallel to the axis of rotation, in this paper we are interested in a general outflow velocity. In order to use L ^p -techniques we introduce a new coordinate system, in which we obtain a non-autonomous partial differential equation with an unbounded drift term. We prove that the linearized problem in \mathbb R^d{\mathbb R^d} is solved by an evolution system on L^p_s(\mathbb R^d){L^p_{\sigma}(\mathbb R^d)} for 1 < p < ∞. For this we use results about time-dependent Ornstein–Uhlenbeck operators. Finally, we prove, for p ≥ d and initial data u₀ ? L^p_s(\mathbb R^d){u_0\in L^p_{\sigma}(\mathbb R^d)}, the existence of a unique mild solution to the full Navier–Stokes system.     

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate ɛ_τ) equations were considered. The emphasis of this paper is focused on the effects of the ɛ_τ-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate of change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters, given by /ɛ, _τ/ɛ_τ, , Sk/ɛ and G/ɛ, becoming constant. Here, and _τ are the production of turbulent kinetic energy k and temperature variance , respectively, ɛ and ɛ_τ are their respective dissipation rates, R is the mixed time scale ratio, G is the buoyant production of k and S is the mean shear gradient. Calculations show that the ɛ_τ-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular ɛ_τ-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the ɛ-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence. Received 21 April 2000 and accepted 21 February 2001     

Stability for the 3D Navier–Stokes Equations with Nonzero far Field Velocity on Exterior Domains

Concerning to the non-stationary Navier–Stokes flow with a nonzero constant velocity at infinity, just a few results have been obtained, while most of the results are for the flow with the zero velocity at infinity. The temporal stability of stationary solutions for the Navier–Stokes flow with a nonzero constant velocity at infinity has been studied by Enomoto and Shibata (J Math Fluid Mech 7:339–367, 2005), in L ^p spaces for p ≥ 3. In this article, we first extend their result to the case \frac32 < p{\frac{3}{2} < p} by modifying the method in Bae and Jin (J Math Fluid Mech 10:423–433, 2008) that was used to obtain weighted estimates for the Navier–Stokes flow with the zero velocity at infinity. Then, by using our generalized temporal estimates we obtain the weighted stability of stationary solutions for the Navier–Stokes flow with a nonzero velocity at infinity.