Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Force-induced globule-coil transition in laminin binding protein and its role for viral-cell membrane fusion.

The specific interactions of the pairs laminin binding protein (LBP)-purified tick-borne encephalitis viral surface protein E and certain recombinant fragments of this protein, as well as West Nile viral surface protein E and certain recombinant fragments of that protein, are studied by combined methods of single-molecule dynamic force spectroscopy (SMDFS), enzyme immunoassay and optical surface waves-based biosensor measurements. The experiments were performed at neutral pH (7.4) and acid pH (5.3) conditions. The data obtained confirm the role of LBP as a cell receptor for two typical viral species of the Flavivirus genus. A comparison of these data with similar data obtained for another cell receptor of this family, namely human αVβ3 integrin, reveals that both these receptors are very important. Studying the specific interaction between the cell receptors in question and specially prepared monoclonal antibodies against them, we could show that both interaction sites involved in the process of virus-cell interaction remain intact at pH 5.3. At the same time, for these acid conditions characteristic for an endosome during flavivirus-cell membrane fusion, SMDFS data reveal the existence of a force-induced (effective already for forces as small as 30-70 pN) sharp globule-coil transition for LBP and LBP-fragments of protein E complexes. We argue that this conformational transformation, being an analog of abrupt first-order phase transition and having similarity with the famous Rayleigh hydrodynamic instability, might be indispensable for the flavivirus-cell membrane fusion process. Copyright © 2014 John Wiley & Sons, Ltd.

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

Analysis of micro- and nano-structures of the corneal surface of Drosophila and its mutants by atomic force microscopy and optical diffraction.

Drosophila melanogaster is a model organism instrumental for numerous biological studies. The compound eye of this insect consists of some eight hundred individual ommatidia or facets, ca. 15 µm in cross-section. Each ommatidium contains eighteen cells including four cone cells secreting the lens material (cornea). High-resolution imaging of the cornea of different insects has demonstrated that each lens is covered by the nipple arrays--small outgrowths of ca. 200 nm in diameter. Here we for the first time utilize atomic force microscopy (AFM) to investigate nipple arrays of the Drosophila lens, achieving an unprecedented visualization of the architecture of these nanostructures. We find by Fourier analysis that the nipple arrays of Drosophila are disordered, and that the seemingly ordered appearance is a consequence of dense packing of the nipples. In contrast, Fourier analysis confirms the visibly ordered nature of the eye microstructures--the individual lenses. This is different in the frizzled mutants of Drosophila, where both Fourier analysis and optical imaging detect disorder in lens packing. AFM reveals intercalations of the lens material between individual lenses in frizzled mutants, providing explanation for this disorder. In contrast, nanostructures of the mutant lens show the same organization as in wild-type flies. Thus, frizzled mutants display abnormal organization of the corneal micro-, but not nano-structures. At the same time, nipples of the mutant flies are shorter than those of the wild-type. We also analyze corneal surface of glossy-appearing eyes overexpressing Wingless--the lipoprotein ligand of Frizzled receptors, and find the catastrophic aberration in nipple arrays, providing experimental evidence in favor of the major anti-reflective function of these insect eye nanostructures. The combination of the easily tractable genetic model organism and robust AFM analysis represents a novel methodology to analyze development and architecture of these surface formations.

Surface reconstruction from microscopic images in optical lithography.

This paper presents a method to reconstruct 3D surfaces of silicon wafers from 2D images of printed circuits taken with a scanning electron microscope. Our reconstruction method combines the physical model of the optical acquisition system with prior knowledge about the shapes of the patterns in the circuit; the result is a shape-from-shading technique with a shape prior. The reconstruction of the surface is formulated as an optimization problem with an objective functional that combines a data-fidelity term on the microscopic image with two prior terms on the surface. The data term models the acquisition system through the irradiance equation characteristic of the microscope; the first prior is a smoothness penalty on the reconstructed surface, and the second prior constrains the shape of the surface to agree with the expected shape of the pattern in the circuit. In order to account for the variability of the manufacturing process, this second prior includes a deformation field that allows a nonlinear elastic deformation between the expected pattern and the reconstructed surface. As a result, the minimization problem has two unknowns, and the reconstruction method provides two outputs: 1) a reconstructed surface and 2) a deformation field. The reconstructed surface is derived from the shading observed in the image and the prior knowledge about the pattern in the circuit, while the deformation field produces a mapping between the expected shape and the reconstructed surface that provides a measure of deviation between the circuit design models and the real manufacturing process.

