Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

Simultaneous steady state and dynamic design of process systems using structural indicators

The simultaneous design of the steady-state and dynamic performance of a process has the ability to satisfy much more demanding dynamic performance criteria than the design of dynamics only by the connection of a control system. A method for designing process dynamics based on the use of a linearised systems' eigenvalues has been developed. The eigenvalues are associated with system states using the unit perturbation spectral resolution (UPSR), characterising the dynamics of each state. The design method uses a homotopy approach to determine a final design which satisfies both steady-state and dynamic performance criteria. A highly interacting single stage forced circulation evaporator system, including control loops, was designed by this method with the goal of reducing the time taken for the liquid composition to reach steady-state. Initially the system was successfully redesigned to speed up the eigenvalue associated with the liquid composition state, but this did not result in an improved startup performance. Further analysis showed that the integral action of the composition controller was the source of the limiting eigenvalue. Design changes made to speed up this eigenvalue did result in an improved startup performance. The proposed approach provides a structured way to address the design-control interface, giving significant insight into the dynamic behaviour of the system such that a systematic design or redesign of an existing system can be undertaken with confidence.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Entanglement in the steady state of a collective-angular-momentum (Dicke) model

We model the behavior of an ion trap with all ions driven simultaneously and coupled collectively to a heat bath. The equations for this system are similar to the irreversible dynamics of a collective angular momentum system known as the Dicke model. We show how the steady state of the ion trap as a dissipative many-body system driven far from equilibrium can exhibit quantum entanglement. We calculate the entanglement of this steady state for two ions in the trap and in the case of more than two ions we calculate the entanglement between two ions by tracing over all the other ions. The entanglement in the steady state is a maximum for the parameter values corresponding roughly to a bifurcation of a fixed point in the corresponding semiclassical dynamics. We conjecture that this is a general mechanism for entanglement creation in driven dissipative quantum systems.

Experimental and model-based analysis of the steady-state and dynamic operating behaviour of the direct methanol fuel cell (DMFC)

Direct methanol fuel cell, DMFC, model, mass transport, Maxwell-Stefan, Flory-Huggins, crossover, polymer electrolyte membrane, Nafion

Second-Order Approximation to the Rotemberg Model around a Distorted Steady State

Less is known about social welfare objectives when it is costly to change prices, as in Rotemberg (1982), compared with Calvo-type models. We derive a quadratic approximate welfare function around a distorted steady state for the costly price adjustment model. We highlight the similarities and differences to the Calvo setup. Both models imply inflation and output stabilization goals. It is explained why the degree of distortion in the economy influences inflation aversion in the Rotemberg framework in a way that differs from the Calvo setup.

Plerixafor and granulocyte colony-stimulating factor for first-line steady-state autologous peripheral blood stem cell mobilization in lymphoma and multiple myeloma: results of the prospective PREDICT trial.

In Europe, the combination of plerixafor + granulocyte colony-stimulating factor is approved for the mobilization of hematopoietic stem cells for autologous transplantation in patients with lymphoma and myeloma whose cells mobilize poorly. The purpose of this study was to further assess the safety and efficacy of plerixafor + granulocyte colony-stimulating factor for front-line mobilization in European patients with lymphoma or myeloma. In this multicenter, open label, single-arm study, patients received granulocyte colony-stimulating factor (10 μg/kg/day) subcutaneously for 4 days; on the evening of day 4 they were given plerixafor (0.24 mg/kg) subcutaneously. Patients underwent apheresis on day 5 after a morning dose of granulocyte colony-stimulating factor. The primary study objective was to confirm the safety of mobilization with plerixafor. Secondary objectives included assessment of efficacy (apheresis yield, time to engraftment). The combination of plerixafor + granulocyte colony-stimulating factor was used to mobilize hematopoietic stem cells in 118 patients (90 with myeloma, 25 with non-Hodgkin's lymphoma, 3 with Hodgkin's disease). Treatment-emergent plerixafor-related adverse events were reported in 24 patients. Most adverse events occurred within 1 hour after injection, were grade 1 or 2 in severity and included gastrointestinal disorders or injection-site reactions. The minimum cell yield (≥ 2 × 10(6) CD34(+) cells/kg) was harvested in 98% of patients with myeloma and in 80% of those with non-Hodgkin's lymphoma in a median of one apheresis. The optimum cell dose (≥ 5 × 10(6) CD34(+) cells/kg for non-Hodgkin's lymphoma or ≥ 6 × 10(6) CD34(+) cells/kg for myeloma) was harvested in 89% of myeloma patients and 48% of non-Hodgkin's lymphoma patients. In this prospective, multicenter European study, mobilization with plerixafor + granulocyte colony-stimulating factor allowed the majority of patients with myeloma or non-Hodgkin's lymphoma to undergo transplantation with minimal toxicity, providing further data supporting the safety and efficacy of plerixafor + granulocyte colony-stimulating factor for front-line mobilization of hematopoietic stem cells in patients with non-Hodgkin's lymphoma or myeloma.

