多孔材料中声波的传播与演化

采用两相多孔介质的拉格朗日模型来描述一种理论流体充填的多孔弹性固体材料,其中孔隙度的变化满足一个附加的平衡方程。     

Propagation of SH waves in an irregular monoclinic crustal layer

The present paper discusses the dispersion equation for SH waves in a monoclinic layer over a semi-infinite elastic medium with an irregularity. In the absence of the irregularity, the dispersion equation reduces to standard dispersion equation for SH waves in a monoclinic layer over an isotropic semi-infinite medium. The dispersion curves for different size of the irregularity are computed and compared for the half-space without any irregularity. It can be seen that the phase velocity is strongly influenced by the wave number and the depth of the irregularity.     

On the dispersion of a solute in Poiseuille flow of a third-grade fluid in a channel

Gill and Sankarasubramanian's analysis of the dispersion of Newtonian fluids in laminar flow between two parallel walls are extended to the flow of non-Newtonian viscoelastic fluid (known as third-grade fluid) using a generalized dispersion model which is valid for all times after the solute injection. The exact expression is obtained for longitudinal convective coefficient K₁(Γ), which shows the effect of the added viscosity coefficient Γ on the convective coefficient. It is seen that the value of the K₁(Γ) for Γ≠0 is always smaller than the corresponding value for a Newtonian fluid. Also, the effect of the added viscosity coefficient on the K₂(t,Γ) (which is a measure of the longitudinal dispersion coefficient of the solute) is explored numerically. Finally, the axial distribution of the average concentration C_m of the solute over the channel cross-section is determined at a fixed instant after the solute injection for several values of the added viscosity coefficient.     

Discrete element modeling of acoustic emission in rock fracture

The acoustic emission (AE) features in rock fracture are simulated numerically with discrete element model (DEM). The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.     

Propagation of longitudinal and transverse waves in a multimodulus elastic medium

A problem of propagation of longitudinal and transverse waves in a multimodulus elastic isotropic medium is considered. In the model used, the medium is described by a potential depending on three invariants of strains, which allows the influence of preliminary deformation of the medium on the longitudinal and transverse velocities to be taken into account. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 176–182, July–August, 2009.     

Non-Equilibrium Mixing of Turbulence Scales Using a Continuous Hybrid RANS/LES Approach: Application to the Shearless Mixing Layer

A new subgrid-scale (SGS) model based on partially integrated transport method (PITM) is applied to the case of a turbulent spectral non-equilibrium flow created by the mixing of two turbulence fields of differing scales: the shearless mixing layer. The method can be viewed as a continuous hybrid RANS/LES approach. In this model the SGS length scale is no longer given by the size of the discretization step, but is dynamically estimated using an additional transport equation for the dissipation rate. The results are compared to those corresponding to the classical model of Smagorinsky and to the experimental data of Veeravalli and Warhaft. A method for creating an anisotropic analytical pseudo-random field for inflow conditions is also proposed. This approach based on subgrid-scale transport modelling combined with anisotropic inlet conditions gives better results for the prediction of the shearless mixing layer.     

Propagation of Harmonic Waves in a Cylindrical Shell (Timoshenko Model)

Equations for studying the axisymmetric deformation of a cylindrical shell are derived on the basis of Timoshenko shell theory. A dispersion equation is set up to study natural harmonic waves in a cylindrical isotropic shell. A numerical approach to plotting dispersion diagrams is proposed. The wave velocities obtained coincide with analytical solutions     

Kinetic modeling of multiple scattering of elastic waves in heterogeneous anisotropic media

In this paper we develop a multiple scattering model for elastic waves in random anisotropic media. It relies on a kinetic approach of wave propagation phenomena pertaining to the situation whereby the wavelength is comparable to the correlation length of the weak random inhomogeneities—the so-called weak coupling limit. The waves are described in terms of their associated energy densities in the phase space position  ×$\times$  wave vector. They satisfy radiative transfer equations in this scaling, characterized by collision operators depending on the correlation structure of the heterogeneities. The derivation is based on a multi-scale asymptotic analysis using spatio-temporal Wigner transforms and their interpretation in terms of semiclassical operators, along the same lines as Bal (2005). The model accounts for all possible polarizations of waves in anisotropic elastic media and their interactions, as well as for the degeneracy directions of propagation when two phase speeds possibly coincide. Thus it embodies isotropic elasticity which was considered in several previous publications. Some particular anisotropic cases of engineering interest are derived in detail.     

Controlling the Turbulent Flow Past a Thick Airfoil by Means of Flow Enhancement in Vortex Cells Using Suction from Central Body Surfaces

On the basis of a numerical simulation of the turbulent steady-state flow past a thick airfoil with vortex cells built into the body contour, an unconventional technique for controlling flow separation by means of distributed suction from central bodies embedded in the cells is analyzed over a wide range of Reynolds numbers and suction velocities.     

Identification of Impact Damage in S-2 Glass Composite Missile Casings Using Complementary Vibration and Wave Propagation Approaches

