The purpose of this study is to present the development and application of residual flow procedure for analysis of three-dimensional (3-D) steady-state and transient seepage. The finite element equations are derived using a pseudo-variational principle which leads to a transient residual flow (load) vector that, in turn, is used to correct the position of the free surface iteratively. The procedure involves a fixed mesh which requires no mesh regeneration during transient analysis and during iterations. The procedure is also capable of handling material nonhomogeneities and anisotropy with relative ease. Several applications are made including verification with respect to closed-form solutions, and with results from a laboratory glass bead model simulating three-dimensional situations. For these glass beads, the coefficients of permeability and specific storage are also evaluated experimentally.

A method used to calculate the fundamental properties that predict the overall wicking performance of a fabric was proposed and executed. The combination of a horizontal and downward wicking test provided detailed measurements of the pertinent properties to wicking performance: capillary pressure and permeability. This method was proposed due to flaws found in standard vertical wicking tests as well as erroneous assumptions made in other wicking tests. Assumptions that capillary pressure and permeability are characteristic constants of porous structures are incorrect and will produce misleading information about that substrate. It was experimentally proven that these properties were a function of the saturation level found within the voids of a fabric. To obtain relevant capillary pressure and permeability data for a given fabric, a range of saturation levels were tested and analyzed. It was shown that saturation levels decreased as the vertical distance traveled by moisture increased. This phenomenon occurs as a result of capillary pressure within the voids dropping below the functional range needed to support flow in those voids at increasing heights. As height is increased, capillary pressure needs to also increase; therefore, only smaller radii pores will fill. Once saturation levels are known at specific heights, capillary pressure and permeability calculations were made using Darcys law and the Lucas-Washburn equation. Although this phenomenon is well known in civil engineering, it has not been widely addressed in the textile sciences, especially in its implications for wicking tests.

It has been shown by Neuman [1990], Di Federico and Neuman [1997, 1998a,b] and Di Federico et al. [1999] that observed multiscale behaviors of subsurface fluid flow and transport variables can be explained within the context of a unified stochastic framework, which views hydraulic conductivity as a random fractal characterized by a power variogram. Any such random fractal field is statistically nonhomogeneous but possesses homogeneous spatial increments. When the field is statistically isotropic, it is associated with a power variogram γ(s) = Cs²ᴴ where C is a constant, s is separation distance, and If is a Hurst coefficient (0 < H< 1). If the field is Gaussian it constitutes fractional Brownian motion (fBm). The authors have shown that the power variogram of a statistically isotropic or anisotropic fractal field can be constructed as a weighted integral from zero to infinity of exponential or Gaussian vario grams of overlapping, homogeneous random fields (modes) having mutually uncorrelated increments and variance proportional to a power 2H of the integral (spatial correlation) scale. Low- and high-frequency cutoffs are related to length scales of the sampling window (domain) and data support (sample volume), respectively. Intermediate cutoffs account for lacunarity due to gaps in the multiscale hierarchy, created by a hiatus of modes associated with discrete ranges of scales. In this dissertation, I investigate the effects of domain and support scales on the multiscale properties of random fractal fields characterized by a power variogram using real and synthetic data. Neuman [1994] and Di Federico and Neuman [1997] have concluded empirically, on the basis of hydraulic conductivity data from many sites, that a finite window of length-scale L filters out (truncates) all modes having integral scales λ larger than λ = μL where μ ≃ 1/3. I confii in their finding computationally by generating truncated fBm realizations on a large grid, using various initial values of μ, and demonstrating that μ ≃ 1/3 for windows smaller than the original grid. My synthetic experiments also show that generating an fl3m realization on a finite grid using a truncated power variogram yields sample variograms that are more consistent with theory than those obtained when the realization is generated using a power variogram. Interpreting sample data from such a realization using wavelet analysis yields more reliable estimates of the Hurst coefficient than those obtained when one employs variogram analysis. Di Federico et al. [1997] developed expressions for the equivalent hydraulic conductivity of a box-shaped support volume, embedded in a log-hydraulic conductivity field characterized by a power variogram, under the action of a mean uniform hydraulic gradient. I demonstrate that their expression and empirically derived value of μ ≃ 1/3 are consistent with a pronounced permeability scale effect observed in unsaturated fractured tuff at the Apache Leap Research Site (ALRS) near Superior, Arizona. I then investigate the compatibility of single-hole air permeability data, obtained at the ALRS on a nominal support scale of about 1 m, with various scaling models including fBm, fGn (fractional Gaussian noise), fLm (fractional Lévy motion), bfLm (bounded fractional Lévy motion) and UM (Universal Multifractals). I find that the data have a Lévy-like distribution at small lags but become Gaussian as the lag increases (corresponding to bfLm). Though this implies multiple scaling, it is not consistent with the UM model, which considers a unique distribution. If one nevertheless applies a UM model to the data, one obtains a very small codimension which suggests that multiple scaling is of minor consequence (applying the UM model to permeability rather than log-permeability data yields a larger codimension but is otherwise not consistent with these data). Variogram and resealed range analyses of the log-permeability data yield comparable estimates of the Hurst coefficient. Resealed range analysis shows that the data are not compatible with an fGn model. I conclude that the data are represented most closely by a truncated fBm model.

