Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Network Scale Ubiquitous Volume Estimation Using Tree-Based Ensemble Learning Methods

Currently real-time, ubiquitous volume data for roadway networks remains the key missing dimension in traffic operations. Most volume data are average annual daily traffic (AADT) measures derived from the Highway Performance Monitoring System (HPMS). Although methods to factor the AADT to hourly averages for typical day of week exist, actual volume data is limited to a sparse collection of locations in which volumes are continuously recorded. This paper explores the use of state-of-art machine learning techniques to estimate accurate real-time volume measures that span the highway network providing ubiquitous coverage in space, and point-in-time measures for a specific date and time. Three tree-based ensemble learning models, random forest (RF), gradient boost machine (GBM), and extreme gradient boost (XGBoost), were tested for volume estimation by learning from combined dataset of commercial probe data provided by TomTom, the FHWA's Travel Monitoring Analysis System (TMAS) data, and other infrastructure attributes such as number of lanes, speed limit, and weather. The methods were tested on major corridors and freeways in the metropolitan area of Denver. All three machine learning methods were able to provide hourly volume estimates 24 hours a day, 7 days a week, and 365 days a year with around 18% mean absolute error to true volume and about 5% of error with respect to roadway capacity. The low error measures allow the potential application by transportation agencies.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Highly cascaded first-order sapphire optical fiber Bragg gratings fabricated by a femtosecond laser

This Letter reports an innovative technique for fabricating large-scale, highly cascaded first-order sapphire optical fiber Bragg gratings (FBGs) using a femtosecond laser-assisted point-by-point inscription method. For the first time, to the best of our knowledge, this study successfully demonstrates a distributed array of 10 FBGs within highly multimode sapphire crystal fiber, made possible by employing a high-power laser technique to generate larger reflectors with a Gaussian intensity profile. These first-order FBGs offer advantages such as enhanced reflectivity, shorter fabrication time, and simplified spectral characteristics, making them easier to interpret compared with high-order FBGs. The FBGs’ resilience and effectiveness are analyzed by subjecting them to temperature tests, proving their capacity for accurate temperature monitoring up to 1500°C—a testament to their suitability for harsh environments. This novel approach broadens the scope for sensing and communication applications in sapphire fibers, particularly under challenging conditions. The novelty of our work lies in successfully overcoming the limitations of previous designs by integrating a cascade of 10 FBGs in sapphire fibers, thereby enhancing multiplexing capabilities, minimizing overlapping of FBG peaks, and ensuring reliable temperature monitoring in industries and applications with thermal gradients.

47 OTHER INSTRUMENTATION↗
NASA NTRS1994

Modeling Thermal Contact Resistance

One difficulty in using cryocoolers is making good thermal contact between the cooler and the instrument being cooled. The connection is often made through a bolted joint. The temperature drop associated with this joint has been the subject of many experimental and theoretical studies. The low temperature behavior of dry joints have shown some anomalous dependence on the surface condition of the mating parts. There is also some doubts on how well one can extrapolate from the test samples to predicting the performance of a real system. Both finite element and analytic models of a simple contact system have been developed. The model assumes (a) the contact is dry (contact limited to a small portion of the total available area and the spaces in-between the actual contact patches are perfect insulators), (b) contacts are clean (conductivity of the actual contact is the same as the bulk), (c) small temperature gradients (the bulk conductance may be assumed to be temperature independent), (d) the absolute temperature is low (thermal radiation effects are ignored), and (e) the dimensions of the nominal contact area are small compared to the thickness of the bulk material (the contact effects are localized near the contact). The models show that in the limit of actual contact area much less than the nominal area (a much less than A), that the excess temperature drop due to a single point of contact scales as a(exp -1/2). This disturbance only extends a distance approx. A(exp 1/2) into the bulk material. A group of identical contacts will result in an excess temperature drop that scales as n(exp -1/2), where n is the number of contacts and n dot a is constant. This implies that flat rough surfaces will have a lower excess temperature drop than flat polished surfaces.

Kittel, Peter↗
NASA NTRS2002

Rayleigh-Taylor Gravity Waves and Quasiperiodic Oscillation Phenomenon in X-ray Binaries

