On the development of the convective boundary layer in a shear-free thermally forced stably stratified fluid setting: a 2D and 3D experimental investigation using image analysis techniques coupled with temperature measurements

The motion of buoyancy driven plumes is, on all scales, the most common heat and momentum transfer mechanism in geophysical flows, well known as Free Convection. Similarly, density stratification due to heating inequalities is also an ordinary scenario in nature. Free Convection phenomenon coupled with a density stratified fluid setting leads to the so-called Penetrative Free Convection (PFC). When a fluid, in static equilibrium, is stably stratified a thermal forcing can produce an unstable configuration ensuing internal waves formation of increasing amplitude. If the perturbation is strong enough, it can definitely erode the initial stratification and cause the motion of turbulent buoyant updrafts, dome-shaped, compensated by denser downdrafts. The entrainment phenomenon occurring at the interface between the turbulent and non-turbulent region justifies the penetrative feature of convection and causes the non linear growth of the Convective Boundary Layer (CBL) of well mixed fluid against the adjacent stably stratified region. In addition to the wide engineering applications, the environmental impact mostly motivates PFC studies. The upper lakes and oceans, under calm conditions, usually exhibit a continuous, moderately stable density distribution. Turbulent convective flow can be generated both by the free-surface cooling and wind shear-stress, eroding the stable stratification on a daily or seasonal time scale. Domes with large downward velocities are generated at the free surface, balanced by updrafts with lower velocity but larger area. Because of the relatively rapid mixing, the density distribution is approximately uniform in the upper layer and it deepens with time as a result of the entrainment and erosion of the underlying denser water. An analogous phenomenon is observed in the lower troposphere when surface heating due to solar radiation results in a growing unstable layer adjacent to the ground which replaces a nocturnal inversion from below. In this case, the initially stable environment near the ground is affected by convection characterized by relatively narrow and fast plumes of rising horizontal surfaces balanced by larger regions of downward slower motion. Resulting internal waves generated within the stable layer take place at or below the Brunt-Väisälä frequency, which is related to the vertical temperature gradient. In nature, the dynamics of the CBL influences the transport and mixing features of a given stratified fluid-body. The amount of materials being mixed due to penetrative convection is a crucial issue both in water or air quality monitoring and forecast with important implications in environmentally-friendly studies. Focusing on the environmental aspect, the pollutant dispersion is a matter of particular importance. On one hand, mixing processes inside the CBL help increasing dispersion with some positive consequences: the concentration of potential harmful pollutants in high risk zones tends to decrease, the turnover and the redistribution of vital substances, like oxygen and nutrients, is guaranteed; the latter plays a major role in large water bodies. These advantages are essential for the safety of populations living close to urban or industrial areas, or for preserving coastal human activities and ecosystems. On the other hand each pollutant, released inside the CBL mainly by human activities, remains confined inside it because of the interface with the non-turbulent region, which acts as a barrier for outward transport. Considering the proximity of the CBL with biosphere, a deeper insight into dispersion and entrainment processes appears mandatory either for sustainable engineering design or for monitoring purposes. Given the great applicability of the topic in several scientific and engineering fields, large amount of experimental, theoretical and numerical investigations on CBL development in a continuously and linearly stratified fluid setting had been conducted in the past since sixties. It appears there is a lack of consensus about the dependence of integral parameters of convective entrainment (in particular the CBL growth rate) on the initial stratification strength and convective phenomenon evolution. Driven by all these premises the main aims of the present research are related to better understanding dispersion of a passive scalar inside the CBL with a lagrangian, non-local approach and a fully three-dimensional (3D 3C) experimental technique, which for the first time has been applied to the topic. When turbulent convection occurs, in fact, dispersion is mostly due to transport by large organized structures while molecular diffusion can be neglected. Given this assumption, a non local approach based on a fully advective-like behaviour of the tracer is necessary and a lagrangian description of flow may be more suitable. Moreover, the knowledge of the horizontal and vertical extension of the structures dominating the flow field appears to be mandatory. Furthermore, turbulence is fully three-dimensional on the scales of motion characterising the phenomena in nature. In order to better understanding and likely describing the evolution of turbulent structures inside the convective layer, a three dimensional experimental technique is strongly required. In the present work shear free convection experiments in a stably stratified environment were performed in a thermally controlled convection chamber. The experimental set up was arranged in different optical configurations to ensure the two components and three components of velocity to be measured in an illuminated plane and volume respectively through different image analysis techniques with sub-pixel accuracy. Temperature measurements were simultaneously acquired using thermocouples of accuracy less than 0.1 °C. The experiments highlighted the time evolution of the convective structure characteristic spatial scales (CBL height and horizontal spacing between thermals) and the non local description of transport and mixing inside the growing CBL. The main novelty of the present contribution covers the improvement of techniques and methods to achieve more reliable, statistically robust and likely results. The experimental effort based on image analysis techniques (mainly Feature Tracking and photogrammetric 3D-PTV) resulted suitable for reconstructing longer trajectories (always more than 103 trajectories longer than 20 consecutive snapshots) and analyzing larger particle density images (reliable results for more than 2000 particles in a volume of 15X15X15 cm3 for 3D-PTV). Velocity statistics are then more robust than those from classical Particle Tracking Velocimetry. Moreover, for the first time, a fully three-dimensional particle tracking technique has been applied to penetrative convection experiments. 3D-PTV allows a more realistic description of the velocity field, which occurs during the evolution of the convective mixed layer, than more traditional 2D techniques. Furthermore, photogrammetric 3D-PTV rather than “scanning” 3D-PTV results in more accuracy when the tracer particle density is high, because particles may be tracked directly in the 3D space rather than through matching of 2D projections. The broader impact of the research mainly refers to the prediction of the CBL growth as a function of initial and boundary conditions with better accuracy than conventional and well established techniques. The experimental study can thus give a positive contribution on real pollutant dispersion studies in urban and natural environments for environmental protection and sustainable design purposes. Field experiments aimed at measuring the turbulence budget of the CBL have shown that the mechanical generation of kinetic energy by wind shear is often confined close to the heat source supporting the validity of laboratory models in which no wind is present. According to this assumption, the similarity proposed by Deardorff (1970) is employed to compute scaling parameters and to make results comparable with real scales. Through normalizing the quantities measured at different stages of the experiment, the phenomenon can be considered as a succession of steady states, according to an evolution of the variables of interest that may be defined quasi-steady state. The experimental apparatus employed to run the experiments is the same as in Cenedese and Querzoli (1994), Querzoli (1996), Cenedese and Querzoli (1997) and Moroni and Cenedese (2006). The spatial resolution of velocity data is largely increased here by means of 2D and 3D image analysis techniques (Feature Tracking, FT, and 3D Particle Tracking Velocimetry, 3D-PTV) used instead of Laser-Doppler Anemometry or 2D Particle Tracking Velocimetry as in Cenedese and Querzoli, 1994; Querzoli, 1996 and Cenedese and Querzoli, 1997. Moreover the photogrammetric 3D-PTV here applied allows fully three-dimensional descriptions of both the Eulerian velocity field and Lagrangian particle trajectories for a more likely understanding of the phenomenon than through the 2D approach used in Moroni and Cenedese (2006). Furthermore, the combined use of thermocouples and flow visualization techniques allows cross-validating different methods to estimate the evolutions of the key parameters and the plume characteristic dimensions. Two different experimental arrangements were set-up, 2D and a 3D models. A large set of data were firstly acquired using a 2D model and employing FT. The preliminary investigation was mainly focused on better understanding the physics of the phenomenon, finding a time scaling law, testing different methods to compute the variables of interest and comparing results with classical methods found in literature. When 2D techniques are employed to detect the velocity field, the flow is illuminated with a thin light sheet and only the velocity components within this sheet can be evaluated. Driven by the idea that only a fully three-dimensional technique can significantly improve our laboratory model in term of a more likely description of free convection structures we run a second set of experiments, by using a stereoscopic arrangement of cameras focused on an illuminated volume. Although some methods do exist for reconstructing 3D velocities in a point (3D laser Doppler Anemometry; Hinsch and Hinrichs, 1996) or plane (3D stereo-PIV; Stuer et al., 1999), only a fully 3D technique based on the illumination of a flow volume rather than a flow sheet will give the information needed to construct the instantaneous 3D velocity fields. A number of imaging-based measurement techniques exist for determining 3D velocity fields in an observation volume. Among these are: scanning, photogrammetric, holographic or photogrammetric techniques dependending on which principle is recalled to reconstruct the third dimension from a 2D image/s: The present study was focused on 3D-PTV which is a 3D extension of the 2D particle-tracking methods. 3D-PTV is based on reconstructing 3D trajectories of reflecting tracer particles through a photogrammetric recording of image sequences. The 3D particle trajectories obtained can be used to calculate the 3D velocity field. The 3D-PTV optical system has been designed with the following capabilities: image a volume far away the boundary walls, lengthen the trajectories, and improve the accuracy of the procedure through a careful test on synthetically generated data. A physically-based photogrammetric calibration of the stereoscopic arrangement was employed and its accuracy tested. The effects of multimedia geometry on calibration parameters were taken into account. The combination of image- and object-space based information was employed to establish the correspondences between particle positions (structure from stereo reconstruction). A particle tracking algorithm was then employed to reconstructed 3D trajectories. Sensitivity tests conducted on the matching algorithm proved that the calibration accuracy is fundamental to obtain the correct matching and particle tracking; small errors in calibration parameters or neglecting water refraction effects reduce matching performance. On the other hand, accuracy less than 1 pixel was reached with our calibration procedure ensuring good results in the matching procedure. Tests on synthetic data demonstrated a multi-choice strategy with a 3 camera arrangement is the best solution for matching data. It is less sensitive to errors in the calibration data set when both the percentage of correctly matched triplets and the number of outputs of the algorithm were considered. Original and cross-validating methods to compute the CBL height and horizontal scale of thermals were applied based on temperature, velocity and fluorescence imaging. The spatial covariance of the velocity field, providing the plume horizontal spacing, allows the spatial extension of the mixed region to be determined. Dome characteristic vertical dimension is of the same order of magnitude as the mixing layer height, while their horizontal dimension becomes similar to the vertical one at the end of the experiment when the structure dimensions are comparable to the test section side and border effects are no longer negligible. The mixing layer growth was computed by exploiting both temperature and velocity data. Outputs were then compared to more classical methods as the zero-order mixed-layer model and the zero-heat flux level method. The accuracy of our methods was computed as well. Present results, normalized accordingly to the Deardorff similarity for free convection, were compared with literature data and LES meteorological models. Outcomes from different experimental configurations, literature and LES models are in fairly good agreement. The comparison with literature data at real scale demonstrates the validity of our experimental task and its applicability for the study of the real atmospheric boundary layer and its monitoring for environmental purposes. On the other hand the agreement with LES models at different boundary conditions and domain aspect ratios proves that both scale and border effects of the experimental model are negligible if data are not processed for too long time

