Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

大气湍流对斜程传输准单色高斯谢尔光束时间相干性的影响

由第一类零阶贝塞尔函数的级数展开推导出波结构函数在任意湍流条件下的近似表达式。由广义惠更斯菲涅耳原理、随高度变化的Hufnagel-Valley湍流廓线模型以及波结构函数在任意湍流条件下的近似表达式,导出了斜程传输时准单色高斯谢尔光束互相干函数的解析式。然后,利用表征光束时间相干性的纵向相干长度(可由互相干函数导出),研究了斜程传输时大气湍流对准单色高斯谢尔光束时间相干性的影响。研究结果表明,准单色高斯谢尔光束的时间相干性在整个斜程传输过程中保持不变。最后,对该结果在物理上给予了定性解释。

大气湍流对斜程传输准单色高斯-谢尔光束空间相干性的影响

由湍流大气中斜程传输时准单色高斯-谢尔(GSM)光束互相干函数的解析式导出了该光束的复相干度.然后,利用表征光束空间相干性的横向相干长度,研究了斜程传输时大气湍流对准单色GSM光束空间相干性的影响.研究结果表明:1)当传输路径偏离水平方向较大(即θ≤88°)时,准单色GSM光束横向相干长度随传输距离均为先迅速增加,后缓慢增加,最后基本保持不变.2)当传输路径接近水平方向(即θ≥89°)时,准单色GSM光束横向相干长度随传输距离均为先增大,达到一个最大值后开始下降并持续减小.3)空间相干性不同的准单色GSM光束在湍流大气中传输一定距离后空间相干性趋于相同.最后,对以上结果在物理上给予了解释.

A study of flow conditions near the substratum in an experimental channel [Translation from: Verhandlungen der Internationalen Vereinigung fur Theoretische und Angewandte Limnologie 18, 718-725, 1972]

In a slow flow, on a smooth uniform substratum, a limited bed allows the existence of currents slow enough for benthic invertebrates. These conditions rarely occur naturally. The investigations carried out in this work aimed, on an intermediary scale, to define the influence of irregularities in the substratum on flow near the bottom. The substrata used were made of glass marbles. The investigations were carried out in a transparent channel of 70 cm in length and a rectangular section 10 x 5 cm. The data was analysed to study the general evolution of flow in terms of average speeds and the appearance of the turbulence near the bottom.

Mean-speed measurements in two-dimensional, incompressible, fully-developed turbulent channel flow

Mean velocity profiles were measured in the 5” x 60” wind channel of the turbulence laboratory at the GALCIT, by the use of a hot-wire anemometer. The repeatability of results was established, and the accuracy of the instrumentation estimated. Scatter of experimental results is a little, if any, beyond this limit, although some effects might be expected to arise from variations in atmospheric humidity, no account of this factor having been taken in the present work. Also, slight unsteadiness in flow conditions will be responsible for some scatter.

Irregularities of a hot-wire in close proximity to a solid boundary at low speeds were observed, as have already been found by others.

That Kármán’s logarithmic law holds reasonably well over the main part of a fully developed turbulent flow was checked, the equation u/u_t = 6.0 + 6.25 log₁₀ y^ut/v being obtained, and, as has been previously the case, the experimental points do not quite form one straight line in the region where viscosity effects are small. The values of the constants for this law for the best over-all agreement were determined and compared with those obtained by others.

The range of Reynolds numbers used (based on half-width of channel) was from 20,000 to 60,000.

Simulations and mechanisms of subtropical low-cloud response to climate change

This thesis focuses on improving the simulation skills and the theoretical understanding of the subtropical low cloud response to climate change.

First, an energetically consistent forcing framework is designed and implemented for the large eddy simulation (LES) of the low-cloud response to climate change. The three representative current-day subtropical low cloud regimes of cumulus (Cu), cumulus-over-stratocumulus, and stratocumulus (Sc) are all well simulated with this framework, and results are comparable to the conventional fixed-SST approach. However, the cumulus response to climate warming subject to energetic constraints differs significantly from the conventional approach with fixed SST. Under the energetic constraint, the subtropics warm less than the tropics, since longwave (LW) cooling is more efficient with the drier subtropical free troposphere. The surface latent heat flux (LHF) also increases only weakly subject to the surface energetic constraint. Both factors contribute to an increased estimated inversion strength (EIS), and decreased inversion height. The decreased Cu-depth contributes to a decrease of liquid water path (LWP) and weak positive cloud feedback. The conventional fixed-SST approach instead simulates a strong increase in LHF and deepening of the Cu layer, leading to a weakly negative cloud feedback. This illustrates the importance of energetic constraints to the simulation and understanding of the sign and magnitude of low-cloud feedback.

