SEARCH · Search NASA

Results for “BALANCE EQUATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗

DOE OSTI2023

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Triboplasma assisted chemical conversion in granular systems: a semi-analytic model

We present a semi-analytic model for a novel plasma-assisted chemical conversion pathway using triboplasmas generated in granular flows. Triboelectric charge relaxation is a well known phenomena where the potential generated from contact charging of particles exceeds the breakdown voltage of the background gas. In this work, we extend the triboelectric charge relaxation theory to include non equilibrium plasma energy and particle balance equations to predict the formation of dissociated and excited species that act as precursors to chemical conversion, for example in plasma-assisted ammonia synthesis. Our example case study with nitrogen background gas and teflon/aluminum tribomaterial system yielded high excited nitrogen species densities per collision that are comparable to current plasma-assisted conversion pathways. We also present a regime diagram for various gases where Paschen breakdown parameters are used to determine whether triboplasmas can be formed for a given effective work-function difference between two materials. Our sensitivity studies indicate particle velocity, particle radius, solids fraction and space charge effects play a critical role in overall plasma densities and excited species production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge

The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Power Flow Solution Methods for Combined Transmission and Distribution Systems

High distributed PV penetration makes the distribution systems more active and increases their impact on the upstream transmission grid. To facilitate the studies on this area, this paper proposes a power flow (PF) model and formulation of combined transmission and distribution (T&D) power systems. The formulation is based on the physical coupling of transmission and distribution grids. The unbalanced nature of distribution feeders and voltage regulators are modeled so that the resulting PF balance equations in distribution systems are similar to those in the interconnected transmission system. The formulated nonlinear PF equations are solved using the proposed unified and sequential methods. Both methods are verified and compared with each other using several combined T&D power systems. A multi-period PF simulation is also included to demonstrate the capability of the proposed formulation and solution methods in capturing the influences of high distributed PV penetration on transmission level.

Nguyen, Quan H.↗

DOE OSTI2020

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Mathematical Model for Li-S Cell with Shuttling-Induced Capacity Loss Approximation

Lithium sulfur (Li-S) batteries have the potential to outperform the current lithium ion batteries and transform the technology of the future. However, dissolution, diffusion, and shuttling of the dissolved polysulfides result in parasitic reactions and substantial capacity loss. To provide a better understanding of the shuttling process, a 1D porous electrode mathematical model has been developed in this paper. An approximation method is used to account for the shuttling-induced capacity loss by adding an extra source/sink term in the material balance equations for the species involved in the parasitic reactions. Shuttling constants used in the source terms can be determined by fitting the model predictions to the experimental measurements. Here, the results showed that by including the approximation method, the model was able to predict the active material loss and the continuous decrease of volume fractions of Li 2 S on the cathode surface. The model sheds light on the capacity loss mechanism occurring inside the cell as a result of the shuttling of polysulfides.

25 ENERGY STORAGE↗

DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗

DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗

DOE OSTI2023

Turbulent Kinetic Energy Budgets over Gentle Topography Covered by Forests

Large-eddy simulations of flow over a “horizontally” uniform model forest are used to investigate the effects of gentle topography on the turbulent kinetic energy (TKE) budget within the canopy roughness sublayer. Despite significant differences between simulations using idealized sinusoidal topography and real topography of the Amazon forest, results indicate that the effects of topography are located predominantly in the upper canopy and above, and are mostly caused by mean advection of TKE. The “horizontally” averaged TKE budget from idealized and real gentle topographies are almost identical to that for flat terrain, including a clear inertial layer above the roughness sublayer in which shear production is balanced by local dissipation. At topography crests, where observational towers are usually located, mean vertical advection of TKE can be as important as horizontal advection. Here, we propose the use on an approximate TKE balance equation to estimate mean advection from single tower measurements, and introduce a new advection index that can be used as a proxy to quantify the importance of the topography on the TKE budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

On the Diffusivity of Moist Static Energy and Implications for the Polar Amplification Response to Climate Warming

Energy balance models (EBMs) have been widely used in a range of climate problems, but the assumption of constant diffusivity in the parameterization of the moist static energy (MSE) flux can be hardly justified. We demonstrate in this study that the diffusive MSE flux can be derived from the basic energy balance equation with a few tolerable assumptions. The estimated diffusivity is both spatially and seasonally dependent, and its midlatitude average is then tested against several scaling theories for the midlatitude eddy diffusivity. The result supports the diffusivity theory of Held and Larichev (1996) modified for the moist atmosphere, affording a dynamics-based parameterization of MSE diffusivity. The implementation of the parameterization in an EBM leads to an interactive MSE diffusivity that accounts for the midlatitude eddy response to climate forcing perturbations. Under a uniform radiative forcing, the EBM with a diffusivity so parameterized produces a weakening of the midlatitude diffusivity and a modestly polar-amplified surface temperature response as an inevitable outcome under the dual constraints of the nonlinear Clausius-Clapeyron relation and the temperature gradient-dependent diffusivity, even in the absence of any poleward amplifying radiative feedbacks. As the consequence of more isothermal temperature and reduced diffusivity, the variance of the midlatitude surface temperature also decreases with warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Plasma Emission, Wave Speeds, and Temperatures of Laser-Supported Plasmas in Air

A study of the generation and subsequent expansion of plasmas using a single beam (532 nm) from the second harmonic generator of a Q-switched Nd:YAG laser is described. The plasmas in air are studied using an intensified charge-coupled device camera over periods of 1000 ns following the laser pulses. The plasma experimental results are interpreted in terms of wave speeds and plasma temperatures. Energy balance equations are considered to estimate the plasma temperatures and the laser-supported radiation (LSR) and the laser-supported detonation (LSD) wave speeds. Within 1000 ns observation time, the circular spectral images turn elliptical because of plasma expansion toward the laser beam. The plasma expansion increases with increasing laser irradiance. A laser irradiance range of [Formula: see text] yields LSR wave speeds between 45 and 170 km/s and LSD wave speeds of 90–130 km/s. The measured wave speeds are compared with the LSR and the LSR wave speed models. The plasma temperatures are between 90,000 and 120,000 K at 15 ns after the laser pulse based on measured emission spectra and are between 640,000 and 760,000 K at 3 ns based on the LSR model.

Engineering↗