A Capillary Microstructure of the Wetting Front

This article reports the experimental results of a study of the wetting-front microscale structure formed only by capillary forces in homogeneous and random etched glass capillary models. In the homogeneous model, water propagates through the capillary system, evenly filling the capillaries across the direction of flow. Air is trapped by the pinch-off mechanism inside the pore bodies in the form of individual bubbles. The experiments specified three consecutive steps of the pinch-off mechanism, film flow, snap-off, and interface movement. In the random model, both the bypass and pinch-off, forming bypass/cut-off mechanism, create residual air structure. Bypass traps air inside large capillary-pore aggregates which are bounded by small-diameter capillaries in where pinch-off traps air in the adjacent pores. An analysis of the residual air distribution versus depth below the surface in the homogeneous and random micromodels made it possible to identify three successive zones, namely a transition zone, a transmission zone, and a wetting-and-front zone. In the transition zone, the residual air content increases with depth from zero to the constant value in the transmission zone where it remains practically constant. The capillary processes within the wetting-and-front combined zone govern air replacement with wetting and formation of the transmission zone.     

On the localness of the spectral energy transfer in turbulence

By utilizing available experimental data for net energy transfer spectra for homogeneous turbulence, contributions P(, ) to the energy transfer at a wavenumber from various other wavenumbers are calculated. This is done by fitting a truncated power-exponential series in and to the experimental data for the net energy transfer T(), and using known properties of P(, ). Although the contributions P(, ) obtained by using this procedure are not unique, the results obtained by using various assumptions do not differ significantly. It seems clear from the results that for a region where the energy entering a wavenumber band dominates that leaving, much of the energy entering the band comes from wavenumbers which are about an order of magnitude smaller. That is, the energy transfer is rather nonlocal. This result is not significantly dependent on Reynolds number (for turbulence Reynolds numbers based on microscale from 3 to 800). For lower wavenumbers, where more energy leaves than enters a wavenumber band, the energy transfer into the band is more local, but much of the energy then leaves at distant wavenumbers.     

LDA measurements in low Reynolds number turbulent boundary layer

LDA measurements of the mean velocity in a low Reynolds number turbulent boundary layer allow a direct estimate of the friction velocity U from the value of /y at the wall. The trend of the Reynolds number dependence of / is similar to the direct numerical simulations of Spalart (1988).     

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

ON THE UNIFICATION OF THE HAMILTON PRINCIPLES IN NONHOLONOMIC SYSTEM AND IN HOLONOMIC SYSTEM

ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬＥＳＩＮＮＯＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭＡＮＤＩＮＨＯＬＯＮＯＭＩＣＳＹＳＴＥＭ（梁立孚）（韦扬）ＯＮＴＨＥＵＮＩＦＩＣＡＴＩＯＮＯＦＴＨＥＨＡＭＩＬＴＯＮＰＲＩＮＣＩＰＬ...     

Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Influence of vorticity on the heat transfer to blunt bodies in hypersonic flight

When blunt bodies are in hypersonic flight, a high-entropy layer of gas with nonzero vorticity is formed near their surface. The transverse gradients of the entropy, density, and gas velocity in the layer are high, which makes it necessary to take into account its absorption by the boundary layer of finite thickness . This vortex interaction is usually accompanied by an increase in the heat flux q and the frictional stress on the wall compared with their values as calculated in accordance with the classical scheme of a thin boundary layer, when the parameters on the outer edge of the boundary layer are set equal to the inviscid parameters on the body. This effect has been investigated on the side surface of slender (with angle 1 to the undisturbed flow) blunt bodies in a hypersonic stream [1–3]. It is shown in the present paper that the effect can have a stronger and even qualitative influence on the flow over blunt bodies with 1 if the radius of curvature R_s of the detached shock wave on the axis is small compared with the midsection radius R of the body. It is shown that the distributions of the heat fluxes with allowance for the vorticity of the inviscid shock layer are similar in the case of slightly blunt (r₀/R 0) cones with half-angles less than a critical _*.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 50–57, March–April, 1981.     

On near-wall turbulence-generating events in a turbulent boundary layer on a riblet surface