Second, an extended eddy-diffusivity mass-flux (EDMF) closure for the unified representation of sub-grid scale (SGS) turbulence and convection processes in general circulation models (GCM) is presented. The inclusion of prognostic terms and the elimination of the infinitesimal updraft fraction assumption makes it more flexible for implementation in models across different scales. This framework can be consistently extended to formulate multiple updrafts and downdrafts, as well as variances and covariances. It has been verified with LES in different boundary layer regimes in the current climate, and further development and implementation of this closure may help to improve our simulation skills and understanding of low-cloud feedback through GCMs.

Estudio de la precisión de varios modelos de turbulencia en la aplicación a perfiles aerodinámicos

[ES]El objetivo principal de este proyecto se centra en modelizar correctamente la capa límite sobre el perfil alar donde se produce la transición del régimen laminar al régimen transitorio. Como objetivo secundario se encuentra el afianzamiento de las bases teóricas de mecánica de fluidos obtenidas en la escuela y la adquisición de más conocimientos relacionados con la aerodinámica, concretamente con la capa límite. En una primera parte se tratarán los conceptos generales de los perfiles alares y se hará una breve introducción a los distintos tipos de mallado existentes. También se explicará el concepto de capa límite y todo lo relacionado con ella. A continuación, se establecerán los criterios de selección del modelo de turbulencia más adecuado y se mostrarán los resultados obtenidos de los distintos tipos de modelos de turbulencia anteriormente mencionados. Una vez seleccionado un modelo de turbulencia se profundizará en su estudio, aplicándolo a varios perfiles NACA. Se analizarán los resultados obtenidos y los errores y se buscarán posibles soluciones. Finalmente, se procederá a sacar las conclusiones del modelo escogido y se comparará con una serie de ensayos experimentales con objeto de poder validarlo.

Turbulent buoyant jets into stratified or flowing ambient fluids

Theoretical and experimental studies were made on two classes of buoyant jet problems, namely:

1) an inclined, round buoyant yet in a stagnant environment with linear density-stratification;

2) a round buoyant jet in a uniform cross stream of homogenous density.

Using the integral technique of analysis, assuming similarity, predictions can be made for jet trajectory, widths, and dilution ratios, in a density-stratified or flowing environment. Such information is of great importance in the design of disposal systems for sewage effluent into the ocean or waste gases into the atmosphere.

The present study of a buoyant jet in a stagnant environment has extended the Morton type of analysis to cover the effect of the initial angle of discharge. Numerical solutions have been presented for a range of initial conditions. Laboratory experiments were conducted for photographic observations of the trajectories of dyed jets. In general the observed jet forms agreed well with the calculated trajectories and nominal half widths when the value of the entrainment coefficient was taken to be α = 0.082, as previously suggested by Morton.

The problem of a buoyant jet in a uniform cross stream was analyzed by assuming an entrainment mechanism based upon the vector difference between the characteristic jet velocity and the ambient velocity. The effect of the unbalanced pressure field on the sides of the jet flow was approximated by a gross drag term. Laboratory flume experiments with sinking jets which are directly analogous to buoyant jets were performed. Salt solutions were injected into fresh water at the free surface in a flume. The jet trajectories, dilution ratios and jet half widths were determined by conductivity measurements. The entrainment coefficient, α, and drag coefficient, C_d, were found from the observed jet trajectories and dilution ratios. In the ten cases studied where jet Froude number ranged from 10 to 80 and velocity ratio (jet: current) K from 4 to 16, α varied from 0.4 to 0.5 and C_d from 1.7 to 0.1. The jet mixing motion for distance within 250D was found to be dominated by the self-generated turbulence, rather than the free-stream turbulence. Similarity of concentration profiles has also been discussed.

Supernovae explosions induced by pair production instability

Stars with a core mass greater than about 30 M_⊙ become dynamically unstable due to electron-positron pair production when their central temperature reaches 1.5-2.0 x 10^{9 0}K. The collapse and subsequent explosion of stars with core masses of 45, 52, and 60 M_⊙ is calculated. The range of the final velocity of expansion (3,400 – 8,500 km/sec) and of the mass ejected (1 – 40 M_⊙) is comparable to that observed for type II supernovae.

An implicit scheme of hydrodynamic difference equations (stable for large time steps) used for the calculation of the evolution is described.

For fast evolution the turbulence caused by convective instability does not produce the zero entropy gradient and perfect mixing found for slower evolution. A dynamical model of the convection is derived from the equations of motion and then incorporated into the difference equations.

The impact of Cow Green Reservoir on invertebrate populations in the River Tees

