Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source term”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Land Water Storage within the Congo Basin Inferred from GRACE Satellite Gravity Data

GRACE satellite gravity data is used to estimate terrestrial (surface plus ground) water storage within the Congo Basin in Africa for the period of April, 2002 - May, 2006. These estimates exhibit significant seasonal (30 +/- 6 mm of equivalent water thickness) and long-term trends, the latter yielding a total loss of approximately 280 km(exp 3) of water over the 50-month span of data. We also combine GRACE and precipitation data set (CMAP, TRMM) to explore the relative contributions of the source term to the seasonal hydrological balance within the Congo Basin. We find that the seasonal water storage tends to saturate for anomalies greater than 30-44 mm of equivalent water thickness. Furthermore, precipitation contributed roughly three times the peak water storage after anomalously rainy seasons, in early 2003 and 2005, implying an approximately 60-70% loss from runoff and evapotranspiration. Finally, a comparison of residual land water storage (monthly estimates minus best-fitting trends) in the Congo and Amazon Basins shows an anticorrelation, in agreement with the 'see-saw' variability inferred by others from runoff data.

Crowley, John W.↗
NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗
NASA NTRS1985

Interpretation of ozone temperature correlations. II - Analysis of SBUV ozone data

In an effort to understand the behavior of ozone in the transition region, the terms in the perturbation continuity equation have been evaluated from solar backscattered ultraviolet ozone data, National Meteorological Center temperature data, and geostrophic winds derived from the temperature data at 2 mbar, winter 1978-1979. A linear parameterization is used to estimate the chemical source term. The perturbation ozone field calculated by integrating the net forcing evaluated from the data compares well with the measured perturbation ozone field. There are occasions during the winter when the sum of the dynamic terms is small due to cancellation, and the ozone and temperature perturbations are anticorrelated. The -1 correlation between the ozone and temperature perturbations does not necessarily imply that chemical forcing is dominant.

Douglass, A. R.↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Integral Velocity Sampling for Unsteady Rotor Models on Cartesian Meshes

An unsteady body-force rotor model utilizing integral velocity sampling is developed for Cartesian meshes with a focus on accurate performance predictions. In this approach, the rotors are modeled as body-force source terms, and integral sampling is used to define the local freestream velocity vector. The method is verified in two dimensions for airfoils, providing accurate angle of attack and freestream Mach number estimates even on relatively coarse meshes, and demonstrating excellent mesh convergence across a variety of source region thicknesses. Extensions of the model to three dimensions for both fixed and rotating wings are considered, with numerical experiments demonstrating the need for a spanwise-continuous representation of the induced velocity or angle of attack. The model provides reasonable spanwise lift distributions without any explicit tip-loss model. Comparisons to high-fidelity simulations and experimental data show excellent agreement for predictions of an isolated DJI Phantom 3 propeller in hover. Finally, simulations of the full aircraft in edgewise forward flight also compare well with wind tunnel data across a variety of pitch angles and propeller speeds.

ARMD↗
NASA NTRS1979

Scattering and distortion of the unsteady motion on transversely sheared mean flows

It is shown that the pressure and velocity fluctuations of the unsteady motion on a transversely sheared mean flow can be expressed entirely in terms of the derivatives of two potential functions. One of these is a convected quantity that can be specified as a boundary condition and is related to a transverse component of the upstream velocity field. The other can be determined by solving an inhomogeneous wave equation whose source term is also a convected quantity that can be specified as a boundary condition in any given problem. The general theory is used to study the interaction of an unsteady flow with a semi-infinite plate embedded in a shear layer. The acoustic field produced by this interaction is calculated in the limits of low and high frequency. The results are compared with experimental one-third octave sound pressure level radiation patterns. The agreement is found to be excellent, especially in the low frequency range, where the mean-flow and convective effects are shown to have a strong influence on the directivity of the sound.

Goldstein, M. E.↗
NASA NTRS1993

Examination of nonequilibrium effects in an ionized nitrogen flow

A thermo-chemical nonequilibrium ionized nitrogen flow in a shock tube is examined. A one-dimensional computational fluid dynamics code has been developed to study the flowfield incorporating a novel treatment of the electron energy equation. In the present approach, the electron pressure is included in the electron energy flux. In contrast, previous work has removed this term from the flux to facilitate the flux splitting. The code uses a five-species gas model for ionized nitrogen characterized by translational-rotational, vibrational, and electron-electronic temperatures. The results give good agreement with experimental data except for an inadequacy in the electronic source terms. A small difference is observed between the new and old splittings which increases with higher levels of ionization.

