Fluid dynamic analysis and modelling of industrial chemical equipment

The research is aimed at contributing to the identification of reliable fully predictive Computational Fluid Dynamics (CFD) methods for the numerical simulation of equipment typically adopted in the chemical and process industries. The apparatuses selected for the investigation, specifically membrane modules, stirred vessels and fluidized beds, were characterized by a different and often complex fluid dynamic behaviour and in some cases the momentum transfer phenomena were coupled with mass transfer or multiphase interactions. Firs of all, a novel modelling approach based on CFD for the prediction of the gas separation process in membrane modules for hydrogen purification is developed. The reliability of the gas velocity field calculated numerically is assessed by comparison of the predictions with experimental velocity data collected by Particle Image Velocimetry, while the applicability of the model to properly predict the separation process under a wide range of operating conditions is assessed through a strict comparison with permeation experimental data. Then, the effect of numerical issues on the RANS-based predictions of single phase stirred tanks is analysed. The homogenisation process of a scalar tracer is also investigated and simulation results are compared to original passive tracer homogenisation curves determined with Planar Laser Induced Fluorescence. The capability of a CFD approach based on the solution of RANS equations is also investigated for describing the fluid dynamic characteristics of the dispersion of organics in water. Finally, an Eulerian-Eulerian fluid-dynamic model is used to simulate mono-disperse suspensions of Geldart A Group particles fluidized by a Newtonian incompressible fluid as well as binary segregating fluidized beds of particles differing in size and density. The results obtained under a number of different operating conditions are compared with literature experimental data and the effect of numerical uncertainties on axial segregation is also discussed.

Muscle ischaemia in patients with orthostatic hypotension assessed by velocity recovery cycles

Patients with orthostatic hypotension may experience neck pain radiating to the occipital region of the skull and the shoulders while standing (so-called coat-hanger ache). This study assessed muscle membrane potential in the trapezius muscle of patients with orthostatic hypotension and healthy subjects during head-up tilt (HUT), by measuring velocity recovery cycles (VRCs) of muscle action potentials as an indicator of muscle membrane potential.

A Stellar Census of the Tucana-Horologium Moving Group

We report the selection and spectroscopic confirmation of 129 new late-type (SpT = K3-M6) members of the Tucana-Horologium moving group, a nearby (d similar to 40 pc), young (tau similar to 40 Myr) population of comoving stars. We also report observations for 13 of the 17 known Tuc-Hor members in this spectral type range, and that 62 additional candidates are likely to be unassociated field stars; the confirmation frequency for new candidates is therefore 129/191 = 67%. We have used radial velocities, Ha emission, and Li-6708 absorption to distinguish between contaminants and bona fide members. Our expanded census of Tuc-Hor increases the known population by a factor of similar to 3 in total and by a factor of similar to 8 for members with SpT >= K3, but even so, the K-M dwarf population of Tuc-Hor is still markedly incomplete. Our expanded census allows for a much more detailed study of Tuc-Hor than was previously feasible. The spatial distribution of members appears to trace a two-dimensional sheet, with a broad distribution in X and Y, but a very narrow distribution (+/- 5 pc) in Z. The corresponding velocity distribution is very small, with a scatter of +/- 1.1 km s(-1) about the mean UVW velocity for stars spanning the entire 50 pc extent of Tuc-Hor. We also show that the isochronal age (tau similar to 20-30 Myr) and the lithium depletion boundary age (tau similar to 40 Myr) disagree, following the trend in other pre-main-sequence populations for isochrones to yield systematically younger ages. The H alpha emission line strength follows a trend of increasing equivalent width with later spectral type, as is seen for young clusters. We find that moving group members have been depleted of measurable lithium for spectral types of K7.0-M4.5. None of our targets have significant infrared excesses in the WISE W3 band, yielding an upper limit on warm debris disks of F < 0.7%. Finally, our purely kinematic and color-magnitude selection procedure allows us to test the efficiency and completeness for activity-based selection of young stars. We find that 60% of K-M dwarfs in Tuc-Hor do not have ROSAT counterparts and would have been omitted in X-ray-selected samples. In contrast, GALEX UV-selected samples using a previously suggested criterion for youth achieve completeness of 77% and purity of 78%, and we suggest new SpT-dependent selection criteria that will yield > 95% completeness for tau similar to 40 Myr populations with GALEX data available.

