Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.4 is a program that obtains a steady-state flow and temperature solution, in geometrical extent, for a nuclear reactor core, or a single fuel assembly, including the effect of manufacturing and modelling uncertainties by the means of hot channel factors. It is based on an evolutionary sequence of codes originally used for calculating plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over the past 30 years. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability for select design and safety parameters is available to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
DOE OSTI2021

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗
DOE OSTI2023

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗
DOE OSTI2026

Oil Filled Flow Shield Performance Evaluation

When hydrophones are used to monitor underwater sound in areas with currents, like rivers or tidal channels, flow noise caused by turbulent flow past a hydrophone can mask low frequency acoustic signals of interest. Flow shields can reduce the impacts of flow noise on hydrophone data. A prior study conducted by PNNL compared three flow shield designs and determined that an oil-filled urethane shield performed better than other designs and was likely sufficiently robust for long-duration marine deployments. However, that study did not assess the high frequency performance of the oil-filled flow shield or conduct long-duration field testing. Further, the prototype flow shield was difficult to assemble. This report details the results of design refinement and testing to address these gaps. Design refinements resulted in a straightforward and repeatable assembly method. Field testing indicated that the flow shields reduce flow noise by over 20 dB below 50 Hz at flow speeds greater than 0.5 m/s and are proven to survive deployments more than 6 months long. Calibration of the same hydrophone with and without a flow shield installed indicated that the shields attenuate propagating sound by no more than 2 dB at frequencies below 30 kHz.

42 ENGINEERING↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2024

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2020

Comparison of loads from HAWC2 and OpenFAST for the IEA Wind 15 MW Reference Wind Turbine

Reference wind turbines (RWTs) that reflect the state-of-the-art of current wind energy technology are necessary in order to properly evaluate innovative methods in wind turbine design and evaluation. The International Energy Agency (IEA) Wind Technology Collaboration Platform (TCP) Task 37 has recently developed a new RWT geared towards offshore floating-foundation applications: the IEA Wind 15 MW. The model has been implemented in two aeroelastic codes, OpenFAST and HAWC2, based on an underlying common ontology. However, these toolchains result in slightly different structural parameters, and the two codes utilise different structural models. Thus, to increase the utility of the model, it is necessary to compare the aeroelastic responses. This paper compares aeroelastic loads calculated using different fidelities of the blade model in OpenFAST (ElastoDyn and BeamDyn) and HAWC2 (prismatic Timoshenko without torsion and Timoshenko with fully populated stiffness matrix), where both codes use the DTU Basic controller and the same turbulence boxes to reduce discrepancies. The aeroelastic responses to steady wind, step wind and turbulent wind (per IEC 61400-1 wind class IB) are considered. The results indicate a generally good agreement between the loads dominated by aerodynamic thrust and force, especially for the no-torsion blade models. Discrepancies were observed in other load channels, partially due to differences in the asymmetric loading of the rotor and partially due to differing closed-loop dynamics, and they will be the subject of future investigations.

17 WIND ENERGY↗
DOE OSTI2024

Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble

Abstract. We explore the cloud system evolution of non-precipitating marine stratocumuli with a focus on the impacts of the diurnal cycle and free-tropospheric (FT) humidity based on an ensemble of 244 large-eddy simulations generated by perturbing initial thermodynamic profiles and aerosol conditions. Cases are categorized based on their degree of decoupling and the cloud liquid water path (LWPc, based on model columns with cloud optical depths greater than one). A budget analysis method is proposed to analyze the evolution of cloud water in both coupled and decoupled boundary layers. More coupled clouds start with a relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of net radiative cooling. During the nighttime, a positive correlation between FT humidity and the LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The LWPc is more likely to decrease during the nighttime for a larger LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation.

Chen, Yao-Sheng (ORCID:0000000208354132)↗
DOE OSTI2025

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗
DOE OSTI2021

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗
DOE OSTI2023

Airborne flux measurements of ammonia over the southern Great Plains using chemical ionization mass spectrometry

