Variable permeability effect on convection in binary mixtures saturating a porous layer

The Darcy Model with the Boussinesq approximation is used to study natural convection in a shallow porous layer, with variable permeability, filled with a binary fluid. The permeability of the medium is assumed to vary exponentially with the depth of the layer. The two horizontal walls of the cavity are subject to constant fluxes of heat and solute while the two vertical ones are impermeable and adiabatic. The governing parameters for the problem are the thermal Rayleigh number, R _T, the Lewis number, Le, the buoyancy ratio, φ, the aspect ratio of the cavity, A, the normalized porosity, ε, the variable permeability constant, c, and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in an infinite layer, an analytical solution of the steady form of the governing equations is obtained on the basis of the parallel flow approximation. The onset of supercritical convection, or subcritical, convection are predicted by the present theory. A linear stability analysis of the parallel flow model is conducted and the critical Rayleigh number for the onset of Hopf’s bifurcation is predicted numerically. Numerical solutions of the full governing equations are found to be in excellent agreement with the analytical predictions.     

Well-posedness and Regularity for the Elasticity Equation with Mixed Boundary Conditions on Polyhedral Domains and Domains with Cracks

We prove a regularity result for the anisotropic linear elasticity equation ${P u := {\rm div} \left( \boldmath\mathsf{C} \cdot \nabla u\right) = f}We prove a regularity result for the anisotropic linear elasticity equationP u : = div ( C ·?u) = f{P u := {\rm div} \left( \boldmath\mathsf{C} \cdot \nabla u\right) = f} , with mixed (displacement and traction) boundary conditions on a curved polyhedral domain W ì \mathbbR³{\Omega \subset \mathbb{R}^3} in weighted Sobolev spaces K^m+1_a+1(W){\mathcal {K}^{m+1}_{a+1}(\Omega)} , for which the weight is given by the distance to the set of edges. In particular, we show that there is no loss of K^m_a{\mathcal {K}^{m}_{a}} -regularity. Our curved polyhedral domains are allowed to have cracks. We establish a well-posedness result when there are no neighboring traction boundary conditions and |a| < η, for some small η > 0 that depends on P, on the boundary conditions, and on the domain Ω. Our results extend to other strongly elliptic systems and higher dimensions.     

Analytical and numerical study of buoyancy-driven convection in a vertical enclosure filled with nanofluids

This paper presents an analytical and numerical study of natural convection of nanofluids contained in a rectangular enclosure subject to uniform heat flux along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number Ra, the Prandtl number Pr, the aspect ratio of the cavity A and the solid volume fraction of nanoparticles, Φ. Three types of nanoparticles are taken into consideration: Cu, Al₂O₃ and TiO₂. Various models are used for calculating the effective viscosity and thermal conductivity of nanofluids. In the first part of the analytical study, a scale analysis is made for the boundary layer regime situation. In the second part, an analytical solution based on the parallel flow approximation is reported for tall enclosures (A ≫ 1). In the boundary layer regime a good agreement is obtained between the predictions of the scale analysis and those of the analytical solution. Solutions for the flow fields, temperature distributions and Nusselt numbers are obtained explicitly in terms of the governing parameters of the problem. A numerical study of the same phenomenon, obtained by solving the complete system of the governing equations, is also conducted. A good agreement is found between the analytical predictions and the numerical simulations.     

On the three-dimensional vibrations of a hollow elastic torus of annular cross-section

This paper offers a three-dimensional elasticity-based variational Ritz procedure to examine the natural vibrations of an elastic hollow torus with annular cross-section. The associated energy functional minimized in the Ritz procedure is formulated using toroidal coordinates (r,q, j)({r,\theta , \varphi}) comprised of the usual polar coordinates (r, θ) originating at each circular cross-sectional center and a circumferential coordinate j{\varphi} around the torus originating at the torus center. As an enhancement to conventional use of algebraic–trigonometric polynomials trial series in related solid body vibration studies in the associated literature, the assumed torus displacement, u, v and w in the r, θ and j{\varphi } toroidal directions, respectively, are approximated in this work as a triplicate product of Chebyshev polynomials in r and the periodic trigonometric functions in the θ and j{\varphi} directions along with a set of generalized coefficients. Upon invoking the stationary condition of the Lagrangian energy functional for the elastic torus with respected to these generalized coefficients, the usual characteristic frequency equations of natural vibrations of the elastic torus are derived. Upper bound convergence of the first seven non-dimensional frequency parameters accurate to at least five significant figures is achieved by using only ten terms of the trial torus displacement functions. Non-dimensional frequencies of elastic hollow tori are examined showing the effects of varying torus radius ratio and cross-sectional radius ratio.     

