A physical perspective of the length scales in gradient elasticity through the prism of wave dispersion

The goal of this study is to understand the physical meaning and evaluate the intrinsic length scale parameters, featured in the theories of gradient elasticity, by deploying the analytical treatment and experimental measurements of the dispersion of elastic waves. The developments are focused on examining the propagation of longitudinal waves in an aluminum rod with periodically varying cross-section. First, the analytical solution for the dispersion relationship, based on the periodic cell analysis of a bi-layered laminate and Bloch theorem, is compared to two competing models of gradient elasticity. It is shown that the customary gradient elastic model with two length-scale parameters is able to capture the dispersion accurately up to the beginning of the first band gap. On the other hand, the gradient elastic model with an additional length scale (affiliated with the fourth-order time derivative in the field equation) is shown to capture not only the first dispersion branch before the band gap, but also the band gap itself and the preponderance of the second branch. Closed form relations between the microstructure parameters and the intrinsic length scales are obtained for both gradient elasticity models. By way of the asymptotic treatment in the limit of a weak contrast between the laminae, a clear physical meaning and scaling of the length-scale parameters was established in terms of: (i) the microstructure (given by the size of the unit cell and the contrast between the laminae), and (ii) thus induced dispersion relationship (characterized by the location and the width of the band gap). The analysis is verified through an experimental observation of wave dispersion, and wave attenuation within the band gap. A comparison between the analytical treatment, the gradient elastic model with three intrinsic length scales, and experimental measurements demonstrates a good agreement over the range of frequencies considered.     

Influence of Internal Structure and Medium Length on Transport and Deposition of Suspended Particles: A Laboratory Study

A laboratory study was undertaken on the transport and the deposition of suspended particles (silt of modal diametre 6 μm) in three columns of different length, filled with glass beads or gravel. Tracer tests were carried out at various flow velocities by short pulses of a mixture of suspended particles/dissolved tracer. The breakthrough curves were competently described with the analytical solution of a convection dispersion equation with a first-order deposition rate and the hydro-dispersive parameters were deduced. For the same experimental conditions, the results showed a difference in the behaviour of the suspended particles transport and deposition rates within the two porous media tested. The internal structure of both media governs the particle-grain collision frequency as well as the particles trapping. The scale effect was highlighted and affects the dispersivity, the size exclusion effect, the recovery rates and the deposition rates. Longitudinal dispersion increases with mean pore velocity and is described with a nonlinear relationship. The dispersivity increases with the column length. The size exclusion effect is more important in the short column. The recovery rate increases with flow velocity and decreases while increasing column length. The deposition rates increases until a critical flow velocity then decreases. This critical velocity is also sensitive to the scale effect, and increases with the column length.     

Single- and Dual-domain Models of Solute Transport in Alluvial Sediments: the Effects of Heterogeneity Structure and Spatial Scale

Fine-scale heterogeneity of alluvial aquifers controls solute transport in groundwater at the scales relevant for practical applications: the architecture of sedimentary structures might create preferential flow paths (PFPs) or hydraulic barriers, which affect the breakthrough curves (BTCs). Objective of this paper was the assessment of the relevance of single- and dual-domain models for different heterogeneity patterns and scale lengths in alluvial sediments. Three case studies have been analysed with a classical single-domain model (SDM) and with three dual-domain models (DDMs): a dual-porosity model (DPorM) and two dual-permeability models (DPerM), which differ for the presence or the absence of solute exchange between the two domains. The first case study includes numerical tracer tests in metre-scale blocks of alluvial sediments; the second is a laboratory experiment of tracer injection in a decimetre-scale column of homogeneous sand; the third is a field tracer test performed at hectometre scale at the Cape Cod site. The relevance of the solute exchange in the DDMs is analysed with the characteristic advection and exchange times and with the Péclet and Damköhler numbers. The SDM is satisfactory for alluvial sediments with unstructured heterogeneity. The uncoupled DPerM is shown to be a better approach than the DPorM in sediments with PFPs; in this case, the coupled DPerM does not improve significantly the results of the uncoupled DPerM. A minor difference between the results of the three DDMs is observed for sediments in which the non-Fickian behaviour is not clearly determined by the presence of PFPs.     

A theory of macrodispersion for the scale-up problem