Accretion onto compact objects in X-ray binaries (black hole, neutron star (NS), white dwarf) is characterized by non-uniform flow density profiles. Such an effect of heterogeneity in presence of gravitational forces and pressure gradients exhibits Rayleigh-Taylor gravity waves (RTGW). They should be seen as quasiperiodic wave oscillations (QPO) of the accretion flow in the transition (boundary) layer between the Keplerian disk and the central object. In this paper the author shows that the main QPO frequency, which is very close to the Keplerian frequency, is split into separate frequencies (hybrid and low branch) under the influence of the gravitational forces in the rotational frame of reference. The RTGWs must be present and the related QPOs should be detected in any system where the gravity, buoyancy and Coriolis force effects cannot be excluded (even in the Earth and solar environments). The observed low and high QPO frequencies are an intrinsic signature of the RTGW. The author elaborates the conditions for the density profile when the RTGW oscillations are stable. A comparison of the inferred QPO frequencies with QPO observations is presented. The author finds that hectohertz frequencies detected from NS binaries can be identified as the RTGW low branch frequencies. The author also predicts that an observer can see the double NS spin frequency during the NS long (super) burst events when the pressure gradients and buoyant forces are suppressed. The Coriolis force is the only force which acts in the rotational frame of reference and its presence causes perfect coherent pulsations with a frequency twice of the NS spin. The QPO observations of neutron binaries have established that the high QPO frequencies do not go beyond of the certain upper limit. The author explains this observational effect as a result of the density profile inversions. Also the author demonstrates that a particular problem of the gravity waves in the rotational frame of reference in the approximation of very small pressure gradients is reduced to the problem of the classical oscillator in the rotational frame of reference which was previously introduced and applied for the interpretation of kHZ QPO observation by Osherovich & Titarchuk.

Titarchuk, Lev↗
DOE OSTI2020

Soil microbial function and fungal community data along an ectomycorrhizal basal area gradient, Vermont, May 2019

Mycorrhizal fungi can affect soil organic matter cycling through several mechanisms including priming, nutrient competition, and direct enzyme production. Differences in nutrient foraging strategies between ectomycorrhizal (EcM) and arbuscular mycorrhizal (AM) fungi produce divergent belowground dynamics: where EcM can take up organic nitrogen and directly break down SOM by producing enzymes, AM fungi are limited to scavenging mineral N. EcM-associated tree species also have leaf litter with relatively higher ratios of carbon to nitrogen (C:N), and belowground saprotrophic communities more dominated by fungi. Consequently, free-living microbes in EcM-dominated soils should experience nitrogen limitation, with subsequent increases in enzyme production and decreased carbon use efficiency (CUE). However, the relative importance of the effects of substrate quality and fungal community composition on enzyme production and CUE are unclear. To assess this distinction, we sampled the organic horizon and 10 cm of the mineral horizon in northern temperate forest soils along a gradient of EcM dominance (see SCxTS_DBH.csv for full tree species information). We characterized fungal community composition by measuring EcM relative abundances from extracted fungal DNA (To run raw ASV data through R code, use SCxTS_ASV_taxonomy.csv, and R code ITS_processing.Rmd; to see processed data with EcM relative abundance, use file SCxTS_sum_afitch_Feb18.20.csv) and the fungal to bacterial ratios (see data file SCxTS_plfa.csv) from phospholipid fatty acid analysis. Soil microbial functions were measured as potential activities of five hydrolytic and two oxidative enzymes ( To run raw data through R code us SCxTS_metadata.csv, hydrolytic_oxidative_raw.csv, MUBraw.csv, Enzyme_soilweight.csv, and incubationtime.csv, and R code SCxTS_enzymes.Rmd; to use Enzyme data ready to run through statistics code, use file Enzyme.data.csv) and microbial carbon use efficiency (see file SCxTS_CUE.csv and for full calculations file 18O_CUEcalcs_SCxTS_2019.xlsx). We assessed soil substrate quality as the soil carbon:nitrogen ratio (See file Organic_CN.csv). Statistical analysis can be run using R code from file SCxTS_analysis.Rmd.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The magnetic gradient scale length explains why certain plasmas require close external magnetic coils

Abstract The separation between the last closed flux surface of a plasma and the external coils that magnetically confine it is a limiting factor in the construction of fusion-capable plasma devices. This plasma-coil separation must be large enough so that components such as a breeding blanket and neutron shielding can fit between the plasma and the coils. Plasma-coil separation affects reactor size, engineering complexity, and particle loss due to field ripple. For some plasmas it can be difficult to produce the desired flux surface shaping with distant coils, and for other plasmas it is infeasible altogether. Here, we seek to understand the underlying physics that limits plasma-coil separation and explain why some configurations require close external coils. In this paper, we explore the hypothesis that the limiting plasma-coil separation is set by the shortest scale length of the magnetic field as expressed by the ∇ B tensor. We tested this hypothesis on a database of > 40 stellarator and tokamak configurations. Within this database, the coil-to-plasma distance compared to the minor radius varies by over an order of magnitude. The magnetic scale length is well correlated to the coil-to-plasma distance of actual coil designs generated using the REGCOIL method (Landreman 2017 Nucl. Fusion 57 046003). Additionally, this correlation reveals a general trend that larger plasma-coil separation is possible with a small number of field periods.