On the development of the convective boundary layer in a shear-free thermally forced stably stratified fluid setting: a 2D and 3D experimental investigation using image analysis techniques coupled with temperature measurements

The motion of buoyancy driven plumes is, on all scales, the most common heat and momentum transfer mechanism in geophysical flows, well known as Free Convection. Similarly, density stratification due to heating inequalities is also an ordinary scenario in nature. Free Convection phenomenon coupled with a density stratified fluid setting leads to the so-called Penetrative Free Convection (PFC). When a fluid, in static equilibrium, is stably stratified a thermal forcing can produce an unstable configuration ensuing internal waves formation of increasing amplitude. If the perturbation is strong enough, it can definitely erode the initial stratification and cause the motion of turbulent buoyant updrafts, dome-shaped, compensated by denser downdrafts. The entrainment phenomenon occurring at the interface between the turbulent and non-turbulent region justifies the penetrative feature of convection and causes the non linear growth of the Convective Boundary Layer (CBL) of well mixed fluid against the adjacent stably stratified region. In addition to the wide engineering applications, the environmental impact mostly motivates PFC studies. The upper lakes and oceans, under calm conditions, usually exhibit a continuous, moderately stable density distribution. Turbulent convective flow can be generated both by the free-surface cooling and wind shear-stress, eroding the stable stratification on a daily or seasonal time scale. Domes with large downward velocities are generated at the free surface, balanced by updrafts with lower velocity but larger area. Because of the relatively rapid mixing, the density distribution is approximately uniform in the upper layer and it deepens with time as a result of the entrainment and erosion of the underlying denser water. An analogous phenomenon is observed in the lower troposphere when surface heating due to solar radiation results in a growing unstable layer adjacent to the ground which replaces a nocturnal inversion from below. In this case, the initially stable environment near the ground is affected by convection characterized by relatively narrow and fast plumes of rising horizontal surfaces balanced by larger regions of downward slower motion. Resulting internal waves generated within the stable layer take place at or below the Brunt-Väisälä frequency, which is related to the vertical temperature gradient. In nature, the dynamics of the CBL influences the transport and mixing features of a given stratified fluid-body. The amount of materials being mixed due to penetrative convection is a crucial issue both in water or air quality monitoring and forecast with important implications in environmentally-friendly studies. Focusing on the environmental aspect, the pollutant dispersion is a matter of particular importance. On one hand, mixing processes inside the CBL help increasing dispersion with some positive consequences: the concentration of potential harmful pollutants in high risk zones tends to decrease, the turnover and the redistribution of vital substances, like oxygen and nutrients, is guaranteed; the latter plays a major role in large water bodies. These advantages are essential for the safety of populations living close to urban or industrial areas, or for preserving coastal human activities and ecosystems. On the other hand each pollutant, released inside the CBL mainly by human activities, remains confined inside it because of the interface with the non-turbulent region, which acts as a barrier for outward transport. Considering the proximity of the CBL with biosphere, a deeper insight into dispersion and entrainment processes appears mandatory either for sustainable engineering design or for monitoring purposes. Given the great applicability of the topic in several scientific and engineering fields, large amount of experimental, theoretical and numerical investigations on CBL development in a continuously and linearly stratified fluid setting had been conducted in the past since sixties. It appears there is a lack of consensus about the dependence of integral parameters of convective entrainment (in particular the CBL growth rate) on the initial stratification strength and convective phenomenon evolution. Driven by all these premises the main aims of the present research are related to better understanding dispersion of a passive scalar inside the CBL with a lagrangian, non-local approach and a fully three-dimensional (3D 3C) experimental technique, which for the first time has been applied to the topic. When turbulent convection occurs, in fact, dispersion is mostly due to transport by large organized structures while molecular diffusion can be neglected. Given this assumption, a non local approach based on a fully advective-like behaviour of the tracer is necessary and a lagrangian description of flow may be more suitable. Moreover, the knowledge of the horizontal and vertical extension of the structures dominating the flow field appears to be mandatory. Furthermore, turbulence is fully three-dimensional on the scales of motion characterising the phenomena in nature. In order to better understanding and likely describing the evolution of turbulent structures inside the convective layer, a three dimensional experimental technique is strongly required. In the present work shear free convection experiments in a stably stratified environment were performed in a thermally controlled convection chamber. The experimental set up was arranged in different optical configurations to ensure the two components and three components of velocity to be measured in an illuminated plane and volume respectively through different image analysis techniques with sub-pixel accuracy. Temperature measurements were simultaneously acquired using thermocouples of accuracy less than 0.1 °C. The experiments highlighted the time evolution of the convective structure characteristic spatial scales (CBL height and horizontal spacing between thermals) and the non local description of transport and mixing inside the growing CBL. The main novelty of the present contribution covers the improvement of techniques and methods to achieve more reliable, statistically robust and likely results. The experimental effort based on image analysis techniques (mainly Feature Tracking and photogrammetric 3D-PTV) resulted suitable for reconstructing longer trajectories (always more than 103 trajectories longer than 20 consecutive snapshots) and analyzing larger particle density images (reliable results for more than 2000 particles in a volume of 15X15X15 cm3 for 3D-PTV). Velocity statistics are then more robust than those from classical Particle Tracking Velocimetry. Moreover, for the first time, a fully three-dimensional particle tracking technique has been applied to penetrative convection experiments. 3D-PTV allows a more realistic description of the velocity field, which occurs during the evolution of the convective mixed layer, than more traditional 2D techniques. Furthermore, photogrammetric 3D-PTV rather than “scanning” 3D-PTV results in more accuracy when the tracer particle density is high, because particles may be tracked directly in the 3D space rather than through matching of 2D projections. The broader impact of the research mainly refers to the prediction of the CBL growth as a function of initial and boundary conditions with better accuracy than conventional and well established techniques. The experimental study can thus give a positive contribution on real pollutant dispersion studies in urban and natural environments for environmental protection and sustainable design purposes. Field experiments aimed at measuring the turbulence budget of the CBL have shown that the mechanical generation of kinetic energy by wind shear is often confined close to the heat source supporting the validity of laboratory models in which no wind is present. According to this assumption, the similarity proposed by Deardorff (1970) is employed to compute scaling parameters and to make results comparable with real scales. Through normalizing the quantities measured at different stages of the experiment, the phenomenon can be considered as a succession of steady states, according to an evolution of the variables of interest that may be defined quasi-steady state. The experimental apparatus employed to run the experiments is the same as in Cenedese and Querzoli (1994), Querzoli (1996), Cenedese and Querzoli (1997) and Moroni and Cenedese (2006). The spatial resolution of velocity data is largely increased here by means of 2D and 3D image analysis techniques (Feature Tracking, FT, and 3D Particle Tracking Velocimetry, 3D-PTV) used instead of Laser-Doppler Anemometry or 2D Particle Tracking Velocimetry as in Cenedese and Querzoli, 1994; Querzoli, 1996 and Cenedese and Querzoli, 1997. Moreover the photogrammetric 3D-PTV here applied allows fully three-dimensional descriptions of both the Eulerian velocity field and Lagrangian particle trajectories for a more likely understanding of the phenomenon than through the 2D approach used in Moroni and Cenedese (2006). Furthermore, the combined use of thermocouples and flow visualization techniques allows cross-validating different methods to estimate the evolutions of the key parameters and the plume characteristic dimensions. Two different experimental arrangements were set-up, 2D and a 3D models. A large set of data were firstly acquired using a 2D model and employing FT. The preliminary investigation was mainly focused on better understanding the physics of the phenomenon, finding a time scaling law, testing different methods to compute the variables of interest and comparing results with classical methods found in literature. When 2D techniques are employed to detect the velocity field, the flow is illuminated with a thin light sheet and only the velocity components within this sheet can be evaluated. Driven by the idea that only a fully three-dimensional technique can significantly improve our laboratory model in term of a more likely description of free convection structures we run a second set of experiments, by using a stereoscopic arrangement of cameras focused on an illuminated volume. Although some methods do exist for reconstructing 3D velocities in a point (3D laser Doppler Anemometry; Hinsch and Hinrichs, 1996) or plane (3D stereo-PIV; Stuer et al., 1999), only a fully 3D technique based on the illumination of a flow volume rather than a flow sheet will give the information needed to construct the instantaneous 3D velocity fields. A number of imaging-based measurement techniques exist for determining 3D velocity fields in an observation volume. Among these are: scanning, photogrammetric, holographic or photogrammetric techniques dependending on which principle is recalled to reconstruct the third dimension from a 2D image/s: The present study was focused on 3D-PTV which is a 3D extension of the 2D particle-tracking methods. 3D-PTV is based on reconstructing 3D trajectories of reflecting tracer particles through a photogrammetric recording of image sequences. The 3D particle trajectories obtained can be used to calculate the 3D velocity field. The 3D-PTV optical system has been designed with the following capabilities: image a volume far away the boundary walls, lengthen the trajectories, and improve the accuracy of the procedure through a careful test on synthetically generated data. A physically-based photogrammetric calibration of the stereoscopic arrangement was employed and its accuracy tested. The effects of multimedia geometry on calibration parameters were taken into account. The combination of image- and object-space based information was employed to establish the correspondences between particle positions (structure from stereo reconstruction). A particle tracking algorithm was then employed to reconstructed 3D trajectories. Sensitivity tests conducted on the matching algorithm proved that the calibration accuracy is fundamental to obtain the correct matching and particle tracking; small errors in calibration parameters or neglecting water refraction effects reduce matching performance. On the other hand, accuracy less than 1 pixel was reached with our calibration procedure ensuring good results in the matching procedure. Tests on synthetic data demonstrated a multi-choice strategy with a 3 camera arrangement is the best solution for matching data. It is less sensitive to errors in the calibration data set when both the percentage of correctly matched triplets and the number of outputs of the algorithm were considered. Original and cross-validating methods to compute the CBL height and horizontal scale of thermals were applied based on temperature, velocity and fluorescence imaging. The spatial covariance of the velocity field, providing the plume horizontal spacing, allows the spatial extension of the mixed region to be determined. Dome characteristic vertical dimension is of the same order of magnitude as the mixing layer height, while their horizontal dimension becomes similar to the vertical one at the end of the experiment when the structure dimensions are comparable to the test section side and border effects are no longer negligible. The mixing layer growth was computed by exploiting both temperature and velocity data. Outputs were then compared to more classical methods as the zero-order mixed-layer model and the zero-heat flux level method. The accuracy of our methods was computed as well. Present results, normalized accordingly to the Deardorff similarity for free convection, were compared with literature data and LES meteorological models. Outcomes from different experimental configurations, literature and LES models are in fairly good agreement. The comparison with literature data at real scale demonstrates the validity of our experimental task and its applicability for the study of the real atmospheric boundary layer and its monitoring for environmental purposes. On the other hand the agreement with LES models at different boundary conditions and domain aspect ratios proves that both scale and border effects of the experimental model are negligible if data are not processed for too long time

