Shock Train/Boundary-Layer Interaction in Rectangular Scramjet Isolators

Numerous studies of the dual-mode scramjet isolator, a critical component in preventing inlet unstart and/or vehicle loss by containing a collection of flow disturbances called a shock train, have been performed since the dual-mode propulsion cycle was introduced in the 1960s. Low momentum corner flow and other three-dimensional effects inherent to rectangular isolators have, however, been largely ignored in experimental studies of the boundary layer separation driven isolator shock train dynamics. Furthermore, the use of two dimensional diagnostic techniques in past works, be it single-perspective line-of-sight schlieren/shadowgraphy or single axis wall pressure measurements, have been unable to resolve the three-dimensional flow features inside the rectangular isolator. These flow characteristics need to be thoroughly understood if robust dual-mode scramjet designs are to be fielded. The work presented in this thesis is focused on experimentally analyzing shock train/boundary layer interactions from multiple perspectives in aspect ratio 1.0, 3.0, and 6.0 rectangular isolators with inflow Mach numbers ranging from 2.4 to 2.7. Secondary steady-state Computational Fluid Dynamics studies are performed to compare to the experimental results and to provide additional perspectives of the flow field. Specific issues that remain unresolved after decades of isolator shock train studies that are addressed in this work include the three-dimensional formation of the isolator shock train front, the spatial and temporal low momentum corner flow separation scales, the transient behavior of shock train/boundary layer interaction at specific coordinates along the isolator's lateral axis, and effects of the rectangular geometry on semi-empirical relations for shock train length prediction. A novel multiplane shadowgraph technique is developed to resolve the structure of the shock train along both the minor and major duct axis simultaneously. It is shown that the shock train front is of a hybrid oblique/normal nature. Initial low momentum corner flow separation spawns the formation of oblique shock planes which interact and proceed toward the center flow region, becoming more normal in the process. The hybrid structure becomes more two-dimensional as aspect ratio is increased but corner flow separation precedes center flow separation on the order of 1 duct height for all aspect ratios considered. Additional instantaneous oil flow surface visualization shows the symmetry of the three-dimensional shock train front around the lower wall centerline. Quantitative synthetic schlieren visualization shows the density gradient magnitude approximately double between the corner oblique and center flow normal structures. Fast response pressure measurements acquired near the corner region of the duct show preliminary separation in the outer regions preceding centerline separation on the order of 2 seconds. Non-intrusive Focusing Schlieren Deflectometry Velocimeter measurements reveal that both shock train oscillation frequency and velocity component decrease as measurements are taken away from centerline and towards the side-wall region, along with confirming the more two dimensional shock train front approximation for higher aspect ratios. An updated modification to Waltrup \& Billig's original semi-empirical shock train length relation for circular ducts based on centerline pressure measurements is introduced to account for rectangular isolator aspect ratio, upstream corner separation length scale, and major- and minor-axis boundary layer momentum thickness asymmetry. The latter is derived both experimentally and computationally and it is shown that the major-axis (side-wall) boundary layer has lower momentum thickness compared to the minor-axis (nozzle bounded) boundary layer, making it more separable. Furthermore, it is shown that the updated correlation drastically improves shock train length prediction capabilities in higher aspect ratio isolators. This thesis suggests that performance analysis of rectangular confined supersonic flow fields can no longer be based on observations and measurements obtained along a single axis alone. Knowledge gained by the work performed in this study will allow for the development of more robust shock train leading edge detection techniques and isolator designs which can greatly mitigate the risk of inlet unstart and/or vehicle loss in flight.

Air Entrainment in the Turbulent Ship Hull Boundary Layer

Turbulent fluctuations in the vicinity of the water free surface along a flat, vertically oriented surface-piercing plate are studied experimentally using a laboratory-scale experiment. In this experiment, a meter-wide stainless steel belt travels horizontally in a loop around two rollers with vertically oriented axes, which are separated by 7.5 meters. This belt device is mounted inside a large water tank with the water level set just below the top edge of the belt. The belt, rollers, and supporting frame are contained within a sheet metal box to keep the device dry except for one 6-meter-long straight test section between rollers. The belt is launched from rest with an acceleration of up to 3-g in order to quickly reach steady state velocity. This creates a temporally evolving boundary layer analogous to the spatially evolving boundary layer created along a flat-sided ship moving at the same velocity, with a length equivalent to the length of belt that has passed the measurement region since the belt motion began. Surface profile measurements in planes normal to the belt surface are conducted using cinematic Laser Induced Fluorescence and quantitative surface profiles are extracted at each instant in time. Using these measurements, free surface fluctuations are examined and the propagation behavior of these free surface ripples is studied. It is found that free surface fluctuations are generated in a region close to the belt surface, where sub-surface velocity fluctuations influence the behavior of these free surface features. These rapidly-changing surface features close to the belt appear to lead to the generation of freely-propagating waves far from the belt, outside the influence of the boundary layer. Sub-surface PIV measurements are performed in order to study the modification of the boundary layer flow field due to the effects of the water free surface. Cinematic planar PIV measurements are performed in horizontal planes parallel to the free surface by imaging the flow from underneath the tank, providing streamwise and wall-normal velocity fields. Additional planar PIV experiments are performed in vertical planes parallel to the belt surface in order to study the bahvior of streamwise and vertical velocity fields. It is found that the boundary layer grows rapidly near the free surface, leading to an overall thicker boundary layer close to the surface. This rapid boundary layer growth appears to be linked to a process of free surface bursting, the sudden onset of free surface fluctuations. Cinematic white light movies are recorded from beneath the water surface in order to determine the onset location of air entrainment. In addition, qualitative observations of these processes are made in order to determine the mechanisms leading to air entrainment present in this flow.

