A high voltage, low current power source for long term space applications
Experimental design of high voltage low current power source for unmanned, long term space mission
Grygiel, M. L.↗
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.
Colloid-Facilitated Actinide Transport in a Cementitious-Impacted SRS Groundwater
The objective of this report is to provide a literature review relevant to the question of whether colloid facilitated transport of plutonium is enhanced in a cementitious groundwater system on the Savannah River Site (SRS). Contaminant transport is commonly described as taking place in a system with a mobile aqueous phase and an immobile solid phase. There has been an increasing awareness of a third phase, a mobile solid phase, also referred to as a mobile colloidal phase. Mobile colloids consist of organic and/or inorganic submicron-particles that move with groundwater flow. When radionuclides are associated with the mobile colloids, the net effect is that radionuclides can move faster through the subsurface system than would be predicted by transport models that do not include mobile colloids. It is important to distinguish between subsurface colloids and subsurface mobile colloids. The subsurface environment includes an enormous reservoir of colloids, but only a tiny fraction, if any, are mobile. Mobile colloid formation is commonly described as involving a three-step process: genesis, stabilization, and transport. It was concluded that there is a strong likelihood that submicron colloids of plutonium exist near the source term, (i.e., the genesis step is likely completed). While such particles have not been directly detected, it is highly likely that plutonium could either attach to submicron particles in the sediment or that submicron plutonium fragments exist at the source. However, the tendency for these plutonium colloids to move in the SRS subsurface is extremely low, especially in engineered cementitious environments where the ionic strength of the solution and the elevated concentrations of divalent cations greatly curtails colloid suspension stability (i.e., the stabilization step is likely not completed). Together these data strongly indicate that plutonium-bearing colloids would likely exist in the near field of a tank closure facility, but the colloids would be unlikely to be mobile, thereby providing a vector for enhanced transport. Under cementitious leachate impacted groundwater conditions, it is reasonable to assume plutonium transport occurs primarily as a two-phase system, a mobile aqueous phase and an immobile solid phase.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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.↗
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.↗
Radionuclide Screening Analysis and Transport Parameters for Pahute Mesa Detonations, Nevada National Security Site
A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
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.↗
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↗
Local wave-number model for inhomogeneous two-fluid mixing
We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
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.↗
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.↗
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.↗
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.↗
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↗
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.↗
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.↗
Turbulence-parameter estimation for current-energy converters using surrogate model optimization
Surrogate models maximize information utility by building predictive models in place of computational or experimentally expensive model runs. Marine hydrokinetic current energy converters require large-domain simulations to estimate array efficiencies and environmental impacts. Meso-scale models typically represent turbines as actuator discs that act as momentum sinks and sources of turbulence and its dissipation. An OpenFOAM model was developed where actuator disc k-ε turbulence was characterized using an approach developed for flows through vegetative canopies. Turbine-wake data from laboratory flume experiments collected at two influent turbulence intensities were used to calibrate parameters in the turbulence-source terms in the k-ε equations. Additionally, parameter influences on longitudinal wake profiles were estimated using Gaussian process regression with subsequent optimization minimizing the objective function within 3.1% of those obtained using the full model representation, but for 74% of the computational cost (far fewer model runs). This framework facilitates more efficient parameterization of the turbulence-source equations using turbine-wake data.
16 TIDAL AND WAVE POWER↗
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.↗