Physics↗
DOE OSTI2014

Band Limited Correlation Estimates for A(ξω/U) and B(ηω/U) Using Beresh et. al. 2013 Data Sets

Here we discuss an improved Corcos (Corcos (1963), (1963)) style cross spectral density utilizing zero pressure gradient, supersonic (Beresh et. al. (2013)) data sets. Using the connection between narrow band measurements with broadband cross-spectral density, i.e. Γ(ξ ,η ,ω )= Φ (ω) A(ωη/U )exp (-i ωξ/U) we focus on estimating coherence expressions of the form: A (ξω nb /U) and B (ηω nb / U) where ω nb denotes the narrow band frequency, i.e. the band center frequency value and ξ and η are sensors spacing in streamwise/longitudinal and cross-stream/lateral directions, respectively. A methodology to estimate the parameters which retains the Corcos exponential functional form, A(ξω/U)=exp(-k lat ηω/U) but identifies new parameters (constants) consistent with the Beresh et. al. data sets is discussed. The Corcos result requires that the data be properly explained by self-similar variable: ξω/U and ηω/U. The longitudinal (streamwise) variable ξω/U tends to provide a better data collapse, while, consistent with the literature the lateral ηω/U is only successful for higher band center frequencies. Assuming the similarity variables provide a useful description of the data, the longitudinal coherence decay constant result using the Beresh et. al. data sets yields a value for the longitudinal constant k long ≈0.36-0.28 that is approximately 3x larger than the “traditional” (low speed, large Reynolds number and zero pressure gradient) of k long ≈0.11. We suggest that the most likely reason that the Beresh et. al. data sets incur increased longitudinal decay which results in reduced coherence lengths is due to wall shear induced compression causing an adverse pressure gradient. Focusing on the higher band center frequency measurements where the frequency dependent similarity variables are applicable, the lateral or transverse coherence decay constant k lat ≈0.7 is consistent with the “traditional” (low speed, large Reynolds number and zero pressure gradient). It should be noted, that the longitudinal/streamwise coherence decay deviates from the value observed by other researchers while the lateral/ cross-stream value is consistent has been observed by other researchers. We believe that while the measurements used to obtain new decay constant estimates are from internal wind tunnel tests, they likely provide a useful estimate expected reentry flow behavior and are therefore recommended for use. These data could also be useful in determining the uncertainty of correlation length for a uncertainty quantification (UQ) analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Crystalline Bilayer Graphene with Preferential Stacking from Ni–Cu Gradient Alloy

We developed a high-yield synthesis of highly crystalline bilayer graphene (BLG) with two preferential stacking modes using a Ni–Cu gradient alloy growth substrate. Previously reported approaches for BLG growth include flat growth substrates of Cu or Ni–Cu uniform alloys and “copper pocket” structures. Use of flat substrates has the advantage of being scalable, but the growth mechanism is either “surface limited” (for Cu) or carbon precipitation (for uniform Ni–Cu), which results in multicrystalline BLG grains. For copper pockets, growth proceeds through a carbon back-diffusion mechanism, which leads to the formation of highly crystalline BLG, but scaling of the copper pocket structure is expected to be difficult. Here we demonstrate a Ni–Cu gradient alloy that combines the advantages of these earlier methods: the substrate is flat, so easy to scale, while growth proceeds by a carbon back-diffusion mechanism leading to high-yield growth of BLG with high crystallinity. Here, the BLG layer stacking was almost exclusively Bernal or twisted with an angle of 30°, consistent with first-principles calculations we conducted. Furthermore, we demonstrated scalable production of transistor arrays based crystalline Bernal-stacked BLG with a band gap that was tunable at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Large-scale anisotropy of TeV-band cosmic rays

The expected anisotropy in the 1 to 10{sup 4} TeV energy range is calculated for Galactic cosmic rays with both anisotropy in the diffusion tensor and source discreteness taken into account. We find that if the sources are distributed radially (but with azimuthal symmetry) in proportion to Galactic pulsars, the expected anisotropy almost always exceeds the observational limits by one order of magnitude in the case of isotropic diffusion. If the radial diffusion is more than an order of magnitude smaller than the azimuthal diffusion rate, the radial gradient of the sources can be accommodated about 5% of the time. If the sources are concentrated in the spiral arms, then the anisotropy depends on our location between them, but in some spatial window, roughly equidistant from adjacent spiral arms, the observational constraints on anisotropy are obeyed roughly 20%-30% of the time for extremely anisotropic diffusion. The solar system is in that window less than 10% of the time, but it may be there now. Under the assumption of isotropic diffusion, nearby supernovae are found to produce a discreteness anisotropy that is nearly two orders of magnitude in excess of the observational limit if all supernovae are assumed to contribute equally with a source rate 1 in every 100 years.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

