Uniform Local Existence for Inhomogeneous Rotating Fluid Equations

We investigate the equations of anisotropic incompressible viscous fluids in , rotating around an inhomogeneous vector B(t, x ₁, x ₂). We prove the global existence of strong solutions in suitable anisotropic Sobolev spaces for small initial data, as well as uniform local existence result with respect to the Rossby number in the same functional spaces under the additional assumption that B = B(t, x ₁) or B = B(t, x ₂). We also obtain the propagation of the isotropic Sobolev regularity using a new refined product law.     

Interactions in general continua with microstructure

This paper studies the behavior of the one dimensional Broadwell model of a discrete three velocity gas on a bounded domain 0 x 1 with specularly reflective boundary condition at x = 0, 1. For smooth initial data we show that the initial boundary value problem possesses a unique smooth solution which tends as t to a free state consisting of traveling waves f _1e (x – ct), f _2e (x + ct), f _3e (x) where each f _ie is 2-periodic. The convergence is in the weak* topology of an appropriate Orlicz-Banach state space. No smallness assumptions are made on the data.In memory of Ronald J. DiPerna     

Analysis of fetal cortical complexity from MR images using 3D entropy based information fractal dimension

The fetal cortical complexity is a significant quantification for assessing the development of fetal brain. This study attempts to quantify the development of fetal cortical complexity using the concept of fractal dimension (FD) analysis. Thirty-two fetal MR images were selected from Taipei Veterans General Hospital at 27–37 weeks of gestational age (GA). To investigate the FD of fetal cortical complexity, the entropy based information fractal dimension method (FD _EBI), which is modified from Box-Counting method, was adopted and extended from 2D to 3D. The FD results from overall whole fetal brains show that the increase of cortical complexity is highly correlated with the gestational age of the fetus. Moreover, the FD values of right hemispheric brain are larger than those of left hemispheric brain, show that the development of right hemispheric fetal cortical complexity earlier than the left. These results are in good agreement with normal fetal brain development and suggest that the FD is an effective means for the quantification of fetal cortical complexity.     

Contact Transformation Group Classification of Nonlinear Wave Equations

Interest in nonlinear wave equations has been stimulated bynumerous physical applications, such as telecommunication (e.g.nonlinear telegrapher equation), gasdynamics, anisotropic plasticity andnonlinear elasticity, etc. Mathematical models of these phenomena canoften be reduced to particular types of the equation u _tt = f(x, u _x ) u _xx + g(x, u _x ). In this paper, the problem ofclassification of the latter equation with respect to admitted contacttransformation groups is reduced to the investigation of pointtransformation groups of the equivalent system of first-orderquasi-linear equations v _t =a(x, v)w _x , w _t = b(x,v)v _x .     

Two‐dimensional dispersion analyses of finite element approximations to the shallow water equations

Dispersion analysis of discrete solutions to the shallow water equations has been extensively used as a tool to define the relationships between frequency and wave number and to determine if an algorithm leads to a dual wave number response and near 2Δx oscillations. In this paper, we explore the application of two‐dimensional dispersion analysis to cluster based and Galerkin finite element‐based discretizations of the primitive shallow water equations and the generalized wave continuity equation (GWCE) reformulation of the harmonic shallow water equations on a number of grid configurations. It is demonstrated that for various algorithms and grid configurations, contradictions exist between the results of one‐dimensional and two‐dimensional dispersion analysis as a result of subtle changes in the mass matrix. Numerical experiments indicate that the two‐dimensional dispersion analysis correctly predicts the existence and onset of near 2Δx noise in the solution. Copyright © 2004 John Wiley & Sons, Ltd.     

Laminar diffusion of suspended particulate matter using a two phase flow model

The present paper envisages laminar mixing of a two‐dimensional jet of particulate suspension in an incompressible carrier fluid with a free stream in direction of the jet axis. Finite difference technique has been employed for finding out solution of governing equations. It is found that the diffusion parameter ε, the ratio of particle diffusion coefficient and kinematic viscosity of the carrier fluid, have significant influence on the concentration of particles. A large value of ε has the effect in increasing the perturbation velocity u_p and perturbation density ρ_p. It is observed that the volume fraction φ, has no significant effect on perturbation velocity u and u_p but has profound effect on perturbation velocity v and v_p. It is also found that the particle phase as well as the carrier fluid velocity attain free stream value for the large ξ, the modified x‐co‐ordinate. Further the magnitude of the perturbation quantities u, u_p, v, v_p decreases as ξ increases i.e. at far away from the nozzle exit. Copyright © 2002 John Wiley & Sons, Ltd.     

