Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Non-existence of polar factorisations and polar inclusion of a vector-valued mapping

R.J. DOUGLAS, Non-existence of polar factorisations and polar inclusion of a vector-valued mapping. Intern. Jour. Of Pure and Appl. Math., (IJPAM) 41, no. 3 (2007).

Specific and unspecific interactions in polar fluids in view of wideband dielectric far-infrared spectra

A simple molecular analytical theory of dielectric relaxation in strongly polar fluids is considered in terms of a semi- phenomenological approach. Theoretical spectra epsilon(v), a(v) of complex permittivity and absorption coefficient are fully determined by a form of intermolecular potential well, in which a dipole reorients. In a recent publication by VI. Gaiduk, O.F. Nielsen, and T.S. Perova [J. Molliq 95 (1002) 1-25] the wideband spectra of liquid H2O and D2O were described in terms of a composite model comprising the rectangular and the cosine squared potential wells. Much better results are achieved in this work, where the rectangular well is replaced by a well with a rounded bottom termed the hat-curved well. The spectrum of the auto-correlation function (ACF) is calculated for such a potential. The proposed theory of a composite model, comprising hat-curved and parabolic wells, is applied for liquid water. This model is capable for describing the Debye relaxation region, the second relaxation region in the submillimeter wavelength range, and the far infra-red (FIR) e(v), a(v) spectra, where an intense librational band and an additional weak band are placed, respectively, near 700 cm(-1) and 200 cm(-1). The latter band reflects the features of so-called specific (viz. directly related to H-bonds) interactions and the former band reflects the features of unspecific interactions. The physical mechanisms connected with these types of interactions are discussed in terms of two relevant types of water structure (types of molecular rotation). The proposed theory is also applied to a non-associated liquid in terms of one hat-curved potential well. (C) 2004 Elsevier B.V. All rights reserved.

Non-genetic mechanisms communicating antibiotic resistance:rethinking strategies for antimicrobial drug design

Introduction: Infections by multidrug-resistant bacteria are of great concern worldwide. In many cases, resistance is not due to the presence of specific antibiotic-modifying enzymes, but rather associated with a general impermeability of the bacterial cell envelope. The molecular bases of this intrinsic resistance are not completely understood. Moreover, horizontal gene transfers cannot solely explain the spread of intrinsic resistance among bacterial strains. Areas covered: This review focuses on the increased intrinsic antibiotic resistance mediated by small molecules. These small molecules can either be secreted from bacterial cells of the same or different species (e.g., indole, polyamines, ammonia, and the Pseudomonas quinolone signal) or be present in the bacterial cell milieu, whether in the environment, such as indole acetic acid and other plant hormones, or in human tissues and body fluids, such as polyamines. These molecules are metabolic byproducts that act as infochemicals and modulate bacterial responses toward antibiotics leading to increasing or decreasing resistance levels. Expert opinion: The non-genetic mechanisms of antibiotic response modulation and communication discussed in this review should reorient our thinking of the mechanisms of intrinsic resistance to antibiotics and its spread across bacterial cell populations. The identification of chemical signals mediating increased intrinsic antibiotic resistance will expose novel critical targets for the development of new antimicrobial strategies.

Tetrabutylammonium methacrylate as a novel receptor for selective extraction of sulphonylurea drugs from biological fluids using molecular imprinting

Glibenclamide (GLIB), an oral antidiabetic medication of the sulphonylurea drugs family, was stoichiometrically imprinted using tetrabutylammonium methacrylate as the functional monomer, for the first time in molecular imprinting, and utilising the sulphonylurea affinity for carboxylate anions. Solution association between the drug and the novel functional monomer was studied by 1H-NMR titrations, whereby evidence of sulphonylurea deprotonation followed by the formation of “narcissistic” GLIB dimers was found when tested in CDCl3, while an affinity constant in excess of 105 L mol-1 was measured in DMSO-d6. Detailed analysis of GLIB binding on the subsequently prepared imprinted and non-imprinted polymers confirmed deactivation of binding sites by exchange of a proton between GLIB and methacrylate, followed by extraction of the tetrabutylammonium counterion from the polymer matrix, resulting in overall reduced binding capacities and affinities by the imprinted material under equilibrium conditions. An optimised MI-SPE protocol, which included a binding site re-activation step, was developed for the extraction of GLIB from blood serum, whereby recoveries of up to 92.4% were obtained with exceptional sample clean-up.

Diffusion coefficients in dense and supercritical fluids

