不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

Coupled acoustic systems

The coupled systems in the acoustic theories of (i) a viscous compressible fluid, (ii) a gassy sediment, (iii) Landau's superfluid, and (iv) elastothermodiffusion form a subclass of coupled quasi-isotropic systems. They are, accordingly, treated through reduced versions of the general principles and general results covering coupled quasi-isotropic systems.     

THE GENERALITY OF THE GENERAL EQUATIONS OF DYNAMICS 1)

理论力学中动力学普遍方程,在分析力学中称为d'Alembert--Lagrange原理。动力学普遍方程之普遍在于,由它不仅可导出动力学普遍定理,可导出完整约束系统和非完整约束系统的运动微分方程,还可导出积分变分原理。     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

Exact asymptotic relations for the effective response of linear viscoelastic heterogeneous media

This article addresses the asymptotic response of viscoelastic heterogeneous media in the frequency domain, at high and low frequencies, for different types of elementary linear viscoelastic constituents. By resorting to stationary principles for complex viscoelasticity and adopting a classification of the viscoelastic behaviours based on the nature of their asymptotic regimes, either elastic or viscous, four exact relations are obtained on the overall viscoelastic complex moduli in each case. Two relations are related to the asymptotic uncoupled heterogeneous problems, while the two remaining ones result from the viscoelastic coupling that manifests itself in the transient regime. These results also provide exact conditions on certain integrals in time of the effective relaxation spectrum. This general setting encompasses the results obtained in preceding studies on mixtures of Maxwell constituents [1], [2].     

Instability of stationary and uniformly moving cylindrical fluid bodies—II. Viscoelastic threads and experimental observations

The instability analysis of Part I is extended to the breakup of viscoelastic threads in fluid media (also possibly viscoelastic). Critical Growth rates and wave-numbers are calculated in terms of the viscosity ratio, the Ohnesorge numbers (continuous and dispersed phases), and elasticity numbers for each of the respective phases. Comparisons with results for Newtonian systems indicate viscoelastic threads to be less stable than Newtonian threads under similar conditions. Also, the critical wave-numbers observed with viscoelastic threads can differ significantly from those observed with Newtonian systems, particularly if the relative magnitudes of elasticity of the dispersed and continuous phases are quite different. Systems with similar magnitudes of elasticity in each phase exhibit wave-numbers similar to Newtonian systems of similar viscosities.Experimental results obtained from observations of fluid thread breakup in a Taylor four-roller device provide a basis for checking the predictions of the lineararized theory for both Newtonian and viscoelastic systems. In general, the agreement is good and the theoretical predictions of Parts I and II seem to be reasonable representations of experimental fact.     

Unified expressions of all differential variational principles

A mathematical expression of the quantitative causal principle is given, using the expression this letter shows the unified expressions of D’Alembert–Lagrange, virtual work, Jourdian, Gauss and general D’Alembert–Lagrange principles of differential style, first finds the intrinsic relations among these variational principles, the conservation quantities of the above principles are first found, finally the Noether conservation charges of the all differential variational principles in the systems with the symmetry of Lie group D_m are first obtained.     

Effective Diffusivities of Point-Like Molecules in Isotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of point-like molecules in random and structured media are used to determine and characterize the effective diffusion coefficients of the molecules in the media. Simulations were carried out in 2D and 3D media. Monte Carlo simulation results in 2D and 3D media are compared with those obtained by analytical techniques. Simulation results indicate that for the structured, isotropic media the effective diffusivities can be characterized according to percolation thresholds in addition to porosity. The effective diffusivities in two isotropic media with the same porosity but different percolation thresholds can differ significantly. The effects of dimensionality on the effective diffusivities can also be significant. It is shown that in general the effective diffusion coefficients obtained from 2D simulation are not a good approximation to those of 3D, especially when the percolation thresholds of the 2D media and the 3D media are very different.     

Electrogas- and electrohydrodynamic control of gas and liquid jets and flows. 2. Physical principles

The paper deals with the electrophysical principles of the method of electrogas- and electrohydrodynamic conversion of electric signals to pneumatic (hydraulic) signals and vice versa for the control of gas and liquid jets and flows in electropneumohydraulic systems, including problems of formulating requirements for the working media and electric circuits of converters, estimating the limiting ranges of pressure, velocity, and temperatures of the working media and limiting dynamic possibilities, and determining conditions for production of potential signals, stable control, and extension of the control ranges. St. Petersburg State Technical University, St. Petersburg 195251. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 34–42, May–June, 2000.     

Overall constitutive description of symmetric layered media based on scattering of oblique SH waves

This papers investigates the scattering of oblique shear horizontal (SH) waves off finite periodic media made of elastic and viscoelastic layers. It further considers whether a Willis-type constitutive matrix (in temporal and spatial Fourier domain) may reproduce the scattering matrix (SM) of such a system. In answering this question the procedure to determine the relevant overall constitutive parameters for such a medium is presented. To do this, first the general form of the dispersion relation and impedances for oblique SH propagation in such coupled Willis-type media are developed. The band structure and scattering of layered media are calculated using the transfer matrix (TM) method. The dispersion relation may be derived based on the eigen-solutions of an infinite periodic domain. The wave impedances associated with the exterior surfaces of a finite thickness slab are extracted from the scattering of such a system. Based on reciprocity and available symmetries of the structure and each constituent layer, the general form of the dispersion and impedances may be simplified. The overall quantities may be extracted by equating the scattering data from TM with those expected from a Willis-type medium. It becomes evident that a Willis-type coupled constitutive tensor with components that are assumed independent of wave vector is unable to reproduce all oblique scattering data. Therefore, non-unique wave vector dependent formulations are introduced, whose SM matches that of the layered media exactly. It is further shown that the dependence of the overall constitutive tensors of such systems on the wave vector is not removable even at very small frequencies and incidence angles and that analytical considerations significantly limit the potential forms of the spatially dispersive constitutive tensors.     

