Existence of the minimal positive solution of some nonlinear elliptic systems when the nonlinearity is the sum of a sublinear and a superlinear term

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｉｓｐａｐｅｒ,ｗｅｃｏｎｓｉｄｅｒｔｈｅｅｌｌｉｐｔｉｃｓｙｓｔｅｍ(1λ)　-Δｕ=ｆ(λ,ｘ,ｕ)-ｖ　　(ｉｎΩ),-Δｖ=δｕ-γｖ(ｉｎΩ),ｕ=ｖ=0(ｏｎΩ),ｗｈｅｒｅΩｉｓａｓｍｏｏｔｈｂｏｕｎｄｅｄｄｏｍａｉｎｉｎＲＮ(Ｎ≥2)ａｎｄλｉｓａｒｅａｌｐａｒａｍｅｔｅｒ.Ｔｈｅｓｏｌｕｔｉｏｎｓ(ｕ,ｖ)ｏｆｔｈｉｓｓｙｓｔｅｍｒｅｐｒｅｓｅｎｔｓｔｅａｄｙｓｔａｔｅｓｏｌｕｔｉｏｎｓｏｆｒｅａｃｔｉｏｎｄｉｆｆｕｓｉｏｎｓｙｓｔｅｍｓｄｅｒｉｖｅｄｆｒｏｍｓｅｖｅｒａｌａｐ…     

Transmission characteristics of EM wave in a finite thickness plasma

One of the key factors for solving the problems of re-entry communication interruption is electromagnetic(EM) wave transmission characteristics in a plasma.Theoretical and experimental studies were carried out on specific transmission characteristics for different plasma sheath characteristic under thin sheath condition in re-entry state.The paper presents systematic studies on the variations of wave attenuation characteristics versus plasma sheath thickness L,collision frequency ν,electron density n e and wave working frequency f in a φ 800 mm high temperature shock tube.In experiments,L is set to 4 cm and 38 cm.ν is 2 GHz and 15 GHz.n e is from 1×10 10 cm(-3) to 1×10 13 cm(-3),and f is set to 2,5,10,14.6 GHz,respectively.Meanwhile,Wentzel-Kramers-Brillouin(WKB) and finite-difference time-domain(FDTD) methods are adopted to carry out theoretical simulation for comparison with experimental results.It is found that when L is much larger than EM wavelength λ(thick sheath) and ν is large,the theoretical result is in good agreement with experimental one,when sheath thickness L is much larger than λ,while ν is relatively small,two theoretical results are obviously different from the experimental ones.It means that the existing theoretical model can not fully describe the contribution of ν.Furthermore,when L and λ are of the same order of magnitude(thin sheath),the experimental result is much smaller than the theoretical values,which indicates that the current model can not properly describe the thin sheath effect on EM attenuation characteristics.     

A Study on Elliptical Gas Flow in Tri-porosity Gas Reservoirs

Elliptical flow is common in the near vertical fracture area and in anisotropic reservoirs. However, the classical radial flow models cannot provide a complete analysis for elliptical flow. This article presents a new mathematic model for gas elliptical flow in tri-porosity gas reservoirs. The differential equation of the new model is written in Mathieu equation, so that the solution can also be expressed by Mathieu functions. The numerical solution of the corresponding Mathieu functions ce _2n (ξ, −q), Ke _2n (ξ, −q) and their derivatives are obtained to derive the dimensionless pseudo pressure drop in Laplace space. The sensitivities of tri-porosity systems, including the parameters related to anisotropies C _De^2S and ξ _w, the storativity ratios ω _f and ω _m, and the interporosity flow coefficients λ_vf and λ_mf, are studied using Laplace numerical inversion. The new solution includes not only the factors considered in classic solutions in previous articles, but also incorporates the effect of reservoir anisotropy.     

Wave Packets, Signaling and Resonances in a Homogeneous Waveguide