Pulmonary venous flow velocity patterns in 404 individuals without cardiovascular disease

OBJECTIVE To determine the pulmonary venous flow velocity (PVFV) values in a large normal population. DESIGN Prospective study in consecutive individuals. SETTING University hospital. METHODS Among 404 normal individuals, the flow velocity pattern in the right upper pulmonary vein was recorded in 315 subjects using transthoracic echocardiography, and in both upper pulmonary veins in 100 subjects using transoesophageal echocardiography. Subjects were divided into five age groups. The PVFV values were compared between transthoracic and transoesophageal echocardiography within the age groups, and intraindividually between the right and left upper pulmonary veins in transoesophageal echocardiography. RESULTS Normal PVFV values for the right upper pulmonary vein in transthoracic and transoesophageal echocardiography are presented. The duration of flow reversal at atrial contraction was overestimated using transthoracic echocardiography (mean (SD): 96 (21) ms in transoesophageal echocardiography, 120 (28) ms in transthoracic echocardiography, p < 0.0001). Systolic to diastolic peak flow velocity ratio (S:D) increased earlier with advancing age with transoesophageal echocardiography than with transthoracic echocardiography. Similar results were found for the corresponding time-velocity integrals. Data from the left and right upper pulmonary veins differed with respect to onset and deceleration of flow velocities, but not for flow durations or peak velocities. CONCLUSIONS Normal PVFV values generally show a wide range. The data presented will be of value in assessing left ventricular diastolic function and mitral regurgitation using the PVFV pattern.

A randomised double-blind, cross-over trial of 4-aminopyridine for downbeat nystagmus--effects on slowphase eye velocity, postural stability, locomotion and symptoms

Objective The effects of 4-aminopyridine (4-AP) on downbeat nystagmus (DBN) were analysed in terms of slow-phase velocity (SPV), stance, locomotion, visual acuity (VA), patient satisfaction and side effects using standardised questionnaires. Methods Twenty-seven patients with DBN received 5 mg 4-AP four times a day or placebo for 3 days and 10 mg 4-AP four times a day or placebo for 4 days. Recordings were done before the first, 60 min after the first and 60 min after the last drug administration. Results SPV decreased from 2.42 deg/s at baseline to 1.38 deg/s with 5 mg 4-AP and to 2.03 deg/s with 10 mg 4-AP (p<0.05; post hoc: 5 mg 4-AP: p=0.04). The rate of responders was 57%. Increasing age correlated with a 4-AP-related decrease in SPV (p<0.05). Patients improved in the ‘get-up-and-go test’ with 4-AP (p<0.001; post hoc: 5 mg: p=0.025; 10 mg: p<0.001). Tandem-walk time (both p<0.01) and tandem-walk error (4-AP: p=0.054; placebo: p=0.059) improved under 4-AP and placebo. Posturography showed that some patients improved with the 5 mg 4-AP dose, particularly older patients. Near VA increased from 0.59 at baseline to 0.66 with 5 mg 4-AP (p<0.05). Patients with idiopathic DBN had the greatest benefit from 4-AP. There were no differences between 4-AP and placebo regarding patient satisfaction and side effects. Conclusions 4-AP reduced SPV of DBN, improved near VA and some locomotor parameters. 4-AP is a useful medication for DBN syndrome, older patients in particular benefit from the effects of 5 mg 4-AP on nystagmus and postural stability.

Velocity of early BCR-ABL transcript elimination as an optimized predictor of outcome in chronic myeloid leukemia (CML) patients in chronic phase on treatment with imatinib

