Incompressible limit and well-posedness of PDE models of tissue growth

Both compressible and incompressible porous medium models are used in the literature to describe the mechanical aspects of living tissues. Using a stiff pressure law, it is possible to build a link between these two different representations. In the incompressible limit, compressible models generate free boundary problems where saturation holds in the moving domain. Our work aims at investigating the stiff pressure limit of reaction-advection-porous medium equations motivated by tumor development. Our first study concerns the analysis and numerical simulation of a model including the effect of nutrients. A coupled system of equations describes the cell density and the nutrient concentration and the derivation of the pressure equation in the stiff limit was an open problem for which the strong compactness of the pressure gradient is needed. To establish it, we use two new ideas: an L3-version of the celebrated Aronson-Bénilan estimate, and a sharp uniform L4-bound on the pressure gradient. We further investigate the sharpness of this bound through a finite difference upwind scheme, which we prove to be stable and asymptotic preserving. Our second study is centered around porous medium equations including convective effects. We are able to extend the techniques developed for the nutrient case, hence finding the complementarity relation on the limit pressure. Moreover, we provide an estimate of the convergence rate at the incompressible limit. Finally, we study a multi-species system. In particular, we account for phenotypic heterogeneity, including a structured variable into the problem. In this case, a cross-(degenerate)-diffusion system describes the evolution of the phenotypic distributions. Adapting methods recently developed in the context of two-species systems, we prove existence of weak solutions and we pass to the incompressible limit. Furthermore, we prove new regularity results on the total pressure, which is related to the total density by a power law of state.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Development of filter media treatments for liquid filtration

Woven monofilament, multifilament, and spun yarn filter media have long been the standard media in liquid filtration equipment. While the energy for a solid-liquid separation process is determined by the engineering work, it is the interface between the slurry and the equipment - the filter media - that greatly affects the performance characteristics of the unit operation. Those skilled in the art are well aware that a poorly designed filter medium may endanger the whole operation, whereas well-performing filter media can make the operation smooth and economical. As the mineral and pulp producers seek to produce ever finer and more refined fractions of their products, it is becoming increasingly important to be able to dewater slurries with average particle sizes around 1 ¿m using conventional, high-capacity filtration equipment. Furthermore, the surface properties of the media must not allow sticky and adhesive particles to adhere to the media. The aim of this thesis was to test how the dirt-repellency, electrical resistance and highpressure filtration performance of selected woven filter media can be improved by modifying the fabric or yarn with coating, chemical treatment and calendering. The results achieved by chemical surface treatments clearly show that the woven media surface properties can be modified to achieve lower electrical resistance and improved dirt-repellency. The main challenge with the chemical treatments is the abrasion resistance and, while the experimental results indicate that the treatment is sufficiently permanent to resist standard weathering conditions, they may still prove to be inadequately strong in terms of actual use.From the pressure filtration studies in this work, it seems obvious that the conventional woven multifilament fabrics still perform surprisingly well against the coated media in terms of filtrate clarity and cake build-up. Especially in cases where the feed slurry concentration was low and the pressures moderate, the conventional media seemed to outperform the coated media. In the cases where thefeed slurry concentration was high, the tightly woven media performed well against the monofilament reference fabrics, but seemed to do worse than some of the coated media. This result is somewhat surprising in that the high initial specific resistance of the coated media would suggest that the media will blind more easily than the plain woven media. The results indicate, however, that it is actually the woven media that gradually clogs during the coarse of filtration. In conclusion, it seems obvious that there is a pressure limit above which the woven media looses its capacity to keep the solid particles from penetrating the structure. This finding suggests that for extreme pressures the only foreseeable solution is the coated fabrics supported by a strong enough woven fabric to hold thestructure together. Having said that, the high pressure filtration process seems to follow somewhat different laws than the more conventional processes. Based on the results, it may well be that the role of the cloth is most of all to support the cake, and the main performance-determining factor is a long life time. Measuring the pore size distribution with a commercially available porometer gives a fairly accurate picture of the pore size distribution of a fabric, but failsto give insight into which of the pore sizes is the most important in determining the flow through the fabric. Historically air, and sometimes water, permeability measures have been the standard in evaluating media filtration performance including particle retention. Permeability, however, is a function of a multitudeof variables and does not directly allow the estimation of the effective pore size. In this study a new method for estimating the effective pore size and open pore area in a densely woven multifilament fabric was developed. The method combines a simplified equation of the electrical resistance of fabric with the Hagen-Poiseuille flow equation to estimate the effective pore size of a fabric and the total open area of pores. The results are validated by comparison to the measured values of the largest pore size (Bubble point) and the average pore size. The results show good correlation with measured values. However, the measured and estimated values tend to diverge in high weft density fabrics. This phenomenon is thought to be a result of a more tortuous flow path of denser fabrics, and could most probably be cured by using another value for the tortuosity factor.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Changes in respiratory mechanics and gas exchange during the acute respiratory distress syndrome