We solve the initial-boundary-value linear stability problem for small localised disturbances in a homogeneous elastic waveguide formally by applying a combined Laplace – Fourier transform. An asymptotic evaluation of the solution, expressed as an inverse Laplace – Fourier integral, is carried out by means of the mathematical formalism of absolute and convective instabilities. Wave packets, triggered by perturbations localised in space and finite in time, as well as responses to sources localised in space, with the time dependence satisfying e^−iωt + O(e^−ɛt ), for t → ∞, where Im ω₀ = 0 and ω > 0 , that is, the signaling problem, are treated. For this purpose, we analyse the dispersion relation of the problem analytically, and by solving numerically the eigenvalue stability problem. It is shown that due to double roots in a wavenumber k of the dispersion relation function D(k, ω), for real frequencies ω, that satisfy a collision criterion, wave packets with an algebraic temporal decay and signaling with an algebraic temporal growth, that is, temporal resonances, are present in a neutrally stable homogeneous waveguide. Moreover, for any admissible combination of the physical parameters, a homogeneous waveguide possesses a countable set of temporally resonant frequencies. Consequences of these results for modelling in seismology are discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Vorticities in a LES Model for 3D Periodic Turbulent Flows

This paper is devoted to the study of a LES model to simulate turbulent 3D periodic flow. We focus our attention on the vorticity equation derived from this LES model for small values of the numerical grid size δ. We obtain entropy inequalities for the sequence of corresponding vorticities and corresponding pressures independent of δ, provided the initial velocity u₀ is in L_x² while the initial vorticity ω₀ = ∇ × u₀ is in L_x¹. When δ tends to zero, we show convergence, in a distributional sense, of the corresponding equations for the vorticities to the classical 3D equation for the vorticity.     

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997     

Conversion of longitudinal waves into electromagnetic radiation in a slowly varying plasma under W.K.B. approximation

This paper investigates the conversion of a dispersive longitudinal oscillation into reflected and transmitted electromagnetic radiation fields in slowly varying unmagnetized warm fluid plasmas, using W.K.B. approximations. The expressions for the power of the transmitted and reflected electromagnetic radiations, generated by electron acoustic waves, have also been obtained. It is shown that this conversion process becomes most efficient under certain conditions.

Nomenclature

In § 2 H magnetic field - H ₁ - u electron fluid velocity - k _t wave number of the transverse wave - k ₁ wave number of the longitudinal wave in electron fluid - m electronic mass - N ₀ number density of electrons in the unperturbed state - N perturbation in the electron number density - p perturbation in the electron fluid pressure - v _e adiabatic sound velocity of the electron fluid - K _t ² c ^{2 2}–e² - K ₁ ² v _e ^{2 2}–e² - wave frequency - _e electron plasma frequency - 1– _e ² / ² - c velocity of light in vacuum In § 3 K ₀ wave number in the 0X direction - K ₁ ² K ₁ ² –K ₀ ² - K ₂ ² K _t ² –K ₀ ² - K ₃ K ₁–K ₂ - K ₄ K ₁+K ₂ - K ₅ (K ₁ K ₂)^1/2 See Appendix A - A ₁ pressure amplitude of the reflected part of the incident wave - B ₁ pressure amplitude of the transmitted part of the incident wave - L characteristic length of variation ofN ₀ - e _x unit vector along 0X - e _z unit vector along 0Z In § 4 S _t Poynting flux of the transverse electromagnetic radiation - S _tZ ^/t Average of the transmitted part of the poynting flux along 0Z over the time period 2/ - S _tZ ^/r Average of the reflected part of the poynting flux along 0Z over the time period 2/ In § 5 S ₁ Energy flux carried by the longitudinal pressure wave - S _1Z ^/t Average of the transmitted part ofS ₁ along 0Z over the time period 2/     

Experiments on the lift of a spinning sphere in a range of intermediate Reynolds numbers

 The lift force experienced by a spinning sphere moving in a viscous fluid, with constant linear and angular velocities, is measured by means of a trajectographic technique. Measurements are performed in the range of dimensionless angular velocities γ=aω/V lying between 1 and 6, and in the range of Reynolds numbers Re=2aV/ν lying between 10 and 140 (a sphere radius, ω angular velocity, V relative velocity of the sphere centre, ν fluid kinematic viscosity). A notable departure from the theoretical relationship at low Reynolds number, C _L =2γ, is obtained, the ratio C _L /γ being found to significantly decrease with increasing γ and increasing Re. The following correlation is finally proposed to estimate the lift coefficient in the range 10<Re<140: C _L ≅0.45+(2γ−0.45) exp (−0.075γ^0.4 Re ^0.7) Received: 20 May 1996/Accepted: 9 November 1997     

Exact analytic solutions for free convection boundary layers on a heated vertical plate with lateral mass flux embedded in a saturated porous medium

 The problem of the self-similar boundary flow of a “Darcy-Boussinesq fluid” on a vertical plate with temperature distribution T _w(x) = T _∞+A·x ^λ and lateral mass flux v _w(x) = a·x ^(λ−1)/2, embedded in a saturated porous medium is revisited. For the parameter values λ = 1,−1/3 and −1/2 exact analytic solutions are written down and the characteristics of the corresponding boundary layers are discussed as functions of the suction/ injection parameter in detail. The results are compared with the numerical findings of previous authors. Received on 8 March 1999     

