Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local mixing”

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 253 records · Page 14

On modeling pressure diffusion in non-homogeneous shear flows

New models are proposed for the 'slow and 'rapid' parts of the pressure diffusive transport based on the examination of DNS databases for plane mixing layers and wakes. The model for the 'slow' part is non-local, but requires the distribution of the triple-velocity correlation as a local source. The latter can be computed accurately for the normal component from standard gradient diffusion models, but such models are inadequate for the cross component. More work is required to remedy this situation.

Demuren, A. O.↗

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that was reinforced with stringers. Shear loading cases the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell element. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. A failure index was calculated by correlating computed mixed-mode failure criterion of the graphite/epoxy material.

Krueger, Ronald↗

Fuel-air mixing apparatus for reducing gas turbine combustor exhaust emissions

A fuel-air mixer for use in a combustion chamber of a gas turbine engine is provided. The fuel air mixing apparatus comprises an annular fuel injector having a plurality of discrete plain jet orifices, a first swirler wherein the first swirler is located upstream from the fuel injector and a second swirler wherein the second swirler is located downstream from the fuel injector. The plurality of discrete plain jet orifices are situated between the highly swirling airstreams generated by the two radial swirlers. The distributed injection of the fuel between two highly swirling airstreams results in rapid and effective mixing to the desired fuel-air ratio and prevents the formation of local hot spots in the combustor primary zone. A combustor and a gas turbine engine comprising the fuel-air mixer of the present invention are also provided as well as a method using the fuel-air mixer of the present invention.

Zupanc, Frank J.↗

Pulling Marbles from a Bag: Deducing the Regional Impact History of the SPA Basin from Impact Melt Rocks

The South Pole-Aitken (SPA) basin is an important target for absolute age-dating. Vertical and lateral impact mixing ensures that regolith within SPA will contain rock fragments from SPA itself, local impact craters, and faraway giant basins. About 20% of the regolith at any given site is foreign [1, 2], but much of this material will be cold ejecta, not impact melt. We calculated the fraction of contributed impact melt using scaling laws to estimate the amount and provenance of impact melt, demonstrating that SPA melt is the dominant impact melt rock (>70%) likely to be present. We also constructed a statistical model to illustrate how many randomly-selected impact-melt fragments would need to be dated, and with what accuracy, to confidently reproduce the impact history of a site. A detailed impact history becomes recognizable after a few hundred to a thousand randomly-selected marbles, however, it will be useful to have more information (e.g. compositional, mineralogical, remote sensing) to group fragments. These exercises show that SPA melt has a high probability of being present in a scoop sample and that dating of a few hundred to a thousand impact-melt fragments will yield the impact history of the SPA basin.

Cohen, Barbara A.↗

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↗

Benchmarking Bayesian Optimization Frameworks and Acquisition Strategies for Materials Discovery and Autonomous Laboratories

Bayesian optimization (BO) can accelerate materials discovery by guiding expensive experiments toward the most promising processing conditions. We systematically compare five BO surrogate and framework combinations (Gaussian processes in Ax, Gaussian processes and Monte-Carlo neural networks in BayBE, random forests in Lolopy, and tree-structured Parzen (TPE) estimators in Hyperopt) on three benchmarks that mimic common materials design tasks (a discrete solid-electrolyte composition space, a hybrid discrete/continuous laminate-composite design problem solved with micromechanics modeling, and the continuous Ishigami analytic function which is a standard optimization benchmark). Each BO surrogate is paired with posterior mean, probability of improvement, and expected improvement acquisition functions and run for 100 trials from randomized initial samples with uniform random search providing a control. Across five random seeds per setting, BayBE’s Gaussian-process surrogate with expected improvement consistently reached ≥95 % of the known optimum in the fewest evaluations, while Lolopy’s random forest matched or exceeded GP performance on purely categorical or mixed spaces at a higher computational cost. Posterior mean alone often stagnated at local optima, underscoring the need for exploration, whereas probability and expected improvement balanced exploration and exploitation leading to better optimization in fewer trials. Execution times ranged from milliseconds for TPE to minutes for neural-network and random-forest surrogates. These results establish baseline expectations for BO in automated materials laboratories and highlight expected improvement with Gaussian processes as a reliable first choice, with random forests offering a strong alternative when categorical variables dominate. The benchmark suite and code are released to facilitate future surrogate, acquisition, and constraint-handling research in data-driven materials optimization.