BACKGROUND: The time course of impairment of respiratory mechanics and gas exchange in the acute respiratory distress syndrome (ARDS) remains poorly defined. We assessed the changes in respiratory mechanics and gas exchange during ARDS. We hypothesized that due to the changes in respiratory mechanics over time, ventilatory strategies based on rigid volume or pressure limits might fail to prevent overdistension throughout the disease process. METHODS: Seventeen severe ARDS patients {PaO2/FiO2 10.1 (9.2-14.3) kPa; 76 (69-107) mmHg [median (25th-75th percentiles)] and bilateral infiltrates} were studied during the acute, intermediate, and late stages of ARDS (at 1-3, 4-6 and 7 days after diagnosis). Severity of lung injury, gas exchange, and hemodynamics were assessed. Pressure-volume (PV) curves of the respiratory system were obtained, and upper and lower inflection points (UIP, LIP) and recruitment were estimated. RESULTS: (1) UIP decreased from early to established (intermediate and late) ARDS [30 (28-30) cmH2O, 27 (25-30) cmH2O and 25 (23-28) cmH2O (P=0.014)]; (2) oxygenation improved in survivors and in patients with non-pulmonary etiology in late ARDS, whereas all patients developed hypercapnia from early to established ARDS; and (3) dead-space ventilation and pulmonary shunt were larger in patients with pulmonary etiology during late ARDS. CONCLUSION: We found a decrease in UIP from acute to established ARDS. If applied to our data, the inspiratory pressure limit advocated by the ARDSnet (30 cmH2O) would produce ventilation over the UIP, with a consequent increased risk of overdistension in 12%, 43% and 65% of our patients during the acute, intermediate and late phases of ARDS, respectively. Lung protective strategies based on fixed tidal volume or pressure limits may thus not fully avoid the risk of lung overdistension throughout ARDS.

Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

Ab initio quantum chemical studies on neutral-radical reactions of ethynyl (C2H) and cyano (CN) with unsaturated hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA's Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth's primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan's atmosphere, the ISM, and cold celestial bodies.<.

Ab Initio Quantum Chemical Studies on Neutral-Radical Reactions of Ethynyl (C2H) and Cyano (CN) with Unsaturated Hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA’s Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth’s primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan’s atmosphere, the ISM, and cold celestial bodies.

Design of pressure vessels using shape optimization: An integrated approach

Previous papers related to the optimization of pressure vessels have considered the optimization of the nozzle independently from the dished end. This approach generates problems such as thickness variation from nozzle to dished end (coupling cylindrical region) and, as a consequence, it reduces the optimality of the final result which may also be influenced by the boundary conditions. Thus, this work discusses shape optimization of axisymmetric pressure vessels considering an integrated approach in which the entire pressure vessel model is used in conjunction with a multi-objective function that aims to minimize the von-Mises mechanical stress from nozzle to head. Representative examples are examined and solutions obtained for the entire vessel considering temperature and pressure loading. It is noteworthy that different shapes from the usual ones are obtained. Even though such different shapes may not be profitable considering present manufacturing processes, they may be competitive for future manufacturing technologies, and contribute to a better understanding of the actual influence of shape in the behavior of pressure vessels. (C) 2011 Elsevier Ltd. All rights reserved.