The boundary layer over a drag reducing riblet surface is investigated using hot-wire anemometry and flow visualisation. The concept of a riblet sublayer is introduced, and a definition is proposed in terms of a region of reduced turbulence energy production formed near the wall by the addition of riblets. The hot wire records are examined using a modified form of quadrant analysis, and results obtained over plain and riblet surfaces are compared. Close to the wall, the addition of riblets produces a marked reduction in the occurrence of ejection (2nd quadrant) events. A corresponding increase in the incidence of sweep (4th quadrant) events is accompanied by the development of a strong tendency toward a preferred event duration, and a preferred interval between events. These changes diminish rapidly with distance from the surface, becoming almost undetectable beyondy ⁺=40. They are discussed in the light of flow visualisation results, and interpreted in terms of mechanisms associated with the interaction between the riblets and the inner boundary layer flow structures. A conceptual model of the flow mechanisms in the riblet sublayer is proposed.     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Approximate relationships for the determination of the pressure on the surface of a sphere or cylinder with an arbitrary mach number in the free stream

Numerical calculations have been made [1–4] of the pressure distribution over the surface of a sphere or cylinder during transverse flow in the range 0 /2, where is the angle reckoned from the stagnation point along the meridional plane, and on the basis of these results simple analytical equations have been proposed in order to determine the pressure for arbitrary Mach numbers M in the free stream. The gas is assumed to be ideal and perfect.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 185–188, March–April, 1985.     

On the vanishing of the additive measures of strain and rotation for finite deformations

Explicit formulae for the finite strain and rotation measures are given, in the cases when either one of the infinitesimal tensors of strain and rotation vanishes. Conversely, when the finite strain or rotation measure vanishes, explicit formulae for the infinitesimal tensors of strain and rotation are also obtained.     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

An experimental study on vortex breakdown in a differentially-rotating cylindrical container

The vortex breakdown phenomenon in a closed cylindrical container with a rotating endwall disk was reproduced. Visualizations were performed to capture the prominent flow characteristics. The locations of the stagnation points of breakdown bubbles and the attendant global flow features were in excellent agreement with the preceding observations. Experiments were also carried out in a differentially-rotating cylindrical container in which the top endwall rotates at a relatively high angular velocity _t, and the bottom endwall and the sidewall rotate at a low angular velocity _sb. For a fixed cylinder aspect ratio, and for a given relative rotational Reynolds number based on the angular velocity difference _t–_sb, the flow behavior is examined as |_sb/_t| increases. For a co-rotation (_sb/_t>0), the breakdown bubble is located closer to the bottom endwall disk. However, for a counter-rotation (_sb/_t<0), the bubble is seen closer to the top endwall disk. For sufficiently large values of _sb, the bubble ceases to exist for both cases.     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

Aspect ratio effect on flow-induced forces on circular cylinders in a cross-flow

Finite-span circular cylinders with two different aspect ratios, placed in a cross-flow, are investigated experimentally at a cylinder Reynolds number of 46,000. Simultaneous measurements of the flow-induced unsteady forces on the cylinders and the stream velocity in the wake are carried out. These results together with mean drag measurements along the span and available literature data are used to evaluate the flow mechanisms responsible for the induced unsteady forces and the effect of aspect ratio on these forces. The coherence of vortex shedding along the span of the cylinder is partially destroyed by the separated flow emanating from the top and by the recirculating flow behind the cylinder. As a result, the fluctuating lift decreases drastically. Based on the data collected, it is conjectured that the fluctuating recirculating flow behind the cylinder is the flow mechanism responsible for the unsteady drag and causes it to increase beyond the fluctuating lift. The fluctuating recirculating flow is a direct consequence of the unsteady separated flow. The unsteady forces vary along the span, with lift increasing and drag decreasing towards the cylinder base. When the cylinder span is large compared to the wall boundary layer thickness, a submerged two-dimensional region exists near the base. As the span decreases, the submerged two-dimensional region becomes smaller and eventually vanishes. Altogether, these results show that fluctuating drag is the dominant unsteady force in finite-span cylinders placed in a cross-flow. Its characteristic frequency is larger than that of the vortex shedding frequency.List of symbols a span of active element on cylinder, = 2.5 cm - C _D local rms drag coefficient, 2D/ U ² da - C _L local rms lift coefficient, 2l/ U ² da - C _D local mean drag coefficient, 2D/ U ² da - C _D spanwise-averaged C _D for finite-span cylinder - (C _D ) _2D spanwise-averaged mean drag coefficient for two-dimensional cylinder - C _p pressured coefficient - -(C _p ) _b pressure coefficient at = - d diameter of cylinder, = 10.2 cm - D fluctuating component of instantaneous drag - D local rms of fluctuating drag - D local mean drag - E _D power spectrum of fluctuating drag, defined as - E _L power spectra of fluctuating lift, defined as - f _D dominant frequency of drag spectrum - f _L dominant frequency of lift spectrum - f _u dominant frequency of velocity spectrum - h span of cylinder - H height of test section, = 30.5 cm - L fluctuating component of instantaneous lift - L local rms of fluctuating lift - R _Du () cross-correlation function of streamwise velocity and local drag, - R _Lu () cross-correlation function of stream wise velocity and local lift, - Re Reynolds number, U d/y - S _L Strouhal number based on f _L ,f _L d/U - S _D Strouhal number based on f _D ,f _D d/U - S _u Strouhal number based on f _u , f _u d/U - t time - u fluctuating component of instantaneous streamwise velocity - U mean streamwise velocity - mean stream velocity upstream of cylinder - x streamwise distance measured from axis of cylinder - y transverse distance measured from axis of test section - z spanwise distance measured from cylinder base - angular position on cylinder circumference measured from forward stagnation - kinematic viscosity of air - density of air - time lag in cross-correlation function - _D normalized spectrum of fluctuating drag - _L normalized spectrum of fluctuating lift     

