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.     

Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

Further Study on Pure Shear

It is known that the Cauchy stress tensor T is a pure shear when trT = 0. An elementary derivation is given for a coordinate system such that, when referred to this coordinate system, the diagonal elements of T vanish while the off-diagonal elements τ ₁, τ ₂, τ ₃, are the pure shears. The structure of τ _i (i = 1, 2, 3) depends on one non-dimensional parameter q = 54(detT)² / [tr(T ²)]³, 0 ≤ q ≤ 1. When q = 0, one of the three τ _i vanishes. A coordinate system can be chosen such that the remaining two have the same magnitude or one of the remaining two also vanishes. When q = 1, all three τ _i have the same magnitude. However, there is a one-parameter family of coordinate systems that gives the same three τ _i . For q ≠ 0 or 1, none of the three τ _i vanishes and the three τ _i in general have different magnitudes. Nevertheless, a coordinate system can be chosen such that two of the three τ _i have the same magnitude. ^★Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

Local and global behavior of solutions of quasilinear equations of Emden-Fowler type

We consider the equation div for p≦N, 0<p−1<q. We study the isolated singularities and the behavior near infinity of nonradial positive solutions when q <N(p −1)/(N − p), and give a complete classification of local and global radial solutions of any sign, for any q.     

Surface action of a heat pulse in a tangential flow

The three-dimensional unsteady problem of the action of a localized radiation pulse on an aluminum plate of finite thickness in a tangential airflow is solved on the basis of the Euler equations in the gas phase and the heat conduction equation in the solid. The process is accompanied by the intense evaporation of plate material and the subsequent (after the action of the radiation pulse ceases) intense condensation of aluminum vapor and condensate particles on the surface. The problem is solved on the assumption that the surface is not screened by its disintegration products and the processes associated with the self-radiation of the gas are disregarded. Heat pulses of intensityq=10⁶−10⁷ W/cm² and durationt _i ≈10⁻⁴ sec are investigated. Various regimes are calculated and the effect of the governing parameters on the impulse transmitted to the surface and ablation of the target material is established. A nontraditional approach, based on the idea of statistical modeling of the motion of an inviscid non-heat-conducting gas at the kinetic level, is used for solving the gas dynamic problem. The application of the approach in question to the solution of problems of this class is justified. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.1, pp. 139–146, January–February, 1994.     

Range-pair-counting fatigue gage

In this paper, the surface roughening phenomenon of an aluminum foil during fatigue process is utilized for monitoring of a variable amplitude stress. It is found that constant mean stress does not affect surface roughness. When the mean stress varies, however, it causes surface roughening. In such a case, the equivalent stress defined by (Σσ _i ^α N _i /ΣN _i )^1/α dominates surface roughness of an aluminum foil, where σ _i andN _i are stress amplitude and the number of cycles counted by the rangepair method, respectively, and α is a gage factor which is peculiar to the foil used. It is concluded that the aluminum foil can be used as the range-pair counting fatigue gage with high accuracy.     

Fatigue gage utilizing slip-initiation phenomenon in electroplated copper foil

Fundamental experiments are carried out to examine the parameter that dominates the slip-band initiation in electroplated copper foil under the condition where the mean stress as well as the stress amplitude varies. In the case of constant-amplitude stressing, the relation between the critical stress for the slip-band initiation σ _p and the number of cyclesN is well represented by σ _p ^α N=C. In other words, the slip bands appear when the total hysteresis energy applied to the copper foil attains a critical value. In the case of variable stresses, the range-pair mainly dominates the occurrence of the slip bands, and Miner's linear cumulative damage rule holds for the accumulation of the fatigue damage for the slip-band initiation. Accordingly, the parameter (Σσ _i ^α n _i/Σn _i)^1/α is equivalent to the critical stress σ _p under constant amplitude stressing, where σ _i andn _i are the stress amplitude and the number of cycles counted by the range-pair method, respectively, and α is the exponent of the σ _p -N relation. Based on these results, the applicability of the copper foil to the fatigue gage that accumulates and indicates a load experience is discussed.     

Analytical and numerical analysis for the FGM thick sphere under combined pressure and temperature loading