Bayesian optimization↗

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modeled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to a single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modelled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

HSRL-2 Observations of Aerosol Variability and Mixing During Boundary Layer Evolution in Houston

The NASA Langley airborne multi-wavelength High Spectral Resolution Lidar (HSRL-2) provides vertical distribution of aerosol optical properties as ?curtains? of aerosol extinction, backscatter and depolarization along the flight track, plus intensive properties that are used to infer aerosol type and external mixing of types. Deployed aboard the NASA Langley King Air on the DISCOVER-AQ field mission in Houston in September 2013, HSRL-2 flew a pattern that included 18 ground sites, repeated four times a day, coordinated with a suite of airborne in situ measurements. The horizontally and vertically resolved curtains of HSRL-2 measurements give an unparalleled view of the spatial and temporal variability of aerosol, which provide broad context for interpreting other measurements and models. In Houston, HSRL-2 generally observed significant variability with distinct layering: boundary layer, residual layer, and frequent upper layers of smoke transported from the Mississippi Valley. The period from Sep. 11-14 is notable for a large aerosol build-up and persistent layers in the free troposphere. We investigate the aerosol properties and evolution using the vertically resolved HSRL-2 measurements, typing and mixture analysis techniques, and boundary layer detection. Between morning and afternoon overpasses, as the boundary layer grows, many distinctions between the layers are lost as the aerosols become mixed. As the boundary layer collapses overnight, the aerosols are cut off and are observed in a distinct residual layer the following day. HSRL-2 measurements of the upper smoke layers suggest slightly different properties each day as new smoke enters the region, while the morning boundary layer indicates more similarity in local emissions day-to-day. HSRL-2 intensive variables (indicators of aerosol type) reflect complex yet predictable mixing. We will present the analysis of aerosol mixtures, and explore the WRF-Chem chemical transport model along the HSRL-2 flight path.

Sharon P Burton↗

Radiometric calibration of Landsat Thematic Mapper Thermal Band

Radiometric calibration of satellite-acquired data is essential for quantitative scientific studies, as well as for a variety of image-processing applications. This paper describes a multiyear, on-orbit radiometric calibration of the Landsat Thematic Mapper (TM) Band 6 conducted at DOE's Pacific Northwest Laboratory. Numerous Landsat TM scenes acquired and analyzed included day and night coverages at several geographical locations over several seasons. Concurrent with Landsat overpasses, thermal field and local meteorological (surface and radiosonde) measurements were collected. At-satellite (uncorrected) radiances and temperatures for water and nonwater land cover were compared to ground truth (GT) measurements after making adjustments for atmospheric (using LOWTRAN), mixed-pixel, and emissivity effects. Results indicate that, for both water and nonwater features, TM Band 6 average corrected temperature determinations using local radiosonde data to adjust for atmospheric effects, and using appropriate emissivities, are within 1.0 C of GT temperature values. Temperatures of water pixels derived from uncorrected TM Band 6 data varied roughly between 1 and 3 C of ground truth values for water temperatures ranging between 4 and 24 C. Moreover, corrections using nonlocal and noncoincident radiosonde data resulted in errors as large as 12 C. Corrections using the U.S. Standard Atmosphere gave temperature values within 1 to 2 C of GT. The average uncertainty for field instruments was + or - 0.2 C; average uncertainty for Landsat TM corrected temperature determinations was + or - 0.4 C.

Wukelic, G. E.↗

Optical Polarization as a Probe of the Local Interstellar Medium