The elliptic paraboloid failure surface for transversely isotropic materials off-axis loaded

The elliptic paraboloid failure surface has been well established as a potential criterion for yielding and failure of transversely isotropic materials, presenting also the strength differential effect [1]. This was done by extending well established criteria for isotropic materials presenting the strength differential effect (SDE), through an introduction process which maintained basic physical principles for the anisotropic materials. All previous literature concerned the special case where the principal axes of the external loading coincided with the principal strength axes of the material. In this paper the most general case where the two systems of frames are arbitrarily oriented relatively to each other is considered. In this situation the simplifications derived from the coincidence of the external principal stress and material principal strength axes are lost and the material should be considered as a general orthotropic one. The general properties for such types of loading of the transversely isotropic material are established by maintaining the general features of the failure locus invariant. Then, this study completes the investigation of yielding and failure mode of such materials considering the most general case of their loading.     

Propagation of electromagnetic disturbances and stability of stationary states in media with a nonlinear ohm's law

The principles of propagation and development of small disturbances in nonlinear conducting media are studied.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 3–15, May–June, 1972.The author thanks G. A. Lyubimov and S. A. Regirer for their evaluation of basic results, and A. G. Kulikovski for his valuable advice.     

Diagnostics of the medium structure by long wave of finite amplitude

The averaged systems of hydrodynamic equations for a structured medium in the Lagrangian and the Eulerian coordinates are discussed. In the general case, the equations cannot be reduced to the average hydrodynamic terms. Under propagation of long waves in media with structure, the non-linear effects appear and they are analyzed in the framework of the asymptotic averaged model. The heterogeneity in a medium structure always increases the non-linear effects for the long-wave perturbations. A new method for diagnostics of the properties of medium components by long non-linear waves is suggested (inverse problem). The mass contents of components in the media can be determined by this diagnostic method.     

Effective Diffusivity Tensors of Point-Like Molecules in Anisotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of random walks in anisotropic structured media are performed to determine the dependence of effective diffusivities on geometrical properties. The anisotropic media used in this study are periodic systems, which are generated by extending primitive, face-centered, and body-centered unit cells indefinitely in all axial directions. Results of simulations compare well with published experimental data and the calculations by the volume averaging method. In addition, these results suggest that if the 2D media with percolation thresholds subtantially differ from those of 3D, 2D approximations of 3D media are not satisfactory. When percolation thresholds are the same, the effective diffusivity tensors depend solely on the porosity. This fact has been suggested for isotropic media and it seems to hold for anisotropic media.     

Layered Thermohaline Convection in Hypersaline Geothermal Systems

Thermohaline convection occurs in hypersaline geothermal systems due to thermal and salinity effects on liquid density. Because of its importance in oceanography, thermohaline convection in viscous liquids has received more attention than thermohaline convection in porous media. The fingered and layered convection patterns observed in viscous liquid thermohaline convection have been hypothesized to occur also in porous media. However, the extension of convective dynamics from viscous liquid systems to porous media systems is complicated by the presence of the solid matrix in porous media. The solid grains cause thermal retardation, hydrodynamic dispersion, and permeability effects. We present simulations of thermohaline convection in model systems based on the Salton Sea Geothermal System, California, that serve to point out the general dynamics of porous media thermohaline convection in the diffusive regime, and the effects of porosity and permeability, in particular. We use the TOUGH2 simulator with residual formulation and fully coupled solution technique for solving the strongly coupled equations governing thermohaline convection in porous media. We incorporate a model for brine density that takes into account the effects of NaCl and CaCl2. Simulations show that in forced convection, the increased pore velocity and thermal retardation in low-porosity regions enhances brine transport relative to heat transport. In thermohaline convection, the heat and brine transport are strongly coupled and enhanced transport of brine over heat cannot occur because buoyancy caused by heat and brine together drive the flow. Random permeability heterogeneity has a limited effect if the scale of flow is much larger than the scale of permeability heterogeneity. For the system studied here, layered thermohaline convection persists for more than one million years for a variety of initial conditions. Our simulations suggest that layered thermohaline convection is possible in hypersaline geothermal systems provided the vertical permeability is smaller than the horizontal permeability, as is likely in sedimentary basins such as the Salton Trough. Layered thermohaline convection can explain many of the observations made at the Salton Sea Geothermal System over the years.     

On tortuosity and areosity tensors for porous media

An average streamwise channel velocity is proposed as a more accurate representation of the actual intrapore velocity than the intrinsic phase average velocity. A relationship is derived between the average streamwise channel velocity and the interstitial velocity and superficial velocity. New definitions of tortuosity and areosity as second-order tensors are proposed for porous media in general. Novel names, semantically in line with the respective physical meanings, are proposed for these quantities. The definitions produce results which conform with several other published results and are applicable to anisotropic media in general.