Thermo-mechanical analysis of functionally graded hollow sphere subjected to mechanical loads and one-dimensional steady-state thermal stresses is carried out in this study. The material properties are assumed to vary non-linearly in the radial direction, and the Poisson’s ratio is assumed constant. The temperature distribution is assumed to be a function of radius, with general thermal and mechanical boundary conditions on the inside and outside surfaces of the sphere. In the analysis presented here, the effect of non-homogeneity in FGM thick sphere was implemented by choosing a dimensionless parameter, named β _i (i = 1, . . . , 3), which could be assigned an arbitrary value affecting the stresses in the sphere. It is observed that the solution process for β _i (i = 3) = −1 are different from those obtained for other values of β _i (i = 1, . . . , 3). Cases of pressure, temperature, and combined loadings were considered separately. It is concluded that by changing the value of β _i (i = 1 . . . 3), the properties of FGM can be so modified that the lowest stress levels are reached. The present results agree well with existing results. Using FEM simulations, the analytical findings for FGM spheres under the influence of internal pressure and temperature gradient were compared to finite element results.     

Detonability of THDCPD-exo–air mixtures

In the frame of industrial risk and propulsive application, the detonability study of JP10–air mixtures was performed. The simulation and measurements of detonation parameters were performed for THDCPD-exo/air mixtures at various initial pressure (1 bar < P ₀ < 3 bar) and equivalence ratio (0.8 < Φ < 1.6) in a heated tube (T ₀ ~ 375 K). Numerical simulations of the detonation were performed with the STANJAN code and a detailed kinetic scheme of the combustion of THDCPD. The experimental study deals with the measurements of detonation velocity and cell size λ. The measured velocity is in a good agreement with the calculated theoretical values. The cell size measurements show a minimum value for Φ ~ 1.2 at every level of initial pressure studied and the calculated induction length L _i corresponds to cell size value with a coefficient k = λ/L _i = 24 at P ₀ = 1 bar. Based on the comparison between the results obtained during this study and those available in the literature on the critical initiation energy E _c, critical tube diameter d _c and deflagration to detonation transition length L _DDT, we can conclude that the detonability of THDCPD–air mixtures corresponds to that of hydrocarbon–air mixtures.

---

This paper is based on the work presented at the 33rd International Pyrotechnics Seminar, IPS 2006, Fort Collins, July 16–21, 2006.     

Predictions of first-ply failure load of laminated composite shafts

The failure strengths of laminated composite shafts subjected to static bending load and/or torque are studied using both theoretical and experimental approaches. We investigate optimal angle-ply orientations of laminated composite shafts with the design of symmetric main fiber lay-ups [θ/−θ/...] _S and anti-symmetric reinforcing fiber lay-ups [α/−α/...]_−S for maximum stiffness. Different analytical methods, together with various failure criteria, are used to predict the first-ply failure strengths of the composite shafts. A hand lay-up and pressure bag technique was used to manufacture the composite shafts, which includes the stacking of pieces of prepreg in the desired orientations and the curing process. The first-ply failure loads of laminated composite shafts were determined using acoustic emission. The failure modes of the shafts are studied and experimental results are used to verify the theoretical predictions. The experimental results have been proven to be efficient and effective in the theoretical prediction of first-ply failure loads of laminated composite shafts.     

Effects of couple stresses on the oscillatory flow past an infinite plate with constant suction

Summary An analysis of a two dimensional oscillatory flow past an infinite porous plate with contant suction is carried out on taking into account the couple stresses. Here the free stream velocity oscillates about a nonzero constant mean. Approximate solutions are derived to coupled linear equations, and the expressions for the mean velocity, the transient velocity, the mean skinfriction, the amplitude and the phase of skin-friction are obtained. The solutions are followed by discussion. the effects of variations of α(ν_r/ν), β(Iν/γ) and λ, the frequency are graphically represented and physically interpreted. It is observed that the reverse type of flow does not occur in the presence of the couple stresses.

---

Sommario In questo lavoro è svolta un'analisi di un flusso oscillatorio bidimensionale sopra una piastra porosa, infinita, con aspirazione costante, tenendo conto delle coppie di sforzo. La velocità della corrente libera oscilla attorno ad un valore medio costante diverso da zero. Si deducono le soluzioni approssimate per le equazioni lineari accoppiate e si ottengono le espressioni per la velocità media, la velocità transitoria, l'attrito superficiale, l'ampiezza e la fase. Si discutono le soluzioni. Si rappresentano graficamente e si interpretano fisicamente gli effetti delle variazioni di α(ν_r/ν), β(Iν/γ) e λ. Si osserva che in presenza delle coppie di sforzo nel fluido non si ha il tipo inverso di flusso.

     

The libration points in photogravitational restricted three-body problem