Determination of chemical composition and its relationship with optical properties of Ti-N and Ti-V-N sputtered thin films

In the investigation of thin films of transition metal nitrides, an essential role is played by the accurate determination of their chemical composition. Actually the chemical composition depends on the deposition parameters and influences the optical properties. These relations are illustrated in thin films of TiNx and (Ti1-yVy)N-x deposited by reactive magnetron sputtering from composite targets of the elements. By variation of the nitrogen partial pressure and the target composition, different samples have been obtained. The chemical composition has been measured by electron probe microanalysis at low irradiation voltages. The optical properties are evaluated by ex-situ ellipsometry. Using the screened Drude model, they are correlated with the differences in composition. Adding vanadium or nitrogen in Ti-N is shown to have the same effect on the optical properties.

Two waves of recombinase gene expression in developing thymocytes.

During T cell development in the thymus, T cell receptor (TCR) alpha, beta, gamma, and delta genes are rearranged and expressed. TCR rearrangement strictly depends upon the coordinate activity of two recombinase activating genes, Rag-1 and Rag-2. In this study we have followed the expression of these genes at different stages of intrathymic development. The results indicate that there are two periods of high Rag-1 and Rag-2 mRNA expression. The first wave peaks early at the CD25+CD4-CD8-CD3- stage of development and coincides with the initial appearance of transcripts derived from fully rearranged TCR beta, gamma, and delta genes, whereas the second wave occurs later at the CD4+CD8+ stage coincident with full-length TCR alpha mRNA expression. Active downregulation of Rag-1 and Rag-2 mRNA expression appears to occur in vivo between the two peaks of recombinase activity. This phenomenon can be mimicked in vitro in response to artificial stimuli such as phorbol myristate acetate and calcium ionophore. Collectively our data suggest that recombinase expression is actively regulated during early thymus development independently of cell surface expression of a mature heterodimeric TCR protein complex.

Surface deterioration of dental materials after simulated toothbrushing in relation to brushing time and load.

OBJECTIVES: (1) To evaluate the changes in surface roughness and gloss after simulated toothbrushing of 9 composite materials and 2 ceramic materials in relation to brushing time and load in vitro; (2) to assess the relationship between surface gloss and surface roughness. METHODS: Eight flat specimens of composite materials (microfilled: Adoro, Filtek Supreme, Heliomolar; microhybrid: Four Seasons, Tetric EvoCeram; hybrid: Compoglass F, Targis, Tetric Ceram; macrohybrid: Grandio), two ceramic materials (IPS d.SIGN and IPS Empress polished) were fabricated according to the manufacturer's instructions and optimally polished with up to 4000 grit SiC. The specimens were subjected to a toothbrushing (TB) simulation device (Willytec) with rotating movements, toothpaste slurry and at three different loads (100g/250g/350g). At hourly intervals from 1h to 10h TB, mean surface roughness Ra was measured with an optical sensor and the surface gloss (Gl) with a glossmeter. Statistical analysis was performed for log-transformed Ra data applying two-way ANOVA to evaluate the interaction between load and material and load and brushing time. RESULTS: There was a significant interaction between material and load as well as between load and brushing time (p<0.0001). The microhybrid and hybrid materials demonstrated more surface deterioration with higher loads, whereas with the microfilled resins Heliomolar and Adoro it was vice versa. For ceramic materials, no or little deterioration was observed over time and independent of the load. The ceramic materials and 3 of the composite materials (roughness) showed no further deterioration after 5h of toothbrushing. Mean surface gloss was the parameter which discriminated best between the materials, followed by mean surface roughness Ra. There was a strong correlation between surface gloss and surface roughness for all the materials except the ceramics. The evaluation of the deterioration curves of individual specimens revealed a more or less synchronous course suspecting hinting specific external conditions and not showing the true variability in relation to the tested material. SIGNIFICANCE: The surface roughness and gloss of dental materials changes with brushing time and load and thus results in different material rankings. Apart from Grandio, the hybrid composite resins were more prone to surface changes than microfilled composites. The deterioration potential of a composite material can be quickly assessed by measuring surface gloss. For this purpose, a brushing time of 10h (=72,000 strokes) is needed. In further comparative studies, specimens of different materials should be tested in one series to estimate the true variability.