Scalar turbulent behavior in the roughness sublayer of an Amazonian forest.

An important current problem in micrometeorology is the characterization of turbulence in the roughness sublayer (RSL), where most of the measurements above tall forests are made. There, scalar turbulent fluctuations display significant departures from the predictions of Monin?Obukhov similarity theory (MOST). In this work, we analyze turbulence data of virtual temperature, carbon dioxide, and water vapor in the RSL above an Amazonian forest (with a canopy height of 40 m), measured at 39.4 and 81.6 m above the ground under unstable conditions. We found that dimensionless statistics related to the rate of dissipation of turbulence kinetic energy (TKE) and the scalar variance display significant departures from MOST as expected, whereas the vertical velocity variance follows MOST much more closely. Much better agreement between the dimensionless statistics with the Obukhov similarity variable, however, was found for the subset of measurements made at a low zenith angle Z, in the range 0°  <  |Z|  <  20°. We conjecture that this improvement is due to the relationship between sunlight incidence and the ?activation?deactivation? of scalar sinks and sources vertically distributed in the forest. Finally, we evaluated the relaxation coefficient of relaxed eddy accumulation: it is also affected by zenith angle, with considerable improvement in the range 0°  <  |Z|  <  20°, and its values fall within the range reported in the literature for the unstable surface layer. In general, our results indicate the possibility of better stability-derived flux estimates for low zenith angle ranges.

Scalar turbulent behavior in the roughness sublayer of an Amazonian forest.

An important current problem in micrometeorology is the characterization of turbulence in the roughness sublayer (RSL), where most of the measurements above tall forests are made. There, scalar turbulent fluctuations display significant departures from the predictions of Monin?Obukhov similarity theory (MOST). In this work, we analyze turbulence data of virtual temperature, carbon dioxide, and water vapor in the RSL above an Amazonian forest (with a canopy height of 40?m), measured at 39.4 and 81.6?m above the ground under unstable conditions. We found that dimensionless statistics related to the rate of dissipation of turbulence kinetic energy (TKE) and the scalar variance display significant departures from MOST as expected, whereas the vertical velocity variance follows MOST much more closely. Much better agreement between the dimensionless statistics with the Obukhov similarity variable, however, was found for the subset of measurements made at a low zenith angle Z, in the range 0°??stable surface layer. In general, our results indicate the possibility of better stability-derived flux estimates for low zenith angle ranges.

A process-based evaluation of dust-emitting winds in the CMIP5 simulation of HadGEM2-ES

Despite the importance of dust aerosol in the Earth system, state-of-the-art models show a large variety for North African dust emission. This study presents a systematic evaluation of dust emitting-winds in 30 years of the historical model simulation with the UK Met Office Earth-system model HadGEM2-ES for the Coupled Model Intercomparison Project Phase 5. Isolating the effect of winds on dust emission and using an automated detection for nocturnal low-level jets (NLLJs) allow an in-depth evaluation of the model performance for dust emission from a meteorological perspective. The findings highlight that NLLJs are a key driver for dust emission in HadGEM2-ES in terms of occurrence frequency and strength. The annually and spatially averaged occurrence frequency of NLLJs is similar in HadGEM2-ES and ERA-Interim from the European Centre for Medium-Range Weather Forecasts. Compared to ERA-Interim, a stronger pressure ridge over northern Africa in winter and the southward displaced heat low in summer result in differences in location and strength of NLLJs. Particularly the larger geostrophic winds associated with the stronger ridge have a strengthening effect on NLLJs over parts of West Africa in winter. Stronger NLLJs in summer may rather result from an artificially increased mixing coefficient under stable stratification that is weaker in HadGEM2-ES. NLLJs in the Bodélé Depression are affected by stronger synoptic-scale pressure gradients in HadGEM2-ES. Wintertime geostrophic winds can even be so strong that the associated vertical wind shear prevents the formation of NLLJs. These results call for further model improvements in the synoptic-scale dynamics and the physical parametrization of the nocturnal stable boundary layer to better represent dust-emitting processes in the atmospheric model. The new approach could be used for identifying systematic behavior in other models with respect to meteorological processes for dust emission. This would help to improve dust emission simulations and contribute to decreasing the currently large uncertainty in climate change projections with respect to dust aerosol.

