An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

Effect of NaCl additive on solute-solvent interactions in aqueous polyethylene glycol-sodium sulfate two-phase systems

Abstract Partition behavior of eight small organic compounds and six proteins was examined in poly(ethylene glycol)-8000-sodium sulfate aqueous two-phase systems containing 0.215 M NaCl and 0.5 M osmolyte (sorbitol, sucrose, TMAO) and poly(ethylene glycol)-10000-sodium sulfate-0.215 M NaCl system, all in 0.01 M sodium phosphate buffer, pH 6.8. The differences between the solvent properties of the coexisting phases (solvent dipolarity/polarizability, hydrogen bond donor acidity, and hydrogen bond acceptor basicity) were characterized with solvatochromic dyes using the solvatochromic comparison method. Differences between the electrostatic properties of the phases were determined by analysis of partitioning of sodium salts of dinitrophenylated (DNP-) amino acids with aliphatic alkyl side-chain. The partition coefficients of all compounds examined (including proteins) were described in terms of solute-solvent interactions. The results obtained in the study show that solute-solvent interactions of nonionic organic compounds and proteins in polyethylene glycol-sodium sulfate aqueous two-phase system change in the presence of NaCl additive.

Effect of sodium chloride on solutesolvent interactions in aqueous polyethylene glycolsodium sulfate two-phase systems

Partition behavior of eight small organic compounds and six proteins was examined in poly(ethylene glycol)-8000sodium sulfate aqueous two-phase systems containing 0.215 M NaCl and 0.5 M osmolyte (sorbitol, sucrose, TMAO) and poly(ethylene glycol)-10000sodium sulfate0.215 M NaCl system, all in 0.01 M sodium phosphate buffer, pH 6.8. The differences between the solvent properties of the coexisting phases (solvent dipolarity/polarizability, hydrogen bond donor acidity, and hydrogen bond acceptor basicity) were characterized with solvatochromic dyes using the solvatochromic comparison method. Differences between the electrostatic properties of the phases were determined by analysis of partitioning of sodium salts of dinitrophenylated (DNP-) amino acids with aliphatic alkyl side-chain. The partition coefficients of all compounds examined (including proteins) were described in terms of solutesolvent interactions. The results obtained in the study show that solutesolvent interactions of nonionic organic compounds and proteins in polyethylene glycolsodium sulfate aqueous two-phase system change in the presence of NaCl additive.

Interrelationship between partition behavior of organic compounds and proteins in aqueous dextran-Polyethylene glycol and Polyethylene glycol-sodium sulfate two-phase systems

Partition behavior of adenosine and guanine mononucleotides was examined in aqueous dextran-polyethylene glycol (PEG) and PEG-sodium sulfate two-phase systems. The partition coefficients for each series of mononucleotides were analyzed as a functions of the number of phosphate groups and found to be dependent on the nature of nucleic base and on the type of \ATPS\ utilized. It was concluded that an average contribution of a phosphate group into logarithm of partition coefficient of a mononucleotide cannot be used to estimate the difference between the electrostatic properties of the coexisting phases of ATPS. The data obtained in this study were considered together with those for other organic compounds and proteins reported previously, and the linear interrelationship between logarithms of partition coefficients in dextran-PEG, PEG-Na2SO4 and PEG-Na2SO4-0.215 M NaCl (all in 0.01 M Na- or K/Na-phosphate buffer, pH 7.4 or 6.8) was established. Similar relationship was found for the previously reported data for proteins in Dex-PEG, PEG-600-Na2SO4, and PEG-8000-Na2SO4 ATPS. It is suggested that the linear relationships of the kind established in \ATPS\ may be observed for biological properties of compounds as well.

Microfluidic based on aqueous two-phase system for plasmidic DNA purification

Dissertação de mestrado em Bioquímica Aplicada (área de especialização em Biotecnologia)

Analytical and numerical investigation of two-phase flows

Multiphase flows, hyperbolic model, Godunov method, nozzle flow, nonstrictly hyperbolic

Optical measurements in disperse two-phase flows : application to rain formation in cumulus clouds

Magdeburg, Univ., Fak. für Verfahrens- und Systemtechnik, Diss., 2012

Simulation of the deformation behaviour of two-phase composites

Two-phase polycrystal, elasto-viscoplastic material model, heterogeneity, anisotropy, flow behaviour, crystallographic texture, strain field

Comparison of the methane oxidation rate in four media

Landfill gas emissions are one of the main sources of anthropogenic methane (CH4), a major greenhouse gas. In this paper, an economically attractive alternative to minimize greenhouse gas emissions from municipal solid waste landfills was sought. This alternative consists in special biofilters as landfill covers with oxidative capacity in the presence of CH4. To improve the quality/cost ratio of the project, compost was chosen as one of the cover substrates and soil (Typic red yellow-silt-clay Podzolic) as the other. The performance of four substrates was studied in laboratory experiments: municipal solid waste (MSW) compost, soil, and two soil-compost at different proportions. This study aimed to evaluate the suitability and environmental compatibility as a means of CH4 oxidation in biofilters. Four biofilters were constructed in 60 cm PVC tubes with an internal diameter of 10 cm. Each filter contained 2.3 L of oxidizing substrate at the beginning of the experiment. The gas used was a mixture of CH4 and air introduced at the bottom of each biofilter, at a flow of 150 mL min-1, by a flow meter. One hundred days after the beginning of the experiment, the best biofilter was the MSW compost with an oxidation rate of 990 g m-3 day-1 , corresponding to an efficiency of 44 %. It can be concluded that the four substrates studied have satisfactory oxidative capacity, and the substrates can be used advantageously as cover substrate of MSW landfills.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

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.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.