Earth pressure on cantilever walls at design retained heights

There are many methods for the analysis and design of embedded cantilever retaining walls. They involve various different simplifications of the pressure distribution to allow calculation of the limiting equilibrium retained height and the bending moment when the retained height is less than the limiting equilibrium value, i.e. the serviceability case. Recently, a new method for determining the serviceability earth pressure and bending moment has been proposed. This method makes an assumption defining the point of zero net pressure. This assumption implies that the passive pressure is not fully mobilised immediately below the excavation level. The finite element analyses presented in this paper examine the net pressure distribution on walls in which the retained height is less, than the limiting equilibrium value. The study shows that for all practical walls, the earth pressure distributions on the front and back of the wall are at their limit values, Kp and K-a respectively, when the lumped factor of safety F-r is less than or equal to2.0. A rectilinear net pressure distribution is proposed that is intuitively logical. It produces good predictions of the complete bending moment diagram for walls in the service configuration and the proposed method gives results that have excellent agreement with centrifuge model tests. The study shows that the method for determining the serviceability bending moment suggested by Padfield and Mair(1) in the CIRIA Report 104 gives excellent predictions of the maximum bending moment in practical cantilever walls. It provides the missing data that have been needed to verify and justify the CIRIA 104 method.

Does shoot water status limit leaf expansion of nitrogen-deprived barley?

The role of shoot water status in mediating the decline in leaf elongation rate of nitrogen (N)-deprived barley plants was assessed. Plants were grown at two levels of N supply, with or without the application of pneumatic pressure to the roots. Applying enough pressure (balancing pressure) to keep xylem sap continuously bleeding from the cut surface of a leaf allowed the plants to remain at full turgor throughout the experiments. Plants from which N was withheld required a greater balancing pressure during both day and night. This difference in balancing pressure was greater at high (2.0 kPa) than low (1.2 kPa) atmospheric vapour pressure deficit (VPD). Pressurizing the roots did not prevent the decline in leaf elongation rate induced by withholding N at either high or low VPD. Thus low shoot water status did not limit leaf growth of N-deprived plants.

Experimental assessment of pressure drop in air stripping columns

One important step in the design of air stripping operations for the removal of VOC is the choice of operating conditions, which are based in the phase ratio. This parameter sets on directly the stripping factor and the efficiency of the operation. Its value has an upper limit determined by the flooding regime, which is previewed using empirical correlations, namely the one developed by Eckert. This type of approach is not suitable for the development of algorithms. Using a pilot scale column and a convenient solution, the pressure drop was determined in different operating conditions and the experimental values were compared with the estimations. This particular research will be incorporated in a global model for simulating the dynamics of air stripping using a multi variable distributed parameter system.

Quantification of clenbuterol at trace level in human urine by ultra-high pressure liquid chromatography-tandem mass spectrometry.

Clenbuterol is a β2 agonist agent with anabolic properties given by the increase in the muscular mass in parallel to the decrease of the body fat. For this reason, the use of clenbuterol is forbidden by the World Anti-Doping Agency (WADA) in the practice of sport. This compound is of particular interest for anti-doping authorities and WADA-accredited laboratories due to the recent reporting of risk of unintentional doping following the eating of meat contaminated with traces of clenbuterol in some countries. In this work, the development and the validation of an ultra-high pressure liquid chromatography coupled to electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) method for the quantification of clenbuterol in human urine is described. The analyte was extracted from urine samples by liquid-liquid extraction (LLE) in basic conditions using tert butyl-methyl ether (TBME) and analyzed by UHPLC-MS/MS with a linear gradient of acetonitrile in 9min only. The simple and rapid method presented here was validated in compliance with authority guidelines and showed a limit of quantification at 5pg/mL and a linearity range from 5pg/mL to 300pg/mL. Good trueness (85.8-105%), repeatability (5.7-10.6% RSD) and intermediate precision (5.9-14.9% RSD) results were obtained. The method was then applied to real samples from eighteen volunteers collecting urines after single oral doses administration (1, 5 and 10μg) of clenbuterol-enriched yogurts.