The future United States military arsenal is transitioning from homogeneous metallic materials to stronger and lighter heterogeneous composite materials. Although these composites offer numerous advantages such as the ability to tune strength to weight ratios for each particular application, composites are susceptible to numerous damage mechanisms and environmental factors. Accidental in-field impacts resulting from mishandling or transportation loads have been identified by the U.S. Army as the primary cause of damage in composite weapon systems. This paper presents a hybrid approach to detect, locate, and quantify damage in filament wound canisters using a complementary set of vibration-based (transmissibility and embedded sensitivity) and wave propagation-based (phased array beamforming) methods. It is shown that this hybrid approach accurately detects, locates, and quantifies the damage imposed by 1.36 and 6.78 N-m impact energies. By measuring the static stiffness of the baseline canister at the mid-point, it was determined that the 1.36 N-m impact resulted in an overall 6% reduction of the tubes bending stiffness, while the 6.78 N-m impact resulted in a 28% reduction.     

Propagation of harmonic waves in longitudinally reinforced cylindrical shells with low shear stiffness

The exact solutions of the equations of motion derived under refined theories of shells and ribs based on the Timoshenko model are used to plot dispersion curves for harmonic waves propagating along a cylindrical shell reinforced with discrete longitudinal ribs __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 5, pp. 37–41, May 2006.     

Kinetic modeling of particles in stratified flow – Evaluation of dispersion tensors in inhomogeneous turbulence

The continuum equations for a dilute particle distribution in inhomogeneous turbulence are tested against results from a Langevin particle tracking simulation. Reeks’ version of the kinetic theory is used to generate the mass, momentum and kinetic stress equations for the particle distribution. The particle tracking data are used to directly evaluate the dispersion tensors λ and μ which serve as closure relations for the continuum equations. These exact forms are compared to approximate, local forms. Even for low Stokes numbers (corresponding to low particle inertia defined by τ/τ_p ? 1), the tensor λ is strongly affected by the inhomogeneity and depends on turbulence parameters in the volume corresponding to the particle path dispersion over the particle Lagrangian integral timescale τ. In contrast, the locally homogeneous form of the velocity dispersion tensor μ is a sufficient approximation, since it depends on the dispersion volume over the much smaller particle relaxation time τ_p. It is demonstrated that the body force due to the dispersion vector γ cannot be neglected. In the limit of passive tracers (zero stopping distance), γ is equal to the gradient of λ, if the physical setting is such that we can invoke constant tracer density in this limit.     

Numerical Investigation of the Temperature Macroparameters in a Shock Wave in a Binary Gas Mixture Using the Kinetic Boltzmann Equation

The effect of the Mach number and the concentration and mass ratios on the behavior of the parallel, radial, and total temperatures of the components in a shock wave in a binary gas mixture is studied. The results obtained are compared with the theoretical and experimental results of other investigators.     

Upscaling and Simulation of Waterflooding in Heterogeneous Reservoirs Using Wavelet Transformations: Application to the SPE-10 Model

We have developed an accurate and highly efficient method for upscaling and simulation of immiscible displacements in three-dimensional (3D) heterogeneous reservoirs, which is an extension of the technique that we developed previously for 2D systems. The method utilizes wavelet transformations (WTs) to upscale the geological model of a reservoir, based on the spatial distribution of the single-phase permeabilities and the locations of the wells in the reservoir. It generates a non-uniform grid in which the resolved structure of the fine grid around the wells, as well as in the high-permeability sectors, are preserved, but the rest of the grid is upscaled. A robust uplayering procedure is used to reduce the number of the layers, and the WTs are used to upscale each layer areally. To demonstrate the method’s accuracy and efficiency, we have applied it to the geological model of a highly heterogeneous reservoir put forward in the tenth Society of Petroleum Engineers comparative solution project (the SPE-10 model), and carried out simulation of waterflooding in the upscaled model. Various upscaling scenarios were examined, and although some of them resulted in efficient simulations and accurate predictions, the results when non-uniform upscaling is used based on the WT technique are in excellent agreement with the solution of the same problem in the fine grid of the SPE-10 model. Most importantly, the speed-up factors that we obtain are several orders of magnitude. Hence, the method renders it unnecessary to use massively parallel computations for such problems.     

The Onset of Convection in an Anisotropic Porous Layer Using a Thermal Non-Equilibrium Model

The stability of a horizontal fluid saturated anisotropic porous layer heated from below and cooled from above is examined analytically when the solid and fluid phases are not in local thermal equilibrium. Darcy model with anisotropic permeability is employed to describe the flow and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is implemented to compute the critical Rayleigh number and the corresponding wavenumber for the onset of convective motion. The effect of thermal non-equilibrium and anisotropy in both mechanical and thermal properties of the porous medium on the onset of convection is discussed. Besides, asymptotic analysis for both very small and large values of the interphase heat transfer coefficient is also presented. An excellent agreement is found between the exact and asymptotic solutions. Some known results, which correspond to thermal equilibrium and isotropic porous medium, are recovered in limiting cases.     

Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.