Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

Influence of wave dispersion on vortex pairing in a jet

The geometry of large-scale structures in the turbulent mixing layer of a moderate Reynolds number jet is deduced from measurements of the fluctuating pressure in the hydrodynamic near field. The structures are rings of concentrated vorticity that distort with downstream distance until statistical axisymmetry disappears. The rings are spaced quasi-periodically and coalesce with each other, producing larger spacings. Statistical and flow-visualization techniques are applied to free and forced jets over a range of Reynolds numbers from 5000 to 50,000 to demonstrate that rings of a given spacing do not coalesce with each other until they are far enough downstream that the local mixing layer has attained some critical thickness which scales with the wavelength of the vortex pair. Wave dispersion is evaluated as a plausible mechanism for localizing the coalescences. The central feature of the model is the observation that a shear layer is dispersive to wavelengths much longer than its thickness and nondispersive to shorter waves.

Petersen, R. A.↗

Vortices in turbulence

Theoretical and experimental results of studies to determine the flow parameters and structures of plane mixing layers are reported. Smoke visualization, combined with hot-wire anemometry, sheets of laser light, and photography were used to gather data from the wall flow. The behavior of vortex rods was examined, noting that the rods persisted only if new vortex energy was supplied from the sublayer. A power spectral density was defined for the velocity fluctuations, as was a hierarchy of velocity scales for geometrically similar vortices. The length scale grows linearly with downstream distance, where the flow structures are fed by longitudinal vortices. A model is developed for vortex pairing in sequential order from the bottom of the mixing layer outward in a repetitive process involving vortex stretching. The model is actually a migration strategy that satisfies the flow self-preservation constraints.

Perry, A. E.↗

Laminar mixing of a compressible fluid

A theoretical investigation of the velocity profiles for laminar mixing of a high-velocity stream with a region of fluid at rest has been made assuming that the Prandtl number is unity. A method which involves only quadratures is presented for calculating the velocity profile in the mixing layer for an arbitrary value of the free-stream Mach number. Detailed velocity profiles have been calculated for free-stream Mach numbers of 0, 1, 2, 3, and 5. For each Mach number, velocity profiles are presented for both a linear and a 0.76-power variation of viscosity with absolute temperature. The calculations for a linear variation are much simpler than those for a 0.76-power variation. It is shown that by selecting the constant of proportionality in the liner approximation such that it gives the correct value for the viscosity in the high-temperature part of the mixing layer, the resulting velocity profiles are in excellent agreement with those calculated by a 0.76-power variation.

Chapman, Dean R↗

Biogenic nonmethane hydrocarbon emissions estimated from tethered balloon observations

A new technique for estimating surface fluxes of trace gases, the mixed-layer gradient technique, is used to calculate isoprene and terpene emissions from forests. The technique is applied to tethered balloon measurements made over the Amazon forest and a pine-oak forest in Alabama at altitudes up to 300 m. The observations were made during the dry season Amazon Boundary Layer Experiment (ABLE 2A) and the Rural Oxidants in the Southern Environment 1990 experiment (ROSE I). Results from large eddy simulations of scalar transport in the clear convective boundary layer are used to infer fluxes from the balloon profiles. Profiles from the Amazon give a mean daytime emission of 3630 +/- 1400 micrograms isoprene sq m/h, where the uncertainty represents the standard deviation of the mean of eight flux estimates. Twenty profiles from Alabama give emissions of 4470 +/- 3300 micrograms isoprene sq m/h, 1740 +/- 1060 micrograms alpha-pinene sq m/h, and 790 +/- 560 micrograms beta-pinene sq m/h, respectively. These results are in agreement with emissions derived from chemical budgets. The emissions may be overestimated because of uncertainty about how to incorporate the effects of the canopy on the mixed-layer gradients. The large variability in these emission estimates is probably due to the relatively short sampling times of the balloon profiles, though spatially heterogeneous emissions may also play a role. Fluxes derived using this technique are representative of an upwind footprint of several kilometers and are independent of hydrocarbon oxidation rate and mean advection.

Davis, K. J.↗

Airborne High Spectral Resolution Lidar Measurements of Aerosol Distributions and Properties During the NASA CAMP2Ex Mission