UNLABELLED Early assessment of response at 3 months of tyrosine kinase inhibitor treatment has become an important tool to predict favorable outcome. We sought to investigate the impact of relative changes of BCR-ABL transcript levels within the initial 3 months of therapy. In order to achieve accurate data for high BCR-ABL levels at diagnosis, beta glucuronidase (GUS) was used as a reference gene. Within the German CML-Study IV, samples of 408 imatinib-treated patients were available in a single laboratory for both times, diagnosis and 3 months on treatment. In total, 301 of these were treatment-naïve at sample collection. RESULTS (i) with regard to absolute transcript levels at diagnosis, no predictive cutoff could be identified; (ii) at 3 months, an individual reduction of BCR-ABL transcripts to the 0.35-fold of baseline level (0.46-log reduction, that is, roughly half-log) separated best (high risk: 16% of patients, 5-year overall survival (OS) 83% vs 98%, hazard ratio (HR) 6.3, P=0.001); (iii) at 3 months, a 6% BCR-ABL(IS) cutoff derived from BCR-ABL/GUS yielded a good and sensitive discrimination (high risk: 22% of patients, 5-year OS 85% vs 98%, HR 6.1, P=0.002). Patients at risk of disease progression can be identified precisely by the lack of a half-log reduction of BCR-ABL transcripts at 3 months.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.

An Intercomparison of Modelling Approaches for Simulating the Atmospheric Dispersion of Ammonia Emitted by Agricultural Sources