Os coeficientes de difusão (D 12) são propriedades fundamentais na investigação e na indústria, mas a falta de dados experimentais e a inexistência de equações que os estimem com precisão e confiança em fases comprimidas ou condensadas constituem limitações importantes. Os objetivos principais deste trabalho compreendem: i) a compilação de uma grande base de dados para valores de D 12 de sistemas gasosos, líquidos e supercríticos; ii) o desenvolvimento e validação de novos modelos de coeficientes de difusão a diluição infinita, aplicáveis em amplas gamas de temperatura e densidade, para sistemas contendo componentes muito distintos em termos de polaridade, tamanho e simetria; iii) a montagem e teste de uma instalação experimental para medir coeficientes de difusão em líquidos e fluidos supercríticos. Relativamente à modelação, uma nova expressão para coeficientes de difusão a diluição infinita de esferas rígidas foi desenvolvida e validada usando dados de dinâmica molecular (desvio relativo absoluto médio, AARD = 4.44%) Foram também estudados os coeficientes de difusão binários de sistemas reais. Para tal, foi compilada uma extensa base de dados de difusividades de sistemas reais em gases e solventes densos (622 sistemas binários num total de 9407 pontos experimentais e 358 moléculas) e a mesma foi usada na validação dos novos modelos desenvolvidos nesta tese. Um conjunto de novos modelos foi proposto para o cálculo de coeficientes de difusão a diluição infinita usando diferentes abordagens: i) dois modelos de base molecular com um parâmetro específico para cada sistema, aplicáveis em sistemas gasosos, líquidos e supercríticos, em que natureza do solvente se encontra limitada a apolar ou fracamente polar (AARDs globais na gama 4.26-4.40%); ii) dois modelos de base molecular biparamétricos, aplicáveis em todos os estados físicos, para qualquer tipo de soluto diluído em qualquer solvente (apolar, fracamente polar e polar). Ambos os modelos dão origem a erros globais entre 2.74% e 3.65%; iii) uma correlação com um parâmetro, específica para coeficientes de difusão em dióxido de carbono supercrítico (SC-CO2) e água líquida (AARD = 3.56%); iv) nove correlações empíricas e semi-empíricas que envolvem dois parâmetros, dependentes apenas da temperatura e/ou densidade do solvente e/ou viscosidade do solvente. Estes últimos modelos são muito simples e exibem excelentes resultados (AARDs entre 2.78% e 4.44%) em sistemas líquidos e supercríticos; e v) duas equações preditivas para difusividades de solutos em SC-CO2, em que os erros globais de ambas são inferiores a 6.80%. No global, deve realçar-se o facto de os novos modelos abrangerem a grande variedade de sistemas e moléculas geralmente encontrados. Os resultados obtidos são consistentemente melhores do que os obtidos com os modelos e abordagens encontrados na literatura. No caso das correlações com um ou dois parâmetros, mostrou-se que estes mesmos parâmetros podem ser ajustados usando um conjunto muito pequeno de dados, e posteriormente serem utilizados na previsão de valores de D 12 longe do conjunto original de pontos. Uma nova instalação experimental para medir coeficientes de difusão binários por técnicas cromatográficas foi montada e testada. O equipamento, o procedimento experimental e os cálculos analíticos necessários à obtenção dos valores de D 12 pelo método de abertura do pico cromatográfico, foram avaliados através da medição de difusividades de tolueno e acetona em SC-CO2. Seguidamente, foram medidos coeficientes de difusão de eucaliptol em SC-CO2 nas gamas de 202 – 252 bar e 313.15 – 333.15 K. Os resultados experimentais foram analisados através de correlações e modelos preditivos para D12.

Multilayered non-invasive temperature estimation from backscattered ultrasound

Roughly one in four breast cancer survivors report some degree of arm oedema. Lymphoedema is a build-up of excess lymph fluids in the tissues. Persistent lymphoedema leads to pain, diminished limb function, increased risk of infection, soft tissue fibrosis, and severe cases can be grossly disfiguring. From a mechanics perspective, the lymphoedemous tissue may be thought of as a two phase composite, consisting of both fluid and solid phases. Here we discuss the use of composites mixture theory to model the mechanics of lymphoedemous tissues. By treating the tissue as a fluid-solid composite, rules-of-mixtures may be used to estimate the effective moduli in terms of the properties of the individual components and their respective volume fractions in these two states.

Efficient non-linear data assimilation in geophysical fluid dynamics

New ways of combining observations with numerical models are discussed in which the size of the state space can be very large, and the model can be highly nonlinear. Also the observations of the system can be related to the model variables in highly nonlinear ways, making this data-assimilation (or inverse) problem highly nonlinear. First we discuss the connection between data assimilation and inverse problems, including regularization. We explore the choice of proposal density in a Particle Filter and show how the ’curse of dimensionality’ might be beaten. In the standard Particle Filter ensembles of model runs are propagated forward in time until observations are encountered, rendering it a pure Monte-Carlo method. In large-dimensional systems this is very inefficient and very large numbers of model runs are needed to solve the data-assimilation problem realistically. In our approach we steer all model runs towards the observations resulting in a much more efficient method. By further ’ensuring almost equal weight’ we avoid performing model runs that are useless in the end. Results are shown for the 40 and 1000 dimensional Lorenz 1995 model.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.