Beyond Mixing-length Theory: a step toward 321D

We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different P\'eclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated.Comment: 22 pages, 4 figures, 2 tables; significantly re-written, critique of Pasetto, et al. model added, accepted for publication by Ap

Stability of parallel, quasi-parallel and stationary flows, The

CER73-74-RK-RNM12.Prepared under Office of Naval Research, project no. NR 062-414/6-6-68 (Code 438), U.S. Department of Defense.Includes bibliographical references.The methods of linear perturbation theory have been used to study the stability of various flows, among them being (i) The stability of boundary layers along concave heated walls; (ii) The stability of boundary layers along concave walls with suction; (iii) The stability of wall jets along concave and convex walls; (iv) The spin up of a two-dimensional cylinder in an infinite medium; (v) The stability of stationary layers of fluid with arbitrary temperature stratification; (vi) The stability of natural convection flow along inclined plates. During the course of this work, three different solution techniques were employed; one of them was an approximate analytic technique, the remaining two were numerical. Three-dimensional spatially and temporally amplifying disturbances were considered in this study. The results indicated that the normal velocity component of the mean flow in a boundary layer, although much smaller than the stream wise component had a profound effect in reducing the stability of the flow. On the other hand, suction at the wall improved the stability characteristics. For the flow of parallel layers of fluid along heated walls with small curvature, it was found that a unique stability curve for neutral disturbances may be obtained if the quantity plotted along the abscissa is Ra + KsNg2 where Ra is the Rayleigh number, Ng is the Goertler number and Ks is a constant which expresses the relative importance of the mean temperature and velocity profiles. It was demonstrated also that wall jets are unstable on concave as well as convex walls. The results obtained for the stability of the spin up of a cylinder in an infinite medium are in qualitative agreement with experiment. The dependence of the onset of convective overturning in an unstable layer of fluid with a nonlinear basic temperature profile and bounded above by fluid of varying stability on Rayleigh number was established. The angle at which the two-dimensional wave instability passes into the three-dimensional mode in natural convection along an inclined plate was calculated. The result was found to be in good agreement with experiment. Other results obtained for this flow were in good qualitative agreement with experiment. Finally, some simple wind tunnel experiments with boundary layers along curved heated walls were performed. Photographic evidence of longitudinal vortices was obtained together with some qualitative data.Under contract no. N00014-68-A-0493-0001

