Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 1,261 records · Page 70

Numerical simulation of underexpanded plumes using upwind algorithms

The solution of the flow characteristics of highly underexpanded plumes is presently sought by means of a high-resolution total variation-diminishing upwind scheme featuring finite-volume metrics, Roe averaging with entropy correction, and monotonic upstream schemes for conservation laws. The existing approximate factorization schemes are used both as preconditioners and in order to apply the conjugate residual algorithm. The plumes treated are axisymmetric and laminar, and for an ideal gas that is allowed to be thermally incomplete. Simulations of the SSME and SRB plumes have been successfully conducted.

Obayashi, Shigeru↗

A computational study on the interaction between a vortex and a shock wave

A computational study of two-dimensional shock vortex interaction is discussed in this paper. A second order upwind finite volume method is used to solve the Euler equations in conservation form. In this method, the shock wave is captured rather than fitted so that the cases where shock vortex interaction may cause secondary shocks can also be investigated. The effects of vortex strength on the computed flow and acoustic field generated by the interaction are qualitatively evaluated.

Meadows, Kristine R.↗

A detailed examination of a X-line region in the distant tail: ISEE-3 observations of jet flow and B(sub z) reversals and a pair of slow shocks

We report an observation of Petschek-type magnetic reconnection at a distant neutral line (X = -230 R(sub e)) with a full set of signatures of the magnetic merging process. These features include a reversal of plasma flows from earthward to tailward, a pair of slow shocks and the magnetic field X-type line. These two slow shocks are shown to satisfy the shock criteria used by Feldman et al. (1987). The spacecraft first crosses a slow shock to enter the earthward flowing plasmasheet with velocity of about 440 km/s. The embedded magnetic field has a positive B(sub z) component. The spacecraft next enters a region of tailward plasma flow with speed approximately 670 km/s and an embedded negative B(sub z), indicating entry into the plasmasheet tailward of the X-line. These observed velocities are comparable to calculated velocities based on Rankine-Hugoniot conservation relationships. The spacecraft subsequently returns into the south tail lobe by crossing another slow shock. Coplanarity analyses shows that the two slow shocks have orientations consistent with that predicted by the Petschek reconnection model. We note that this event occurs during northward interplanetary magnetic fields. Thus, a magnetic stress built-up in the distant tail may be responsible for this reconnection process.

Ho, C. M.↗

deal.II Implementation of a Weak Galerkin Finite Element Solver for Darcy Flow

This paper presents a weak Galerkin (WG) finite element solver for Darcy flow and its implementation on the deal.II platform. The solver works for quadrilateral and hexahedral meshes in a unified way. It approximates pressure by Q-type degree $k(≥ 0)$ polynomials separately defined in element interiors and on edges/faces. In this work, numerical velocity is obtained in the unmapped Raviart-Thomas space $RT_{[k]}$ via postprocessing based on the novel concepts of discrete weak gradients. The solver is locally mass-conservative and produces continuous normal fluxes. The implementation in deal.II allows polynomial degrees up to 5. Numerical experiments show that our new WG solver performs better than the classical mixed finite element methods.

54 ENVIRONMENTAL SCIENCES↗

Development of a Coupled Hydrological/Sediment Yield Model for a Watershed at Regional Level

Development of a hydrologic model for the study of environmental conservation requires a comprehensive understanding of individual-storm affecting hydrologic and sedimentologic processes. The hydrologic models that we are currently coupling are the Simulator for Hydrology and Energy Exchange at the Land Surface (SHEELS) and the Distributed Runoff Model (DRUM). SHEELS runs continuously to estimate surface energy fluxes and sub-surface soil water fluxes, while DRUM operates during and following precipitation events to predict surface runoff and peak flow through channel routing. The lateral re-distribution of surface water determined by DRUM is passed to SHEELS, which then adjusts soil water contents throughout the profile. The model SHEELS is well documented in Smith et al. (1993) and Laymen and Crosson (1995). The model DRUM is well documented in Vieux et al. (1990) and Vieux and Gauer (1994). The coupled hydrologic model, SHEELS/DRUM, does not simulate sedimentologic processes. The simulation of the sedimentologic process is important for environmental conservation planning and management. Therefore, we attempted to develop a conceptual frame work for coupling a sediment yield model with SHEELS/DRUM to estimate individual-storm sediment yield from a watershed at a regional level. The sediment yield model that will be used for this study is the Universal Soil Loss Equation (USLE) with some modifications to enable the model to predict individual-storm sediment yield. The predicted sediment yield does not include wind erosion and erosion caused by irrigation and snow melt. Units used for this study are those given by Foster et al. (1981) for SI units.

Rajbhandaril, Narayan↗