Abstract. Ammonia (NH3) is an abundant trace gas in the atmosphere and an important player in atmospheric chemistry, aerosol formation and the atmosphere–surface exchange of nitrogen. The accurate determination of NH3 emission rates remains a challenge, partly due to the propensity of NH3 to interact with instrument surfaces, leading to high detection limits and slow response times. In this paper, we present a new method for quantifying ambient NH3, using chemical ionization mass spectrometry (CIMS) with deuterated benzene cations as reagents. The setup aimed at limiting sample–surface interactions and achieved a 1σ precision of 10–20 pptv and an immediate 1/e response rate of < 0.4 s, which compares favorably to the existing state of the art. The sensitivity exhibited an inverse humidity dependence, in particular in relatively dry conditions. Background of up to 10 % of the total signal required consideration as well, as it responded on the order of a few minutes. To showcase the method's capabilities, we quantified NH3 mixing ratios from measurements obtained during deployment on a Gulfstream I aircraft during the HI-SCALE (Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems) field campaign in rural Oklahoma during May 2016. Typical mixing ratios were 1–10 parts per billion by volume (ppbv) for the boundary layer and 0.1–1 ppbv in the lower free troposphere. Sharp plumes of up to tens of ppbv of NH3 were encountered as well. We identified two of their sources as a large fertilizer plant and a cattle farm, and our mixing ratio measurements yielded upper bounds of 350 ± 50 and 0.6 kg NH3 h−1 for their respective momentary source rates. The fast response of the CIMS also allowed us to derive vertical NH3 fluxes within the turbulent boundary layer via eddy covariance, for which we chiefly used the continuous wavelet transform technique. As expected for a region dominated by agriculture, we observed predominantly upward fluxes, implying net NH3 emissions from the surface. The corresponding analysis focused on the most suitable flight, which contained two straight-and-level legs at ∼ 300 m above ground. We derived NH3 fluxes between 1 and 11 mol km−2 h−1 for these legs, at an effective spatial resolution of 1–2 km. The analysis demonstrated how flux measurements benefit from suitably arranged flight tracks with sufficiently long straight-and-level legs, and it explores the detrimental effect of measurement discontinuities. Following flux footprint estimations, comparison to the NH3 area emissions inventory provided by the U.S. Environmental Protection Agency indicated overall agreement but also the absence of some sources, for instance the identified cattle farm. Our study concludes that high-precision CIMS measurements are a powerful tool for in situ measurements of ambient NH3 mixing ratios, and even allow for the airborne mapping of the air–surface exchange of NH3.

58 GEOSCIENCES↗
DOE OSTI2022

Airborne flux measurements of ammonia over the Southern Great Plains using chemical ionization mass spectrometry

Ammonia (NH 3 ) is an abundant trace gas in the atmosphere and an important player in atmospheric chemistry, aerosol formation and the atmosphere-surface exchange of nitrogen. It is recognized as a major source of aerosol pollution, and it may limit the formation of cloud nuclei in remote or cold parts of the atmosphere. For soil and plants, NH 3 -mediated nitrogen can act as a harmful pollutant or as a desirable nutrient, mostly in natural and agricultural settings, respectively. Agriculture is also the main source of atmospheric NH 3 via volatilization from fertilizers and manure processing in livestock farming. The accurate determination of NH 3 emission rates remains a challenge, partly due to the propensity of NH 3 to interact with instrument surfaces leading to high detection limits and slow response times. In this paper, we present a new method for quantifying ambient NH 3 , using chemical ionization mass spectrometry (CIMS) with deuterated benzene cations as reagents. The setup aimed at limiting sample-surface interactions and achieved a 1-σ precision of 10–20 pptv and an immediate 1/e response rate < 0.4 s, which compares favorably to the existing state of the art. The sensitivity exhibited an inverse humidity dependence, in particular in relatively dry conditions. Background of up to 10 % of the total signal required consideration as well, as it responded on the order of a few minutes. To showcase the method’s capabilities, we quantified NH3 mixing ratios from measurements obtained during deployment on a Gulfstream I aircraft during the HI-SCALE (Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems) field campaign in rural Oklahoma during May 2016. Typical mixing ratios were 1–10 parts per billion by volume (ppbv) for the boundary layer and 0.1–1 ppbv in the lower free troposphere. Sharp plumes of up to 10s of ppbv of NH3 were encountered as well. We identified two of their sources as a large fertilizer plant and a cattle farm, and our mixing ratio measurements yielded upper bounds of 350 ± 50 and 0.6 kg NH 3 h –1 for their respective momentary source rates. The fast response of the CIMS also allowed us to derive vertical NH 3 fluxes within the turbulent boundary layer via eddy covariance, for which we chiefly used the continuous wavelet transform technique. As expected for a region dominated by agriculture, we observed predominantly upward fluxes, implying net NH 3 emissions from surface. The corresponding analysis focused on the most suitable flight, which contained two straight-and-level legs at ~300 m above ground. We derived NH 3 fluxes between –4 and 18 mol km –2 h –1 for these legs, at an effective spatial resolution of 1–2 km. The analysis demonstrated how flux measurements benefit from suitably arranged flight tracks with sufficiently long straight-and-level legs, and explores the detrimental effect of measurement discontinuities. Following flux footprint estimations, comparison to the NH 3 area emissions inventory provided by the US Environmental Protection Agency indicated overall agreement, but also the absence of some sources, for instance the identified cattle farm. Our study concludes that high-precision CIMS measurements are a powerful tool for in-situ measurements of ambient NH 3 mixing ratios, and even allow for the airborne mapping of the air-surface exchange of NH 3 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