Comparison of the power spectral changes of the voluntary surface electromyogram and M wave during intermittent maximal voluntary contractions.

INTRODUCTION: To compare the power spectral changes of the voluntary surface electromyogram (sEMG) and of the compound action potential (M wave) in the vastus medialis and vastus lateralis muscles during fatiguing contractions. METHODS: Interference sEMG and force were recorded during 48 intermittent 3-s isometric maximal voluntary contractions (MVC) from 13 young, healthy subjects. M waves and twitches were evoked using supramaximal femoral nerve stimulation between the successive MVCs. Mean frequency (F mean), and median frequency were calculated from the sEMG and M waves. Muscle fiber conduction velocity (MFCV) was computed by cross-correlation. RESULTS: The power spectral shift to lower frequencies was significantly greater for the voluntary sEMG than for the M waves (P < 0.05). Over the fatiguing protocol, the overall average decrease in MFCV (~25 %) was comparable to that of sEMG F mean (~22 %), but significantly greater than that of M-wave F mean (~9 %) (P < 0.001). The mean decline in MFCV was highly correlated with the mean decreases in both sEMG and M-wave F mean. CONCLUSIONS: The present findings indicated that, as fatigue progressed, central mechanisms could enhance the relative weight of the low-frequency components of the voluntary sEMG power spectrum, and/or the end-of-fiber (non-propagating) components could reduce the sensitivity of the M-wave spectrum to changes in conduction velocity.

Investigation of high-resolution functional magnetic resonance imaging by means of surface and array radiofrequency coils at 7 T.

In this investigation, high-resolution, 1x1x1-mm(3) functional magnetic resonance imaging (fMRI) at 7 T is performed using a multichannel array head coil and a surface coil approach. Scan geometry was optimized for each coil separately to exploit the strengths of both coils. Acquisitions with the surface coil focused on partial brain coverage, while whole-brain coverage fMRI experiments were performed with the array head coil. BOLD sensitivity in the occipital lobe was found to be higher with the surface coil than with the head array, suggesting that restriction of signal detection to the area of interest may be beneficial for localized activation studies. Performing independent component analysis (ICA) decomposition of the fMRI data, we consistently detected BOLD signal changes and resting state networks. In the surface coil data, a small negative BOLD response could be detected in these resting state network areas. Also in the data acquired with the surface coil, two distinct components of the positive BOLD signal were consistently observed. These two components were tentatively assigned to tissue and venous signal changes.

Off-axis low-coherence interferometry for surface topology measurement.

In this communication we introduce a low or reduced coherence interferometry technique that can be used to retrieve surface topology on samples with high roughness. Moreover, we will show that the approach enables surface topology measurement also at the interface of so-called turbid media, where multiple scattering inside tissues can be a major issue, preventing accurate measurements.

Membrane insertion and antibody recognition of a glycosylphosphatidylinositol-anchored protein: an optical study.

The kinetics of binding of a glycolipid-anchored protein (the promastigote surface protease, PSP) to planar lecithin bilayers is studied by an integrated optics technique, in which the bilayer membrane is supported on an optical wave guide and the phase velocities of guided light modes in the wave guide are measured. From these velocities, the optical parameters of the membrane and PSP layers deposited on the waveguide are determined, yielding in particular the mass of PSP bound to the membrane, which is followed in real time. From a comparison of the binding rates of PSP and PSP from which the lipid moiety has been removed, it is shown that the lipid moiety plays a key role in anchoring the protein to the membrane. Specific and nonspecific binding of antibodies to membrane-anchored PSP is also investigated. As little as a fifth of a monolayer of PSP is sufficient to suppress the appreciable nonspecific binding of antibodies to the membrane.