Dispersion is the result, observable on large length scales, of events which are random on small length scales. When the length scale on which the randomness operates is not small, relative to the observations, then classical dispersion theory fails. The scale up problem refers to situations in which randomness occurs on all length scales, and for which classical dispersion theory necessarily fails. The purpose of this article is to present non-Fickian, theories of dispersion, which do not assume a scale separation between the randomness and the observed consequences, and which do not assume a single length scale.Porous media flow properties are heterogeneous on all length scales. The geological variation on length scales below the observational length scale can be regarded as unknown and unknowable, and thus as a random variable.We develop a systematic theory relating scaling behavior of the geological heterogeneity to the scaling behavior of the fluid dispersivity. Three qualitatively distinct regimes (Fickian, non-Fickian and nonrenormalizable) are found. The theory gives consistent answers within several distinct analytic approximations, and with numerical simulation of the equations of porous media flow.Comparison to field data is made. The use of Kriging to generate constrained ensembles for conditional simulation is discussed.     

Periodic structure of the dispersed phase in a forced jet and their effects on the particle dispersion

This paper presents an experimental study on an acoustically forced particle-laden jet. The acoustic disturbances cause a train of strong periodical vortex due to the selected frequency and its high excitation level. The jet Reynolds number is not small (Re = 11,900) and the particle Stokes number is about one, responding partly to the forcing. The flow was tested using a Phase Doppler Anemometer (PDA). The paper includes measurements of both gas and dispersed phase over the whole forcing cycle. An external post-processing (developed by the authors) carefully corrects the bias inherent to the operation principles of the PDA in all supplied averages (including the phase-averaged values). This post-processing gives also some variables which were defined ad-hoc to characterize the periodic structure of the flow. Such information is never given in the previous literature. This work continues a previous study done by the authors.Measurements detect three axial zones. The strong periodic gas vortices control the flow in the area close to the nozzle exit. They generate highly concentrated clusters of particles as well as tongue-shaped structures of radially ejected particles (or radial streaks). Downstream, the gas vortices vanish and inertia plays a central role in the development of the dispersed phase. The particle clustering ends here. Finally, all periodic motion disappears and flow degenerates into an unforced two-phase jet. Radially, the inertial zone of the particulate phase covers the outermost layers. The influence of the particle size is also discussed.The radial dispersion of particles across certain section is quantified by means of a suitably defined parameter. This dispersion radius was measured at the end of the area disturbed by forcing for both the forced and unforced jet. Thus, the comparison assesses the total effect of forcing on the transversal dispersion. The dispersion of the whole size distribution and of each particle size is quantified. Results show that forcing enhances the dispersion and it is controlled mainly by the periodic streaks while turbulence has a secondary role. The streak shape is accurately computed from the measurements and its extension has been successfully related with the particle's history and it size by means of a suitably defined Stokes number.Finally, this study supplies a set of high quality data useful to validate inherently unsteady numerical models. As stated by other authors, there is a lack of periodic well-characterized experiments for validation purposes which mimic the interaction between the particles and the large scales of turbulence.     

Concentration fluctuations and dilution in two-dimensionally periodic heterogeneous porous media

Dilution of solute in two-dimensionally periodic heterogeneous porous media is assessed by numerically simulating advection-dispersion. The concentration fluctuations, created by advective heterogeneity, are destroyed by local dispersion, over a characteristic variance residence time (VRT). For an impulse introduction of solute, initially, plumes become increasingly irregular with time—the coefficient of variation (CV) of concentration grows with time. A priori, the spatial second moment and mean concentrations are insufficient measures of dilution, because concentration fluctuations can be large. At large times (t > VRT) the relative concentration fluctuations weaken—the concentration CV was observed to slowly decrease with time. At the center of mass the general trend of the decreasing CV follows VRT/t (predicted by Kapoor and Gelhar). The VRT is found to be an increasing function of the log hydraulic conductivity microscale. In employing effective dispersion coefficents to model the mean concentration, it needs to be recognized that excursions of concentrations around the mean are singularly determined by local dispersion.     

An experimental investigation of dispersion in layered porous media

Experiments were run in three linear, homogeneous, nonuniform porous media constructed in lucite columns using spherical glass beads. The columns were also joined end to end to create an in series layered heterogeneous porous media. Each column, all combinations of columns and several permutations were studied with a factorial experimental design to determine the effects of porosity, permeability, velocity, length, and column order upon dispersion. Attempts to predict the heterogeneous results from the homogeneous results were made, and a statistical regression based on the factorial design was calculated. Results showed that no simple averaging procedure accurately predicted the heterogeneous results. The statistical regression showed permeability, velocity, viscosity, length and column order to be significant.     