The stability of coronal loops - Finite-length and pressure-profile limits

Results are described from a quickly converging, necessary-and-sufficient, MHD-stability test for coronal-loop models. The primary stabilizing influence arises from magnetic line tying at the photosphere, and this end conditions requires a series expansion of possible loop excitations. The stability boundary is shown to quickly approach a limit as the number of terms increases, providing a critical length for the loop in proportion to its transverse magnetic scale. Several models of force-free-field profiles are tested and the stability behavior of a localized current channel, embedded in an external current-free region, is shown to be superior to that of other, broader, current profiles. Pressure-gradient effects, leading to increased or decreased stability, are shown to be amplified by line tying. Long loops must either conduct low net current, or exhibit an axial-field reversal coexisting with a low-pressure core. The limits on stability depend on the magnetic aspect ratio, the plasma-to-magnetic pressure ratio, and the field orientation at the loop edge. Applications of these results to the structure of coronal loops are described.

Einaudi, G.↗
DOE OSTI2020

Rewilding of Fukushima's human evacuation zone

There is substantial interest in understanding the ecological impacts of the nuclear accidents at the Chernobyl and Fukushima Daiichi nuclear power plants. However, population‐level data for large mammals have been limited, and there remains much speculation regarding the status of wildlife species in these areas. Using a network of remote cameras placed along a gradient of radiological contamination and human presence, we collected data on population‐level impacts to wildlife (that is, abundance and occupancy patterns) following the 2011 Fukushima Daiichi nuclear accident. We found no evidence of population‐level impacts in mid‐ to large‐sized mammals or gallinaceous birds, and show several species were most abundant in human‐evacuated areas, despite the presence of radiological contamination. These data provide unique evidence of the natural rewilding of the Fukushima landscape following human abandonment, and suggest that if any effects of radiological exposure in mid‐ to large‐sized mammals in the Fukushima Exclusion Zone exist, they occur at individual or molecular scales, and do not appear to manifest in population‐level responses.

Lyons, Phillip C.↗
DOE OSTI2024

Speeding up and reducing memory usage for scientific machine learning via mixed precision

Scientific machine learning (SciML) has emerged as a versatile approach to address complex computational science and engineering problems. Within this field, physics-informed neural networks (PINNs) and deep operator networks (DeepONets) stand out as the leading techniques for solving partial differential equations by incorporating both physical equations and experimental data. However, training PINNs and DeepONets require significant computational resources, including long computational times and large amounts of memory. In search of computational efficiency, training neural networks using half precision (float16) rather than the conventional single (float32) or double (float64) precision has gained substantial interest, given the inherent benefits of reduced computational time and memory consumed. However, we find that float16 cannot be applied to SciML methods, because of gradient divergence at the start of training, weight updates going to zero, and the inability to converge to a local minima. To overcome these limitations, we explore mixed precision, which is an approach that combines the float16 and float32 numerical formats to reduce memory usage and increase computational speed. Our experiments showcase that mixed precision training not only substantially decreases training times and memory demands but also maintains model accuracy. Here, we also reinforce our empirical observations with a theoretical analysis. The research has broad implications for SciML in various computational applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An analysis of the structure of an n-dodecane spray flame using TPDF modelling