Scale-Up of Electrode Coating and Flow-Field for Commercial Hydrogen Peroxide Electrolyzer: Cooperative Research and Development Final Report, CRADA Number CRD-17-00687

Hydrogen peroxide is currently produced at central chemical plants via the anthraquinone oxidation process. This process produces environmental pollutants that are costly to remediate, requires hazardous long distance shipping of highly concentrated peroxide (50% or 70%), and necessitates extra handling costs related to storage and dilution. Peroxygen Systems, Inc. (PSi) is developing breakthrough technology for on-site hydrogen peroxide production. PSi’s on-site on-demand electrolyzer can reduce the cost of producing hydrogen peroxide by 50%, while also completely eliminating the cost and safety issues associated with shipping and handling of high concentration hydrogen peroxide. The challenge for PSi is scaling. To support the next step toward commercialization (customer pilot tests), scaling the prototype into larger single cells and 20-40 cell stacks is required. In addition to internal hardware and flow-field design efforts at PSi, NREL will address three critical problems for this scale-up effort: (1) demonstrating a large scale roll-to-roll (R2R) process to coat uniform electrode materials for 100 cm2 and 500 cm2 stack testing, (2) demonstrating an in-line diagnostic to achieve better electrode quality control, and (3) performing in situ cell/stack testing to better understand and optimize the performance of the flow field design.

28 EE - Advanced Manufacturing Office (EE-5A)↗

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optimal Power Dispatch of DGs in Radial and Mesh AC Grids: A Hybrid Solution Methodology between the Salps Swarm Algorithm and Successive Approximation Power Flow Method

In this paper, we address the problem of the optimal power dispatch of Distributed Generators (DGs) in Alternating Current (AC) networks, better known as the Optimal Power Flow (OPF) problem. We used, as the objective function, the minimization of power losses (P loss ) associated with energy transport, which are subject to the set of constraints that compose AC networks in an environment of distributed generation. To validate the effectiveness of the proposed methodology in solving the OPF problem in any network topology, we employed one 10-node mesh test system and three radial text systems: 10, 33, and 69 nodes. In each test system, DGs were allowed to inject 20%, 40%, and 60% of the power supplied by the slack generator in the base case. To solve the OPF problem, we used a master–slave methodology that integrates the optimization method Salps Swarm Algorithm (SSA) and the load flow technique based on the Successive Approximation (SA) method. Moreover, for comparison purposes, we employed some of the algorithms reported in the specialized literature to solve the OPF problem (the continuous genetic algorithm, the particle swarm optimization algorithm, the black hole algorithm, the antlion optimization algorithm, and the Multi-Verse Optimizer algorithm), which were selected because of their excellent results in solving such problems. The results obtained by the proposed solution methodology demonstrate its superiority and convergence capacity in terms of minimization of P loss in both radial and mesh systems. It provided the best reduction in minimum P loss in short processing times and showed excellent repeatability in each test system and scenario under analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The attractor flow for AdS 5 black holes in $\mathcal{N}$ = 2 gauged supergravity

We study the flow equations for BPS black holes in $\mathcal{N}$ = 2 five-dimensional gauged supergravity coupled to any number of vector multiplets via FI couplings. We develop the Noether-Wald procedure in this context and exhibit the conserved charges as explicit integrals of motion, in the sense that they can be computed at any radius on the rotating spacetime. The boundary conditions needed to solve the first order differential equations are discussed in great detail. We extremize the entropy function that controls the near horizon geometry and give explicit formulae for all geometric variables at their supersymmetric extrema. We have also considered a complexification of the near-horizon variables that elucidates some features of the theory from the near-horizon perspective.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gauge invariance of the helicity continuity equation

The derivation of the helicity continuity equation in electromagnetic theory is performed without specifying a gauge. In contrast with previous proposals, the form of this equation is shown to be gauge invariant without invoking a Helmholtz decomposition. The helicity and its flow, the latter associated with the spin in quantized fields, involve two sets of a vector and a scalar potential, where each set can independently undergo a gauge transformation. There are alternative definitions of the helicity and flow densities that arise from different grouping of terms in the continuity differential equation. The various definitions acquire an unambiguous meaning, depending on the gauge and the physical context. The helicity density, defined as ϱ{sub A{sub C}{sup (2)}} : =με(A·B-C·E) and flow density J{sub AC}{sup (2)} : =με(E-∇ϕ{sub A})×A+(B-∇ϕ{sub C})×C, include all the rotational content of the free fields regardless of the gauge. In free space, these quantities satisfy a gauge invariant conservation equation without gauge-fixing source terms. A further asset of the present formulation is that charge and current source terms can be readily incorporated. The helicity source terms are of the form μB·∫Jdt-μ∫Bdt·J. A helicity continuity equation in terms solely of transverse fields is derived in the Coulomb gauge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessing Nutrient and Carbon Responses to Agricultural Conservation Practices in Two Midwest Watersheds