The use of interstellar polarization as a tool for measuring interstellar dust is discussed. Problems resulting from dust and magnetic field configurations becoming mixed up are discussed, as is the availability of sufficiently bright stars to obtain the photons needed for precision measurements. It is proposed that: (1) on the scale of several hundred parsec, there is a preferential magnetic field direction, as evidenced by observations at the Galactic poles and selected longitudes in the Galactic plane; (2) the local (r 50 pc) region is devoid of dust, as evidenced by the mean square degree of polarization as a function of distance; and, less certainly, that (3) at a distance of less than 5 pc, there is a patch of dust which may be of interest in connection with cloud models.

Tinbergen, J.↗

Tidal Energy Available for Deep Ocean Mixing: Bounds from Altimetry Data

Maintenance of the large-scale thermohaline circulation has long presented an interesting problem. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Ray, Richard D.↗

Tidal Energy Available for Deep Ocean Mixing: Bounds From Altimetry Data

Maintenance of the large-scale thermohaline circulation has long presented a problem to oceanographers. Observed mixing rates in the pelagic ocean are an order of magnitude too small to balance the rate at which dense bottom water is created at high latitudes. Recent observational and theoretical work suggests that much of this mixing may occur in hot spots near areas of rough topography (e.g., mid-ocean ridges and island arcs). Barotropic tidal currents provide a very plausible source of energy to maintain these mixing processes. Topex/Poseidon satellite altimetry data have made precise mapping of open ocean tidal elevations possible for the first time. We can thus obtain empirical, spatially localized, estimates of barotropic tidal dissipation. These provide an upper bound on the amount of tidal energy that is dissipated in the deep ocean, and hence is available for deep mixing. We will present and compare maps of open ocean tidal energy flux divergence, and estimates of tidal energy flux into shallow seas, derived from T/P altimetry data using both formal data assimilation methods and empirical approaches. With the data assimilation methods we can place formal error bars on the fluxes. Our results show that 20-25% of tidal energy dissipation occurs outside of the shallow seas, the traditional sink for tidal energy. This suggests that up to 1 TW of energy may be available from the tides (lunar and solar) for mixing the deep ocean. The dissipation indeed appears to be concentrated over areas of rough topography.

Egbert, Gary D.↗

Thermonuclear runaways in nova outbursts. 2: Effect of strong, instantaneous, local fluctuations

In an attempt to understand the manner in which nova outbursts are initiated on the surface of a white dwarf, we investigate the effects fluctuations have on the evolution of a thermonuclear runaway. Fluctuations in temperature density, or the composition of material in the burning shell may arise due to the chaotic flow field generated by convection when it occurs, or by the accretion process itself. With the aid of two-dimensional reactive flow calculations, we consider cases where a strong fluctutation in temperature arises during the early, quiescent accretion phase or during the later, more dynamic, explosion phase. In all cases we find that an instantaneous, local temperature fluctuation causes the affected material to become Rayleigh-Taylor unstable. The rapid rise and subsequent expansion of matter immediately cools the hot blob, which prevents the lateral propagation of burning. This suggests that local temperature fluctuations do not play a significant role in directly initiating the runaway, especially during the early stages. However, they may provide an efficient mechanism of mixing core material into the envelope (thereby pre-enriching the fuel for subsequent episodes of explosive hydrogen burning) and of mixing substantial amounts of the radioactive nucleus N-13 into the surface layers, making novae potential gamma-ray sources. This suggests that it is the global not the local, evolution of the core-envelope interface to high temperatures which dominates the development of the runaway. We also present a possible new scenario for the initiation of nova outbursts based on our results.

Shankar, Anurag↗

Turbulent compressible convection in a deep atmosphere. I - Preliminary two-dimensional results

The turbulent convection of a compressible fluid in a deep atmosphere is simulated by two-dimensional numerical computations, displaying coexisting 'cells' whose sizes range from the total depth of the convection zone to the smallest scale height at the top. While the largest cells traverse the entire zone, smaller ones are clustered near the top. The vertical correlation length of the vertical velocity is proportional to the local pressure or density scale height, in a way reminiscent of the concept of mixing length.

Chan, K. L.↗