On the nonlinear slewing dynamics and control of the Space Station based Mobile Servicing System

A relatively general Lagrangian formulation for studying the nonlinear dynamics and control of space-craft with interconnected flexible members in a tree-type topology is developed. Versatility of the formulation is illustrated through a dynamical study of the Space Station based two-link Mobile Servicing System (MSS). The performance of the MSS undergoing inplane and out-of-plane slewing maneuvers is compared. Results indicate that, in absence of control, the maneuvers induce undesirable librational motion of the Space Station as well as vibration of the links. Nonlinear control, based on the Feedback Linearization Technique (FLT), appears promising. Quasi-Closed Loop Control (QCLC), a variation of the FLT, is applied to control the libration of the Space Station. Once the attitude of the Space Station is controlled, the performance of the MSS improves significantly. For a 5-minute maneuver of the MSS, the maximum control torque required is only 34.5 Nm.Nomenclature f _i ¹ , f _i,j ¹ fundamental frequency of bodies B _i and B _i,j, respectively - l _c, l _i, l _i,j length of bodies B _c, B _i, and B _i,j, respectively - m _c, m _i, m _i,j mass of bodies B _c, B _i, and B _i,j, respectively - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbhv2BYDwAHbacfiGab8xCayaaraqefavySfgDP52BGWuAU9gD% 5bxzaGGbciaa+zgacaWFSaGaa8hiaiqa-fhagaqeaiaa-jhaaaa!4B1F!\[\bar qf, \bar qr\] vector representing flexible and rigid generalized coordinates - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbhv2BYDwAHbacfiGaa8hkaiqa-fhagaqeaiaa-jhacaWFPaqe% favySfgDP52BGWuAU9gD5bxzaGGbciaa+rgaaaa!4A18!\[(\bar qr)d\] vector representing the desired rigid generalized coordinates - (I _xx)k, (I _yy)k, (I _zz)k principal inertia of body B _k about X _k, Y _k, and Z _k axes, respectively; ksc, i or i, j - K _p, K _v displacement and velocity gain matrices - N _q total number of generalized coordinates - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbwvMCKfMBHbacfiGab8xuayaaraqefavySfgDP52BGWuAU9gD% 5bxzaGGbciaa+zgaieaacaqFSaGaa0hiaiqa-ffagaqeaGqaciaa8j% haaaa!4AEF!\[\bar Qf, \bar Qr\] control effort vectors for flexible and rigid coordinates, respectively - Q , Q , Q control effort for pitch, roll and yaw degree of freedom, respectively - _k ^y , _k ^z tip deflection of a beam type appendage (B _k) in the Y _k and Z _k directions, respectively.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

[0, k i]1 m -Factorizations Orthogonal to a Subgraph