Natural Convection Above Heated Inclined Surfaces.

Motivated by a need to assess pollutant transport by upslope flows, an investigation has been conducted into the fundamentals of natural convection flow over inclined surfaces. Particular attention was focused on the influence of ambient fluid stability. Field studies were performed using tracer gas releases into the upslope flow over a Southern California mountain range. The field studies served to reveal the presence of a split slope flow recirculation and demonstrated the impact of this recirculation on the transport of pollutants from a valley. In order to pursue a controlled investigation of the phenomena found in the field work, a laboratory model was developed using water as a working fluid. Extensive dye studies demonstrated the presence of this recirculation to varying degrees in nearly every configuration with a stable layer present. Heat transfer experiments were conducted with the laboratory model to refine and validate the experimental techniques used. Comparisons are made with existing theoretical and empirical predictions where available. Existing correlations for inclined surfaces are extended two orders of magnitude lower in Rayleigh number. A modification to vertical theory based on simply replacing g with gcos$\theta$ is shown to be useful for inclinations down to 75\sp\circ from vertical. Transition ranges and empirical correlations are expressed for inclinations of 0\sp\circ, 15\sp\circ, 30\sp\circ, 45\sp\circ, 60\sp\circ, 75\sp\circ, and 90\sp\circ. Overall correlations are also reported with apparently far less scatter of data than for any previously reported research with inclined surfaces. Experimentation with stratified ambient fluid resulted in the observation that, with turbulent flow, stratification could apparently be disregarded and heat transfer simply calculated from local conditions. The observed heat transfer coefficients are essentially independent of position along the slope suggesting that an approximate analytical model of upslope flows developed by L. Prandtl in 1942 is applicable. However, quantitative laboratory results showed that Prandtl\u27s one dimensional theory underpredicted the observed boundary layer depths. The results are in good agreement with observed characteristics of atmospheric slope flows