A Hierarchy of Plate Models Derived from Nonlinear Elasticity by Gamma-Convergence

We derive a hierarchy of plate models from three-dimensional nonlinear elasticity by Γ-convergence. What distinguishes the different limit models is the scaling of the elastic energy per unit volume ∼h^β, where h is the thickness of the plate. This is in turn related to the strength of the applied force ∼h^α. Membrane theory, derived earlier by Le Dret and Raoult, corresponds to α=β=0, nonlinear bending theory to α=β=2, von Kármán theory to α=3, β=4 and linearized vK theory to α>3. Intermediate values of α lead to certain theories with constraints. A key ingredient in the proof is a generalization to higher derivatives of our rigidity result [29] which states that for maps v:(0,1)³→ℝ³, the L² distance of ∇v from a single rotation is bounded by a multiple of the L² distance from the set SO(3) of all rotations.     

Optimization of the deflagration to detonation transition: reduction of length and time of transition

The aim of this experimental investigation is the study of Deflagration to Detonation Transition (DDT) in tubes in order to (i) reduce both run-up distance and time of transition (L _DDT and t _DDT) in connection with Pulsed Detonation Engine applications and to (ii) attempt to scale L _DDT with λ_CJ (the detonation cellular structure width). In DDT, the production of turbulence during the long flame run-up can lead to L _DDT values of several meters. To shorten L _DDT, an experimental set-up is designed to quickly induce highly turbulent initial flow. It consists of a double chamber terminated with a perforated plate of high Blockage Ratio (BR) positioned at the beginning of a 26 mm inner diameter tube containing a “Shchelkin spiral” of BR ≈ 0.5. The study involves stoichiometric reactive mixtures of H₂, CH₄, C₃H₈, and C₂H₄ with oxygen and diluted with N₂ in order to obtain the same cell width λ_CJ≈10 mm at standard conditions. The results show that a shock-flame system propagating with nearly the isobaric speed of sound of combustion products, called the choking regime, is rapidly obtained. This experimental set-up allows a L _DDT below 40 cm for the mixtures used and a ratio L _DDT/λ_CJ ranging from 23 to 37. The transition distance seems to depend on the reduced activation energy (E _a/RT _c) and on the normalized heat of reaction (Q/a ₀ ²). The higher these quantities are, the shorter the ratio L _DDT/λ_CJ is. PACS 47.40.Rs · 47.60.+i · 47.70.Pq · 47.80.CbThis paper was based on the work that was presented at the 19th International Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27–August 1, 2003.     

The Decay Rate and Higher Approximation of Mild Solutions to the Navier�CStokes Equations

In this paper, we consider v(t) = u(t) − e ^tΔ u ₀, where u(t) is the mild solution of the Navier–Stokes equations with the initial data u₀ ? L²(\mathbb Rⁿ)?Lⁿ(\mathbb Rⁿ){u_0\in L^2({\mathbb R}^n)\cap L^n({\mathbb R}^n)} . We shall show that the L ² norm of D ^β v(t) decays like t^{-\frac |b|-1 2-\frac n4}{t^{-\frac {|\beta|-1} {2}-\frac n4}} for |β| ≥ 0. Moreover, we will find the asymptotic profile u ₁(t) such that the L ² norm of D ^β (v(t) − u ₁(t)) decays faster for 3 ≤ n ≤ 5 and |β| ≥ 0. Besides, higher-order asymptotics of v(t) are deduced under some assumptions.     

Reynolds Number Dependence of Velocity Structure Functions in a Turbulent Pipe Flow

Low-order moments of the increments δu andδv where u and v are the axial and radial velocity fluctuations respectively, have been obtained using single and X-hot wires mainly on the axis of a fully developed pipe flow for different values of the Taylor microscale Reynolds numberR _λ. The mean energy dissipation rate〉ε〈 was inferred from the uspectrum after the latter was corrected for the spatial resolution of the hot-wire probes. The corrected Kolmogorov-normalized second-order structure functions show a continuous evolution withR _λ. In particular, the scaling exponentζ _v , corresponding to the v structure function, continues to increase with R _λ in contrast to the nearly unchanged value of ζ _u . The Kolmogorov constant for δu shows a smaller rate of increase with R _λ than that forδv. The level of agreement with local isotropy is examined in the context of the competing influences ofR _λ and the mean shear. There is close but not perfect agreement between the present results on the pipe axis and those on the centreline of a fully developed channel flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Elasticity and dynamics of particle gels in non-Newtonian melts