The NASA Clouds, Aerosols and Monsoon Processes-Philippines Experiment (CAMP2Ex) occurred during August-October, 2019,primarily over the ocean near the Philippines, including the ocean west of Luzon. CAMP2Ex attempts to characterize the role of anthropogenic and natural aerosol particles in modulating warm and mixed-phase precipitation in the maritime continent during the Southwest Monsoon. This mission is also examining relationships between spatial heterogeneity in the aerosol and cloud fields and radiative heating rates and studying the effects of aerosol concentration and composition on the microphysical properties of shallow cumulus and cumulus congestus clouds. Critical to the successful investigation of these objectives is profiling of aerosol, cloud and meteorological properties throughout the atmospheric column. For the CAMP2Ex experiment, we deployed the second generation NASA Langley airborne High Spectral Resolution Lidar (HSRL-2) on the NASA P-3 aircraft to obtain profile measurements of aerosol and thin cloud extinction and optical thickness via the HSRL technique at 355 and 532 nm, and of backscatter and depolarization at 355, 532, and 1064 nm. The HSRL technique allows for the independent measurement of extinction and backscatter without the need to assume modeled lidar ratios. Additional HSRL-2 data products include aerosol type, mixed layer depth, and range-resolved aerosol microphysical parameters. The HSRL measurements reveal the temporal, spatial, and vertical variability of aerosol optical properties over these locations. The HSRL-2 measured these backscatter, extinction, and depolarization profiles during all nineteen science flights and the transit flights to and from the aircraft operations base at Clark International Airport.In this presentation, we use these HSRL-2 measurements to show the horizontal and vertical distributions of aerosols in the CAMP2Ex operational area as well as provide insight into the optical and physical characteristics of these aerosols. We also discuss how these aerosols are distributed within and above the mixed layer.

aerosol↗

A Supersonic Argon/Air Coaxial Jet Experiment for Computational Fluid Dynamics Code Validation

A non-reacting experiment is described in which data has been acquired for the validation of CFD codes used to design high-speed air-breathing engines. A coaxial jet-nozzle has been designed to produce pressure-matched exit flows of Mach 1.8 at 1 atm in both a center jet of argon and a coflow jet of air, creating a supersonic, incompressible mixing layer. The flowfield was surveyed using total temperature, gas composition, and Pitot probes. The data set was compared to CFD code predictions made using Vulcan, a structured grid Navier-Stokes code, as well as to data from a previous experiment in which a He-O2 mixture was used instead of argon in the center jet of the same coaxial jet assembly. Comparison of experimental data from the argon flowfield and its computational prediction shows that the CFD produces an accurate solution for most of the measured flowfield. However, the CFD prediction deviates from the experimental data in the region downstream of x/D = 4, underpredicting the mixing-layer growth rate.

Clifton, Chandler W.↗

Continuous wavelet analysis of coherent structures

We perform an analysis of planar cuts through three dimensional turbulent fields (planar channel flow and mixing layer) using the 2D continuous wavelet transform. We propose two new diagnostics: (1) a measure of intermittency I(r, vector x), which is the ratio of local energy and average energy at a given scale r; and (2) a local Reynolds number, defined on the local velocity contribution at a given scale, computed from the wavelet transform of the three velocity components, the scale of the transform, and molecular viscosity; this gives a representation of the local non-linearity of the flow viewed in both space and scale. We find, for the analyzed flows, strong small-scale intermittency located in the ejection regions for the channel flow and in the vortex core of the mixing layer.

Farge, M.↗

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗

SST Variation Due to Interactive Convective-Radiative Processes

The recent linking of Cloud-Resolving Models (CRMs) to Ocean-Mixed Layer (OML) models has provided a powerful new means of quantifying the role of cloud systems in ocean-atmosphere coupling. This is due to the fact that the CRM can better resolve clouds and cloud systems and allow for explicit cloud-radiation interaction. For example, Anderson (1997) applied an atmospheric forcing associated with a CRM simulated squall line to a 3-D OML model (one way or passive interaction). His results suggested that the spatial variability resulting from the squall forcing can last at least 24 hours when forced with otherwise spatially uniform fluxes. In addition, the sea surface salinity (SSS) variability continuously decreased following the forcing, while some of the SST variability remained when a diurnal mixed layer capped off the surface structure. The forcing used in the OML model, however, focused on shorter time (8 h) and smaller spatial scales (100-120 km). In this study, the 3-D Goddard Cumulus Ensemble Model (GCE; 512 x 512 x 23 cu km, 2-km horizontal resolution) is used to simulate convective active episodes occurring in the Western Pacific warm pool and Eastern Atlantic regions. The model is integrated for seven days, and the simulated results are coupled to an OML model to better understand the impact of precipitation and changes in the planetary boundary layer upon SST variation. We will specifically examine and compare the results of linking the OML model with various spatially-averaged outputs from GCE simulations (i.e., 2 km vs. 10-50 km horizontal resolutions), in order to help understand the SST sensitivity to multi-scale influences. This will allow us to assess the importance of explicitly simulated deep and shallow clouds, as well as the subgrid-scale effects (in coarse-model runs) upon SST variation. Results using both 1-D and 2-D OML models will be evaluated to assess the effects of horizontal advection.

Tao, W.-K.↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

North Atlantic Simulations in Coordinated Ocean-Ice Reference Experiments Phase II (CORE-II) : Inter-Annual to Decadal Variability - Part II

Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958−2007 period from twenty global ocean - sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958−2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.

Ocean model comparisons↗

Automated Identification and Differentiation of Spectrally Similar Hydrothermal Minerals on Mars

