An introduction to flow and transport in fractal models of porous media: Part I

This special issue gathers together a number of recent papers on fractal geometry and its applications to the modeling of flow and transport in porous media. The aim is to provide a systematic approach for analyzing the statics and dynamics of fluids in fractal porous media by means of theory, modeling and experimentation. The topics covered include lacunarity analyses of multifractal and natural grayscale patterns, random packing's of self-similar pore/particle size distributions, Darcian and non-Darcian hydraulic flows, diffusion within fractals, models for the permeability and thermal conductivity of fractal porous media and hydrophobicity and surface erosion properties of fractal structures.

Are power laws that estimate fractal dimension a good descriptor of soil structure and its link to soil biological properties?

The study of interrelationships between soil structure and its functional properties is complicated by the fact that the quantitative description of soil structure is challenging. Soil scientists have tackled this challenge by taking advantage of approaches such as fractal geometry, which describes soil architectural complexity through a scaling exponent (D) relating mass and numbers of particles/aggregates to particle/aggregate size. Typically, soil biologists use empirical indices such as mean weight diameters (MWD) and percent of water stable aggregates (WSA), or the entire size distribution, and they have successfully related these indices to key soil features such as C and N dynamics and biological promoters of soil structure. Here, we focused on D, WSA and MWD and we tested whether: D estimated by the exponent of the power law of number-size distributions is a good and consistent correlate of MWD and WSA; D carries information that differs from MWD and WSA; the fraction of variation in D that is uncorrelated with MWD and WSA is related to soil chemical and biological properties that are thought to establish interdependence with soil structure (e.g., organic C, N, arbuscular mycorrhizal fungi). We analysed observational data from a broad scale field study and results from a greenhouse experiment where arbuscular mycorrhizal fungi (AMF) and collembola altered soil structure. We were able to develop empirical models that account for a highly significant and large portion of the correlation observed between WSA and MWD but we did not uncover the mechanisms that underlie this correlation. We conclude that most of the covariance between D and soil biotic (AMF, plant roots) and abiotic (C. N) properties can be accounted for by WSA and MWD. This result implies that the ecological effects of the fragmentation properties described by D and generally discussed under the framework of fractal models can be interpreted under the intuitive perspective of simpler indices and we suggest that the biotic components mostly impacted the largest size fractions, which dominate MWD, WSA and the scaling exponent ruling number-size distributions. (C) 2010 Elsevier Ltd. All rights reserved.

Uses and abuses of fractal methodology in ecology

Fractals have found widespread application in a range of scientific fields, including ecology. This rapid growth has produced substantial new insights, but has also spawned confusion and a host of methodological problems. In this paper, we review the value of fractal methods, in particular for applications to spatial ecology, and outline potential pitfalls. Methods for measuring fractals in nature and generating fractal patterns for use in modelling are surveyed. We stress the limitations and the strengths of fractal models. Strictly speaking, no ecological pattern can be truly fractal, but fractal methods may nonetheless provide the most efficient tool available for describing and predicting ecological patterns at multiple scales.

Robust regression applied to fractal/multifractal analysis.

Fractal and multifractal are concepts that have grown increasingly popular in recent years in the soil analysis, along with the development of fractal models. One of the common steps is to calculate the slope of a linear fit commonly using least squares method. This shouldn?t be a special problem, however, in many situations using experimental data the researcher has to select the range of scales at which is going to work neglecting the rest of points to achieve the best linearity that in this type of analysis is necessary. Robust regression is a form of regression analysis designed to circumvent some limitations of traditional parametric and non-parametric methods. In this method we don?t have to assume that the outlier point is simply an extreme observation drawn from the tail of a normal distribution not compromising the validity of the regression results. In this work we have evaluated the capacity of robust regression to select the points in the experimental data used trying to avoid subjective choices. Based on this analysis we have developed a new work methodology that implies two basic steps: ? Evaluation of the improvement of linear fitting when consecutive points are eliminated based on R pvalue. In this way we consider the implications of reducing the number of points. ? Evaluation of the significance of slope difference between fitting with the two extremes points and fitted with the available points. We compare the results applying this methodology and the common used least squares one. The data selected for these comparisons are coming from experimental soil roughness transect and simulated based on middle point displacement method adding tendencies and noise. The results are discussed indicating the advantages and disadvantages of each methodology.

花岗岩风化土微观结构及其工程地质特性研究

