Pressure-dependent scaling scenarios in experiments of spontaneous imbibition

The scaling properties of the rough liquid-air interface formed in the spontaneous imbibition of a viscous liquid by a model porous medium are found to be very sensitive to the magnitude of the pressure difference applied at the liquid inlet. Interface fluctuations change from obeying intrinsic anomalous scaling at large negative pressure differences, to being super-rough with the same dynamic exponent z¿3 at less negative pressure differences, to finally obeying ordinary Family-Vicsek scaling with z¿2 at large positive pressure differences. This rich scenario reflects the relative importance on different length scales of capillary and permeability disorder, and the role of surface tension and viscous pressure in damping interface fluctuations.

Experiments of periodic forcing of Saffman-Taylor fingers

We report on an experimental study of long normal Saffman-Taylor fingers subject to periodic forcing. The sides of the finger develop a low amplitude, long wavelength instability. We discuss the finger response in stationary and nonstationary situations, as well as the dynamics towards the stationary states. The response frequency of the instability increases with forcing frequency at low forcing frequencies, while, remarkably, it becomes independent of forcing frequency at large forcing frequencies. This implies a process of wavelength selection. These observations are in good agreement with previous numerical results reported in [Ledesma-Aguilar et al., Phys. Rev. E 71, 016312 (2005)]. We also study the average value of the finger width, and its fluctuations, as a function of forcing frequency. The average finger width is always smaller than the width of the steady-state finger. Fluctuations have a nonmonotonic behavior with a maximum at a particular frequency.

Experiments in a rotating Hele-Shaw cell

We introduce a modification to Hele-Shaw flows consisting of a rotating cell. A viscous fluid (oil) is injected at the rotation axis of the cell, which is open to air. The morphological instability of the oil-air interface is thus driven by centrifugal force and is controlled by the density (not viscosity) difference. We derive the linear dispersion relation and verify the maximum growth rate selection of initial patterns within experimental uncertainty. The nonlinear growth regime is studied in the case of vanishing injection rate. Several characteristic lengths are studied to quantify the patterns obtained. Experimental data exhibit good collapse for two characteristic lengths, namely, the radius of gyration and the radial finger length, which in the nonlinear regime appear to grow linearly in time.

A simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Monte Carlo simulation of kilovolt electron transport in solids

A Monte Carlo procedure to simulate the penetration and energy loss of low¿energy electron beams through solids is presented. Elastic collisions are described by using the method of partial waves for the screened Coulomb field of the nucleus. The atomic charge density is approximated by an analytical expression with parameters determined from the Dirac¿Hartree¿Fock¿Slater self¿consistent density obtained under Wigner¿Seitz boundary conditions in order to account for solid¿state effects; exchange effects are also accounted for by an energy¿dependent local correction. Elastic differential cross sections are then easily computed by combining the WKB and Born approximations to evaluate the phase shifts. Inelastic collisions are treated on the basis of a generalized oscillator strength model which gives inelastic mean free paths and stopping powers in good agreement with experimental data. This scattering model is accurate in the energy range from a few hundred eV up to about 50 keV. The reliability of the simulation method is analyzed by comparing simulation results and experimental data from backscattering and transmission measurements.

Simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Prediction of Infant Drug Exposure Through Breastfeeding: Population PK Modeling and Simulation of Fluoxetine Exposure.

The likelihood of significant exposure to drugs in infants through breast milk is poorly defined, given the difficulties of conducting pharmacokinetics (PK) studies. Using fluoxetine (FX) as an example, we conducted a proof-of-principle study applying population PK (popPK) modeling and simulation to estimate drug exposure in infants through breast milk. We simulated data for 1,000 mother-infant pairs, assuming conservatively that the FX clearance in an infant is 20% of the allometrically adjusted value in adults. The model-generated estimate of the milk-to-plasma ratio for FX (mean: 0.59) was consistent with those reported in other studies. The median infant-to-mother ratio of FX steady-state plasma concentrations predicted by the simulation was 8.5%. Although the disposition of the active metabolite, norfluoxetine, could not be modeled, popPK-informed simulation may be valid for other drugs, particularly those without active metabolites, thereby providing a practical alternative to conventional PK studies for exposure risk assessment in this population.

Arrangement of experiments for simulating the effects of elevated temperatures and elevated CO2 levels on field-sown crops in Finland

Selostus: Koejärjestelyt kohonneen lämpötilan ja CO2-tason vaikutusten simuloimiseksi peltokasveilla Suomessa

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Magnetic microstructures from magnetic force microscopy and Monte Carlo simulation in CoFe-Ag-Cu granular films

CoFe-Ag-Cu granular films, prepared by rf sputtering, displayed magnetic domain microstructures for ferromagnetic concentrations above about 32% at, and below the percolation threshold. All samples have a fcc structure with an (111) texture perpendicular to the film plane. Magnetic force microscopy (MFM) showed a variety of magnetic domain microstructures, extremely sensitive to the magnetic history of the sample, which arise from the balance of the ferromagnetic exchange, the dipolar interactions and perpendicular magnetocrystalline anisotropy, MFM images indicate that in virgin samples, magnetic bubble domains with an out-of-plane component of the magnetization are surrounded by a quasicontinuous background of opposite magnetization domains. The application of a magnetic field in different geometries drastically modifies the microstructure of the system in the remanent state: i) for an in-plane field, the MFM images show that most of the magnetic moments are aligned along the film plane, ii) for an out-of-plane field, the MFM signal increases about one order of magnitude, and out-of-plane striped domains with alternating up and down magnetization are stabilized. Numerical simulations show that a variety of metastable domain structures (similar to those observed experimentally) can be reached, depending on magnetic history, in systems with competing perpendicular anisotropy, exchange and dipolar interactions.