Hatfield, John A.↗
NASA NTRS2003

Various Numerical Applications on Tropical Convective Systems Using a Cloud Resolving Model

In recent years, increasing attention has been given to cloud resolving models (CRMs or cloud ensemble models-CEMs) for their ability to simulate the radiative-convective system, which plays a significant role in determining the regional heat and moisture budgets in the Tropics. The growing popularity of CRM usage can be credited to its inclusion of crucial and physically relatively realistic features such as explicit cloud-scale dynamics, sophisticated microphysical processes, and explicit cloud-radiation interaction. On the other hand, impacts of the environmental conditions (for example, the large-scale wind fields, heat and moisture advections as well as sea surface temperature) on the convective system can also be plausibly investigated using the CRMs with imposed explicit forcing. In this paper, by basically using a Goddard Cumulus Ensemble (GCE) model, three different studies on tropical convective systems are briefly presented. Each of these studies serves a different goal as well as uses a different approach. In the first study, which uses more of an idealized approach, the respective impacts of the large-scale horizontal wind shear and surface fluxes on the modeled tropical quasi-equilibrium states of temperature and water vapor are examined. In this 2-D study, the imposed large-scale horizontal wind shear is ideally either nudged (wind shear maintained strong) or mixed (wind shear weakened), while the minimum surface wind speed used for computing surface fluxes varies among various numerical experiments. For the second study, a handful of real tropical episodes (TRMM Kwajalein Experiment - KWAJEX, 1999; TRMM South China Sea Monsoon Experiment - SCSMEX, 1998) have been simulated such that several major atmospheric characteristics such as the rainfall amount and its associated stratiform contribution, the Qlheat and Q2/moisture budgets are investigated. In this study, the observed large-scale heat and moisture advections are continuously applied to the 2-D model. The modeled cloud generated from such an approach is termed continuously forced convection or continuous large-scale forced convection. A third study, which focuses on the respective impact of atmospheric components on upper Ocean heat and salt budgets, will be presented in the end. Unlike the two previous 2-D studies, this study employs the 3-D GCE-simulated diabatic source terms (using TOGA COARE observations) - radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the Ocean mixed-layer (OML) model.

Shie, C.-L.↗
NASA NTRS2017

Identification of Spurious Signals from Permeable Ffowcs Williams and Hawkings Surfaces

Integral forms of the permeable surface formulation of the Ffowcs Williams and Hawkings (FW-H) equation often require an input in the form of a near field Computational Fluid Dynamics (CFD) solution to predict noise in the near or far field from various types of geometries. The FW-H equation involves three source terms; two surface terms (monopole and dipole) and a volume term (quadrupole). Many solutions to the FW-H equation, such as several of Farassat's formulations, neglect the quadrupole term. Neglecting the quadrupole term in permeable surface formulations leads to inaccuracies called spurious signals. This paper explores the concept of spurious signals, explains how they are generated by specifying the acoustic and hydrodynamic surface properties individually, and provides methods to determine their presence, regardless of whether a correction algorithm is employed. A potential approach based on the equivalent sources method (ESM) and the sensitivity of Formulation 1A (Formulation S1A) is also discussed for the removal of spurious signals.

Lopes, Leonard V.↗
DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗
NASA NTRS2004

High-Resolution Simulation of Hurricane Bonnie (1998): Water Budget - Part 2

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, Scott A.↗
NASA NTRS2005

The Water Budget of a Simulated Hurricane

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weaker updrafts, indicative of a larger role of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where cooling from melting enhances condensation. Unlike in squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance in situ stratiform precipitation development through additional growth by vapor deposition and aggregation. We examine artificial source terms for cloud and precipitation mass associated with setting to zero negative mixing ratios that arise from numerical advection errors. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation and 15-20% of the precipitation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, S.↗
NASA NTRS2001

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of surface evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote geographic sources of surface evaporation for precipitation, based on the implementation of three-dimensional constituent tracers of regional water vapor sources (termed water vapor tracers, WVT) in a general circulation model. The major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In the WVT approach, each tracer is associated with an evaporative source region for a prognostic three-dimensional variable that represents a partial amount of the total atmospheric water vapor. The physical processes that act on a WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be predicted within the model simulation, and can be validated against the model's prognostic water vapor. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional sources, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In most North American continental regions, the local source of precipitation is correlated with total precipitation. There is a general positive correlation between local evaporation and local precipitation, but it can be weaker because large evaporation can occur when precipitation is inhibited. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.↗
DOE OSTI2026