RESUMEN La dispersión del amoniaco (NH3) emitido por fuentes agrícolas en medias distancias, y su posterior deposición en el suelo y la vegetación, pueden llevar a la degradación de ecosistemas vulnerables y a la acidificación de los suelos. La deposición de NH3 suele ser mayor junto a la fuente emisora, por lo que los impactos negativos de dichas emisiones son generalmente mayores en esas zonas. Bajo la legislación comunitaria, varios estados miembros emplean modelos de dispersión inversa para estimar los impactos de las emisiones en las proximidades de las zonas naturales de especial conservación. Una revisión reciente de métodos para evaluar impactos de NH3 en distancias medias recomendaba la comparación de diferentes modelos para identificar diferencias importantes entre los métodos empleados por los distintos países de la UE. En base a esta recomendación, esta tesis doctoral compara y evalúa las predicciones de las concentraciones atmosféricas de NH3 de varios modelos bajo condiciones, tanto reales como hipotéticas, que plantean un potencial impacto sobre ecosistemas (incluidos aquellos bajo condiciones de clima Mediterráneo). En este sentido, se procedió además a la comparación y evaluación de varias técnicas de modelización inversa para inferir emisiones de NH3. Finalmente, se ha desarrollado un modelo matemático simple para calcular las concentraciones de NH3 y la velocidad de deposición de NH3 en ecosistemas vulnerables cercanos a una fuente emisora. La comparativa de modelos supuso la evaluación de cuatro modelos de dispersión (ADMS 4.1; AERMOD v07026; OPS-st v3.0.3 y LADD v2010) en un amplio rango de casos hipotéticos (dispersión de NH3 procedente de distintos tipos de fuentes agrícolas de emisión). La menor diferencia entre las concentraciones medias estimadas por los distintos modelos se obtuvo para escenarios simples. La convergencia entre las predicciones de los modelos fue mínima para el escenario relativo a la dispersión de NH3 procedente de un establo ventilado mecánicamente. En este caso, el modelo ADMS predijo concentraciones significativamente menores que los otros modelos. Una explicación de estas diferencias podríamos encontrarla en la interacción de diferentes “penachos” y “capas límite” durante el proceso de parametrización. Los cuatro modelos de dispersión fueron empleados para dos casos reales de dispersión de NH3: una granja de cerdos en Falster (Dinamarca) y otra en Carolina del Norte (EEUU). Las concentraciones medias anuales estimadas por los modelos fueron similares para el caso americano (emisión de granjas ventiladas de forma natural y balsa de purines). La comparación de las predicciones de los modelos con concentraciones medias anuales medidas in situ, así como la aplicación de los criterios establecidos para la aceptación estadística de los modelos, permitió concluir que los cuatro modelos se comportaron aceptablemente para este escenario. No ocurrió lo mismo en el caso danés (nave ventilada mecánicamente), en donde el modelo LADD no dio buenos resultados debido a la ausencia de procesos de “sobreelevacion de penacho” (plume-rise). Los modelos de dispersión dan a menudo pobres resultados en condiciones de baja velocidad de viento debido a que la teoría de dispersión en la que se basan no es aplicable en estas condiciones. En situaciones de frecuente descenso en la velocidad del viento, la actual guía de modelización propone usar un modelo que sea eficaz bajo dichas condiciones, máxime cuando se realice una valoración que tenga como objeto establecer una política de regularización. Esto puede no ser siempre posible debido a datos meteorológicos insuficientes, en cuyo caso la única opción sería utilizar un modelo más común, como la versión avanzada de los modelos Gausianos ADMS o AERMOD. Con el objetivo de evaluar la idoneidad de estos modelos para condiciones de bajas velocidades de viento, ambos modelos fueron utilizados en un caso con condiciones Mediterráneas. Lo que supone sucesivos periodos de baja velocidad del viento. El estudio se centró en la dispersión de NH3 procedente de una granja de cerdos en Segovia (España central). Para ello la concentración de NH3 media mensual fue medida en 21 localizaciones en torno a la granja. Se realizaron también medidas de concentración de alta resolución en una única localización durante una campaña de una semana. En este caso, se evaluaron dos estrategias para mejorar la respuesta del modelo ante bajas velocidades del viento. La primera se basó en “no zero wind” (NZW), que sustituyó periodos de calma con el mínimo límite de velocidad del viento y “accumulated calm emissions” (ACE), que forzaban al modelo a calcular las emisiones totales en un periodo de calma y la siguiente hora de no-calma. Debido a las importantes incertidumbres en los datos de entrada del modelo (inputs) (tasa de emisión de NH3, velocidad de salida de la fuente, parámetros de la capa límite, etc.), se utilizó el mismo caso para evaluar la incertidumbre en la predicción del modelo y valorar como dicha incertidumbre puede ser considerada en evaluaciones del modelo. Un modelo dinámico de emisión, modificado para el caso de clima Mediterráneo, fue empleado para estimar la variabilidad temporal en las emisiones de NH3. Así mismo, se realizó una comparativa utilizando las emisiones dinámicas y la tasa constante de emisión. La incertidumbre predicha asociada a la incertidumbre de los inputs fue de 67-98% del valor medio para el modelo ADMS y entre 53-83% del valor medio para AERMOD. La mayoría de esta incertidumbre se debió a la incertidumbre del ratio de emisión en la fuente (50%), seguida por la de las condiciones meteorológicas (10-20%) y aquella asociada a las velocidades de salida (5-10%). El modelo AERMOD predijo mayores concentraciones que ADMS y existieron más simulaciones que alcanzaron los criterios de aceptabilidad cuando se compararon las predicciones con las concentraciones medias anuales medidas. Sin embargo, las predicciones del modelo ADMS se correlacionaron espacialmente mejor con las mediciones. El uso de valores dinámicos de emisión estimados mejoró el comportamiento de ADMS, haciendo empeorar el de AERMOD. La aplicación de estrategias destinadas a mejorar el comportamiento de este último tuvo efectos contradictorios similares. Con el objeto de comparar distintas técnicas de modelización inversa, varios modelos (ADMS, LADD y WindTrax) fueron empleados para un caso no agrícola, una colonia de pingüinos en la Antártida. Este caso fue empleado para el estudio debido a que suponía la oportunidad de obtener el primer factor de emisión experimental para una colonia de pingüinos antárticos. Además las condiciones eran propicias desde el punto de vista de la casi total ausencia de concentraciones ambiente (background). Tras el trabajo de modelización existió una concordancia suficiente entre las estimaciones obtenidas por los tres modelos. De este modo se pudo definir un factor de emisión de para la colonia de 1.23 g NH3 por pareja criadora por día (con un rango de incertidumbre de 0.8-2.54 g NH3 por pareja criadora por día). Posteriores aplicaciones de técnicas de modelización inversa para casos agrícolas mostraron también un buen compromiso estadístico entre las emisiones estimadas por los distintos modelos. Con todo ello, es posible concluir que la modelización inversa es una técnica robusta para estimar tasas de emisión de NH3. Modelos de selección (screening) permiten obtener una rápida y aproximada estimación de los impactos medioambientales, siendo una herramienta útil para evaluaciones de impactos en tanto que permite eliminar casos que presentan un riesgo potencial de daño bajo. De esta forma, lo recursos del modelo pueden Resumen (Castellano) destinarse a casos en donde la posibilidad de daño es mayor. El modelo de Cálculo Simple de los Límites de Impacto de Amoniaco (SCAIL) se desarrolló para obtener una estimación de la concentración media de NH3 y de la tasa de deposición seca asociadas a una fuente agrícola. Está técnica de selección, basada en el modelo LADD, fue evaluada y calibrada con diferentes bases de datos y, finalmente, validada utilizando medidas independientes de concentraciones realizadas cerca de las fuentes. En general SCAIL dio buenos resultados de acuerdo a los criterios estadísticos establecidos. Este trabajo ha permitido definir situaciones en las que las concentraciones predichas por modelos de dispersión son similares, frente a otras en las que las predicciones difieren notablemente entre modelos. Algunos modelos nos están diseñados para simular determinados escenarios en tanto que no incluyen procesos relevantes o están más allá de los límites de su aplicabilidad. Un ejemplo es el modelo LADD que no es aplicable en fuentes con velocidad de salida significativa debido a que no incluye una parametrización de sobreelevacion del penacho. La evaluación de un esquema simple combinando la sobreelevacion del penacho y una turbulencia aumentada en la fuente mejoró el comportamiento del modelo. Sin embargo más pruebas son necesarias para avanzar en este sentido. Incluso modelos que son aplicables y que incluyen los procesos relevantes no siempre dan similares predicciones. Siendo las razones de esto aún desconocidas. Por ejemplo, AERMOD predice mayores concentraciones que ADMS para dispersión de NH3 procedente de naves de ganado ventiladas mecánicamente. Existe evidencia que sugiere que el modelo ADMS infraestima concentraciones en estas situaciones debido a un elevado límite de velocidad de viento. Por el contrario, existen evidencias de que AERMOD sobreestima concentraciones debido a sobreestimaciones a bajas Resumen (Castellano) velocidades de viento. Sin embrago, una modificación simple del pre-procesador meteorológico parece mejorar notablemente el comportamiento del modelo. Es de gran importancia que estas diferencias entre las predicciones de los modelos sean consideradas en los procesos de evaluación regulada por los organismos competentes. Esto puede ser realizado mediante la aplicación del modelo más útil para cada caso o, mejor aún, mediante modelos múltiples o híbridos. ABSTRACT Short-range atmospheric dispersion of ammonia (NH3) emitted by agricultural sources and its subsequent deposition to soil and vegetation can lead to the degradation of sensitive ecosystems and acidification of the soil. Atmospheric concentrations and dry deposition rates of NH3 are generally highest near the emission source and so environmental impacts to sensitive ecosystems are often largest at these locations. Under European legislation, several member states use short-range atmospheric dispersion models to estimate the impact of ammonia emissions on nearby designated nature conservation sites. A recent review of assessment methods for short-range impacts of NH3 recommended an intercomparison of the different models to identify whether there are notable differences to the assessment approaches used in different European countries. Based on this recommendation, this thesis compares and evaluates the atmospheric concentration predictions of several models used in these impact assessments for various real and hypothetical scenarios, including Mediterranean meteorological conditions. In addition, various inverse dispersion modelling techniques for the estimation of NH3 emissions rates are also compared and evaluated and a simple screening model to calculate the NH3 concentration and dry deposition rate at a sensitive ecosystem located close to an NH3 source was developed. The model intercomparison evaluated four atmospheric dispersion models (ADMS 4.1; AERMOD v07026; OPS-st v3.0.3 and LADD v2010) for a range of hypothetical case studies representing the atmospheric dispersion from several agricultural NH3 source types. The best agreement between the mean annual concentration predictions of the models was found for simple scenarios with area and volume sources. The agreement between the predictions of the models was worst for the scenario representing the dispersion from a mechanically ventilated livestock house, for which ADMS predicted significantly smaller concentrations than the other models. The reason for these differences appears to be due to the interaction of different plume-rise and boundary layer parameterisations. All four dispersion models were applied to two real case studies of dispersion of NH3 from pig farms in Falster (Denmark) and North Carolina (USA). The mean annual concentration predictions of the models were similar for the USA case study (emissions from naturally ventilated pig houses and a slurry lagoon). The comparison of model predictions with mean annual measured concentrations and the application of established statistical model acceptability criteria concluded that all four models performed acceptably for this case study. This was not the case for the Danish case study (mechanically ventilated pig house) for which the LADD model did not perform acceptably due to the lack of plume-rise processes in the model. Regulatory dispersion models often perform poorly in low wind speed conditions due to the model dispersion theory being inapplicable at low wind speeds. For situations with frequent low wind speed periods, current modelling guidance for regulatory assessments is to use a model that can handle these conditions in an acceptable way. This may not always be possible due to insufficient meteorological data and so the only option may be to carry out the assessment using a more common regulatory model, such as the advanced Gaussian models ADMS or AERMOD. In order to assess the suitability of these models for low wind conditions, they were applied to a Mediterranean case study that included many periods of low wind speed. The case study was the dispersion of NH3 emitted by a pig farm in Segovia, Central Spain, for which mean monthly atmospheric NH3 concentration measurements were made at 21 locations surrounding the farm as well as high-temporal-resolution concentration measurements at one location during a one-week campaign. Two strategies to improve the model performance for low wind speed conditions were tested. These were ‘no zero wind’ (NZW), which replaced calm periods with the minimum threshold wind speed of the model and ‘accumulated calm emissions’ (ACE), which forced the model to emit the total emissions during a calm period during the first subsequent non-calm hour. Due to large uncertainties in the model input data (NH3 emission rates, source exit velocities, boundary layer parameters), the case study was also used to assess model prediction uncertainty and assess how this uncertainty can be taken into account in model evaluations. A dynamic emission model modified for the Mediterranean climate was used to estimate the temporal variability in NH3 emission rates and a comparison was made between the simulations using the dynamic emissions and a constant emission rate. Prediction uncertainty due to model input uncertainty was 67-98% of the mean value for ADMS and between 53-83% of the mean value for AERMOD. Most of this uncertainty was due to source emission rate uncertainty (~50%), followed by uncertainty in the meteorological conditions (~10-20%) and uncertainty in exit velocities (~5-10%). AERMOD predicted higher concentrations than ADMS and more of the simulations met the model acceptability criteria when compared with the annual mean measured concentrations. However, the ADMS predictions were better correlated spatially with the measurements. The use of dynamic emission estimates improved the performance of ADMS but worsened the performance of AERMOD and the application of strategies to improved model performance had similar contradictory effects. In order to compare different inverse modelling techniques, several models (ADMS, LADD and WindTrax) were applied to a non-agricultural case study of a penguin colony in Antarctica. This case study was used since it gave the opportunity to provide the first experimentally-derived emission factor for an Antarctic penguin colony and also had the advantage of negligible background concentrations. There was sufficient agreement between the emission estimates obtained from the three models to define an emission factor for the penguin colony (1.23 g NH3 per breeding pair per day with an uncertainty range of 0.8-2.54 g NH3 per breeding pair per day). This emission estimate compared favourably to the value obtained using a simple micrometeorological technique (aerodynamic gradient) of 0.98 g ammonia per breeding pair per day (95% confidence interval: 0.2-2.4 g ammonia per breeding pair per day). Further application of the inverse modelling techniques for a range of agricultural case studies also demonstrated good agreement between the emission estimates. It is concluded, therefore, that inverse dispersion modelling is a robust technique for estimating NH3 emission rates. Screening models that can provide a quick and approximate estimate of environmental impacts are a useful tool for impact assessments because they can be used to filter out cases that potentially have a minimal environmental impact allowing resources to be focussed on more potentially damaging cases. The Simple Calculation of Ammonia Impact Limits (SCAIL) model was developed as a screening model to provide an estimate of the mean NH3 concentration and dry deposition rate downwind of an agricultural source. This screening tool, based on the LADD model, was evaluated and calibrated with several experimental datasets and then validated using independent concentration measurements made near sources. Overall SCAIL performed acceptably according to established statistical criteria. This work has identified situations where the concentration predictions of dispersion models are similar and other situations where the predictions are significantly different. Some models are simply not designed to simulate certain scenarios since they do not include the relevant processes or are beyond the limits of their applicability. An example is the LADD model that is not applicable to sources with significant exit velocity since the model does not include a plume-rise parameterisation. The testing of a simple scheme combining a momentum-driven plume rise and increased turbulence at the source improved model performance, but more testing is required. Even models that are applicable and include the relevant process do not always give similar predictions and the reasons for this need to be investigated. AERMOD for example predicts higher concentrations than ADMS for dispersion from mechanically ventilated livestock housing. There is evidence to suggest that ADMS underestimates concentrations in these situations due to a high wind speed threshold. Conversely, there is also evidence that AERMOD overestimates concentrations in these situations due to overestimation at low wind speeds. However, a simple modification to the meteorological pre-processor appears to improve the performance of the model. It is important that these differences between the predictions of these models are taken into account in regulatory assessments. This can be done by applying the most suitable model for the assessment in question or, better still, using multiple or hybrid models.

