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.

Dynamic, auxin-responsive plasma membrane-to-nucleus movement of Arabidopsis BRX.

In Arabidopsis, interplay between nuclear auxin perception and trans-cellular polar auxin transport determines the transcriptional auxin response. In brevis radix (brx) mutants, this response is impaired, probably indirectly because of disturbed crosstalk between the auxin and brassinosteroid pathways. Here we provide evidence that BRX protein is plasma membrane-associated, but translocates to the nucleus upon auxin treatment to modulate cellular growth, possibly in conjunction with NGATHA class B3 domain-type transcription factors. Application of the polar auxin transport inhibitor naphthalene phthalamic acid (NPA) resulted in increased BRX abundance at the plasma membrane. Thus, nuclear translocation of BRX could depend on cellular auxin concentration or on auxin flux. Supporting this idea, NPA treatment of wild-type roots phenocopied the brx root meristem phenotype. Moreover, BRX is constitutively turned over by the proteasome pathway in the nucleus. However, a stabilized C-terminal BRX fragment significantly rescued the brx root growth phenotype and triggered a hypocotyl gain-of-function phenotype, similar to strong overexpressors of full length BRX. Therefore, although BRX activity is required in the nucleus, excess activity interferes with normal development. Finally, similar to the PIN-FORMED 1 (PIN1) auxin efflux carrier, BRX is polarly localized in vascular cells and subject to endocytic recycling. Expression of BRX under control of the PIN1 promoter fully rescued the brx short root phenotype, suggesting that the two genes act in the same tissues. Collectively, our results suggest that BRX might provide a contextual readout to synchronize cellular growth with the auxin concentration gradient across the root tip.

NPAS2 as a transcriptional regulator of non-rapid eye movement sleep: genotype and sex interactions.

Because the transcription factor neuronal Per-Arnt-Sim-type signal-sensor protein-domain protein 2 (NPAS2) acts both as a sensor and an effector of intracellular energy balance, and because sleep is thought to correct an energy imbalance incurred during waking, we examined NPAS2's role in sleep homeostasis using npas2 knockout (npas2-/-) mice. We found that, under conditions of increased sleep need, i.e., at the end of the active period or after sleep deprivation (SD), NPAS2 allows for sleep to occur at times when mice are normally awake. Lack of npas2 affected electroencephalogram activity of thalamocortical origin; during non-rapid eye movement sleep (NREMS), activity in the spindle range (10-15 Hz) was reduced, and within the delta range (1-4 Hz), activity shifted toward faster frequencies. In addition, the increase in the cortical expression of the NPAS2 target gene period2 (per2) after SD was attenuated in npas2-/- mice. This implies that NPAS2 importantly contributes to the previously documented wake-dependent increase in cortical per2 expression. The data also revealed numerous sex differences in sleep; in females, sleep need accumulated at a slower rate, and REMS loss was not recovered after SD. In contrast, the rebound in NREMS time after SD was compromised only in npas2-/- males. We conclude that NPAS2 plays a role in sleep homeostasis, most likely at the level of the thalamus and cortex, where NPAS2 is abundantly expressed.

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.

Computer simulation study of liquid lithium at 470 and 843 K

Both structural and dynamical properties of 7Li at 470 and 843 K are studied by molecular dynamics simulation and the results are comapred with the available experimental data. Two effective interatomic potentials are used, i.e., a potential derived from the Ashcroft pseudopotential [Phys. Lett. 23, 48 (1966)] and a recently proposed potential deduced from the neutral pseudoatom method [J. Phys.: Condens. Matter 5, 4283 (1993)]. Although the shape of the two potential functions is very different, the majority of the properties calculated from them are very similar. The differences among the results using the two interaction models are carefully discussed.

Self-dynamic structure factor of dense liquids: Theory and simulation

The self-intermediate dynamic structure factor Fs(k,t) of liquid lithium near the melting temperature is calculated by molecular dynamics. The results are compared with the predictions of several theoretical approaches, paying special attention to the Lovesey model and the Wahnstrm and Sjgren mode-coupling theory. To this end the results for the Fs(k,t) second memory function predicted by both models are compared with the ones calculated from the simulations.

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Time correlation functions between the velocity of a tagged particle and velocities of particles within specified ranges of initial separations have been obtained by molecular dynamics simulation. These correlation functions have allowed us to analyze the momentum transfer between particles in different coordination shells. Two simple liquids at very different densities and two purely repulsive potentials with very different softnesses have been considered. The longitudinal correlations, which are the velocity cross-correlations along the initial direction defined by the centers of two given particles, have been calculated separately. It has been proven that these correlations should be attributed to particles both in front of and behind the central one. As with propagating longitudinal modes, they are strongly dependent on the softness of the potential core. Some characteristic features of the velocity correlation functions after the initial rise should be related to nonlongitudinal correlations. It has been shown that velocity cross-correlations between distinct particles cannot only be attributed to the direct interactions among particles, but also to the motions induced by the movement of a tagged particle on their neighbors.

Magnetic relaxation of a one-dimensional model for small particle systems with dipolar interactions: Monte Carlo simulation

Atribution as a function of the time are analyzed and this study leads to a deeper knowledge of the microscopic processes involved in the magnetic relaxation

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Monte Carlo simulation of a whole-body counter using IGOR phantoms.

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.