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.

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.

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.