Dynamics of Bora wind over the Adriatic sea: atmospheric water balance and role of air-sea fluxes and orography

The Bora wind is a mesoscale phenomenon which typically affects the Adriatic Sea basin for several days each year, especially during winter. The Bora wind has been studied for its intense outbreak across the Dinaric Alps. The properties of the Bora wind are widely discussed in the literature and scientific papers usually focus on the eastern Adriatic coast where strong turbulence and severe gust intensity are more pronounced. However, the impact of the Bora wind can be significant also over Italy, not only in terms of wind speed instensity. Depending on the synoptic pressure pattern (cyclonic or anticyclonic Bora) and on the season, heavy snowfall, severe storms, storm surges and floods can occur along the Adriatic coast and on the windward flanks of the Apennines. In the present work five Bora cases that occurred in recent years have been selected and their evolution has been simulated with the BOLAM-MOLOCH model set, developed at ISAC-CNR in Bologna. Each case study has been addressed by a control run and by several sensitivity tests, performed with the purpose of better understanding the role played by air-sea latent and sensible heat fluxes. The tests show that the removal of the fluxes induces modifications in the wind approching the coast and a decrease of the total precipitation amount predicted over Italy. In order to assess the role of heat fluxes, further analysis has been carried out: column integrated water vapour fluxes have been computed along the Italian coastline and an atmospheric water balance has been evaluated inside a box volume over the Adriatic Sea. The balance computation shows that, although latent heat flux produces a significant impact on the precipitation field, its contribution to the balance is relatively minor. The most significant and lasting case study, that of February 2012, has been studied in more detail in order to explain the impressive drop in the total precipitation amount simulated in the sensitivity tests with removed heat fluxes with respect to the CNTRL run. In these experiments relative humidity and potential temperature distribution over different cross-sections have been examined. With respect to the CNTRL run a drier and more stable boundary layer, characterised by a more pronounced wind shear at the lower levels, has been observed to establish above the Adriatic Sea. Finally, in order to demonstrate that also the interaction of the Bora flow with the Apennines plays a crucial role, sensitivity tests varying the orography height have been considered. The results of such sensitivity tests indicate that the propagation of the Bora wind over the Adriatic Sea, and in turn its meteorological impact over Italy, is influenced by both the large air-sea heat fluxes and the interaction with the Apennines that decelerate the upstream flow.

Use of Land Surface Temperature to Estimate Surface Energy Fluxes: Contributions of Wilfried Brutsaert and Collaborators

Land surface temperature (LST) plays a key role in governing the land surface energy budget, and measurements or estimates of LST are an integral part of many land surface models and methods to estimate land surface sensible heat (H) and latent heat fluxes. In particular, the LST anchors the potential temperature profile in Monin-Obukhov similarity theory, from which H can be derived. Brutsaert has made important contributions to our understanding the nature of surface temperature measurements as well as the practical but theoretically sound use of LST in this framework. His work has coincided with the wide-spread availability of remotely sensed LST measurements. Use of remotely sensed LST estimates inevitably involves complicating factors, such as: varying spatial and temporal scales in measurements, theory, and models; spatial variability of LST and H; the relationship between measurements of LST and the temperature felt by the atmosphere; and the need to correct satellite-based radiometric LST measurements for the radiative effects of the atmosphere. This paper reviews the progress made in research in these areas by tracing and commenting on Brutsaert's contributions.

Table 4. Trace gas concentrations from samples taken at Ginninderra

We assess the performance of an inverse Lagrangian dispersion technique for its suitability to quantify leakages from geological storage of CO2. We find the technique is accurate ((QbLS/Q)=0.99, sigma=0.29) when strict meteorological filtering is applied to ensure that Monin-Obukhov Similarity Theory is valid for the periods analysed and when downwind enrichments in tracer gas concentration are 1% or more above background concentration. Because of their respective baseline atmospheric concentrations, this enrichment criterion is less onerous for CH4 than for CO2. Therefore for geologically sequestered gas reservoirs with a significant CH4 component, monitoring CH4 as a surrogate for CO2 leakage could be as much as 10 times more sensitive than monitoring CO2 alone. Additional recommendations for designing a robust atmospheric monitoring strategy for geosequestration include: continuous concentration data; exact inter-calibration of up- and downwind concentration measurements; use of an array of point concentration sensors to maximise the use of spatial information about the leakage plume; and precise isotope ratio measurement to confirm the source of any concentration elevations detected.