We investigate the relation between the structure and the viscoelastic behavior of a model polymer nanocomposite system based on a mixture of titanium dioxide (TiO₂) nanoparticles and polypropylene. Above a critical volume fraction, Φ _c, the elasticity of the hybrids dramatically increases, and the frequency dependence of the elastic and viscous moduli reflects the superposition of the independent responses of the suspending polymer melt and of an elastic particle network. In addition, the elasticity of the hybrids shows critical behavior around Φ _c. We interpret these observations by hypothesizing the formation of a transient network, which forms due to crowding of particle clusters. Consistent with this interpretation, we find a long-time, Φ-dependent, structural relaxation, which emphasizes the transient character of the structure formed by the particle clusters. For times below this characteristic relaxation time, the elasticity of the network is Φ-independent and reminiscent of glassy behavior, with the elastic modulus, G^′, scaling with frequency, ω, as G^′∼ω ^0.3. We expect that our analysis will be useful for understanding the behavior of other complex fluids where the elasticity of the components could be superimposed.     

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     

Dynamic Contact Behavior of a Golf Ball during an Oblique Impact

The oblique impact between a golf ball and a rigid steel target was studied using a high-speed video camera. Video images recorded before and after the impact were used to determine the inbound velocity v _i, rebound velocity v _r, inbound angle θ_i, rebound angle θ_r, and the coefficient of restitution e. The results showed that θ_r and e decreased as v _i increased. The maximum compression ratio η_c, contact time t _c, average angular velocity , and tangential velocity , along the target were determined from images obtained during the impact. The images demonstrated that η_c increased with v _i while t _c decreased. In addition, and increased almost linearly as v _i increased. A rigid body model was used to estimate the final angular velocity ω* and tangential velocity ν_t* at the end of the impact; these results were then compared with experimental data.     

Single bubble deformation and breakup in simple shear flow

Experiments in a parallel band apparatus and a transparent concentric cylinder device allow the observation of bubble deformation (shape and orientation) and breakup as a function of the viscosity ratio λ and the Capillary number Ca. For viscosity ratios between 3.1 × 10⁻⁷ and 6.7 × 10⁻⁸, critical Capillary numbers Ca _c for bubble breakup between 29 and 45 are found. It is furthermore shown that in the given parameter space no clear distinction between tip breakup and fracture can be made for bubbles. An erratum to this article can be found at     

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     

Dependence of the zero shear-rate viscosity and the viscosity function of linear high-density polyethylenes on the mass-average molar mass and polydispersity

Linear high-density polyethylenes with molar masses M _w between 240 and 1,000,000 g/mol, obtained by metallocene catalysts, were characterized in shear using oscillatory and creep tests. The polydispersities of the molar mass distributions (MMDs) lay between 1 and 16. The resulting zero shear-rate viscosities η₀ covered a range from 2.5×10⁻³ to around 10⁸ Pas. Above a critical molar mass of M _c≈2,900 g/mol, the experimental results can be described by the relation η₀ ∼ M _w^3.6, independently of the MMD. The oscillatory data were fitted with a Carreau–Yasuda equation. The resulting parameters were correlated to molecular structure. The parameter a, being a quantity for the width of the transition between the Newtonian and the non-Newtonian regime, showed a dependence on the molar mass M _w but not on M _w/M _n. The parameter λ of the Carreau-Yasuda equation was found to be the reciprocal crossover frequency for all samples with a log-Gaussian MMD. λ depends on the molar mass M _w and also on M _w/M _n.

Mach reflection of a plane shock wave passing a mountain of 45° inclination

The transition from regular reflection (RR) to Mach reflection (MR) as a plane shock wave diffracts around a triangular mountain of 45° inclination is analysed in this paper, both by optical measurement in a shock tube and by numerical simulation the numerical method developed by Li Yingfan^[1] is of the FLIC type with triangular mesh. The dependence of the critical transition point Lk ofRR→MR on shock Mach numberM _i is analyzed and the variations of the incidence angle ω _i of the impinging shock and the reflection angle ω _r with the distanceL ^* are investigated. Our experimental and numerical results agree well with the theoretical results of Iton and Italya.