Saprolite is the residual soil resulted from completely weathered or highly weathered granite and with corestones of parent rock. It is widely distributed in Hong Kong. Slope instability usually happens in this layer of residual soil and thus it is very important to study the engineering geological properties of Saprolite. Due to the relic granitic texture, the deformation and strength characteristics of Saprolite are very different from normal residual soils. In order to investigate the effects of the special microstructure on soil deformation and strength, a series of physical, chemical and mechanical tests were conducted on Saprolite at Kowloon, Hong Kong. The tests include chemical analysis, particle size analysis, mineral composition analysis, mercury injection, consolidation test, direct shear test, triaxial shear test, optical analysis, SEM & TEM analysis, and triaxial shear tests under real-time CT monitoring.Based on the testing results, intensity and degree of weathering were classified, factors affecting and controlling the deformation and strength of Saprolite were identified, and the interaction between those factors were analyzed.The major parameters describing soil microstructure were introduced mainly based on optical thin section analysis results. These parameters are of importance and physical meaning to describe particle shape, particle size distribution (PSD), and for numerical modeling of soil microstructure. A few parameters to depict particle geometry were proposed or improved. These parameters can be used to regenerate the particle shape and its distribution. Fractal dimension of particle shape was proposed to describe irregularity of particle shapes and capacity of space filling quantitatively. And the effect of fractal dimension of particle shape on soil strength was analyzed. At the same time, structural coefficient - a combined parameter which can quantify the overall microstructure of rock or soil was introduced to study Saprolite and the results are very positive. The study emphasized on the fractal characteristics of PSD and pore structure by applying fractal theory and method. With the results from thin section analysis and mercury injection, it was shown that at least two fractal dimensions Dfl(DB) and Df2 (Dw), exist for both PSD and pore structure. The reasons and physical meanings behind multi-fractal dimensions were analyzed. The fractal dimensions were used to calculate the formation depth and weathering rate of granite at Kowloon. As practical applications, correlations and mathematical models for fractal dimensions and engineering properties of soil were established. The correlation between fractal dimensions and mechanical properties of soil shows that the internal friction angle is mainly governed by Dfl 9 corresponding to coarse grain components, while the cohesion depends on Df2 , corresponding to fine grain components. The correlations between the fractal dimension, friction angle and cohesion are positive linear.Fractal models of PSD and pore size distribution were derived theoretically. Fragmentation mechanism of grains was also analyzed from the viewpoint of fractal. A simple function was derived to define the theoretical relationship between the water characteristic curve (WCC) and fractal dimension, based on a number of classical WCC models. This relationship provides a new analytical tool and research method for hydraulic properties in porous media and solute transportation. It also endues fractal dimensions with new physical meanings and facilitates applications of fractal dimensions in water retention characteristics, ground water movement, and environmental engineering.Based on the conclusions from the fractal characteristics of Saprolite, size effect on strength was expressed by fractal dimension. This function is in complete agreement with classical Weibull model and a simple function was derived to represent the relationship between them.In this thesis, the phenomenon of multi-fractal dimensions was theoretically analyzed and verified with WCC and saprolite PSD results, it was then concluded that multi-fractal can describe the characteristics of one object more accurately, compared to single fractal dimension. The multi-fractal of saprolite reflects its structural heterogeneity and changeable stress environment during the evolution history.

Compositional analysis for an unbiased measure of soil aggregation

Soil aggregation is an index of soil structure measured by mean weight diameter (MWD) or scaling factors often interpreted as fragmentation fractal dimensions (D-f). However, the MWD provides a biased estimate of soil aggregation due to spurious correlations among aggregate-size fractions and scale-dependency. The scale-invariant D-f is based on weak assumptions to allow particle counts and sensitive to the selection of the fractal domain, and may frequently exceed a value of 3, implying that D-f is a biased estimate of aggregation. Aggregation indices based on mass may be computed without bias using compositional analysis techniques. Our objective was to elaborate compositional indices of soil aggregation and to compare them to MWD and D-f using a published dataset describing the effect of 7 cropping systems on aggregation. Six aggregate-size fractions were arranged into a sequence of D-1 balances of building blocks that portray the process of soil aggregation. Isometric log-ratios (ilrs) are scale-invariant and orthogonal log contrasts or balances that possess the Euclidean geometry necessary to compute a distance between any two aggregation states, known as the Aitchison distance (A(x,y)). Close correlations (r>0.98) were observed between MWD, D-f, and the ilr when contrasting large and small aggregate sizes. Several unbiased embedded ilrs can characterize the heterogeneous nature of soil aggregates and be related to soil properties or functions. Soil bulk density and penetrater resistance were closely related to A(x,y) with reference to bare fallow. The A(x,y) is easy to implement as unbiased index of soil aggregation using standard sieving methods and may allow comparisons between studies. (C) 2012 Elsevier B.V. All rights reserved.

A study of habitat fragmentation in Southeastern Brazil using remote sensing and geographic information systems (GIS)

The purpose of this work was to study fragmentation of forest formations (mesophytic forest, riparian woodland and savannah vegetation (cerrado)) in a 15,774-ha study area located in the Municipal District of Botucatu in Southeastern Brazil (São Paulo State). A land use and land cover map was made from a color composition of a Landsat-5 thematic mapper (TM) image. The edge effect caused by habitat fragmentation was assessed by overlaying, on a geographic information system (GIS), the land use and land cover data with the spectral ratio. The degree of habitat fragmentation was analyzed by deriving: 1. mean patch area and perimeter; 2. patch number and density; 3. perimeter-area ratio, fractal dimension (D), and shape diversity index (SI); and 4. distance between patches and dispersion index (R). In addition, the following relationships were modeled: 1. distribution of natural vegetation patch sizes; 2. perimeter-area relationship and the number and area of natural vegetation patches; 3. edge effect caused by habitat fragmentation, the values of R indicated that savannah patches (R = 0.86) were aggregated while patches of natural vegetation as a whole (R = 1.02) were randomly dispersed in the landscape. There was a high frequency of small patches in the landscape whereas large patches were rare. In the perimeter-area relationship, there was no sign of scale distinction in the patch shapes, In the patch number-landscape area relationship, D, though apparently scale-dependent, tends to be constant as area increases. This phenomenon was correlated with the tendency to reach a constant density as the working scale was increased, on the edge effect analysis, the edge-center distance was properly estimated by a model in which the edge-center distance was considered a function of the to;al patch area and the SI. (C) 1997 Elsevier B.V. B.V.