The United States is undertaking efforts to transition to a low-carbon economy in response to the heightened impacts of climate change, which are largely attributed to decades of carbon-intensive development. A concerted effort by the energy sector to achieve net-zero carbon emissions is aimed at strategically reducing the carbon footprint of the energy supply chain, particularly in the production of feedstocks for biofuels. The production of biofuels relies heavily on land use and management practices in biomass cultivation. Among other factors, soil organic carbon (SOC) plays a crucial role in the biofuel carbon cycle, impacting land productivity, greenhouse gas emissions, and water quality. The process of carbon drawdown during plant growth and storage helps to mitigate the release of carbon dioxide (CO 2 ) into the atmosphere. Another key parameter is the release of nitrous oxide (N 2 O) — a potent greenhouse gas with a global warming potential 273 times that of CO 2 — from soil. While prioritizing low-carbon production, sustainable bioenergy also necessitates improved water quality and ecosystem services. Agricultural conservation practices can contribute to enhanced soil health and reduce nutrient and sediment loss into streams. However, studies that explore the broader impact of these practices on soil carbon storage and N 2 O emissions at a watershed scale are limited. Specifically, our understanding of how low-carbon feedstock production and management affect nutrient cycle dynamics is incomplete. This study assesses watershed responses to land management practices in two agriculturally dominant watersheds: the Raccoon River watershed and the Southfork of Iowa River watershed. The study focuses on carbon and nutrient dynamics, characterizing spatial and temporal variations in SOC, N 2 O, and nutrient loadings to examine the relationship between nutrient and carbon responses and agricultural conservation practices. The study employs the newly developed Soil Water Analysis Tool for Carbon (SWAT-C) model, calibrated using 20 years’ of climate and water monitoring data, to simulate and evaluate two agricultural conservation practices: no-till and crop residue harvest with cover crop planting. We compared the calibrated SWAT-C model with a historical baseline model, which allowed us to analyze various aspects of the watershed, including stream flow, suspended sediments, nitrogen, phosphorus, organic carbon, SOC at different depths, and N 2 O emissions from topsoil. Our SWAT-C modeling results indicate that, compared with results from the historical baseline model, no-till practices are positively correlated with SOC accumulation, reduced N 2 O emissions, and decreased soil loss. We also observed no change or slightly increased nitrogen and phosphorus loss to water bodies in the watersheds compared with the baseline. Conversely, crop residue harvest with cover crop planting improved water quality by reducing nutrient and soil losses but also increased N 2 O emissions. This SWAT-C modeling study establishes the groundwork for further smaller-scale and/or sub-basin-level assessments of nitrogen, phosphorus, and carbon cycles. It provides valuable science-based insights to inform policy decisions — particularly in the context of transitioning to biofuel feedstock production in these watersheds — by addressing concerns about nutrient exports to local water bodies and, ultimately, the Mississippi River.

54 ENVIRONMENTAL SCIENCES↗

Three-dimensional, shock-on-shock interaction problem

The unsteady, three-dimensional flowfield resulting from the interaction of a plane shock with a cone-shaped vehicle traveling supersonically is determined, using a second-order, shock-capturing, finite-difference approach. The time-dependent, inviscid gasdynamic equations are transformed to include the self-similar property of the flow, to align various coordinate surfaces with known shock waves, and to cluster points in the vicinity of the intersection of the transmitted incident shock and the surface of the vehicle. The governing partial differential equations in conservation-law form are then solved iteratively using MacCormack's (1969) algorithm.

Kutler, P.↗

Two Algorithms For Hypersonic Computations

Two new algorithms developed containing several desirable features for computation of viscous, hypersonic flows about three-dimensional geometries. Employ upwind-differencing methods, having unsurpassed abiltiy to capture flow-field discontinuities without need for any user-specified smoothing terms. Particularly important in hypersonic computations where many complicated wave structures exist. Both algorithms based on finite-volume formulation to ensure schemes, including boundary conditions, fully conservative.

Molvik, G. A.↗

Sensing and actuation in a liquid propellant rocket engine

Equations describing the conservation of mass, momentum, and energy, and multicomponent transport and chemical reaction in a liquid propellant rocket (LPR) engine are given, along with appropriate boundary conditions. The formulation includes a detailed description of the two-phase (gas/liquid) flow in terms of a moving boundary problem, with droplets and vapor occupying time-dependent subdomains of the LPR chamber. Several sensor and actuator types are considered in detail, and their mathematical descriptions are incorporated into the differential equations and boundary conditions.

Bentsman, Joseph↗