On Aerodynamic Forces for Viscous Compressible Flow

The paper is aimed at reviewing and adding some new results to our recent work on a force theory for viscous compressible flows around a finite body. It has been proposed to analyze aerodynamic forces directly in terms of fluid elements of nonzero vorticity and density gradient. Let ρ denote the density, u the velocity, and ω the vorticity. It is demonstrated that for largely separated flows about bluff bodies, there are two major source elements: R _e(x) =−?u ²∇ρ·∇ϕ and V _e(x) =−u×ω·∇ϕ, where ϕ is an acyclic potential, generated by the solid body moving with unit velocity in the negative direction of the force considered. In particular, under mild conditions, the (unique) choice of ϕ enforces that the elements R _e(x) and V _e(x) decay rapidly away from the body. Four kinds of finite body are considered: a circular cylinder, a sphere, a hemi sphere-cylinder, and a delta wing of elliptic section—all in transonic-to-supersonic regimes. From an extensive numerical study carried out for these bodies, it is found that these two elements contribute to 95% or more of the total drag or lift for all the cases under consideration. Moreover, R _e(x) due to density gradient becomes progressively important relative to V _e(x) due to vorticity as the Mach number increases. The present method of force analysis enables effective analysis and assessment of relative importance of aerodynamics forces, contributed from individual flow structures. The analysis could therefore be very much useful in view of the rapid growth in numerical fluid dynamics; detailed (either local or global) flow information is often available. The paper is dedicated to Sir James Lighthill in honor of his great contributions to aeronautics on the occasion of the publication of his collected works. Received 3 January 1997 and accepted 11 April 1997     

Percolation Cluster Statistics and Conductivity Semi-Variograms

Cluster statistics of percolation theory have been shown to generate expressions for the distribution of hydraulic conductivity values in accord with field studies. Percolation theory yields directly the smallest possible generalized resistance value, R _c, for which a continuous path through an infinite heterogeneous system can avoid all larger resistances. R _c, defines an infinite system hydraulic conductivity. Cluster statistics generate the number of clusters of resistors of a given size with a given R, for which a continuous path through the cluster can avoid resistances larger than R. The probability that a volume of size x ³ falls on a particular cluster gives the probability that volume has a characteristic resistance, R. Determining the semi-variogram of the hydraulic conductivity is now elementary; it is necessary only to determine whether translation h of the center of the volume x ³ removes it from the cluster in question. If the cluster is larger than (x+h)³, then, on the average, the same cluster resistance R will control K. Otherwise, the value of K at x+h will be uncorrelated with its value at x. The condition is then expressed as an integral related to the one, which gives the distribution of K. Then an integral over the derived distribution of K gives the variogram. Results obtained are that the variogram should be similar to either the exponential or Gaussian forms typically in use, if K is a power law function of random variables (as in Poiseuilles Law), or more closely related to the spherical approximation if K is an exponential function of random variables.     

Auto-correlation measurements and integral time scales in three-dimensional turbulent boundary layers

Auto-correlation, time and length scales of the three components of turbulence and power spectra in a three-dimensional turbulent boundary layer developing on a yawed flat plate have been obtained. The measurements indicate that close to the wall, in the region of turbulence production, there is a marked disparity among the time scales but as the outer edge of the boundary layer is approached, the scales become comparable to one another. Also, the behaviour of the length scales and the power spectra across the boundary layer is presented.Nomenclature Boundary layer thickness where Q/Q _e=0.995 - E _u(f) one dimensional frequency spectra - f frequency in Hz - k ₁ wave number defined as k ₁=2f/Q - L length scale defined as: time scale times local mean velocity - Q local mean velocity - Q _e free stream velocity - R _u, R _v, R _w Auto-correlation coefficients of u, v and w respectively as defined in equation (1) - T _u, T _v, T _w the time scales of u, v and w fluctuations as defined in equation (2) - delay time - u fluctuating velocity component in x-direction - v fluctuation velocity component in y-direction - w fluctuation velocity component in z-direction - x coordinate axis in the streamwise direction - y coordinate axis normal to the surface - z coordinate axis normal to the x-direction and parallel to the wall     

