Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

Implications of specimen preparation and of surface contamination for the measurement of the grain boundary carbon concentration of steels using x-ray microanalysis in an UHVFESTEM

The purpose of the present investigation was to gain an understanding of the nature of the carbon contamination on the surface of standard steel transmission electron spectroscopy (TEM) specimens, the effect of exposure of a clean specimen to normal laboratory air, and the efficacy of plasma-cleaning treatments. This knowledge is a necessary prerequisite to the development of appropriate specimen preparation and/or specimen cleaning methods. X-ray photoelectron spectroscopy in combination with argon ion beam profiling was used to characterize the specimen surfaces of X65 steel and 316 stainless steel. The only clean carbon-free surface obtained was that during argon etching of the sample in the surface analysis chamber. Any exposure of a previously cleaned sample to laboratory air resulted in a rapid carbon (hydrocarbon) contamination of the sample surface and the development of surface oxidation, Plasma cleaning with subsequent exposure of the specimen to the laboratory air also resulted in a carbon-contaminated surface. This suggests that procedures of preparation of TEM specimens of steels outside an ultrahigh vacuum chamber are unlikely to result in the lowering of contamination rates on specimens to levels where measurements for carbon in the grain boundaries are possible. What is needed is a cleaning system as an integral part of the specimen insertion system into the field-emission scanning transmission electron microscope. This cleaning could be carried out by argon ion etching. Copyright (C) 2000 John Wiley & Sons, Ltd.

Beach water table fluctuations due to spring-neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection-diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring-neap tides (a bichromatic tidal system with the fundamental frequencies wt and wt) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency omega(1) - omega(2), i.e., the spring-neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations. (C) 2000 Elsevier Science Ltd. All rights reserved.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

The q-deformed supersymmetric t-J model with a boundary

The q-deformed supersymmetric t-J model on a semi-infinite lattice is diagonalized by using the level-one vertex operators of the quantum affine superalgebra U-q[sl(2\1)]. We. give the bosonization of the boundary states. We give an integral expression for the correlation functions of the boundary model, and derive the difference equations which they satisfy.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Algebraic Bethe ansatz for integrable one-dimensional extended Hubbard models with open boundary conditions

Nine classes of integrable open boundary conditions, further extending the one-dimensional U-q (gl (212)) extended Hubbard model, have been constructed previously by means of the boundary Z(2)-graded quantum inverse scattering method. The boundary systems are now solved by using the algebraic Bethe ansatz method, and the Bethe ansatz equations are obtained for all nine cases.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Apparent strain localization and shear wave dispersion in elastic fault gouge with microrotations

Shear deformation of fault gouge or other particulate materials often results in observed strain localization, or more precisely, the localization of measured deformation gradients. In conventional elastic materials the strain localization cannot take place therefore this phenomenon is attributed to special types of non-elastic constitutive behaviour. For particulate materials however the Cosserat continuum which takes care of microrotations independent of displacements is a more appropriate model. In elastic Cosserat continuum the localization in displacement gradients is possible under some combinations of the generalized Cosserat elastic moduli. The same combinations of parameters also correspond to a considerable dispersion in shear wave propagation which can be used for independent experimental verification of the proposed mechanism of apparent strain localization in fault gouge.

Evaluation of elastic modulus and hardness of thin films by nanoindentation

Simple equations are proposed for determining elastic modulus and hardness properties of thin films on substrates from nanoindentation experiments. An empirical formulation relates the modulus E and hardness H of the film/substrate bilayer to corresponding material properties of the constituent materials via a power-law relation. Geometrical dependence of E and H is wholly contained in the power-law exponents, expressed here as sigmoidal functions of indenter penetration relative to film thickness. The formulation may be inverted to enable deconvolution of film properties from data on the film/substrate bilayers. Berkovich nanoindentation data for dense oxide and nitride films on silicon substrates are used to validate the equations and to demonstrate the film property deconvolution. Additional data for less dense nitride films are used to illustrate the extent to which film properties may depend on the method of fabrication.

Dibucaine effects on structural and elastic properties of lipid bilayers