Steady state plasma levels of the enantiomers of trimipramine and of its metabolites in CYP2D6-, CYP2C19- and CYP3A4/5-phenotyped patients.

Steady state plasma concentrations of the (L)- and (D)-enantiomers of trimipramine (TRI), desmethyltrimipramine (DTRI), 2-hydroxytrimipramine (TRIOH) and 2-hydroxydesmethyl-trimipramine (DTRIOH) were measured in 27 patients receiving between 300 and 400 mg/day racemic TRI. The patients were phenotyped with dextromethorphan and mephenytoin, and the 8-hour urinary ratios of dextromethorphan/dextrorphan, dextromethorphan/3-methoxymorphinan, and (S)-mephenytoin/(R)mephenytoin were used as markers of cytochrome P-450IID6 (CYP2D6), CYP3A4/5 and CYP2C19 activities, respectively. One patient was a CYP2D6 and one was a CYP2C19 poor metabolizer. A stereoselectivity in the metabolism of TRI has been found, with a preferential N-demethylation of (D)-TRI and a preferential hydroxylation of (L)-TRI. CYP2D6 appears to be involved in the 2-hydroxylation of (L)-TRI, (L)DTRI and (D)-DTRI, but not of (D)-TRI, as significant correlations were measured between the dextromethorphan/dextrorphan ratios and the (L)-TRI/(L)-TRIOH (r = 0.45, p = 0.019), the (L)-DTRI/(L)-DTRIOH (r = 0.47, p = 0.014), and the (D)-DTRI/(D)-DTRIOH (r = 0.51, p = 0.006), but not with the (D)-TRI/(D)-TRIOH ratios (r = 0.29, NS). CYP2C19, but not CYP2D6, appears to be involved in the demethylation pathway, with a stereoselectivity toward the (D)-enantiomer of TRI, as a significant positive correlation was calculated between the mephenytoin (S)/(R) ratios and the concentrations to dose-to-weight ratios of (D)-TRI (r = 0.69, p = 0.00006). CYP3A4/5 appears to be involved in the metabolism of (L)-TRI to a presently not determined metabolite. The CYP2D6 poor metabolizer had the highest (L)-DTRI and (D)-DTRI concentrations to dose-to-weight ratios, and the CYP2C19 poor metabolizer had the highest (L)-TRI and (D)-TRI concentrations to dose-to-weight ratios of the group.

Deep thiopental anesthesia alters steady-state glucose homeostasis but not the neurochemical profile of rat cortex.

Barbiturates are regularly used as an anesthetic for animal experimentation and clinical procedures and are frequently provided with solubilizing compounds, such as ethanol and propylene glycol, which have been reported to affect brain function and, in the case of (1)H NMR experiments, originate undesired resonances in spectra affecting the quantification. As an alternative, thiopental can be administrated without any solubilizing agents. The aim of the study was to investigate the effect of deep thiopental anesthesia on the neurochemical profile consisting of 19 metabolites and on glucose transport kinetics in vivo in rat cortex compared with alpha-chloralose using localized (1)H NMR spectroscopy. Thiopental was devoid of effects on the neurochemical profile, except for the elevated glucose at a given plasma glucose level resulting from thiopental-induced depression of glucose consumption at isoelectrical condition. Over the entire range of plasma glucose levels, steady-state glucose concentrations were increased on average by 48% +/- 8%, implying that an effect of deep thiopental anesthesia on the transport rate relative to cerebral glucose consumption ratio was increased by 47% +/- 8% compared with light alpha-chloralose-anesthetized rats. We conclude that the thiopental-induced isoelectrical condition in rat cortex significantly affected glucose contents by depressing brain metabolism, which remained substantial at isoelectricity.

Deep thiopental anesthesia alters steady-state glucose homeostasis but not the neurochemical profile of rat cortex.