Toward mass production of microtextured microdevices: linking rapid prototyping with microinjection molding

The possibility of manufacturing textured materials and devices, with surface properties controlled from the design stage, instead of being the result of machining processes or chemical attacks, is a key factor for the incorporation of advanced functionalities to a wide set of micro and nanosystems. Recently developed high-precision additive manufacturing technologies, together with the use of fractal models linked to computer-aided design tools, allow for a precise definition and control of final surface properties for a wide set of applications, although the production of larger series based on these resources is still an unsolved challenge. However, rapid prototypes, with controlled surface topography, can be used as original masters for obtaining micromold inserts for final large-scale series manufacture of replicas using microinjection molding. In this study, an original procedure is presented, aimed at connecting rapid prototyping with microinjection molding, for the mass production of two different microtextured microsystems, linked to tissue engineering tasks, using different thermoplastics as ultimate materials.

On the generative equations of fractal self-similarity in granular media and the related PSD models

The physical appearance of granular media suggests the existence of geometrical scale invariance. The paper discuss how this physico-empirical property can be mathematically encoded leading to different generative models: a smooth one encoded by a differential equation and another encoded by an equation coming from a measure theoretical property.

Porosity imaged by a vector projection algorithm correlates with fractal dimension measured on 3D models obtained by microCT

International audience

A quasi-likelihood method for fractal-dimension estimation

We propose a simple method of constructing quasi-likelihood functions for dependent data based on conditional-mean-variance relationships, and apply the method to estimating the fractal dimension from box-counting data. Simulation studies were carried out to compare this method with the traditional methods. We also applied this technique to real data from fishing grounds in the Gulf of Carpentaria, Australia

Ordering near the percolation threshold in models of two-dimensional interacting bosons with quenched dilution

Randomly diluted quantum boson and spin models in two dimensions combine the physics of classical percolation with the well-known dimensionality dependence of ordering in quantum lattice models. This combination is rather subtle for models that order in two dimensions but have no true order in one dimension, as the percolation cluster near threshold is a fractal of dimension between 1 and 2: two experimentally relevant examples are the O(2) quantum rotor and the Heisenberg antiferromagnet. We study two analytic descriptions of the O(2) quantum rotor near the percolation threshold. First a spin-wave expansion is shown to predict long-ranged order, but there are statistically rare points on the cluster that violate the standard assumptions of spin-wave theory. A real-space renormalization group (RSRG) approach is then used to understand how these rare points modify ordering of the O(2) rotor. A new class of fixed points of the RSRG equations for disordered one-dimensional bosons is identified and shown to support the existence of long-range order on the percolation backbone in two dimensions. These results are relevant to experiments on bosons in optical lattices and superconducting arrays, and also (qualitatively) for the diluted Heisenberg antiferromagnet La-2(Zn,Mg)(x)Cu1-xO4.

Fractal boundaries of basin of attraction of Newton-Raphson method in helicopter trim

A key problem in helicopter aeroelastic analysis is the enormous computational time required for a numerical solution of the nonlinear system of algebraic equations required for trim, particularly when free wake models are used. Trim requires calculation of the main rotor and tail rotor controls and the vehicle attitude which leads to the six steady forces and moments about the helicopter center of gravity to be zero. An appropriate initial estimate of the trim state is needed for successful helicopter trim. This study aims to determine the control inputs that can have considerable effect on the convergence of trim solution in the aeroelastic analysis of helicopter rotors by investigating the basin of attraction of the nonlinear equations (set of initial guess points from which the nonlinear equations converge). It is illustrated that the three main rotor pitch controls of collective pitch, longitudinal cyclic pitch and lateral cyclic pitch have a significant contribution to the convergence of the trim solution. Trajectories of the Newton iterates are shown and some ideas for accelerating the convergence of a trim solution in the aeroelastic analysis of helicopters are proposed. It is found that the basins of attraction can have fractal boundaries. (C) 2010 Elsevier Ltd. All rights reserved.

Fractal dimension method for delamination detection in composite beams

A new delaminated composite beam element is formulated for Timoshenko as well as Euler-Bernoulli beam models. Shape functions are derived from Timoshenko functions; this provides a unified formulation for slender to moderately deep beam analyses. The element is simple and easy to implement, results are on par with those from free mode delamination models. Katz fractal dimension method is applied on the mode shapes obtained from finite element models, to detect the delamination in the beam. The effect of finite element size on fractal dimension method of delamination detection is quantified.