Large eddy simulation of free surface shallow‐water flow

Shallow‐water flow with free surface frequently occurs in ambient water bodies, in which the horizontal scale of motion is generally two orders of magnitude greater than the water depth. To accurately predict this flow phenomenon in more detail, a three‐dimensional numerical model incorporating the method of large eddy simulation (LES) has been developed and assessed. The governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by the QUICKEST scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly using the Gauss–Seidel iteration method. The model has been applied to free surface channel flow for which ample experimental data are available for verification. The inflow boundary condition for turbulence is generated by a spectral line processor. The computed results compare favourably with the experimental data and those results obtained by using a periodic boundary condition. The performance of the model is also assessed for the case in which anisotropic grids and filters with horizontal grid size of the order of the water depth are used for computational efficiency. The coarse horizontal grid was found to cause a significant reduction in the large‐scale turbulent motion generated by the bottom turbulence, and the turbulent motion is predominately described by the sub‐grid scale (SGS) terms. The use of the Smagorinsky model for SGS turbulence in this situation is found inappropriate. A parabolic mixing length model, which accounts for the filtered turbulence, is then proposed. The new model can reproduce more accurately the flow quantities. Copyright © 2000 John Wiley & Sons, Ltd.     

On the Parametric Resonance of Columns Carrying Concentrated Masses

Abstract

The effect of various parameters upon the region of dynamic instability of a uniform simply supported column carrying n concentrated masses and subjected to an axial periodic force at one end is presented. The effect of axial inertia of the mass per unit length of the column is also included. This problem is reduced to a coupled system of two homogeneous partial differential equations of second order with periodic coefficients which can be written in the form of a matric differential equation of the Mathieu-Hill type. A solution methodology is developed and successfully demonstrated through numerical examples.     

A size-dependent Reddy–Levinson beam model based on a strain gradient elasticity theory

A size-dependent Reddy–Levinson beam model is developed based on a strain gradient elasticity theory. Governing equations and boundary conditions are derived by using Hamilton’s principle. The model contains three material length scale parameters, which may effectively capture the size effect in micron or sub-micron. This model can degenerate into the modified couple stress model or even the classical model if two or all material length scale parameters are taken to be zero respectively. In addition, the present model recovers the micro scale Timoshenko and Bernoulli–Euler beam models based on the same strain gradient elasticity theory. To illustrate the new model, the static bending and free vibration problems of a simply supported micro scale Reddy–Levinson beam are solved respectively; the results are compared with the reduced models. Numerical results reveal that the differences in the deflection, rotation and natural frequency predicted by the present model and the other two reduced Reddy–Levinson models are getting larger as the beam thickness is comparable to the material length scale parameters. These differences, however, are decreasing or even diminishing with the increase of the beam thickness. This study may be helpful to characterize the mechanical properties of small scale beam-like structures for a wide range of potential applications.     

Size effect on mechanical behavior of random fiber networks

Bonded random fiber networks are heterogeneous on multiple scales. This leads to a pronounced size effect on their mechanical behavior. In this study we quantify the size effect and determine the minimum model size required to eliminate the size effect for given set of system parameters. These include the network density, the fiber length and the fiber bending and axial stiffness. The results may guide the definition of models and the selection of the size of representative volume elements in sequential multiscale models of fiber networks. To underline the origins of the size effect, we characterize the network heterogeneity by analyzing the geometry of the network (density distribution), the strain field and the strain energy distribution. The dependence of the heterogeneity on the scale of observation and system parameters is discussed.     

A numerical approximation to the elastic properties of sphere-reinforced composites

Three-dimensional cubic unit cells containing 30 non-overlapping identical spheres randomly distributed were generated using a new, modified random sequential adsortion algorithm suitable for particle volume fractions of up to 50%. The elastic constants of the ensemble of spheres embedded in a continuous and isotropic elastic matrix were computed through the finite element analysis of the three-dimensional periodic unit cells, whose size was chosen as a compromise between the minimum size required to obtain accurate results in the statistical sense and the maximum one imposed by the computational cost. Three types of materials were studied: rigid spheres and spherical voids in an elastic matrix and a typical composite made up of glass spheres in an epoxy resin. The moduli obtained for different unit cells showed very little scatter, and the average values obtained from the analysis of four unit cells could be considered very close to the “exact” solution to the problem, in agreement with the results of Drugan and Willis (J. Mech. Phys. Solids 44 (1996) 497) referring to the size of the representative volume element for elastic composites. They were used to assess the accuracy of three classical analytical models: the Mori-Tanaka mean-field analysis, the generalized self-consistent method, and Torquato's third-order approximation.     