Sensitivity of mesh spacing on simulating macrosegregation during directional solidification of a superalloy

The sensitivity of mesh spacing on simulations of macrosegregation, particularly ‘freckles’, during vertical directional solidification of a superalloy in a rectangular mold was systematically analyzed to achieve accurate predictions in finite element calculations. It was observed that a coarser mesh spacing in the x‐direction horizontal tends to minimize the simulated macrosegregation, whereas a coarser mesh spacing in the y‐direction vertical artificially tends to make the system appear to have more macrosegregation. When solidification conditions either lead to a well‐established freckling case or to a well‐established non‐freckling case, the simulated results are not sensitive to the mesh spacing provided the elements are no larger than about 2d₁ by 2 D/V and 3d₁ by 4 D/V respectively, where d₁ is the primary dendrite arm spacing, D is the diffusivity of the alloy solute with the smallest diffusivity in the liquid, and V is the growth rate. However, when solidification conditions are very close to the transition between freckling and no freckling, the simulated results are sensitive to the mesh spacing, especially in the y‐direction. Based on the mesh sensitivity analysis from the two‐dimensional simulations of rectangular castings of René N5, the mesh with element dimensions no larger than 2d₁ in the x‐direction and 1.5 D/V in the y‐direction are recommended as the most stringent element size. Copyright © 2001 John Wiley & Sons, Ltd.     

Explicit expressions of the Barnett-Lothe tensors for anisotropic materials

The three Barnett-Lothe tensorsS, H, andL appear frequently in the real form solutions of two-dimensional anisotropic elasticity problems. Explicit expressions for the components of these tensors are presented for general anisotropic materials. The special cases of monoclinic materials with the plane of material symmetry at x₃=0, x₂=0, and x₁=0 are then deduced. For monoclinic materials with the symmetry plane at x₂=0 or x₁=0, the locations of image singularities for the Green's functions for a half-space have a special geometry.     

Asymptotic Solutions of Problems of One-Dimensional Time-Dependent Combustible Gas Flows in the Presence of Thermal Action

In a development of studies [1, 2], asymptotic solutions of the Navier-Stokes equations are found for one-dimensional combustible gas flows in the presence of various forms of thermal action on a moving surface (x=x _w(t)). In the problems considered, the temperature or the heat flux q _w(t) is specified on the surface or the surface is the interface between a combustible gas and a moving heated piston or another gas (for example, in a shock tube). Use is made of the fact that, as t , in many cases the values of v _w=(dx/dt)_w and q _w are bounded. This leads to a steady-state flow in the flame zone in the coordinate system moving with its front and homogeneous uniform flow ahead of and behind it. Solutions of all these problems are given for the burnt-gas boundary layer region adjacent to the surface. The numerical calculations performed confirm the results obtained. A velocity law leading to time invariability of the flow pattern obtained with allowance for the interaction between the boundary layer and the burnt-gas homogeneous flow is found, including in the problem of the breakdown of an arbitrary discontinuity. The results are generalized to include the case of motion at an angle of incidence with an additional velocity component aligned with the surface.     

Symmetry of Mountain Pass Solutions of Some Vector Field Equations

We prove radial symmetry (or axial symmetry) of the mountain pass solution of variational elliptic systems − AΔu(x) + ∇ F(u(x)) = 0 (or − ∇.(A(r) ∇ u(x)) + ∇ F(r,u(x)) = 0,) u(x) = (u ₁(x),...,u _N (x)), where A (or A(r)) is a symmetric positive definite matrix. The solutions are defined in a domain Ω which can be , a ball, an annulus or the exterior of a ball. The boundary conditions are either Dirichlet or Neumann (or any one which is invariant under rotation). The mountain pass solutions studied here are given by constrained minimization on the Nehari manifold. We prove symmetry using the reflection method introduced in Lopes [(1996), J. Diff. Eq. 124, 378–388; (1996), Eletron. J. Diff. Eq. 3, 1–14].     