X-ray phase contrast and absorption imaging for the quantification of transient cavitation in high-speed nozzle flows

In this work, high-flux synchrotron radiation has been employed in a time-resolved manner to characterise the distinct topology features and dynamics of different cavitation regimes arising in a throttle orifice with an abrupt flow-entry contraction. Radiographs obtained though both X-ray phase- contrast and absorption imaging have been captured at 67,890 frames per second. The flow lied in the turbulent regime (Re=35,500), while moderate (CN=2.0) to well-established (CN=6.0) cavitation conditions were examined encompassing the cloud and vortical cavitation regimes with pertinent transient features, such as cloud-cavity shedding. X-ray Phase-Contrast Imaging (XPCI), exploiting the shift in the X-ray wave phase during interactions with matter, offers sharp-refractive index gradients in the interface region. Hence, it is suitable for capturing fine morphological fluctuations of transient cavitation structures. Nevertheless, the technique cannot provide information on the quantity of vapour within the orifice. Such data have been obtained utilising absorption imaging, where beam attenuation is not associated with scattering and refraction events, and hence can be explicitly correlated with the projected vapour thickness in line-of-sight measurements. A combination of the two methods is proposed, as it has been found is capable of quantifying the vapour content arising in the complex nozzle flow, while also faithfully illustrating the dynamics of the highly-transient cavitation features.

42 ENGINEERING↗
DOE OSTI2025

Vertical column dual-comb spectroscopy to a TBS

Open-path dual-frequency-comb spectroscopy (DCS) is a broadband, high spectral resolution, and high precision method for measuring gas concentrations over kilometer-scale paths. It has been used for detection and quantification of emissions of pollutants, hazardous gasses, and greenhouse gasses (GHGs). DCS has been shown to measure trace-gas mixing ratios with 0.14%-0.4% agreement between instruments. To achieve high signal-to-noise ratios (SNRs) with DCS, comb light is targeted onto a retroreflector at the end of the measurement path, which returns the signal to a detector co-located with the launch signal. Installing the retroreflector on a mobile platform such as a balloon or unmanned aerial vehicle (UAV) extends the capabilities of DCS by enabling variable path lengths, greater mobility, and access to higher altitudes. Mobile-target DCS has many uses in plume and leak detection, emissions modeling, and planetary boundary layer (PBL) studies. It is a promising method for observing vertical distributions of GHGs and mixing processes in the PBL, which are difficult to measure but important for pollution and climate monitoring as well as for understanding transport of gasses through the atmosphere. Tethered balloons are an intriguing platform because they enable longer flight durations and higher altitudes than easily obtainable with a UAV and thus allow for column measurements up to and above the PBL. New measurements completed in October/November 2024 reach the highest altitudes above ground level yet achieved by mobile-target DCS. For these measurements, the retroreflector was mounted on a 7 m diameter tethered helium balloon. An actively tracking gimbal on the ground holds the DCS launch telescope and keeps the 5 cm beam pointed onto the retroreflector while the balloon is lifted, lowered, and moved by wind and turbulence.

atmosphere↗
DOE OSTI2022

Surface Heater Fabrication Using Micro-Lithography for Transpiration Cooling Heat Transfer Coefficient Measurements

Abstract Transpiration cooling is able to provide more uniform coolant coverage than film cooling to effectively protect the component surface from contacting the hot gas. Due to numerous coolant ejection outlets within a small area at the target surface, the experimental thermo-fluid investigation on transpiration cooling becomes a significant challenge. Two classic methods to investigate film cooling, the steady-state foil heater method and the transient thermography technique, fail for transpiration cooling because the foil heater would block numerous coolant outlets, and the semi-infinite solid conduction model no longer holds for porous plates. In this study, a micro-lithography method to fabricate a silver coil pattern on top of the additively manufactured polymer porous media as the surface heater was proposed. The circuit was deliberately designed to cover the solid surface in a combination of series connection, and parallel connection to ensure the power in each unit cell area at the target surface was identical. With uniform heat flux generation, the steady-state tests were conducted to obtain distributions of a pair of parameters, adiabatic cooling effectiveness, and heat transfer coefficient (HTC). The results showed that the adiabatic cooling effectiveness could reach 0.65 with a blowing ratio lower than 0.5. Meanwhile, the heat transfer coefficient ratio (hf/h0) of transpiration cooling was close to 1 with a small blowing ratio at 0.125. A higher HTC ratio was observed for smaller pitch-to-diameter cases due to more turbulence intensity generated at the target surface.

Engineering↗