Barbiturates are regularly used as an anesthetic for animal experimentation and clinical procedures and are frequently provided with solubilizing compounds, such as ethanol and propylene glycol, which have been reported to affect brain function and, in the case of (1)H NMR experiments, originate undesired resonances in spectra affecting the quantification. As an alternative, thiopental can be administrated without any solubilizing agents. The aim of the study was to investigate the effect of deep thiopental anesthesia on the neurochemical profile consisting of 19 metabolites and on glucose transport kinetics in vivo in rat cortex compared with alpha-chloralose using localized (1)H NMR spectroscopy. Thiopental was devoid of effects on the neurochemical profile, except for the elevated glucose at a given plasma glucose level resulting from thiopental-induced depression of glucose consumption at isoelectrical condition. Over the entire range of plasma glucose levels, steady-state glucose concentrations were increased on average by 48% +/- 8%, implying that an effect of deep thiopental anesthesia on the transport rate relative to cerebral glucose consumption ratio was increased by 47% +/- 8% compared with light alpha-chloralose-anesthetized rats. We conclude that the thiopental-induced isoelectrical condition in rat cortex significantly affected glucose contents by depressing brain metabolism, which remained substantial at isoelectricity.

Test of the fluctuation theorem for stochastic entropy production in a nonequilibrium steady state

We derive a simple closed analytical expression for the total entropy production along a single stochastic trajectory of a Brownian particle diffusing on a periodic potential under an external constant force. By numerical simulations we compute the probability distribution functions of the entropy and satisfactorily test many of the predictions based on Seiferts integral fluctuation theorem. The results presented for this simple model clearly illustrate the practical features and implications derived from such a result of nonequilibrium statistical mechanics.

Steady State Concentrations of the Enantiomers of Mianserin and Desmethylmianserin in Poor and in Homozygous and Heterozygous Extensive Metabolizers of Debrisoquine.

AB Summary: Steady state concentrations of (S)- and (R)-mianserin and desmethylmianserin were measured in 21 homozygous extensive metabolizers (as determined by genotyping for mutations 3 [or A] and 4 [or B]), in seven heterozygous extensive metabolizers and in one poor metabolizer of debrisoquine, as well as in one patient receiving very high doses of mianserin (360 mg/day) and fluoxetine (160 mg/day), a strong cytochrome P450IID6 inhibitor. The mean dose of mianserin was (mean +/- SD, range: 67 +/- 63, 10 to 360 mg/day). High dispersions of the (S)/(R)-mianserin and desmethylmianserin ratios were observed (mean +/- SD, range: 2.10+/- 1.01, 0.64 to 4.76, and 0.29 +/- 0.14, 0.08 to 0.57, respectively). The highest (S)/(R)-mianserin ratio was calculated for the poor metabolizer (4.76) agreeing with those results of a single-dose study with poor and extensive metabolizers of debrisoquine, in that the cytochrome P450IID6 is probably involved in the metabolism of mianserin with an enantioselectivity for the (S)-enantiomer. Nevertheless, the mean concentration-to-dose ratios for (S)- or (R)-mianserin or desmethylmianserin were not significantly different between homozygous and heterozygous extensive metabolizers, and no particular values were measured in the poor metabolizer nor in the patient receiving fluoxetine. Furthermore, the(S)/(R)-mianserin ratio measured in the PM was only slightly higher than the second highest ratio (3.85) of an homozygous extensive metabolizer, whereas no particular value (2.92) was calculated for the patient taking fluoxetine. Finally, no significant differences in (S)/(R)-mianserin or(S)/(R)-desmethylmianserin were calculated between homozygous and heterozygous extensive metabolizers. Although the number of patients included in this study is too low to allow definite conclusions, the results suggest that the debrisoquine genotype has only a moderate influence on the steady state concentrations of the enantiomers of mianserin and desmethylmianserin. (C) Lippincott-Raven Publishers

Sulfuric acid - glucose separation using mixed-recycle steady state recycling chromatography

Batch chromatography is a widely used separation technique in a variety of fields meeting difficult separations. Several technologies for improving the performance of chromatography have been studied, including mixed-recycle steady state recycling (MR-SSR) chromatography. Design of MR-SSR has been commonly limited on 100 % purity constraint cases and empirical work. In this study a predictive design method was used to optimize feed pulse size and design a number of experimental MR-SSR separations for a solution of 20 % sulfuric acid and 100 g/L glucose. The design was under target product fraction purities of 98.7 % for H2SO4 and 95 % for glucose. The experiments indicate a maximum of 59 % increase in sulfuric acid productivity and 82 % increase for glucose when compared to corresponding batch separation. Eluent consumption was lowered by approximately 50 % using recycling chromatography. Within this study the target purities and yields set in design were not completely met, and further optimization of the process is deemed necessary.

Simulation of Steady-State Nonlinear Heat Transfer Problems Through the Minimization of Quadratic Functionals

In this work it is presented a systematic procedure for constructing the solution of a large class of nonlinear conduction heat transfer problems through the minimization of quadratic functionals like the ones usually employed for linear descriptions. The proposed procedure gives rise to an efficient and easy way for carrying out numerical simulations of nonlinear heat transfer problems by means of finite elements. To illustrate the procedure a particular problem is simulated by means of a finite element approximation.