A non-perturbative renormalization group study of the stochastic NavierStokes equation

We study the renormalization group flow of the average action of the stochastic Navier-Stokes equation with power-law forcing. Using Galilean invariance, we introduce a nonperturbative approximation adapted to the zero-frequency sector of the theory in the parametric range of the Hölder exponent 4−2 ɛ of the forcing where real-space local interactions are relevant. In any spatial dimension d, we observe the convergence of the resulting renormalization group flow to a unique fixed point which yields a kinetic energy spectrum scaling in agreement with canonical dimension analysis. Kolmogorov's −5/3 law is, thus, recovered for ɛ=2 as also predicted by perturbative renormalization. At variance with the perturbative prediction, the −5/3 law emerges in the presence of a saturation in the ɛ dependence of the scaling dimension of the eddy diffusivity at ɛ=3/2 when, according to perturbative renormalization, the velocity field becomes infrared relevant.

Comparison of particle dispersion obtained with simple stochastic models and les in pipe flow

A hybrid Eulerian-Lagrangian approach is employed to simulate heavy particle dispersion in turbulent pipe flow. The mean flow is provided by the Eulerian simulations developed by mean of JetCode, whereas the fluid fluctuations seen by particles are prescribed by a stochastic differential equation based on normalized Langevin. The statistics of particle velocity are compared to LES data which contain detailed statistics of velocity for particles with diameter equal to 20.4 µm. The model is in good agreement with the LES data for axial mean velocity whereas rms of axial and radial velocities should be adjusted.