Early telescopic observations corroborated hydration related absorptions on Mars in the infrared. Images from the Viking missions led to speculation of hydrothermal alteration and were followed by two missions which mapped the spatial variability of the ~ 3 μm hydration feature. Since then, the Compact Reconnaissance Imager for Mars (CRISM) has provided high spatial resolution (up to 18m) spectral identification of a suite of hydrothermal and diagenetic minerals which have illuminated a range of formation mechanisms. Presence/absence and spatial segregation or mixing of minerals like prehnite, epidote, chlorite amphiboles, and mixed-layer Fe/Mg smectite-chlorite provide valuable evidence for the geologic setting of deposits on Earth, and these phases are often used as temperature and aqueous chemistry indicators in terrestrial systems. Mapping the distribution of these phases will help to answer whether Mars had widespread conditions favorable for low-grade metamorphism and diagenesis, or only focused hydrothermal systems in areas of high heat flow. Further characterizing the chemistry and structure of these phases will then help to answer how most of the widespread Fe/Mg phyllosilicates formed, further defining early geochemical cycling and climate. A fully automated approach for accurate mapping of important hydrothermal mineral phases on Mars has been a challenge. Due to overlapping features in the M-OH region (~2.2-2.4 μm), the strongest absorption features of chlorite, prehnite, and epidote in the short-wave infrared are difficult to distinguish from one another and from the most commonly occurring hydrated silicates on Mars, Fe/Mg smectites. Weaker absorptions are present in both prehnite and epidote which help to distinguish them from chlorite and smectites, but their relative strength in the presence of noise and spatial mixing is often too low to confidently identify them without the noise suppression and feature enhancement methods described here. The spectral signatures of mixed-layer Fe/Mg smectite-chlorite and partially chloritized Fe/Mg smectites have not yet been adequately assessed. Here we evaluate the effectiveness of two empirical and statistical methods for identifying and differentiating these phases using CRISM data.

Rasmussen, B. P.↗

Mixing of the lunar regolith

A probabilistic model for mixing and turnover rates for the lunar regolith due to meteoritic impact is presented and evaluated using results from laboratory impact experiments and estimated meteoritic fluxes. The upper millimeter of the lunar surface is shown to be the primary mixing zone in the regolith and an important source for impact melts and vapors. Below this 'mixing layer' the rate of mixing and turnover decreases very rapidly with increasing depth, consistent with well-preserved stratigraphy and resident times deduced from deep drill core tube samples.

Gault, D. E.↗

Simulations of the stratocumulus-topped boundary layer with a third-order closure model

A third order closure model is proposed by Andre et al. (1982), in which the time rate of change terms, the relaxation and rapid effects for the pressure related terms, and the clipping approximation are included along with the quasi-normal closure, to study turbulence in a cloudy layer which is cooled radiatively from above. A spurious oscillation which is strongest near the inversion occurs. An analysis of the problem shows that the oscillation arises from the mean gradient and buoyancy terms of the triple moment equations; these terms are largest near the cloud top. The oscillation is physical, rather than computational. In nature the oscillation is effectively damped, by a mechanism which apparently is not included in our model. In the stably stratified layer just above the mixed layer top, turbulence can excite gravity waves, whose energy is radiated away. Because the closure assumption for the pressure terms does not take into account the transport of wave energy, the model generates spurious oscillations. Damping of the oscillations is possible by introducing diffusion terms into the triple moment equations. With a large enough choice for the diffusion coefficient, the oscillation is effectively eliminated. The results are quite sensitive to the ad hoc eddy coefficient.

Moeng, C. H.↗

Doppler lidar measurement of profiles of turbulence and momentum flux

A short-pulse CO2 Doppler lidar with 150-m range resolution measured vertical profiles of turbulence and momentum flux. Example measurements are reported of a daytime mixed layer with strong mechanical mixing caused by a wind speed of 15 m/sec, which exceeded the speed above the capping inversion. The lidar adapted an azimuth scanning technique previously demonstrated by radar. Scans alternating between two elevation angles allow determination of mean U-squared, V-squared, and W-squared. Expressions were derived to estimate the uncertainty in the turbulence parameters. A new processing method, partial Fourier decomposition, has less uncertainty than the filtering used earlier.

Eberhard, Wynn L.↗

Budget analysis of the boundary layer grid flights during FIFE 1987

An aircraft grid pattern was flown by the Canadian Twin Otter to map the low-level fluxes and structure over the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) research area in 1987. The time dependence and horizontal advection of heat and moisture were extracted from these flights, combined with surface flux measurements and boundary layer top measurements from radiosondes, to analyze the boundary layer budget using a mixed layer model. The results confirm the suggestion of an earlier study that the boundary layer top entrainment (when parameterized using the buoyancy flux) is nearly double the value used by many modeling studies. Both surface and aircraft data have been revised, and it now appears that the direct measurements of the sensible and latent heat fluxes by the aircraft underestimated these fluxes by about 20 percent, because of filtering and undersampling of long wavelength contributions.

Betts, A. K.↗

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