With a view to understanding ignition and combustion behaviours in diesel engines, this study investigates several aspects of ignition and combustion of an n-dodecane spray in a high pressure, high temperature chamber, known as Spray A, using data resulting from modelling using the transported probability density function (TPDF) method. The model has been validated comprehensively with good to excellent agreement in our previous work against all available experimental data including for mixture-fraction and velocity fields in non-reacting cases, and flame lift-off length and ignition delay in reacting cases. This good agreement encourages further investigation of the numerical model results to help understand the structure of this flame, which serves to complement the experimental information that is available, which is very limited due to the difficult experimental conditions in which this flame exists. For example, quantitative experimental measurements of local mixture-fraction, temperature, velocity gradients, etc. are not yet possible in reacting cases. Analysis of the model results shows that two-stage ignition is found to occur across the ambient temperature conditions considered: the first stage is rapidly initiated on the lean side where temperatures are high and sequentially moves to richer, cooler conditions. The first stage is extremely resilient to turbulence, occurring in a region of very low Damkohler number. The second stage of ignition occurs first in rich mixtures in a region behind the head of the fuel jet where mixture gradients are low, and appears to be influenced strongly by turbulence. Relative to a homogeneous reactor, it is delayed on the lean side but advanced on the rich side, suggesting entrainment and mixing from the early igniting lean regions into richer mixtures is an important moderator of the ignition process. The second-stage ignition front propagates at very high velocities initially, suggesting it is a sequential ignition moving according to gradients of ignition delay and/or residence time. The flame stabilises however on the lean side in a region of much lower velocity, where turbulent velocity fluctuations are sufficiently high such that turbulent transport influences the propagation. It stabilises in a region of low Damkohler number which implies that a competition of chemistry versus micro-mixing might also be involved in stabilisation. The stabilisation mechanism is investigated by an analysis of the transport budgets, showing the flame is stabilised by autoignition but moderated by turbulent diffusion. Further analysis of the flame index supports this stabilisation mechanism, and demonstrates the simultaneous existence of non-premixed and premixed combustion modes in the same flame. Analysis of the flow fields also reveals that local entrainment and dilatation are important flow features near the flame base. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Diesel↗
DOE OSTI2022

Machine learning moment closure models for the radiative transfer equation I: Directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for the radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. Finally, a variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

Computer science↗
DOE OSTI2016

Layer-by-Layer Sorting of Rhenium Disulfide via High-Density Isopycnic Density Gradient Ultracentrifugation

Isopycnic density gradient ultracentrifugation (iDGU) has been widely applied to sort nanomaterials by their physical and electronic structure. Yet, the commonly used density-gradient medium iodixanol has a finite maximum buoyant density that prevents the use of iDGU for high-density nanomaterials. Here, we overcome this limit by adding cesium chloride (CsCl) to iodixanol, thus increasing its maximum buoyant density to the point where the high-density two-dimensional nanomaterial rhenium disulfide (ReS 2 ) can be sorted in a layer-by-layer manner with iDGU. The resulting aqueous ReS 2 dispersions show photoluminescence at ~1.5 eV, which is consistent with its direct bandgap semiconductor electronic structure. Moreover, photocurrent measurements on thin films formed from solution-processed ReS 2 show a spectral response that is consistent with optical absorbance and photoluminescence data. In addition to offering a pathway for effective solution processing of ReS 2 , this work establishes a general methodology for sorting high-density nanomaterials via iDGU.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impacts of Latent Energy and Snow Fall Speed on a Wintertime Midlatitude Cyclone

This study aimed to understand the impacts of latent energy as well as snow fall speeds on precipitation properties and synoptic-scale storm characteristics of a wintertime midlatitude cyclone. Simulations of a potent winter storm that impacted the Rocky Mountains and northern Great Plains of the United States in February 2017 were performed using the Weather Research and Forecasting model with latent heating or cooling from individual microphysical processes systematically turned off and fall speeds adjusted by ±50%. Results indicated substantial impacts on the microphysical characteristics of the simulated storm to fall speed, cooling from sublimation, and warming from deposition. The impacts of cooling and warming were manifested as differences in accumulated snowfall. Increased (decreased) fall speeds led to smaller (larger) ice crystals and total mass, resulting in offsetting effects on the precipitation flux, and thus minimal impacts on snowfall and large-scale characteristics of the storm. Warming and cooling associated with deposition and sublimation, respectively, impacted the synoptic-scale dynamics, whereby removing warming from deposition resulted in an increased meridional temperature gradient near the jet stream, thus increasing the jet strength and causing it to be more westerly with less curvature aloft. This in turn limited upper-level divergence, creating a weaker surface low and shifting the precipitation shield southward. The opposite occurred with the removal of latent cooling due to sublimation. This study highlights the potential importance of latent energy associated with ice sublimation and deposition and fall speeds in the larger-scale characteristics of winter storms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Photon acceleration in plasma wake wave

The photon acceleration effect in a laser wake field is investigated based on photon Hamiltonian dynamics. A test laser pulse is injected into a plasma wave at an incident angle θ{sub i}, which could slow down the photon velocity along the propagating direction of the wake wave so as to increase the acceleration distance for the photons. The photon trapping condition is analyzed in detail, and the maximum frequency shift of the trapped photon is obtained. The acceleration gradient and dephasing length are emphatically studied. The compression of the test laser pulse is examined and used to interpret the acceleration process. The limit of finite transverse width of the wake wave on photon acceleration is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