Numerical simulation of group combustion of pulverized coal

A mathematical model for the group combustion of pulverized coal particles was developed in a previous work. It includes the Lagrangian description of the dehumidification, devolatilization and char gasification reactions of the coal particles in the homogenized gaseous environment resulting from the three fuels, CO, H2 and volatiles, supplied by the gasification of the particles and their simultaneous group combustion by the gas phase oxidation reactions, which are considered to be very fast. This model is complemented here with an analysis of the particle dynamics, determined principally by the effects of aerodynamic drag and gravity, and its dispersion based on a stochastic model. It is also extended to include two other simpler models for the gasification of the particles: the first one for particles small enough to extinguish the surrounding diffusion flames, and a second one for particles with small ash content when the porous shell of ashes remaining after gasification of the char, non structurally stable, is disrupted. As an example of the applicability of the models, they are used in the numerical simulation of an experiment of a non-swirling pulverized coal jet with a nearly stagnant air at ambient temperature, with an initial region of interaction with a small annular methane flame. Computational algorithms for solving the different stages undergone by a coal particle during its combustion are proposed. For the partial differential equations modeling the gas phase, a second order finite element method combined with a semi-Lagrangian characteristics method are used. The results obtained with the three versions of the model are compared among them and show how the first of the simpler models fits better the experimental results.