Cow Green is a new reservoir situated in Pennine moorland. It has an area of 312 ha, a capacity of 40 . 9 x 10 SUP-6 m SUP-3 and a maximum depth of 22 . 8 m. The function of the reservoir is to regulate flow in the River Tees to provide industrial Teesside with sufficient water during the dry spells. Invertebrate studies were carried out in the Tees to monitor changes resulting from the construction of the reservoir both in the flooded basin and below the dam. The overall effect of the reservoir on the Tees has been to increase the numbers and biomass of certain taxa, but generally not at the expense of previous fauna. Some of the positive effects, ie. increase in number and biomass, and maintenance of faunal diversity, may in part be attributable to the presence of the rapids and waterfall. Turbulence resulting from this rapid flow over heterogeneous bottom is sufficient to prevent clogging of interstitial spaces by silt and to maintain the variety of ecological niches necessary for a diverse fauna.

Changes in the coherence of quasi-monochromatic electromagnetic Gaussian Schell-model beams propagating through turbulent atmosphere

Based on the extended Huygens-Fresnel principle, the mutual coherence function of quasi-monochromatic electromagnetic Gaussian Schell-model (EGSM) beams propagating through turbulent atmosphere is derived analytically. By employing the lateral and the longitudinal coherence length of EGSM beams to characterize the spatial and the temporal coherence of the beams, the behavior of changes in the spatial and the temporal coherence of those beams is studied. The results show that with a fixed set of beam parameters and under particular atmospheric turbulence model, the lateral coherence of an EGSM beam reaches its maximum value as the beam propagates a certain distance in the turbulent atmosphere, then it begins degrading and keeps decreasing along with the further distance. However, the longitudinal coherence length of an EGSM beam keeps unchanging in this propagation. Lastly, a qualitative explanation is given to these results. (c) 2007 Elsevier B.V. All rights reserved.

Physiological ecology of the ciliated protozoon Loxodes

Loxodes faces special problems in living close to the oxic-anoxic boundary. In tightly-stratified ponds like Priest Pot its optimum environment may be quite narrow and it can be displaced by the slightest turbulence. Loxodes cannot sense an O sub(2) gradient directly but its ability to perceive gravity allows it to make relatively long vertical migrations. It is also sensitive to light and oxygen and it uses these environmental cues to modulate the parameters of its random motility: in the dark, it aggregates at a low O sub(2) tension and in bright light it aggregates in anoxic water. The oxic-anoxic boundary is also a zone where O sub(2) may be a scarce and transient resource, but Loxodes) can switch to nitrate respiration and exploit the pool of nitrate that often exists close to the base of the oxycline.

Observação do transporte de sedimentos em suspensão ao longo do Canal Norte do rio Amazonas durante condições de baixa descarga (outubro 2008)

The Amazon river, located in northernBrazil, discharges between 80,000 and 250,000 m3s-1 of water onto the adjacent shelf, creating a plume of brackish water that extends hundreds of kilometers away from the river mouth. This river also carries a large amount of fine sediments to the ocean where fluid mud has been found in the topset and upper foreset layers of the subaqueous delta formed on the mid-shelf. One of the main goals of this dissertation is to describe how turbulence and suspended sediment concentration vary along the Northern Channel of the Amazon river. Water column measurements were carried out in October 2008 at six anchor stations (P1, P3, P5, P6, P8 e P9) located seaward of the river mouth; P1 and P9 were 125 km apart. Each station was occupied during 13 hours during which current speed and direction were continuously sampled with a 600 kHz Teledyne-RDI ADCP; hourly profiles of temperature, salinity, turbidity and depth were also obtained. Water samples were collected for determination of Suspended Particulate Matter (SPM) concentration and calibration of the turbidity sensor. Current speed reached values above 1.5 m s1 in the along-channel direction (NE-SW); a remarkable ebb-flood asymmetry was observed and flows were strongly ebb-dominated. Throughout the water column, SPM concentration at stations P1 and P3 varied between 100 and 300 mg L1 in association with the presence of freshwater. In contrast, a strong salinity gradient was observed between stations P6 and P9, coinciding with the occurrence of concentrations of SPM above 10 g L-1 (fluid mud). At stations P3, P5 and P6, interface between freshwater from the Amazon river and salt water from the continental shelf, shear stresses wereestimated through four diferents methods: Reynolds, Turbulent Kinetic Energy (TKE), modified TKE and Quadratic Law; in the nearbed region (3 mab) the computed values varied between 0 and 3 Pa. At the three stations (P3, P5 and P6) the lowest and the highest shear stress values were obtained through, respectively, the Reynolds and the TKE methods. Over the whole water column turbulence intensity was estimated through the standard deviation of the turbulent component of the along-channel current velocity (root-mean square of u); from these values, it was estimated the turbulent dissipation of energy (G), whose values at 3 mab varied between zero and 20 s1.