Scaling Transformations for Boundary Layer Flow near the Stagnation-Point on a Heated Permeable Stretching Surface in a Porous Medium Saturated with a Nanofluid and Heat Generation/Absorption Effects

In this article, a similarity solution of the steady boundary layer flow near the stagnation-point flow on a permeable stretching sheet in a porous medium saturated with a nanofluid and in the presence of internal heat generation/absorption is theoretically studied. The governing partial differential equations with the corresponding boundary conditions are reduced to a set of ordinary differential equations with the appropriate boundary conditions via Lie-group analysis. Copper (Cu) with water as its base fluid has been considered and representative results have been obtained for the nanoparticle volume fraction parameter f{\phi} in the range 0 ￡ f ￡ 0.2{0\leq \phi \leq 0.2} with the Prandtl number of Pr = 6.8 for the water working fluid. Velocity and temperature profiles as well as the skin friction coefficient and the local Nusselt number are determined numerically. The influence of pertinent parameters such as nanofluid volume fraction parameter, the ratio of free stream velocity and stretching velocity parameter, the permeability parameter, suction/blowing parameter, and heat source/sink parameter on the flow and heat transfer characteristics is discussed. Comparisons with published results are also presented. It is shown that the inclusion of a nanoparticle into the base fluid of this problem is capable to change the flow pattern.     

Double-diffusive parallel flow induced in a horizontal Brinkman porous layer subjected to constant heat and mass fluxes: analytical and numerical studies

Fluid flow and heat and mass transfer induced by double-diffusive natural convection in a horizontal porous layer subjected to vertical gradients of temperature and concentration are studied analytically and numerically using the Brinkman-extended Darcy model. Both cases of rigid and free horizontal boundaries are examined in the present study. The parameters governing the problem are the Rayleigh number R_T, the Lewis number Le, the buoyancy ratio N, the Darcy number Da and the aspect ratio Ar. The analytical solution is based on the parallel flow approximation. The critical Rayleigh number corresponding to the onset of the parallel flow in this system is determined analytically as a function of Le, N and Da. For sufficiently small Da, both free and rigid boundaries yield results which are identical to those predicted by the Darcy model. The present investigation shows that there exists a region in the plane (N, Le) where the convective flow is not possible in the layer regardless of the Rayleigh and Darcy numbers considered. Received on 21 December 1998     

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.     

Positivity and Almost Positivity of Biharmonic Green’s Functions under Dirichlet Boundary Conditions

In general, for higher order elliptic equations and boundary value problems like the biharmonic equation and the linear clamped plate boundary value problem, neither a maximum principle nor a comparison principle or—equivalently—a positivity preserving property is available. The problem is rather involved since the clamped boundary conditions prevent the boundary value problem from being reasonably written as a system of second order boundary value problems. It is shown that, on the other hand, for bounded smooth domains W ì \mathbbRⁿ{\Omega \subset\mathbb{R}^n} , the negative part of the corresponding Green’s function is “small” when compared with its singular positive part, provided n\geqq 3{n\geqq 3} . Moreover, the biharmonic Green’s function in balls B ì \mathbbRⁿ{B\subset\mathbb{R}^n} under Dirichlet (that is, clamped) boundary conditions is known explicitly and is positive. It has been known for some time that positivity is preserved under small regular perturbations of the domain, if n = 2. In the present paper, such a stability result is proved for n\geqq 3{n\geqq 3} .     

Maximal L p – L q -Estimates for the Stokes Equation: a Short Proof of Solonnikov’s Theorem

Introducing a new localization method involving Bogovskiĭ's operator we give a short and new proof for maximal L^p – L^q-estimates for the solution of the Stokes equation. Moreover, it is shown that, up to constants, the Stokes operator is an R{\mathcal{R}}-sectorial operator in L^p_s(W)L^{p}_{\sigma}(\Omega), 1 < p < ￥1 < p < \infty, of R{\mathcal{R}}-angle 0, for bounded or exterior domains of Ω.     

Hydrodynamic Limit of Gradient Exclusion Processes with Conductances