Finite element simulation of dispersion in the Bay of Santander

Two mathematical models are used to simulate pollution in the Bay of Santander. The first is the hydrodynamic model that provides the velocity field and height of the water. The second gives the pollutant concentration field as a resultant. Both models are formulated in two-dimensional equations. Linear triangular finite elements are used in the Galerkin procedure for spatial discretization. A finite difference scheme is used for the time integration. At each time step the calculated results of the first model are input to the second model as field data. The efficiency and accuracy of the models are tested by their application to a simple illustrative example. Finally a case study in simulation of pollution evolution in the Bay of Santander is presented

Caracterização física e química da pólvora

O presente trabalho consiste na caracterização física e química da pólvora. Esta caracterização foi realizada para alguns tipos de pólvora, com o objetivo de se poder empregar este tipo de material noutros âmbitos, que não seja só nas Forças Armadas, nomeadamente, no armamento. A caracterização física abrangeu, essencialmente, a caracterização morfológica das amostras tal-qual, nomeadamente, as pólvoras multi-perfurada, tubular, tubular (rocket), cilíndrica, esférica, lamelar, em fita e a pólvora negra. A técnica utilizada foi a observação através de uma lupa estereoscópica. Após a combustão, foi utilizado o microscópio eletrónico de varrimento. A caracterização química foi realizada no âmbito da análise química elementar, e também no âmbito da combustão, às condições atmosféricas. Na análise química elementar, foram estudadas a pólvora multiperfurada, tubular, tubular rocket e em fita, por intermédio da espectrometria de fluorescência de raios-X - dispersão de energia e por espectrometria de absorção atómica de chama. No âmbito da combustão, foram estudas a taxa de queima e a velocidade de propagação de chama, nas amostras de pólvora multi-perfurada e a pólvora de fita, através de técnicas de medição de massa e de visualização de chama. Na discussão de resultados constatou-se que a maioria dos tipos de pólvora estudados pertencem ao grupo dos propelentes de base dupla. Apesar da variabilidade entre amostras, verificou-se que o principal elemento comum é o chumbo. Quanto à taxa de queima, esta apresenta uma evolução aproximadamente linear em todas as amostras. Foi, ainda, apresentada uma velocidade de propagação de chama característica para os dois tipos de pólvora estudados, tendo sido estabelecida para a pólvora multi-perfurada uma velocidade SR = 1,1 mm/s, para em fita, SR = 6,0 mm/s.

Performance characterization and relative dispersion of drogued drifters deployed in the fjords of Nootka Sound

Senior thesis written for Oceanography 445

Comparison of Lagrangian Coherent Structures and Relative Dispersion in a Mixing Layer

Thesis (Master's)--University of Washington, 2016-06