JOYCE: Jülich Observatory for cloud evolution

The Jülich Observatory for Cloud Evolution (JOYCE), located at Forschungszentrum Jülich in the most western part of Germany, is a recently established platform for cloud research. The main objective of JOYCE is to provide observations, which improve our understanding of the cloudy boundary layer in a midlatitude environment. Continuous and temporally highly resolved measurements that are specifically suited to characterize the diurnal cycle of water vapor, stability, and turbulence in the lower troposphere are performed with a special focus on atmosphere–surface interaction. In addition, instruments are set up to measure the micro- and macrophysical properties of clouds in detail and how they interact with different boundary layer processes and the large-scale synoptic situation. For this, JOYCE is equipped with an array of state-of-the-art active and passive remote sensing and in situ instruments, which are briefly described in this scientific overview. As an example, a 24-h time series of the evolution of a typical cumulus cloud-topped boundary layer is analyzed with respect to stability, turbulence, and cloud properties. Additionally, we present longer-term statistics, which can be used to elucidate the diurnal cycle of water vapor, drizzle formation through autoconversion, and warm versus cold rain precipitation formation. Both case studies and long-term observations are important for improving the representation of clouds in climate and numerical weather prediction models.

SCG postnatal remodelling - hypertrophy and neuron number stability - in Spix's Yellow-toothed Cavies (Galea spixii)

Whilst a fall in neuron numbers seems a common pattern during postnatal development, several authors have nonetheless reported an increase in neuron number, which may be associated with any one of a number of possible processes encapsulating either neurogenesis or late maturation and incomplete differentiation. Recent publications have thus added further fuel to the notion that a postnatal neurogenesis may indeed exist in sympathetic ganglia. In the light of these uncertainties surrounding the effects exerted by postnatal development on the number of superior cervical ganglion (SCG) neurons, we have used state-of-the-art design-based stereology to investigate the quantitative structure of SCG at four distinct timepoints after birth, viz., 1-3 days, 1 month, 12 months and 36 months. The main effects exerted by ageing on the SCG structure were: (i) a 77% increase in ganglion volume; (ii) stability in the total number of the whole population of SCG nerve cells (no change - either increase or decrease) during post-natal development; (iii) a higher proportion of uninucleate neurons to binucleate neurons only in newborn animals; (iv) a 130% increase in the volume of uninucleate cell bodies; and (v) the presence of BrdU positive neurons in animals at all ages. At the time of writing our results support the idea that neurogenesis takes place in the SCG of preas, albeit it warrants confirmation by further markers. We also hypothesise that a portfolio of other mechanisms: cell repair, maturation, differentiation and death may be equally intertwined and implicated in the numerical stability of SCG neurons during postnatal development. (C) 2011 ISDN. Published by Elsevier Ltd. All rights reserved.

A comparison of dry and wet season aerosol number fluxes over the Amazon rain forest

Vertical number fluxes of aerosol particles and vertical fluxes of CO(2) were measured with the eddy covariance method at the top of a 53 m high tower in the Amazon rain forest as part of the LBA (The Large Scale Biosphere Atmosphere Experiment in Amazonia) experiment. The observed aerosol number fluxes included particles with sizes down to 10 nm in diameter. The measurements were carried out during the wet and dry season in 2008. In this study focus is on the dry season aerosol fluxes, with significant influence from biomass burning, and these are compared with aerosol fluxes measured during the wet season. Net particle deposition fluxes dominated in daytime in both seasons and the deposition flux was considerably larger in the dry season due to the much higher dry season particle concentration. The particle transfer velocity increased linearly with increasing friction velocity in both seasons. The difference in transfer velocity between the two seasons was small, indicating that the seasonal change in aerosol number size distribution is not enough for causing any significant change in deposition velocity. In general, particle transfer velocities in this study are low compared to studies over boreal forests. The reasons are probably the high percentage of accumulation mode particles and the low percentage of nucleation mode particles in the Amazon boundary layer, both in the dry and wet season, and low wind speeds in the tropics compared to the midlatitudes. In the dry season, nocturnal particle fluxes behaved very similar to the nocturnal CO(2) fluxes. Throughout the night, the measured particle flux at the top of the tower was close to zero, but early in the morning there was an upward particle flux peak that is not likely a result of entrainment or local pollution. It is possible that these morning upward particle fluxes are associated with emission of primary biogenic particles from the rain forest. Emitted particles may be stored within the canopy during stable conditions at nighttime, similarly to CO(2), and being released from the canopy when conditions become more turbulent in the morning.