Let G be a graph,k₁,…,k_m be positive integers. If the edges of graph G can be decomposed into some edge disjoint. [0,k₁]-factor. F₁,…,[0,k_m]-factor F_m, then we can say F={F₁,…,F_m, is a [0,k_i]_m¹-factorization of G. If H is a subgraph with m edges in graph G and |E(H)∩E(F_i)|=1 for all 1≤i≤m, then we can call that F is orthogonal to H. It is proved that if G is a..[0,k₁+…+k_m-m+1]-graph, H is a subgraph with m edges in G, then graph G has a. [0,k_i]₁^m-factorization orthogonal to H.     

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.     

Effect of hydrodynamic forces on the pressure-difference equation

Because of the influence of hydrodynamic forces, the difference in macroscopic pressure which exists, at static equilibrium, between two immiscible phases located in a porous medium may be different from that which pertains during flow. In this paper, the concept of relative pressure difference, together with a new pressure-difference equation, is used to investigate the impact that the hydrodynamic forces have on the difference in macroscopic pressure which pertains when two immiscible fluids flow simultaneously through a homogeneous, water-wet porous medium. This investigation reveals that, in general, the equation defining the difference in pressure between two flowing phases must include a term which takes proper account of the hydrodynamic effects. Moreover, it is pointed out that, while neglect of the hydrodynamic effects introduces only a small amount of error when the two fluids are flowing cocurrently, such neglect is not permissible during steady-state, countercurrent flow. This is because failure to include the impact of the hydrodynamic effects in the latter case makes it impossible to explain the pressure behaviour observed in steady-state, countercurrent flow. Finally, the results of this investigation are used as a basis for arguing that, during steady-state, countercurrent flow, saturation is uniform, as is the case of steady-state, cocurrent flow.Roman Letters a parameter in Equation (18) - k absolute permeability, m² - k _i effective permeability to phasei;i=1, 2, m² - k _ij generalized effective permeability for phasei;i, j=1, 2, m² - p _d p ₂–p ₁=difference in macroscopic pressure between two flowing phases, N/m² - p _i pressure for phasei;i=1, 2, N/m² - p _h hydrodynamic contribution to difference in macroscopic pressure which exists during flow, N/m² - P _c macroscopic static capillary pressure, N/m² - R ₁₂ function defined by Equation (18) - S _i saturation of phasei;i=1, 2 - S _n normalized saturation of phase 1 - t time, s - u _i flux of phasei;i=1, 2m³/m²/s - x distance in direction of flow, m Greek Letters _R relative pressure difference - _i k _i / _i =mobility of phasei;i=1, 2m²/Pa·s - _ij k _ij / _j =generalized mobility of phasei;i, j=1, 2m²/Pa·s - _i viscosity of phasei;i=1, 2, Pa·s - porosity     

Crack repair using an elastic filler

The effect of repairing a crack in an elastic body using an elastic filler is examined in terms of the stress intensity levels generated at the crack tip. The effect of the filler is to change the stress field singularity from order 1/r^1/2 to 1/r^(1-λ) where r is the distance from the crack tip, and λ is the solution to a simple transcendental equation. The singularity power (1-λ) varies from (the unfilled crack limit) to 1 (the fully repaired crack), depending primarily on the scaled shear modulus ratio γ_r defined by G₂/G₁=γ_rε, where 2πε is the (small) crack angle, and the indices (1, 2) refer to base and filler material properties, respectively. The fully repaired limit is effectively reached for γ_r≈10, so that fillers with surprisingly small shear modulus ratios can be effectively used to repair cracks. This fits in with observations in the mining industry, where materials with G₂/G₁ of the order of 10^-3 have been found to be effective for stabilizing the walls of tunnels. The results are also relevant for the repair of cracks in thin elastic sheets.     

Effects of the bifurcation angle on the steady flow structure in model saccular aneurysms

The effects of the bifurcation angle on the steady flow structure in a straight terminal aneurysm model with asymmetric outflow through the branches have been characterized quantitatively in terms of laser-Doppler velocimetry (LDV)-measured mean velocity and fluctuating intensity distributions. The bifurcation angles investigated were 60°, 90°, and 140° and the Reynolds number based on the bulk average velocity and diameter of the afferent vessel was 500. It is found that the size of the recirculating zones in the afferent vessel, the flow activity (both mean and fluctuating motions) inside the aneurysm, and the shear stresses acting on the aneurysmal wall increase with increasing bifurcation angle. More importantly, both LDV-measured and flow-visualized results of the present study suggest the presence of a critical bifurcation angle below which the aneurysm is susceptible to thrombosis, whereas above this the aneurysm is prone to progression or rupture.List of symbols a aneurysm height - b distance from orifice to fundus - c orifice diameter - D afferent conduit diameter - d fundus diameter - Hz frequency unit = cycle/second - L length of bifurcation zone - Re Reynolds number = U · D/v - U streamwise mean velocity - U _m streamwise bulk mean velocity - u streamwise fluctuating component - X ^* normalized streamwise coordinate: X ^* 0: X ^* = X/a; X ^* <0: X ^* = X/L - Y ^* normalized transverse coordinate: Y ^* = Y/D - Z ^* normalized spanwise coordinate: Z ^* = Z/D - kinematic viscosity - _b angle of bifurcation - _c critical bifurcation angle     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.     

Static Solutions of the Vlasov-Einstein System

Static solutions of the SO(3)-symmetric Vlasov-Einstein system are studied via a variational approach. For the constitutive relation of the Emden-Fowler type φ(E,F)≡E ^{σ+ 1} F ^k we prove the existence of such solutions of sufficiently small mass-energy, provided 0<σ < k+3/2. These solutions are local minimizers of the energy-Casimir functional, subjected to a variational barrier. Accepted July 16, 2000?Published online January 22, 2001     

Determination of the orthotropic plate parameters of stiffened plates and grillages in free vibration

Summary It is well known that the analysis of vibration of orthogonally stiffened rectangular plates and grillages may be simplified by replacing the actual structure by an orthotropic plate. This needs a suitable determination of the four elastic rigidity constants D _x, D _y, D _xy, D ₁ and the mass of the orthotropic plate. A method is developed here for determining these parameters in terms of the sectional properties of the original plate-stiffener combination or the system of interconnected beams. Results of experimental work conducted on aluminium plates agree well with the results of the theory developed here.Notation a, b span and width of stiffened plates, grillages and orthotropic plates - c, spacing of beams in the y and x directions - D rigidity of an unstiffened isotropic plate - D _x, D _y rigidities in the x and y directions of an orthotropic plate - D _xy torsional rigidity of an orthotropic plate - D ₁ a parameter associated with a poisson type ratio - e depth of the common neutral axis of the beam-plate combination below the middle surface of the plate - E, G modulus of elasticity and modulus of shear - H =D ₁+2D _xy - i, j integers - I, , I _i, I _j moment of inertia of beams - J, J _i, J _j polar moment of inertia of beams - m, n integers referring to mode numbers in the x and y directions - p _{m n} natural frequency of an orthotropic plate - r, s number of beams in the transverse and longitudinal direction - T _g, T _o, T _s kinetic energy of grillages, orthotropic plates and stiffened plates - u frequency parameter of grillages - V _g, V _o, V _s potential energy of grillages, orthotropic plates and stiffened plates - W=W(x, y, t) deflection of an orthotropic plate - w=w(x, y) amplitude of vibration of orthotropic plates, stiffened plates and grillages - x, y, x _i, y _i cartesian coordinates - X _m, Y _n beam eigen functions - X _mj, Y _ni X _mand Y _nat x=x _jand y=y _i - X _m, Y _n; X_m, Y_n first and second derivatives of X _mand Y _nrespectively - Y _mn plate eigen functions - _n, _n parameters occurring in the expression for the beam functions - , , _i, _j mass per unit length of beams - _mn = frequency parameter of beam and slab bridges - = - Poisson's ratio - , mass per unit area of unstiffened plates and orthotropic plates - _mn natural frequency of stiffened plates and grillages     

Some properties of the solutions of wave equations

Some properties of solutions of initial value problems and mixed initial-boundary value problems of a class of wave equations are discussed. Wave modes are defined and it is shown that for the given class of wave equations there is a one to one correspondence with the roots _i (k) or k _j () of the dispersion relation W(, k)=0. It is shown that solutions of initial value problems cannot consist of single wave modes if the initial values belong to W ₂ ¹ (–, ); generally such solutions must contain all possible modes. Similar results hold for solutions of mixed initial-boundary value problems. It is found that such solutions are stable, even if some of the singularities of the functions k _j () lie in the upper half of the plane. The implications of this result for the Kramers-Kronig relations are discussed.     

Product inhibition in packed-bed immobilized enzyme reactors for complex processes: Modelling of two-substrate processes

An analytical model was developed for describing the performance of packed-bed enzymic reactors operating with two cosubstrates, and when one of the reaction products is inhibitory to the enzyme. To this aim, the compartmental analysis technique was used. The relevant equations obtained were solved numerically, and the effect of the main operational parameters on the reactor characteristics were studied.Notation C _infa,i ^sup* local concentration of products in the pores of stage i - C _j,i concentration of substrate j in the pores of stage i - D _infa ^sup* internal (pore) diffusion coefficient for the reaction product a - D _j internal (pore) diffusion coefficient of substrate j - J _infa,i ^sup* net flux of product a, taking place from the pores of stage i into the corresponding bulk phase - J _j,i net flux of substrate j, taking place from the bulk phase of stage i into the corresponding pores - K _b inhibition constant - K _m,1, K _m,2 Michaelis constants for substrate 1 and 2, respectively - K _q inhibition constant - n total number of elementary stages in the reactor - Q volumetric flow rate throughout the reactor - R _j,i, R _infa,i ^sup* local reaction rates in pores of stage i, in terms of concentration of substrate j and product a respectively - S _infa,i ^sup* , S _infa,i-1 ^sup* bulk concentration of the reaction product a, in the stages i and i — 1, respectively - S _j,0 concentration of substrate j in the reactor feed - S _j,i-1, S _j,i concentration of substrate j in the bulk phase leaving stages i — 1 and i, respectively - V total volume of the reactor - V _m maximal reaction rate in terms of volumetric units - y axial coordinate of the pores - y ₀ depth of the pores - ^* dimensionless parameter, defined in Equation (22) - ₁ dimensionless parameter, defined in Equation (6) - ₂ dimensionless parameter, defined in Equation (6) - ₁ dimensionless parameter, defined in Equation (6) - ₂ dimensionless parameter, defined in Equation (6) - ^* dimensionless parameter, defined in Equation (22) - ₁ dimensionless parameter, defined in Equation (6) - ₂ dimensionless parameter, defined in Equation (6) - ^* dimensionless parameter, defined in Equation (22) - ^* dimensionless parameter, defined in Equation (22) - volumetric packing density of catalytic particles (dimensionless) - porosity of the catalytic particles (dimensionless) - V _infi ^sup* dimensionless concentration of reaction product in pores of stage i, defined in Equation (17) - j,i dimensionless concentration of substrate j in pores of stage i; defined in Equation (6) - _j,i-1, _j.i dimensionless concentration of substrate j in the bulk phase of stage i; defined in Equation (6) - dimensionless position along the pore; defined in Equation (6)     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

Concentration-dependent diffusion in a composite slab with and without chemical reactions

Concentration-dependent diffusion of solute in a composite slab is investigated. The complex diffusion problem can be described by a set of nonlinear diffusion equations which is coupled to each other through the nonlinear interfacial boundary conditions. A two-layer diffusion is illustrated and the coupled nonlinear diffusion equations are conveniently solved by the orthogonal collocation method. Numerical simulation of the example reveals many interesting diffusion characteristics which are quite different from those in a single slab diffusion system.Nomenclature a _j expansion coefficient - A _i,j element of collocation matrix - B _i,j element of collocation matrix - C _a , C _b surface concentration - C _i concentration in the ith layer - D _i diffusion coefficient in the ith layer - D _i0 diffusion coefficient at very low concentration - k _i reaction rate in the ith layer - K _i dimensionless reaction rate, k _i l _i ² c _a ^m–1 /D ₁₀ - l _i thickness of the ith layer - m order of chemical reaction - n order of the orthogonal polynomial approximation - P _j–1(x _i ) orthogonal polynomial of order j - t time - x _i coordinate of the ith layer - X _i dimensionless coordinate of the ith layer, x _i/l _i - ratio of diffusion coefficient at low concentration, D ₂₀/D ₁₀ - ratio of thicknesses of layer, l ₁/l ₂ - _i dimensionless parameter in the concentration-dependent function of the ith layer - ratio of surface concentration, C _b /C _a - dimensionless time, tD ₁₀/l ₁ ² - _i dimensionless concentration in the ith layer, C _i /C _a      

Observations on the theoretical and experimental foundations of thermoplasticity

The thermodynamics of materials with internal state variables has been employed to study the properties of a class of thermoplastic materials in which the evolution equation for the internal variables is given by equation k⁽ⁱ⁾ = g⁽ⁱ⁾(σ_kl, '_kl, k⁽ⁱ⁾, θ, · '_kl) where g⁽ⁱ⁾ is homogeneous of order one in · '_kl. The most general form of the Helmholtz potential consistent with the assumption of insensitivity of the elastic relations to inelastic deformation has been derived and a geometric interpretation of the Clausius-Duhem restriction has been made employing the concept of a thermodynamic reference stress. Experimental results of one of the authors have been correlated with the theory.     

Measurement of the turbulent diffusivity in the near field of variable density jets using conditional velocimetry

Pulsed laser Mie scattering and laser Doppler velocimetry (LDV), both conditioned on the origin of the seed particles, have been successively performed in turbulent jets with variable density. In the early stages of the jet developments, significant differences are measured between the ensemble average LDV data obtained by jet seeding and those obtained by seeding the ambient air. Careful analysis of the marker statistics shows that this difference is a quantitative measure of the turbulent mixing. The good agreement with gradient–diffusion modelling suggests the validity of a general diffusion equation where the velocities involved are expressed in terms of ensemble conditional Favre averages. This operator accounts for all events (including intermittent ones) and for variations in the density of the marked fluid whose velocity is still specified by the binary origin of the marker.List of symbols D_L laminar diffusivity, m²/s - D_T turbulent diffusivity, m²/s - d diameter of the jet nozzle, m - F_r Froude number - J diffusion vector, m/s - k global sensitivity of the detection system for one particle (signal level) - N_P number of seed particles in the probe volume - N_P,i number of seed particles in sample i - N_P⁽ⁱ⁾ value of N_P in channel i - N_B number of Doppler bursts - count rate of bursts, s^–1 - N_v number of validated Doppler bursts - count rate of validated bursts, s^–1 - N_id number of ideal particles - N_id* number of marked ideal particles - P* probability that an ideal particle be marked by a seed particle - P(z) probability density function for z, m³/kg - probability to have k seed particles in the probe volume - probability of having k seed particle conditioned on a given value of z - r radial coordinate, m - R =⁽¹⁾/⁽²⁾, density ratio - S₁ local signal level with jet seeding - S₁⁽¹⁾ reference signal level in pure stream 1 with jet seeding - s₁ = S_1/S₁⁽¹⁾, normalized signal - v_c volumic capacity of the probe volume, m³ - V velocity vector, m/s - V_x axial velocity component, m/s - V_r radial velocity component, m/s - V_P particulate velocity vector, m/s - V_Pj velocity vector of particle j, m/s - V_Pij velocity vector of the jth particle in sample i, m/s - V_i velocity vector of the marked flow for realization i, m/s - V_1,i velocity vector of the flow such it is marked in realization i by particles issuing only from stream 1, m/s - x axial coordinate, m - Y_i local mass fraction of species i - Z mixture fraction:local mass fraction of jet fluid - Z_i mixture fraction for realization iGreek local density, kg/m³ - _i local density for realization i, kg/m³ - ⁽¹⁾ density in stream 1 (density of the jet fluid), kg/m³ - ₁ time of flight of jet seed particles to reach the probe volume, s - _B duration of a Doppler burst, sAverages <A> ensemble average of A - Ā time average of A - Favre average, , ( ) the present notation is only due to printing problems - A Favre fluctuation,      

Estimation of three-phase relative permeabilities

Starting from the statistical structural model of Alemánet al. (1988), we have developed an alternative to Stone's (1970, 1973; Aziz and Settari, 1979) methods for estimating steady-state, three-phase relative permeabilities from two sets of steady-state, two-phase relative permeabilities. Our result reduces to Stone's (1970; Aziz and Settari, 1979) first method, when the steady-state, two-phase relative permeability of the intermediate-wetting phase with respect to either the wetting phase or the nonwetting phase is a linear function of the saturation of the intermediate-wetting phase. As the curvature of either of these relative permeability functions increases, the deviation of our result from Stone's (1970; Aziz and Settari, 1979) first method increases. Currently, there are no data available that are sufficiently complete to form the basis of a comparison between our result and either of the methods of Stone (1970, 1973; Aziz and Settari, 1979).Notation a free parameter in Equation (19) - B(m, n) Beta function defined by Equation (17) - F ^(w), F^(nw) defined by Equations (31) and (27), respectively - G ⁽ⁱ⁾ defined by Equations (37) and (39) - H ⁽ⁱ⁾ defined by Equations (38) and (40) - k ⁽ⁱ⁾ three-phase relative permeability fo phasei - k ^(i)* defined by Equations (34) through (36) - k ^(i,j) relative permeability to phasei during a two-phase flow with phasej, possibly in the presence of an immobile phase - k ^(i,j)* defined by analogy with Equations (41) and (42) - k ^(i,j)** defined by Equations (49), (50), (53), and (54) - k _max ⁽ⁱ⁾ defined by Equation (11) - k ₁₉₇₀ ^(iw) defined by Equation (10) - k ₁₉₇₃ ^(iw)* defined by Equation (58) - k ₁₉₇₃ ^(iw) defined by Equation (13) - L length and diameter of cylindrical averaging surfaceS - L _t length of an individual capillary tube enclosed byS - L _t ^* defined by Equation (19) - L _t,min length of pore whose radius isR _max - N total number of pores contained within the averaging surfaceS - p ₁ ⁽ⁱ⁾ ,p ₂ ⁽ⁱ⁾ pressure of phasei at entrance and exit of averaging surfaceS, respectively - p defined by Equation (21) - p _c ^(i,j) capillary pressure function - p _c ^(i,j)* defined by Equations (23), (29), and (32) - p ⁽ⁱ⁾ intrinsic average of pressure within phasei defined by Alemánet al. (1988) - R pore radius - R ^* defined by Equation (18) - R _max maximum pore radius that occurs withinS - s ⁽ⁱ⁾ local saturation of phasei - s ^(i)* defined by Equation (7) - s _min ⁽ⁱ⁾ minimum or immobile saturation of phasei - S averaging surface introduced in local volume averaging - V ⁽ⁱ⁾ volume of phasei occupying the pore space enclosed byS Greek Letters , parameters in the Beta distribution defined by Equation (16) - ^(w), ^(nw) functions of only the wetting phase saturation and the non-wetting phase saturation, respectively. Introduced in Equation (6) - ^(i,j) interfacial tension between phasesi andj - (x) Gamma function - defined by Equation (57) - , spherical coordinates in system centered upon the axis of the averaging surfaceS - _max maximum value of , 45 °, in view of assumption (9) - ^(i,j) contact angle between phasesi andj measured through the displacing phase - ^(w),^(nw) functions of only the wetting phase saturation and the non-wetting phase saturation, respectively. Introduced in Equation (12) Other gradient operator Amoco Production Company, PO Box 591 Tulsa, OK 74102, U.S.A.     

Critical blowup and global existence numbers for a weakly coupled system of reaction-diffusion equations

Let D ⊂ R ^N be either all of R ⁿ or else a cone in R ^N whose vertex we may take to be at the origin, without loss of generality. Let p _i, q_j, i = 1, 2, be nonnegative with 0<p ₁+q ₁≦p ₂+q ₂. We consider the long-time behavior of nonnegative solutions of the system