A size-dependent Kirchhoff micro-plate model based on strain gradient elasticity theory

A size-dependent Kirchhoff micro-plate model is developed based on the strain gradient elasticity theory. The model contains three material length scale parameters, which may effectively capture the size effect. The model can also degenerate into the modified couple stress plate model or the classical plate model, if two or all of the material length scale parameters are taken to be zero. The static bending, instability and free vibration problems of a rectangular micro-plate with all edges simple supported are carried out to illustrate the applicability of the present size-dependent model. The results are compared with the reduced models. The present model can predict prominent size-dependent normalized stiffness, buckling load, and natural frequency with the reduction of structural size, especially when the plate thickness is on the same order of the material length scale parameter.     

Contribution to the Modelling of Coupled Heat and Mass Transfers on 3D Real Structure of Heterogeneous Building Materials: Application to Hemp Concrete

Hydraulic flow, electrical flow and the passage of elastic waves through porous media are all linked by electrokinetic processes. In its simplest form, the passage of elastic waves through the porous medium causes fluid to flow through that medium and that flow gives rise to an electrical streaming potential and electrical counter-current. These processes are frequency-dependent and governed by coupling coefficients which are themselves frequency-dependent. The link between fluid pressure and fluid flow is described by dynamic permeability, which is characterised by the hydraulic coupling coefficient (C_hp). The link between fluid pressure and electrical streaming potential is characterised by the streaming potential coefficient (C_sp). While the steady-state values of such coefficients are well studied and understood, their frequency dependence is not. Previous work has been confined to unconsolidated and disaggregated materials such as sands, gravels and soils. In this work, we present an apparatus for measuring the hydraulic and streaming potential coefficients of high porosity, high permeability consolidated porous media as a function of frequency. The apparatus operates in the range 1 Hz to 2 kHz with a sample of 10 mm diameter and 5–30 mm in length. The full design and validation of the apparatus are described together with the experimental protocol it uses. Initial data are presented for three samples of Boise sandstone, which present as dispersive media with the critical transition frequency of 918.3?±?99.4 Hz. The in-phase and in-quadrature components of the measured hydraulic and streaming potential coefficients have been compared to the Debye-type dispersion model as well as theoretical models based on bundles of capillary tubes and porous media. Initial results indicate that the dynamic permeability data present an extremely good fit to the capillary bundle and Debye-type dispersion models, while the streaming potential coefficient presents an extremely good fit to all of the models up to the critical transition frequency, but diverges at higher frequencies. The streaming potential coefficient data are best fitted by the Pride model and its Walker and Glover simplification. Characteristic pore size values calculated from the measured critical transition frequency fell within 1.73% of independent measures of this parameter, while the values calculated directly from the Packard model showed an underestimation by about 12%.

     

Multi-Stage Upscaling: Selection of Suitable Methods

Reservoirs are often composed of an assortment of rock types giving rise to permeability heterogeneities at a variety of length-scales. To predict fluid flow at the full-field scale, it is necessary to be aware of these different types of heterogeneity, to recognise which are likely to have important effects on fluid flow, and to capture them by upscaling. In fact, we may require a series of stages of upscaling to go from small-scales (mm or cm) to a full-field model. When there are two (or more) phases present, we also need to know how these heterogeneities interact with fluid forces (capillary, viscous and gravity). We discuss how these effects may be taken into account by upscaling. This study focusses on the effects of steady-state upscaling for viscous-dominated floods and tests carried out on a range of 2D models are described. Upscaling errors are shown to be reduced slightly by the increase in numerical dispersion at the coarse scale. We select a combination of three different upscaling methods, and apply this approach to a model of a North Sea oil reservoir in a deep marine environment. Six different genetic units (rock types) were identified, including channel sandstone and inter-bedded sandstone and mudstone. These units were modelled using different approaches, depending on the nature of the heterogeneities. Our results show that the importance of small-scale heterogeneity depends on the large-scale distribution of the rock types. Upscaling may not be worthwhile in sparsely distributed genetic units. However, it is important in the dominant rock type, especially if there is good connectivity through the unit between the injector wells (or aquifer) and the producer wells.This revised version was published online in May 2005. In the previous version one of the authors name was missing.     

