Upscaling in Subsurface Transport Using Cluster Statistics of Percolation

Transport/flow problems in soils have been treated in random resistor network representations (RRNs). Two lines of argument can be used to justify such a representation. Solute transport at the pore-space level may probably be treated using a system of linear, first-order differential equations describing inter-pore probability fluxes. This equation is equivalent to a random impedance network representation. Alternatively, Darcys law with spatially variable hydraulic conductivity is equivalent to an RRN. Darcys law for the hydraulic conductivity is applicable at sufficiently low pressure head in saturated soils, but only for steady-state flow in unsaturated soils. The result given here will have two contributions, one of which is universal to any linear conductance problem, i.e., requires only the applicability of Darcys (or Ohms) law. The second contribution depends on the actual distribution of linear conductances appropriate. Although nonlinear effects in RRNs (including changes in resistance values resulting from current, analogous to changes in matric potential resulting from flow) have been treated within the framework of percolation theory, the theoretical development lags the corresponding development of the linear theory, which is, in principle, on a solid foundation. In practice, calculations of the nonlinear conductivity in relatively (compared with soils) well characterized solid-state systems such as amorphous or impure semiconductors, do not agree with each other or with experiment. In semiconductors, however, experiments do at least appear consistent with each other.In the limit of infinite system size the transport properties of a sufficiently inhomogeneous medium are best calculated through application of critical rate analysis with the system resistivity related to the critical (percolating) resistance value, R_c. Here well-known cluster statistics of percolation theory are used to derive the variability, W (R,x) in the smallest maximal resistance, R of a path spanning a volume x³ as well as to find the dependence of the mean value of the conductivity, (x). The functional form of the cluster statistics is a product of a power of cluster size, and a scaling function, either exponential or Gaussian, but which, in either case, cuts off cluster sizes at a finite value for any maximal resistance other than R_c. Either form leads to a maximum in W (R,x) at R=R_c. When the exponential form of the cluster statistics is used, and when individual resistors are exponential functions of random variables (as in stochastic treatments of the unsaturated zone by the McLaughlin group [see Graham and MacLaughlin (1991), or the series of papers by Yeh et al. (1985, 1995), etc.], or as is known for hopping conduction in condensed matter physics), then W (R,x) has a power law decay in R/R_c (or R_c/R, the power being an increasing function of x. If the statistics of the individual resistors are given by power law functions of random variables (as in Poiseiulles Law), then an exponential decay in R for W (R,x) is obtained with decay constant an increasing function of x. Use, instead, of the Gaussian cluster statistics alters the case of power law decay in R to an approximate power, with the value of the power a function of both R and x.     

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

This paper studies the generalized Lorenz canonical form of dynamical systems introduced by elikovský and Chen [International Journal of Bifurcation and Chaos 12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The ilnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have ilnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of ilnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.     

On Ericksen’s Theorem for Unconstrained Hyperelastic Materials

After reviewing the proof of Ericksens theorem for the set of unconstrained, homogeneous, isotropic, hyperelastic materials, Ericksens result is sharpened to cover the subset of materials that possess a natural configuration and satisfy the empirical inequalities. Mathematics Subject Classifications (2000) 74B20.     

Symbolic Sequences Generated by a Special Class of Discrete Systems With Discontinuity and Input

In this paper, we study the symbolic sequences generated by aclass of discrete systems. This class contains the double-modulators, a typical example of discretetime electronic systems with discontinuity and input. First we develop ageneral theory and then we apply it to some examples in order to obtainsets of inadmissible sequences.     

The Experimental Approach to Effective Stress Law of Coal Mass by Effect of Methane

Effective stress law of all kinds of coal samples, including steam coal, fat coal, corking coal, thin coal and anthracite, under pore pressure of gas, is experimentally studied using a newly developed test machine. These samples are taken from Coal Mines in Wuda, Hebi, Yanzhou, Yangquan, Qingshui, and Gujiao in China. The experiment results show that, under pore pressure of gas, the tested coal samples comply with Biots effective stress law, where the Biots coefficient is not a constant, and is bilinear function of volumetric stress () and pore pressure (p), that is, We define four areas according to the numerical feature of , that is, functionless area of pore pressure, normal function area, fracturing function area, and quasi-soil function area. The effective stress law of coal mass introduced by this paper is a constitutive equation in the study of coupled solid and fluid. This has significance in the drainage and outburst of methane in coal seam.     

Permeability Effects of Turbulent Flow Through a Porous Insert in a Backward-Facing-Step Channel

In the present work, a k– model, based on the work of Lee and Howell (Proceedings of the ASME-JSME Thermal Engineering Hawaii, 1987), is rigorously derived based on time average of spatially averaged Navier–Stokes equations. The model is then employed to solve for a flow in a backward-facing step channel with a porous insert. The numerical solver is modified from the STREAM code (Lien and Leschziner, Comput. Meth. Appl. Mech. Eng. 114 (1994a) 123–148), and it has been validated against the experimental data of Seegmiller and Driver (AIAA Journal 23 (1985) 163–171). The code is then used to perform simulation for cases with a porous insert. The resistance of the porous insert can be altered by changing its permeability (), Forchheimers constant (F), or thickness (b). The goal is to examine the influence of each parameter on the resulting flow and turbulent kinetic energy (k) distributions. It is discovered that, by increasing the resistance of the insert, flow eventually enters a transitional regime towards relaminarization. This is due to the contribution of Darcys and Forchheimers terms in the governing equations, and modifying these two terms changes the levels of P_k and, hence, k and . Generally speaking, lowering or raising F results in a greater suppression of P_k than , causing the flow to relaminarize. Meanwhile, if the pore size is reasonably large to sustain turbulence within the porous media, increasing b reduces but does not eliminate the turbulent activity in the porous insert.     

GREEN＇S FUNCTION AND EFFECTIVE ELASTIC STIFFNESS TENSOR FOR ARBITRARY AGGREGATES OF CUBIC CRYSTALS

A closed but approximate formula of Green‘s function for an arbitrary aggregate of cubic crystallites is given to derive the effective elastic stiffness tensor of the polycrystal. This formula, which includes three elastic constants of single cubic crystal and five texture coefficients,accounts for the effects of the orientation distribution function (ODF) up to terms linear in the texture coefficients. Thus it is expected that our formula would be applicable to arbitrary aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy.Three examples are presented to compare predictions from our formula with those from Nishioka and Lothe‘s formula and Synge‘s contour integral through numerical integration. As an application of Green‘s function, we briefly describe the procedure of deriving the effective elastic stiffness tensor for an orthorhombic aggregate of cubic crystallites. The comparison of the computational results given by the finite element method and our effective elastic stiffness tensor is made by an example.     

The stress state of an elastic orthotropic medium with a penny-shaped crack

The static-equilibrium problem for an elastic orthotropic space with a circular (penny-shaped) crack is solved. The stress state of an elastic medium is represented as a superposition of the principal and perturbed states. To solve the problem, Willis approach is used, which is based on the triple Fourier transform in spatial variables, the Fourier-transformed Greens function for an anisotropic material, and Cauchys residue theorem. The contour integrals obtained are evaluated using Gauss quadrature formulas. The results for particular cases are compared with those obtained by other authors. The influence of anisotropy on the stress intensity factors is studied.__________Translated from Prikladnaya Mekhanika, Vol. 40, No. 12, pp. 76–83, December 2004.     

Numerische Untersuchung des elastisch-plastischen Verhaltens eines Fadenverbundwerkstoffes mit metallischer Matrix Eine Fallstudie

Übersicht Die Fehlervorhersage eines einachsigen Fadenverbundwerkstoffes mit metallischer Matrix erfordert die Kenntnis des nichtlinearen elastisch-plastischen Verhaltens in mikroskopischer Abmessung. Als Fallstudie wurde ein spezieller Verbundwerkstoff mit einem FEM-Programm untersucht. Im plastischen Bereich wurde das klassische, v. Mises Potential' mit dem neuen Übergangsfließpotential unter Berücksichtigung der plastischen Volumendehnung verglichen. Unter transversaler Normalbelastung zeigte der Verbund deutliche Unterschiede in der Zunahme und der Ausdehnung der lokalen plastischen Zonen. Alle kritischen Verformungsstadien wurden von dem Übergangsfließpotential bei geringerer Belastung erreicht.

---

Numerical investigation of the elastic-plastic behaviour of a fibre-reinforced composite with a metallic matrix

Summary Failure prediction of unidirectional fibre-reinforced composite with a metallic matrix needs knowledge about the nonlinear elastic-plastic behaviour in a microscopic scale. A specific composite was investigated using a FEM-program as a case study. In the plastic range the classical v. Mises Potential was compared with the new Transition Flow Potential, taking into account the plastic volume dilatation. Subjected to transverse normal loading the composite showed evident differences in the increase and the spread of locally plastic regions. All critical deformation states were reached by the Transition Flow Potential at lower loading.

     

The modified Casson＇s equation and its application to pipe flows of shear-thickening fluid

A few additional data from our previous experiments were plotted to emphasize the shear-thickening behavior of deoxy sickle erythrocyte (SS) suspension. A constitutive equation (named as FX equation) was developed and applied to a cylindrical pipe flow of a shear-thickening fluid. A blunt velocity profile and its volume flow rate were calculated. The flow was non-viscous (potential) in the central part of the pipe (i.e. the central core or the central plug-flow), and became more and more viscous towards the wall of the pipe after a specific radial distance, which was determined by a critical shear rate of (named as Fungs shear rate). Furthermore, combining the FX equation with the original Cassons equation, the author obtained a modified Cassons equation by introducing .The English text was polished by Yunming Chen.     

Analytic Solution to a Problem of Seepage in a Chequer-Board Porous Massif

A study is made of steady two-dimensional seepage in a porous massif composed by a double-periodic system of white and black chequers of arbitrary conductivity. Rigorous matching of Darcy's flows in zones of different conductivity is accomplished. Using the methods of complex analysis, explicit formulae for specific discharge are derived. Stream lines, travel times, and effective conductivity are evaluated. Deflection of marked particles from the natural direction of imposed gradient and stretching of prescribed composition of these particles enables the elucidation of the phenomena of transversal and longitudinal dispersion. A model of pure advection is related with the classical one-dimensional vective dispersion equation by selection of dispersivity which minimizes the difference between the breakthrough curves calculated from the two models.     

Design of a Beam with Controllable Force Elements

A method for solving the problem of design of an intellectual structure formulated for the pair optimal position of actuators, optimal control of actuators is developed. In the method proposed, physical and logical objects are treated as equivalent.     

Über singuläre Lastfälle in einer linearen Schalentheorie und ihre finite Behandlung

Übersicht Ausgehend von bekannten Fundamentallösungen für Platten bzw. Scheiben wird die Erstellung singulärer Ansatzfunktionen gezeigt, wie sie für finite Näherungsverfahren benötigt werden, die von den Funktionalen der totalen Energie bzw. der komplementären Energie ausgehen. Das Vorgehen wird eingehend an Kreiszylinderschalen erläutert.

---

Summary Starting from known fundamental solutions of plates the construction of singular basic estimate functions is shown. These are necessary in finite approximation methods basing on the functionals of total energy or complementary energy. The proceeding is explained in detail in a cylindrical shell analysis.

     

A Direct Method for the Identification of the Permeability Field Based on Flux Assimilation by a Discrete Analog of Darcy's Law

Inverse models to determine the permeability are generally based on existing forward models for the pressure. The permeabilities are adapted in such a way that the calculated pressures match the specified pressures in a number of points. To assimilate a priori knowledge about the flux, we introduce the flux assimilation method, which is based on the vector potential–pressure formulation of Darcy's law. Thanks to an unconventional discretization technique – the edge-based face element method – not only the specified pressures, but also specified information about the flux density can easily be assimilated. A relatively simple, but insightful analytical example illustrates the potential of this method.     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

First- and second-order sensitivities of beams with respect to cross-sectional cracks

This paper deals with a variational formulation for the sensitivity problem of beam systems in the context of deformable solids with cracks. Natural frequencies are defined as state variables involved in the energy functional of the system, while the cracks length and position are treated as system parameters. The hierarchical equation system is formed and solved for the first and second derivatives of the natural frequency functions with respect to the cracks length and position. An analytical procedure for calculations of the second-order sensitivities of natural frequencies is proposed for the non-symmetrical equation system operator. Numerical algorithms are worked out and implemented computationally. Analytical and numerical aspects of the problem are discussed in detail through a number of illustrative results.The support of this work by the State Committee for Scientific Research (KBN) under Grant No. 4-050-0148/17-98-00 is gratefully acknowledged.     

Brittle to Quasi-Brittle Transition

Abstract. The present study focuses on the kinetic and non-deterministic aspects of the brittle to quasi-brittle transition. A solid is approximated by a lattice formed by the interacting continuum particles and the evolution of damage is estimated using particle dynamics. The onset of transition is measured by the rate of the change of correlation length. The proposed method is illustrated on the examples of creep rupture, strain localization and dynamic expansion of a circular void in a brittle plate.Sommario. Viene posta l'attenzione sugli aspetti cinetici e non deterministici della transizione dal comportamento fragile a quello quasi-fragile. Un solido viene approssimato da un reticolo formato da particelle interagenti e l'evoluzione del danno viene stimata tramite la dinamica delle particelle. L'inizio della transizione viene misurato tramite la variazione della lunghezza di correlazione. Il metodo proposto viene illustrato su esempi di rottura per creep, localizzazione della deformazione e l'espansione di un foro circolare in una piastra fragile.     

Non-gaussian diffusion modeling in composite porous media by homogenization: Tail effect

In the paper anomalous diffusion appearing in a porous medium composed of two porous components of considerably different diffusion characteristics is examined. The differences in diffusivities are supposed to result either from two medium types being present or from variations in pore size (double porosity media). The long-tail effect is predicted using the homogenization approach based on the application of multiple scale asymptotic developments. It is shown that, if the ratio of effective diffusion coefficients of two porous media is of the order of magnitude smaller or equal O( ²), where is a homogenization parameter, then the macroscopic behaviour of the composite may be affected by the presence of tail-effect. The results of the theoretical analysis were applied to a problem of diffusion in a bilaminate composite. Analytical calculations were performed to show the presence of the long-tail effect in two particular cases.Notations c _i the concentration of chemical species in water within the medium i - D _i the effective diffusion coefficient for the medium i - D _ij ^eff the macroscopic (or effective) diffusion tensor in the composite - ERV the elementary representative volume - h the thickness of the period - l a chracteristic length of the ERV or the periodic cell - L a characteristic macroscopic length - n the volumetric fraction of the material 2 - 1–n the volumetric fraction of the material 1 - N the unit vector normal to - t the time variable - x the macroscopic (or slow) space variable - y the microscopic (or fast) space variable - c _1c ,C _2c ,D _1c ,D _2c the characteristic quantities - T,T _1L ,T _2L ,T _1l ,T _2l the characteristic times - c ₁ ^* ,c ₂ ^* ,D ₁ ^* ,D ₂ ^* ,t ^* the non-dimensional variables - the homogenization parameter - ₁ the domain occupied by the material 1 - ₂ the domain occupied by the material 2 - the interface between the domains 1 and 2 - the total volume of the periodic cell - /x_i the gradient operator - the gradient operator     

Point-source spark shadowgraphy at the historic birthplace of supersonic transportation

About a decade after its foundation, the most advanced ballistics laboratory in Germany at the time, the Kruppsche Schiessplatz, was utilized by Ernst Mach in 1888. His intent was to validate his pioneering shock wave research using military shells as supersonic vehicles. The 125th anniversary of the Schiessplatz was celebrated in 2002. Along with Machs research, it served to initiate the field of supersonic transportation technology. The specific subject of this paper is the application of point-source spark shadowgraphy at the same laboratory in the 1970s to visualize gas flow over aeroballistic projectiles. However, different from Machs original interest, the new purpose of spark photography at that time in the ballistic ranges of the German Bundeswehr was to find technical solutions to aeroballistic problems when field-testing gave incomplete answers. Both a qualitative and a quantitative understanding of the principles of aeroballistics were sought in this research.Received: 3 February 2003, Accepted: 30 June 2003, Published online: 2 September 2003An abridged version of this paper was presented at the 13th American Physical Society Topical Conference on Shock Compression of Condensed Matter at Portland, Oregon, from July 20 to 25, 2003