In this work we report the interaction effects of the local anesthetic dibucaine (DBC) with lipid patches in model membranes by Atomic Force Microscopy (AFM). Supported lipid bilayers (egg phosphatidylcholine, EPC and dimyristoylphosphatidylcholine, DMPQ were prepared by fusion of unilamellar vesicles on mica and imaged in aqueous media. The AFM images show irregularly distributed and sized EPC patches on mica. On the other hand DMPC formation presents extensive bilayer regions on top of which multibilayer patches are formed. In the presence of DBC we observed a progressive disruption of these patches, but for DMPC bilayers this process occurred more slowly than for EPC. In both cases, phase images show the formation of small structures on the bilayer surface suggesting an effect on the elastic properties of the bilayers when DBC is present. Dynamic surface tension and dilatational surface elasticity measurements of EPC and DMPC monolayers in the presence of DBC by the pendant drop technique were also performed, in order to elucidate these results. The curve of lipid monolayer elasticity versus DBC concentration, for both EPC and DMPC cases, shows a maximum for the surface elasticity modulus at the same concentration where we observed the disruption of the bilayer by AFM. Our results suggest that changes in the local curvature of the bilayer induced by DBC could explain the anesthetic action in membranes. (C) 2008 Elsevier B.V. All rights reserved.

Turbulent spots in boundary layers in a free-piston shock-tunnel flow

In this paper, experiments to detect turbulent spots in the transitional boundary layers, formed on a flat plate in a free-piston shock tunnel how, are reported. Experiments indicate that thin-film heat-transfer gauges are suitable for identifying turbulent-spot activity and can be used to identify parameters such as the convection rate of spots and the intermittency of turbulence.

Collagen and Elastic Content of Abdominal Skin After Surgical Weight Loss

Collapsed skin folds after bariatric weight loss are often managed by plastic procedures, but changes in dermal composition and architecture have rarely been documented. Given the potential consequences on surgical outcome, a prospective histochemical study was designed. The hypothesis was that a deranged dermal fiber pattern would accompany major changes in adipose tissue. Female surgical candidates undergoing postbariatric abdominoplasty (n = 40) and never obese women submitted to control procedures (n = 40) were submitted to double abdominal biopsy, respectively in the epigastrium and hypogastrium. Histomorphometric assessment of collagen and elastic fibers was executed by the Image Analyzer System (Kontron Electronic 300, Zeiss, Germany). Depletion of collagen, but not of elastic fibers, in cases with massive weight loss was confirmed. Changes were somewhat more severe in epigastrium (P = 0.001) than hypogastrium (P = 0.007). Correlation with age did not occur. (1) Patients displayed lax, soft skin lacking sufficient collagen fiber network. (2) Elastic fiber content was not damaged, and was even moderately increased in epigastrium; (3) Preoperative obesity negatively correlated with hypogastric collagen concentration; (4) Future studies should pinpoint the roles of obesity, and especially of massive weight loss, on dermal architecture and response to surgery.

Increased elastic microfibrils and thickening of fibroblastic nuclear lamina in canine cutaneous asthenia

Cutaneous asthenia is a hereditary connective tissue disease, primarily of dogs and cats, resembling Ehlers-Danlos syndrome in man. Collagen dysplasia results in skin hyperextensibility, skin and vessel fragility, and poor wound healing. The purpose of this study was to describe the histological findings in a dog with a collagenopathy consistent with cutaneous asthenia. An 8-month-old crossbreed female dog presented with lacerations and numerous atrophic and irregular scars. The skin was hyperextensible and easily torn by the slightest trauma. Ultrastructurally, the dermis was comprised of elaunin and oxytalan microfibrils. These are immature fibres in which the fibrillar component is increased but elastin is reduced. Collagen fibres were profoundly disorganized. The fibrils had a highly irregular outline and a corroded appearance when viewed in cross-section, and were spiralled and fragmented in a longitudinal view. Dermal fibroblasts displayed a conspicuous thickening of the nuclear lamina. Nuclear lamins form a fibrous nucleoskeletal network of intermediate-sized filaments underlying the inner nuclear membrane. Mutations in lamins or lamin-associated proteins cause a myriad of genetic diseases collectively called laminopathies. Disruption of the nuclear lamina seems to affect chromatin organization and transcriptional regulation of gene expression. A common link among all laminopathies may be a failure of stem cells to regenerate mesenchymal tissue. This could contribute to the connective tissue dysplasia seen in cutaneous asthenia.