A critical review of dispersion in packed beds

The phenomenon of dispersion (transverse and longitudinal) in packed beds is summarized and reviewed for a great deal of information from the literature. Dispersion plays an important part, for example, in contaminant transport in ground water flows, in miscible displacement of oil and gas and in reactant and product transport in packed bed reactors. There are several variables that must be considered, in the analysis of dispersion in packed beds, like the length of the packed column, viscosity and density of the fluid, ratio of column diameter to particle diameter, ratio of column length to particle diameter, particle size distribution, particle shape, effect of fluid velocity and effect of temperature (or Schmidt number). Empirical correlations are presented for the prediction of the dispersion coefficients (D _T and D _L) over the entire range of practical values of Sc and Pe_m, and works on transverse and longitudinal dispersion of non-Newtonian fluids in packed beds are also considered.     

Long wavelength inner-resonance cut-off frequencies in elastic composite materials

We revisit an ancient paper (Auriault and Bonnet, 1985) which points out the existence of cut-off frequencies for long acoustic wavelength in high-contrast elastic composite materials, i.e. when the wavelength is large with respect to the characteristic heterogeneity length. The separation of scales enables the use of the method of multiple scale expansions for periodic structures, a powerful upscaling technique from the heterogeneity scale to the wavelength scale. However, the results remain valid for non-periodic composite materials which show a Representative Elementary Volume (REV). The paper extends the previous investigations to three-component composite materials made of hard inclusions, coated with a soft material, both of arbitrary geometry, and embedded in a connected stiff material. The equivalent macroscopic models are rigorously established as well as their domains of validity. Provided that the stiffness contrast within the soft and the connected stiff materials is of the order of the squared separation of scales parameter, it is demonstrated (i) that the propagation of long wave may coincide with the resonance frequencies of the hard inclusions/soft material system and (ii) that the macroscopic model presents a series of cut-off frequencies given by an eigenvalue problem for the resonating domain in the cell. These results are illustrated in the case of stratified composites and the possible microstructures of heterogeneous media in which the inner dynamics phenomena may occur are discussed.     

Multiscale asymptotic homogenization for multiphysics problems with multiple spatial and temporal scales: a coupled thermo-viscoelastic example problem

A systematic approach for analyzing multiple physical processes interacting at multiple spatial and temporal scales is developed. The proposed computational framework is applied to the coupled thermo-viscoelastic composites with microscopically periodic mechanical and thermal properties. A rapidly varying spatial and temporal scales are introduced to capture the effects of spatial and temporal fluctuations induced by spatial heterogeneities at diverse time scales. The initial-boundary value problem on the macroscale is derived by using the double scale asymptotic analysis in space and time. It is shown that an extra history-dependent long-term memory term introduced by the homogenization process in space and time can be obtained by solving a first order initial value problem. This is in contrast to the long-term memory term obtained by the classical spatial homogenization, which requires solutions of the initial-boundary value problem in the unit cell domain. The validity limits of the proposed spatial–temporal homogenized solution are established. Numerical example shows a good agreement between the proposed model and the reference solution obtained by using a finite element mesh with element size comparable to that of material heterogeneity.     

Effect of sample size on the response of DEM samples with a realistic grading

This paper shows that for DEM simulations of triaxial tests using samples with a grading that is repre- sentative of a real soil, the sample size significantly influences the observed material response. Four DEM samples with identical initial states were produced： three cylindrical samples bounded by rigid wails and one bounded by a cubical periodic cell, When subjected to triaxial loading, the samples with rigid boundaries were more dilative, stiffer and reached a higher peak stress ratio than the sample enclosed by periodic boundaries. For the rigid-wall samples, dilatancy increased and stiffness decreased with increasing sample size, The periodic sample was effectively homogeneous, The void ratio increased and the contact density decreased close to the rigid walls, This heterogeneity reduced with increasing sample size. The positions of the critical state lines （CSLs） of the overall response in e-log p＇ space were sensitive to the sample size, although no difference was observed between their slopes. The critical states of the interior regions of the rigid-wall-bounded samples approached that of the homogeneous periodic sample with increasing sample size. The ultimate strength of the material at the critical state is independent of sample size.