Heat transfer characteristics of boundary-layer flows induced by continuous surfaces stretched with prescribed skin friction

A continuous surface stretched with velocity u _w=u _w (x) and having the temperature distribution T _w=T _w (x) interacts with the viscous fluid in which it is immersed both mechanically and thermally. The thermal interaction is characterized by the surface heat flux q _w=q _w (x) and the mechanical one by the skin friction τ _w=τ _w (x). In the whole previous theoretical research concerned with such processes, either (u _w and T _w) or (u _w and q _w) have been prescribed as known boundary conditions. The goal of the present paper is to initiate the investigation of the boundary layer flows induced by stretching processes for which either (τ _w and T _w ) or (τ _w and q _w) are the prescribed quantities. The case of an isothermal surface stretched with constant skin friction, (τ _w=const., T _w=const. ≠ T _∞) is worked out in detail. The corresponding flow and heat transfer characteristics are compared to those obtained for the (well known) case of a uniformly moving isothermal surface (u _w=const., T _w=const. ≠ T _∞).     

Asymptotic Variational Wave Equations

We investigate the equation (u _t +(f(u)) _x ) _x =f ^{′ ′}(u) (u _x )²/2 where f(u) is a given smooth function. Typically f(u)=u ²/2 or u ³/3. This equation models unidirectional and weakly nonlinear waves for the variational wave equation u _tt − c(u) (c(u)u _x ) _x =0 which models some liquid crystals with a natural sinusoidal c. The equation itself is also the Euler–Lagrange equation of a variational problem. Two natural classes of solutions can be associated with this equation. A conservative solution will preserve its energy in time, while a dissipative weak solution loses energy at the time when singularities appear. Conservative solutions are globally defined, forward and backward in time, and preserve interesting geometric features, such as the Hamiltonian structure. On the other hand, dissipative solutions appear to be more natural from the physical point of view.We establish the well-posedness of the Cauchy problem within the class of conservative solutions, for initial data having finite energy and assuming that the flux function f has a Lipschitz continuous second-order derivative. In the case where f is convex, the Cauchy problem is well posed also within the class of dissipative solutions. However, when f is not convex, we show that the dissipative solutions do not depend continuously on the initial data.     

The Evolution of Stress Intensity Factors and the Propagation of Cracks in Elastic Media

When a crack Γ _s propagates in an elastic medium the stress intensity factors evolve with the tip x(s) of Γ _s . In this paper we derive formulae which describe the evolution of these stress intensity factors for a homogeneous isotropic elastic medium under plane strain conditions. Denoting by ψ=ψ(x,s) the stress potential (ψ is biharmonic and has zero traction along the crack Γ _s ) and by κ(s) the curvature of the crack at the tip x(s), we prove that the stress intensity factors A ₁(s), A ₂(s), as functions of s, satisfy:

On the structure of continua and the mathematical properties of algebraic elastodynamic of a triclinic structural system

This paper is neither laudatory nor derogatory but it simply contrasts with what might be called elastostatic (or static topology), a proposition of the famous six equations. The extension strains and the shearing strains which were derived by A.L. Cauchy, are linearly expressed in terms of nine partial derivatives of the displacement function (u _i ,u _j ,u _h )=u(x ⁱ ,x ^j ,x ^k ) and it is impossible for the inverse proposition to sep up a system of the above six equations in expressing the nine components of matrix (∂(u _i ,u _j ,u _h )/∂(x ⁱ ,x ^j ,x ^k )). This is due to the fact that our geometrical representations of deformation at a given point are as yet incomplete^[1]. On the other hand, in more geometrical language this theorem is not true to any triangle, except orthogonal, for “squared length” in space^[2]. The purpose of this paper is to describe some mathematic laws of algebraic elastodynamics and the relationships between the above-mentioned important questions.     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

Third-order systems: Periodicity conditions

Recently a third-order existence theorem has been proven to establish the sufficient conditions of periodicity for the most general third-order ordinary differential equation

x+f(t,x,x^′,x^″)=0