A Proof of the Asymptotic Variance of Path Length Estimators for Single-Collision Monte Carlo Source Iteration in the Thick Diffusion Limit

Here, we prove a theorem relating the variance of path length estimators for single-collision Monte Carlo source iteration to a parameter that becomes infinitesimally small in an important physical regime arising in radiative transfer. In our usage, “single-collision Monte Carlo source iteration” refers to Monte Carlo Boltzmann transport methods in which each Monte Carlo particle history includes no more than a single collision, and the physics of multiple scattering is modeled by lagging the scattering source term and iterating until this term converges. Our theorem can be used to construct variance reduction techniques which improve the order of the estimator variance. This enables calculations that would otherwise require impractically large sample sizes to achieve practical estimator uncertainties. We believe this is the first postulation of a theorem relating estimator variance to a limiting case parameter for single-collision Monte Carlo source iteration, and the first proof of such a theorem. We illustrate the theorem’s value with an example in which the authors of a transport method used the theorem to design a variance reduction technique that improved the uncertainty of their solution by a factor of about 500 for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material.

Mathematics and Computing↗
DOE OSTI2025

Effect of Time Window and Spectral Measurement Options on Empirical Green’s Function Analysis Using DAS Array and Seismic Stations

The recorded seismic waveform is a convolution of event source term, path term, and station term. Removing high-frequency attenuation due to path effect is a challenging problem. Empirical Green’s function (EGF) method uses nearly collocated small earthquakes to correct the path and station terms for larger events recorded at the same station. However, this method is subject to variability due to many factors. Here, we focus on three events that were well recorded by the seismic network and a rapid response distributed acoustic sensing (DAS) array. Using a suite of high-quality EGF events, we assess the influence of time window, spectral measurement options, and types of data on the spectral ratio and relative source time function (RSTF) results. Increased number of tapers (from 2 to 16) tends to increase the measured corner frequency and reduce the source complexity. Extended long time window (e.g., 30 s) tends to produce larger variability of corner frequency. The multitaper algorithm that simultaneously optimizes both target and EGF spectra produces the most stable corner-frequency measurements. The stacked spectral ratio and RSTF from the DAS array are more stable than two nearby seismic stations, and are comparable to stacked results from the seismic network, suggesting that DAS array has strong potential in source characterization.

58 GEOSCIENCES↗
NASA NTRS2021

Preventing Anomalous Torques in Circumbinary Accretion Simulations

Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms that prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how improved sink methods reduce some of these biases, using simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. We show that sink terms that do not reduce the angular momentum of gas relative to the accreting object (1) reduce the dependence on the sink rate of quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semimajor axis; (2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; (3) alter the variability of accretion onto the binary, regularizing it in time. We also investigate other potential sources of systematic error, such as gravitational softening and simplifications to the viscous stress tensor. Because of the strong dependence of binary orbital evolution on both the torque and the distribution of mass between binary components, the sink method influences the orbital evolution of the binary at all mass ratios, with mass ratios below ∼0.3 most affected.

Alexander Dittmann↗
NASA NTRS1996

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗
DOE OSTI2024

Parameterization of vacancy production rate in phase-field models of fission gas bubble evolution in nuclear fuel

Phase-field modeling has increasingly been used to study microstructural evolution in fission gas bubbles in nuclear fuel to improve understanding of fission gas release. To improve computational efficiency, often only vacancies and gas atoms are included as defect species. In this case, the net effects of vacancy and interstitial production, recombination, and biased sink absorption are included as a net vacancy source, or net vacancy source combined with an effective sink. However, there has been a lack of clarity on what parameter values should be used for these approaches to best match the more complete physical picture that includes interstitials and vacancies. Here, we compare a phase-field model of void growth to analytical models for the source-only and source plus sink approach to gain insight into how the phase-field models can be parameterized effectively. The source-only approach provides greater flexibility to match growth rates determined from the full vacancy-interstitial picture. A strategy was developed for determining the value of the net vacancy source term by comparing to an analytical model that includes vacancy and interstitial production, recombination, and biased sink absorption. Finally, this strategy can be used to parameterize phase-field models of fission gas bubble growth.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