Thesis (MScEng (Mathematical Sciences. Applied Mathematics))--University of Stellenbosch, 2006. / This work entails the improvement of an existing three dimensional pore-scale model. Stagnant zones are included, the closure of the volume averaged pressure gradient is improved and an improved calculation of pore-scale averages, using the RUC, is done for the model to be a more realistic representative of the REV and thus of the foamlike material. Both the Darcy and the Forchheimer regimes are modelled and a general momentum transport equation is derived by means of an asymptotic matching technique. The RUC model is also extended to cover non-Newtonian flow. Since metallic foams are generally of porosities greater than 90%, emphasis is put on the accurate prediction of permeability for these porosities. In order to improve permeability predictions for these high porosity cases an adaptation to the RUC model was considered, whereby rectangular prisms were replaced by cylinders. Although this adaptation appears to give more accurate permeabilities at very high porosities, its implementation in a generalised model seems impractical. The prediction of the characteristic RUC side length is discussed and results of both the cylindrical strand model and the square strand model are compared to experimental work.

Not available / text

Not available / text

Not available / text

Under consideration is the finite-size scaling of two thermomechanical responses of random heterogeneous materials. Stochastic mechanics is applied here to the modeling of heterogeneous materials in order to construct the constitutive relations. Such relations (e.g. Hooke's Law in elasticity or Fourier's Law in heat transfer) are well-established under spatial homogeneity assumption of continuum mechanics, where the Representative Volume Element (RVE) is the fundamental concept. The key question is what is the size L of RVE? According to the separation of scales assumption, L must be bounded according to d&lt;L&lt;&lt;LMacro where d is the microscale (or average size of heterogeneity), and LMacro is the macroscale of a continuum mechanics problem. Statistically, for spatially ergodic heterogeneous materials, when the mesoscale is equal to or bigger than the scale of the RVE, the elements of the material can be considered homogenized. In order to attain the said homogenization, two conditions must be satisfied: (a) the microstructure's statistics must be spatially homogeneous and ergodic; and (b) the material's effective constitutive response must be the same under uniform boundary conditions of essential (Dirichlet) and natural (Neumann) types. / In the first part of this work, the finite-size scaling trend to RVE of the Darcy law for Stokesian flow is studied for the case of random porous media, without invoking any periodic structure assumptions, but only assuming the microstructure's statistics to be spatially homogeneous and ergodic. By analogy to the existing methodology in thermomechanics of solid random media, the Hill-Mandel condition for the Darcy flow velocity and pressure gradient fields was first formulated. Under uniform essential and natural boundary conditions, two variational principles are developed based on minimum potential energy and complementary energy. Then, the partitioning method was applied, leading to scale dependent hierarchies on effective (RVE level) permeability. The proof shows that the ensemble average of permeability has an upper bound under essential boundary conditions and a lower bound under uniform natural boundary conditions. / To quantitatively assess the scaling convergence towards the RVE, these hierarchical trends were numerically obtained for various porosities of random disk systems, where the disk centers were generated by a planar Poisson process with inhibition. Overall, the results showed that the higher the density of random disks---or, equivalently, the narrower the micro-channels in the system---the smaller the size of RVE pertaining to the Darcy law. / In the second part of this work, the finite-size scaling of effective thermoelastic properties of random microstructures were considered from Statistical to Representative Volume Element (RVE). Similarly, under the assumption that the microstructure's statistics are spatially homogeneous and ergodic, the SVE is set-up on a mesoscale, i.e. any scale finite relative to the microstructural length scale. The Hill condition generalized to thermoelasticity dictates uniform essential and natural boundary conditions, which, with the help of two variational principles, led to scale dependent hierarchies of mesoscale bounds on effective (RVE level) properties: thermal expansion strain coefficient and stress coefficient, effective stiffness, and specific heats. Due to the presence of a non-quadratic term in the energy formulas, the mesoscale bounds for the thermal expansion are more complicated than those for the stiffness tensor and the heat capacity. To quantitatively assess the scaling trend towards the RVE, the hierarchies are computed for a planar matrix-inclusion composite, with inclusions (of circular disk shape) located at points of a planar, hard-core Poisson point field. Overall, while the RVE is attained exactly on scales infinitely large relative to microscale, depending on the microstructural parameters, the random fluctuations in the SVE response become very weak on scales an order of magnitude larger than the microscale, thus already approximating the RVE. / Based on the above studies, further work on homogenization of heterogeneous materials is outlined at the end of the thesis. / Keywords: Representative Volume Element (RVE), heterogeneous media, permeability, thermal expansion, mesoscale, microstructure.

Experimental studies on infiltration/soil-water movement processes are vital to better understanding movement of soil-water in the vadose zone. The objective of this experimental research was to investigate infiltration/soil-water movement processes utilizing laboratory experiments and computer modeling. Small scale laboratory soil box infiltration experiments were conducted and utilized for the improved parameterization of the Green-Ampt (GA) saturated moisture content parameter to produce an effective moisture content parameter (Be) for utilization in a modified GA model. By incorporating ⊖e values into GA modeling, modeling results showed greatly improved wetting front prediction across different soil conditions. A new soil packing method was proposed for replicating complex microtopographical surfaces with uniform bulk densities in laboratory soil box experiments which proved efficient and effective at accomplishing both objectives. A rainfall simulator and an instantaneous-profile laser scanner were used to simulate rainfall and quantify surface microtopography for experiments. The results clearly show the effect of microtopography on infiltration and soil-water movement characteristics. This offers valuable insight into infiltration/soil-water movement processes as affected by different soil and surface microtopographic conditions. / National Science Foundation (Grant No. EAR-0907588)

Thesis (M.S.)--Massachusetts Institute of Technology, Dept. of Civil Engineering, 1980. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND ENGINEERING. / Bibliography: leaves 150-155. / by Paul Christopher Damian Milly. / M.S.