ＴＨＥＬＩＢＲＡＴＩＯＮＰＯＩＮＴＳＩＮＰＨＯＴＯＧＲＡＶＩＴＡＴＩＯＮＡＬＲＥＳＴＲＩＣＴＥＤＴＨＲＥＥ－ＢＯＤＹＰＲＯＢＬＥＭＺｈｅｎｇＸｕｅ－ｔａｎｇ（郑学塘）ＹｕＬｉ－ｚｈｏｎｇ（郁丽忠）（Ｄｅｐｔ．ｏｆＡｐｐｌｉｅｄＰｈｙｓｉｃｓ,Ｎａｎｊ...     

Crack tip field in functionally gradient material with exponential variation of elastic constants in two directions

This paper presents an exact solution of the crack tip field in functionally gradient material with exponential variation of elastic constants. The dimensionless Poisson's ratios ν ₀ of the engineering materials (iron, glass . . .) are far less than one; therefore, neglecting them, one can simplify the basic equation and the exact solution is easy to obtain. Although the exact solution for the case ν ₀≠0 is also obtained, it is very complicated and the main result is the same with the case ν ₀=0 (it will be dealt with in Appendix VII). It has been found that the exponential term exp (ax+by) in the constitutive equations becomes exp (ax/2+by/2−kr/2) in the exact solution. The English text was polished by Keren Wang.     

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.     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Billiards and Two-Dimensional Problems of Optimal Resistance

A body moves in a medium composed of noninteracting point particles; the interaction of the particles with the body is completely elastic. The problem is: find the body’s shape that minimizes or maximizes resistance of the medium to its motion. This is the general setting of the optimal resistance problem going back to Newton. Here, we restrict ourselves to the two-dimensional problems for rotating (generally non-convex) bodies. The main results of the paper are the following. First, to any compact connected set with piecewise smooth boundary B ì \mathbbR²{B \subset \mathbb{R}^2} we assign a measure ν _B on ∂(conv B)×[ − π/2, π/2] generated by the billiard in \mathbbR² \B{\mathbb{R}^2 \setminus B} and characterize the set of measures {ν _B }. Second, using this characterization, we solve various problems of minimal and maximal resistance of rotating bodies by reducing them to special Monge–Kantorovich problems.     

The Damped-Driven 2D Navier–Stokes System on Large Elongated Domains

The Navier–Stokes system with damping, which is motivated by Stommel–Charney model of ocean circulation, is considered in a large elongated periodic rectangular domain with area of the order α⁻¹, as α → 0. We obtain estimates for the dimension of the global attractor that are sharp as both α → 0 and ν → 0, where ν is the viscosity coefficient. This work was supported in part by the US Civilian Research and Development Foundation, grant no. RUM1-2654-MO-05 (A.A.I. and E.S.T.). The work of A.A.I. was supported in part by the Russian Foundation for Fundamental Research, grants no. 06-001-0096 and no. 05-01-429, and by the RAS Programme no. 1 ‘Modern problems of theoretical mathematics’. The work of E.S.T. was supported in part by the NSF, grant no. DMS-0204794, the MAOF Fellowship of the Israeli Council of Higher Education, and by the BSF, grant no. 200423.     

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.     

Aerodynamic forces and flow fields of a two-dimensional hovering wing

This paper reports the results of an experimental investigation on a two-dimensional (2-D) wing undergoing symmetric simple harmonic flapping motion. The purpose of this investigation is to study how flapping frequency (or Reynolds number) and angular amplitude affect aerodynamic force generation and the associated flow field during flapping for Reynolds number (Re) ranging from 663 to 2652, and angular amplitudes (α _A) of 30°, 45° and 60°. Our results support the findings of earlier studies that fluid inertia and leading edge vortices play dominant roles in the generation of aerodynamic forces. More importantly, time-resolved force coefficients during flapping are found to be more sensitive to changes in α _A than in Re. In fact, a subtle change in α _A may lead to considerable changes in the lift and drag coefficients, and there appears to be an optimal mean lift coefficient around α _A = 45°, at least for the range of flow parameters considered here. This optimal condition coincides with the development a reverse Karman Vortex street in the wake, which has a higher jet stream than a vortex dipole at α _A = 30° and a neutral wake structure at α _A = 60°. Although Re has less effect on temporal force coefficients and the associated wake structures, increasing Re tends to equalize mean lift coefficients (and also mean drag coefficients) during downstroke and upstroke, thus suggesting an increasing symmetry in the mean force generation between these strokes. Although the current study deals with a 2-D hovering motion only, the unique force characteristics observed here, particularly their strong dependence on α _A, may also occur in a three-dimensional hovering motion, and flying insects may well have taken advantage of these characteristics to help them to stay aloft and maneuver. An erratum to this article can be found at     

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