Investigation of the flow past wings of infinite span with a blunt leading edge in the vorticity interaction regime

An asymptotic analysis of the Navier-Stokes equations is carried out for the case of hypersonic flow past wings of infinite span with a blunt leading edge when 0, Re , and M . Analytic solutions are obtained for an inviscid shock layer and inviscid boundary layer. The results of a numerical solution of the problems of vorticity interaction at the blunt edge and on the lateral surface of the wing are presented. These solutions are compared with the solution of the equations of a thin viscous shock layer and on the basis of this comparison the boundaries of the asymptotic regions are estimated.deceasedTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 120–127, November–December, 1987.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters A interfacial area of the - interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the - interface contained within the averaging volume, m² - A ^* interfacial area of the - interface contained within a unit cell, m² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - p pressure in the -phase, N/m² - p intrinsic phase average pressure for the -phase, N/m² - p –p , spatial deviation of the pressure in the -phase, N/m² - r ₀ radius of the averaging volume and radius of a capillary tube, m - v velocity vector for the -phase, m/s - v phase average velocity vector for the -phase, m/s - v intrinsic phase average velocity vector for the -phase, m/s - v –v , spatial deviation of the velocity vector for the -phase, m/s - V averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ Greek Letters V/V, volume fraction of the -phase - mass density of the -phase, kg/m³ - viscosity of the -phase, Nt/m² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - a _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Harmonic solutions of micropolar elastodynamics

The equations of micropolar elastodynamics are considered for an unbounded continuum subjected to a body force and a body couple. These act harmonically with the same real frequency , but with individual arbitrary spatial distributions. Over a harmonic state, the displacement and microrotation are related to two radiation conditioned harmonic vectors, each acquiring three eigenvalue contributions, assuming a noncritical -frequency. Altogether, four distinct eigenvalues are admissible. If ²<2² ₀, ₀ being a frequency parameter of the continuum, two of these are real while two are purely imaginary. But if ²<2² ₀, then all admissible eigenvalues are real. Each eigenvalue contribution resolves into a series of Hankel and Bessel functions coupled to Hankel type transforms of: (i) spherical integrals which, in turn, can be expanded via spherical harmonics for the 3-dimensional problem, (ii) circular integrals for the 2-dimensional problem. Axisymmetric and spherically symmetric results are deduced in 3-dimensions. Asymptotic solutions are also established; they disclose long-range formation of radially attenuated spherical (or circular) waves propagating with, generally, anisotropic amplitudes but, invariably, isotropic eikonals.If, in the absence of a body couple, a body force acts radially in 3-dimensions with a spherically symmetric strength, then the elastic displacement behaves likewise while the microrotation vanishes identically. Another application is made to a 2-dimensional problem for a 1 × 3 source system of body force plus body couple without longitudinal variation but with magnitudes symmetric about a longitudinal axis.As approaches a certain critical frequency , dependent solely on the continuum, at least two eigenvalues approach the same value. The phenomenon is explored for a continuum consistent with ²<2² ₀ and under the hypothesis ²<2² ₀. All admissible eigenvalues are then real throughout an -neighbourhood of . Here, two associated eigenvalue contributions behave singularly. Nevertheless, their essential singularities cancel out within the relevant combination. Examination of a far-field suggests that critical frequency attainment sets off a slow instability in the 2-dimensional configuration. In the 3-dimensional configuration, however, it preserves stability and eliminates radial attenuation; an exact solution is formulated for this case.