The stability of stratified layers within ventilated enclosures

The project consists of experimental and numerical investigation of buoyancy-dominated flow leading to thermal stratification in ventilated spaces for the range of Re from 4,000 to 50,000 and Rio below 40. In the evolution of turbulence in stratified shear flows, solving for the primitive variables both in time and space is most important. Steady-state Reynolds Averaged Navier Stokes (RANS) simulations have shown in the past to model adequately the flows that become eventually steady-state. Therefore, three-dimensional steady-state CFD models were used to simulate the flow. The experimental set-up used an Environmental Test Chamber to investigate the stratification in buildings. Temperature differential was created by introducing cold air through a terminal at floor level and hot air through a diffuser pointing vertically downwards at the ceiling level. An extract point was located opposite to the inlets at the rear wall. Different combinations of inlet velocities/flow rates, temperatures as well as different arrangements of the exhaust height were studied to evaluate its effect on velocity and temperature distribution and the effect of jet flow interaction with the stratified layers. The experimental data obtained in this work were also used to validate the predictions from the numerical simulations so that CFD was used to perform more parametric simulations. To model turbulence in the flow field, the appropriateness of several eddy-viscosity based turbulence models was evaluated. The turbulence models used were the standard k-c model, the k-e Chen-Kim modified model, a Low- Re modification and the RNG k -c model. All models are in close agreement with each other. The buoyancy extended standard k-e model gives a sharper prediction of mixing in the interface. The physics of the flow are well predicted. The hot jet from the ceiling is buoyed up when it reaches the height of the interface where buoyancy forces dominate gravitational forces marking its appearance on the temperature gradient. From the simulations made at this stage, the characteristics of stratified layers are revealed. It is also shown that the interface height is proportional to the exhaust height for certain room temperatures and velocities. In a comparison between experimental results and numerical simulation, the differences are attributed on the non-adiabatic effects such as heat losses and radiation. A correlation is found between inlet parameters and stratified flow in buildings obtained both experimentally and computationally