Fix a strictly increasing right continuous with left limits function ${W: \mathbb{R} \to \mathbb{R}}Fix a strictly increasing right continuous with left limits function W: \mathbbR ? \mathbbR{W: \mathbb{R} \to \mathbb{R}} and a smooth function F: [l,r] ? \mathbb R{\Phi : [l,r] \to \mathbb R}, defined on some interval [l, r] of \mathbb R{\mathbb R}, such that 0 < b\leqq F￠\leqq b^-1{0 < b\leqq \Phi'\leqq b^{-1}}. On the diffusive time scale, the evolution of the empirical density of exclusion processes with conductances given by W is described by the unique weak solution of the non-linear differential equation ?_t r = (d/dx)(d/dW) F(r){\partial_t \rho = ({\rm d}/{\rm d}x)({\rm d}/{\rm d}W) \Phi(\rho)}. We also present some properties of the operator (d/dx)(d/dW).     

Generalized Solution for 1-D Non-Newtonian Flow in a Porous Domain due to an Instantaneous Mass Injection

Non-Newtonian fluid flow through porous media is of considerable interest in several fields, ranging from environmental sciences to chemical and petroleum engineering. In this article, we consider an infinite porous domain of uniform permeability k and porosity f{\phi} , saturated by a weakly compressible non-Newtonian fluid, and analyze the dynamics of the pressure variation generated within the domain by an instantaneous mass injection in its origin. The pressure is taken initially to be constant in the porous domain. The fluid is described by a rheological power-law model of given consistency index H and flow behavior index n; n, < 1 describes shear-thinning behavior, n > 1 shear-thickening behavior; for n = 1, the Newtonian case is recovered. The law of motion for the fluid is a modified Darcy’s law based on the effective viscosity μ _ef , in turn a function of f, H, n{\phi, H, n} . Coupling the flow law with the mass balance equation yields the nonlinear partial differential equation governing the pressure field; an analytical solution is then derived as a function of a self-similar variable η = rt ^β (the exponent β being a suitable function of n), combining spatial coordinate r and time t. We revisit and expand the work in previous papers by providing a dimensionless general formulation and solution to the problem depending on a geometrical parameter d, valid for plane (d = 1), cylindrical (d = 2), and semi-spherical (d = 3) geometry. When a shear-thinning fluid is considered, the analytical solution exhibits traveling wave characteristics, in variance with Newtonian fluids; the front velocity is proportional to t ^(n-2)/2 in plane geometry, t ^{(2n-3)/(3−n)} in cylindrical geometry, and t ^{(3n-4)/[2(2−n)]} in semi-spherical geometry. To reflect the uncertainty inherent in the value of the problem parameters, we consider selected properties of fluid and matrix as independent random variables with an associated probability distribution. The influence of the uncertain parameters on the front position and the pressure field is investigated via a global sensitivity analysis evaluating the associated Sobol’ indices. The analysis reveals that compressibility coefficient and flow behavior index are the most influential variables affecting the front position; when the excess pressure is considered, compressibility and permeability coefficients contribute most to the total response variance. For both output variables the influence of the uncertainty in the porosity is decidedly lower.     

Well-posedness for the heat flow of harmonic maps and the liquid crystal flow with rough initial data

We investigate the well-posedness of (1) the heat flow of harmonic maps from \mathbb Rⁿ{\mathbb R^n} to a compact Riemannian manifold N without boundary for initial data in BMO; and (2) the hydrodynamic flow (u, d) of nematic liquid crystals on \mathbb Rⁿ{\mathbb R^n} for initial data in BMO⁻¹ × BMO.     

Boundary Regularity in Variational Problems

We prove that, if ${u : \Omega \subset \mathbb{R}^n \to \mathbb{R}^N}We prove that, if u : W ì \mathbbRⁿ ? \mathbbR^N{u : \Omega \subset \mathbb{R}^n \to \mathbb{R}^N} is a solution to the Dirichlet variational problem

minwòW F(x, w, Dw) dx     subject  to     w o u0  on  ?W,\mathop {\rm min}\limits_{w}\int_{\Omega} F(x, w, Dw)\,{\rm d}x \quad {\rm subject \, to} \quad w \equiv u_0\; {\rm on}\;\partial \Omega,      14. Large Time Well-posedness of the Three-Dimensional Capillary-Gravity Waves in the Long Wave Regime      Mei?Ming  Ping?Zhang  Zhifei?ZhangEmail author 《Archive for Rational Mechanics and Analysis》2012,204(2):387-444 In the regime of weakly transverse long waves, given long-wave initial data, we prove that a non-dimensionalized water wave system in an infinite strip under the influence of gravity and surface tension on the upper free interface has a unique solution on [0, T/ e{0, T/ \varepsilon} ] for some e{\varepsilon} independent of constant T. We shall prove in the subsequent paper (Ming et al., The long wave approximation to the three-dimensional capillary gravity waves, 2011) that on the same time interval, these solutions can be accurately approximated by sums of solutions of two decoupled Kadomtsev–Petviashvili (KP) equations.      15. Centers and Isochronous Centers for Two Classes of Generalized Seventh and Ninth Systems      Jaume Llibre  Clàudia Valls 《Journal of Dynamics and Differential Equations》2010,22(4):657-675 We classify new classes of centers and of isochronous centers for polynomial differential systems in \mathbb R2{\mathbb R^2} of arbitrary odd degree d ≥ 7 that in complex notation z = x + i y can be written as [(z)\dot] = (l+i) z + (z[`(z)])\fracd-7-2j2 (A z5+j[`(z)]2+j + B z4+j[`(z)]3+j + C z3+j[`(z)]4+j+D[`(z)]7+2j ),\dot z = (\lambda+i) z + (z \overline z)^{\frac{d-7-2j}2} \left(A z^{5+j} \overline z^{2+j} + B z^{4+j} \overline z^{3+j} + C z^{3+j} \overline z^{4+j}+D \overline z^{7+2j} \right),      16. On the Symmetry of Energy-Minimising Deformations in Nonlinear Elasticity I: Incompressible Materials      Jeyabal Sivaloganathan  Scott J. Spector 《Archive for Rational Mechanics and Analysis》2010,196(2):363-394 Let ${A=\{{\bf x} \in \mathbb{R}^n : a < |{\bf x}| < b\}, n \geqq 2, a > 0}Consider a homogeneous, isotropic, hyperelastic body occupying the region A = {x ? \mathbbRn :  a <  |x |  < b }{A = \{{\bf x}\in\mathbb{R}^{n}\, : \,a <\,|{\bf x} |\,< b \}} in its reference state and subject to radially symmetric displacement, or mixed displacement/traction, boundary conditions. In Part I (Sivaloganathan and Spector in Arch Ration Mech Anal, 2009, in press) the authors restricted their attention to incompressible materials. For a large-class of polyconvex constitutive relations that grow sufficiently rapidly at infinity it was shown that to each nonradial isochoric deformation of A there corresponds a radial isochoric deformation that has strictly less elastic energy than the given deformation. In this paper that analysis is further developed and extended to the compressible case. The key ingredient is a new radial-symmetrisation procedure that is appropriate for problems where the symmetrised mapping must be one-to-one in order to prevent interpenetration of matter. For the pure displacement boundary-value problem, the radial symmetrisation of an orientation-preserving diffeomorphism u : A → A* between spherical shells A and A* is the deformation urad(x)=\fracr(R)Rx,     R=|x|,                        (0.1){\bf u}^{\rm rad}({\bf x})=\frac{r(R)}{R}{\bf x}, \quad R=|{\bf x}|,\qquad\qquad\qquad\qquad(0.1)      17. Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature   总被引：1，自引：0，他引：1   M. Kumari  G. Nath 《Transport in Porous Media》2009,77(3):463-474 Unsteady natural convection flow in a two-dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left-hand vertical wall has temperature T h and the right-hand vertical wall is maintained at temperature T c (T h > T c) and the horizontal walls are insulated. At time t > 0, the left-hand vertical wall temperature is suddenly raised to which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.      18. Stability of Degenerate Stationary Waves for Viscous Gases      Yoshihiro Ueda  Tohru Nakamura  Shuichi Kawashima 《Archive for Rational Mechanics and Analysis》2010,198(3):735-762 This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t −α/4 as t → ∞, provided that the initial perturbation is in the weighted space L2a=L2(\mathbb R+; (1+x)a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t −α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α*(q), where α*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L2a{L^2_\alpha} for α > α*(q) with another critical value α*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.      19. Lower Bound for the Energy of Bloch Walls in Micromagnetics      Radu Ignat  Benoît Merlet 《Archive for Rational Mechanics and Analysis》2011,199(2):369-406 We study a two-dimensional nonconvex and nonlocal energy in micromagnetics defined over S 2-valued vector fields. This energy depends on two small parameters, β and e{\varepsilon} , penalizing the divergence of the vector field and its vertical component, respectively. Our objective is to analyze the asymptotic regime b << e << 1{\beta \ll \varepsilon \ll 1} through the method of Γ-convergence. Finite energy configurations tend to become divergence-free and in-plane in the magnetic sample except in some small regions of typical width e{\varepsilon} (called Bloch walls) where the magnetization connects two directions on S 2. We are interested in quantifying the limit energy of the transition layers in terms of the jump size between these directions. For one-dimensional transition layers, we show by Γ-convergence analysis that the exact line density of the energy is quadratic in the jump size. We expect the same behaviour for the two-dimensional model. In order to prove that, we investigate the concept of entropies. In the prototype case of a periodic strip, we establish a quadratic lower bound for the energy with a non-optimal constant. Then we introduce and study a special class of Lipschitz entropies and obtain lower bounds coinciding with the one-dimensional Γ-limit in some particular cases. Finally, we show that entropies are not appropriate in general for proving the expected sharp lower bound.      20. Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall      M. Pavlidis  Y. Dimakopoulos  J. Tsamopoulos 《Rheologica Acta》2009,48(9):1031-1048 The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]s/([(h)\tilde]s +[(h)\tilde]p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]p\widetilde{\eta}_p and [(h)\tilde]s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]2/([(h)\tilde]p +[(h)\tilde]s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.      设为首页 | 免责声明 | 关于勤云 | 加入收藏 Copyright©北京勤云科技发展有限公司  京ICP备09084417号

Large Time Well-posedness of the Three-Dimensional Capillary-Gravity Waves in the Long Wave Regime

In the regime of weakly transverse long waves, given long-wave initial data, we prove that a non-dimensionalized water wave system in an infinite strip under the influence of gravity and surface tension on the upper free interface has a unique solution on [0, T/ e{0, T/ \varepsilon} ] for some e{\varepsilon} independent of constant T. We shall prove in the subsequent paper (Ming et al., The long wave approximation to the three-dimensional capillary gravity waves, 2011) that on the same time interval, these solutions can be accurately approximated by sums of solutions of two decoupled Kadomtsev–Petviashvili (KP) equations.     

Centers and Isochronous Centers for Two Classes of Generalized Seventh and Ninth Systems

We classify new classes of centers and of isochronous centers for polynomial differential systems in \mathbb R²{\mathbb R^2} of arbitrary odd degree d ≥ 7 that in complex notation z = x + i y can be written as

On the Symmetry of Energy-Minimising Deformations in Nonlinear Elasticity I: Incompressible Materials

Let ${A=\{{\bf x} \in \mathbb{R}^n : a < |{\bf x}| < b\}, n \geqq 2, a > 0}Consider a homogeneous, isotropic, hyperelastic body occupying the region A = {x ? \mathbbRⁿ :  a <  |x |  < b }{A = \{{\bf x}\in\mathbb{R}^{n}\, : \,a <\,|{\bf x} |\,< b \}} in its reference state and subject to radially symmetric displacement, or mixed displacement/traction, boundary conditions. In Part I (Sivaloganathan and Spector in Arch Ration Mech Anal, 2009, in press) the authors restricted their attention to incompressible materials. For a large-class of polyconvex constitutive relations that grow sufficiently rapidly at infinity it was shown that to each nonradial isochoric deformation of A there corresponds a radial isochoric deformation that has strictly less elastic energy than the given deformation. In this paper that analysis is further developed and extended to the compressible case. The key ingredient is a new radial-symmetrisation procedure that is appropriate for problems where the symmetrised mapping must be one-to-one in order to prevent interpenetration of matter. For the pure displacement boundary-value problem, the radial symmetrisation of an orientation-preserving diffeomorphism u : A → A* between spherical shells A and A* is the deformation

urad(x)=\fracr(R)Rx,     R=|x|,                        (0.1){\bf u}^{\rm rad}({\bf x})=\frac{r(R)}{R}{\bf x}, \quad R=|{\bf x}|,\qquad\qquad\qquad\qquad(0.1)      17. Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature   总被引：1，自引：0，他引：1   M. Kumari  G. Nath 《Transport in Porous Media》2009,77(3):463-474 Unsteady natural convection flow in a two-dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left-hand vertical wall has temperature T h and the right-hand vertical wall is maintained at temperature T c (T h > T c) and the horizontal walls are insulated. At time t > 0, the left-hand vertical wall temperature is suddenly raised to which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.      18. Stability of Degenerate Stationary Waves for Viscous Gases      Yoshihiro Ueda  Tohru Nakamura  Shuichi Kawashima 《Archive for Rational Mechanics and Analysis》2010,198(3):735-762 This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t −α/4 as t → ∞, provided that the initial perturbation is in the weighted space L2a=L2(\mathbb R+; (1+x)a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t −α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α*(q), where α*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L2a{L^2_\alpha} for α > α*(q) with another critical value α*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.      19. Lower Bound for the Energy of Bloch Walls in Micromagnetics      Radu Ignat  Benoît Merlet 《Archive for Rational Mechanics and Analysis》2011,199(2):369-406 We study a two-dimensional nonconvex and nonlocal energy in micromagnetics defined over S 2-valued vector fields. This energy depends on two small parameters, β and e{\varepsilon} , penalizing the divergence of the vector field and its vertical component, respectively. Our objective is to analyze the asymptotic regime b << e << 1{\beta \ll \varepsilon \ll 1} through the method of Γ-convergence. Finite energy configurations tend to become divergence-free and in-plane in the magnetic sample except in some small regions of typical width e{\varepsilon} (called Bloch walls) where the magnetization connects two directions on S 2. We are interested in quantifying the limit energy of the transition layers in terms of the jump size between these directions. For one-dimensional transition layers, we show by Γ-convergence analysis that the exact line density of the energy is quadratic in the jump size. We expect the same behaviour for the two-dimensional model. In order to prove that, we investigate the concept of entropies. In the prototype case of a periodic strip, we establish a quadratic lower bound for the energy with a non-optimal constant. Then we introduce and study a special class of Lipschitz entropies and obtain lower bounds coinciding with the one-dimensional Γ-limit in some particular cases. Finally, we show that entropies are not appropriate in general for proving the expected sharp lower bound.      20. Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall      M. Pavlidis  Y. Dimakopoulos  J. Tsamopoulos 《Rheologica Acta》2009,48(9):1031-1048 The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]s/([(h)\tilde]s +[(h)\tilde]p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]p\widetilde{\eta}_p and [(h)\tilde]s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]2/([(h)\tilde]p +[(h)\tilde]s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two-dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left-hand vertical wall has temperature T _h and the right-hand vertical wall is maintained at temperature T _c (T _h > T _c) and the horizontal walls are insulated. At time t > 0, the left-hand vertical wall temperature is suddenly raised to which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.     

Stability of Degenerate Stationary Waves for Viscous Gases

This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t ^−α/4 as t → ∞, provided that the initial perturbation is in the weighted space L²_a=L²(\mathbb R₊; (1+x)^a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t ^−α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α_*(q), where α_*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L²_a{L^2_\alpha} for α > α^*(q) with another critical value α^*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.     

Lower Bound for the Energy of Bloch Walls in Micromagnetics

We study a two-dimensional nonconvex and nonlocal energy in micromagnetics defined over S ²-valued vector fields. This energy depends on two small parameters, β and e{\varepsilon} , penalizing the divergence of the vector field and its vertical component, respectively. Our objective is to analyze the asymptotic regime b << e << 1{\beta \ll \varepsilon \ll 1} through the method of Γ-convergence. Finite energy configurations tend to become divergence-free and in-plane in the magnetic sample except in some small regions of typical width e{\varepsilon} (called Bloch walls) where the magnetization connects two directions on S ². We are interested in quantifying the limit energy of the transition layers in terms of the jump size between these directions. For one-dimensional transition layers, we show by Γ-convergence analysis that the exact line density of the energy is quadratic in the jump size. We expect the same behaviour for the two-dimensional model. In order to prove that, we investigate the concept of entropies. In the prototype case of a periodic strip, we establish a quadratic lower bound for the energy with a non-optimal constant. Then we introduce and study a special class of Lipschitz entropies and obtain lower bounds coinciding with the one-dimensional Γ-limit in some particular cases. Finally, we show that entropies are not appropriate in general for proving the expected sharp lower bound.     

Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall

The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]_s/([(h)\tilde]_s +[(h)\tilde]_p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]_p\widetilde{\eta}_p and [(h)\tilde]_s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]²/([(h)\tilde]_p +[(h)\tilde]_s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.