Summary of Research 1998, Department of Oceangraphy

The views expressed in this report are those of the authors and do not reflect the official policy or position of the Department of Defense or the U.S. Government.This report contains summaries of research projects in the Department of Oceanography. A list of recent publications is also included which consists of conference presentations and publications, books, contributions to books, published journal papers, technical reports, and thesis abstracts

Experimental evidence for the boundary zonal flow in rotating Rayleigh-Benard convection

We report on the presence of the boundary zonal flow in rotating Rayleigh-Benard convection evidenced by two-dimensional particle image velocimetry. Experiments were conducted in a cylindrical cell of aspect ratio between its diameter (D) and height (H) of G=D/H=1. As the working fluid we used various mixtures of water and glycerol, leading to Prandtl numbers in the range 6.6<Pr<75. The horizontal velocity components were measured at a horizontal cross-section at half height. The Rayleigh numbers were in the range 10e8 < Ra < 6.3e9. The effect of rotation is quantified by the Ekman number which was between 1.5e-5<Ek< 1.2e-3 in our experiment. With our results we show the first direct measurements of the boundary zonal flow - BZF that develops near the sidewall and was recently discovered in numerical simulations as well as in sparse and localized temperature measurements. We analyse the thickness of the BZF as well as its maximal velocity as function of Pr, Ra, and Ek, and compare these results with previous DNS results