Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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,045 records · Page 58

Coupling reduced‐order modeling and coarse‐grained CFD‐DEM to accelerate coal gasifier simulation and optimization

Abstract High fidelity three‐dimensional (3‐D) numerical simulations of commercial‐scale coal gasifiers are very time consuming and expensive. This article proposes a reduced‐order modeling approach that uses quasi one‐dimensional (1‐D) CFD–DEM simulation results to serve as an accurate initial condition for the 3‐D simulation, which significantly shortens the physical time from several hours to minutes to achieve steady states. The 3‐D simulation employs a previously validated coarse‐grained CFD‐DEM approach, which further accelerates the simulation. Both 1‐D and 3‐D simulations adopt a particle shrinkage scheme to correctly account for the particle volume change due to chemical reaction and thus provides an accurate bed height. Comparison with available experimental measurements of the bed temperature and synthesis gas composition is used to validate the simulation. The final syngas composition and flowrates are strongly affected by the fuel or gasification agent rate and process operating conditions as expected. This approach is then used to optimize two important operating conditions (air and steam flows) to achieve specific gasifier performance goals for the production of Fischer–Tropsch synthesis gas and production of fuel gas for gas turbines.

Yu, Jia↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗

Inhibition of conductive heat flow by magnetic constriction in the corona and transition region - Dependence on the shape of the constriction

Attention is given to the process by which the shape of the constriction of field lines inhibiting the conduction of heat down from the corona acts in conjunction with the degree of constriction to inhibit the heat flow. The analytical model for the heat flow involves tapered flux tubes in which the plasma properties are constant on cross sections, while the plasma is static and the only energy transfer is by thermal conduction. It is determined that only two specific solutions to the model are applicable to the solar atmosphere: the steady state, which is appropriate for quiet regions and for active ones that are not flaring, and the time-dependent case, in which no heat enters the hot end, which is appropriate for conductive cooling of flare loops.

Moore, R. L.↗

Characterization of Swirl-Venturi Lean Direct Injection Designs for Aviation Gas-Turbine Combustion

Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.

Combustion Chamber↗

Human Cargo Resupply Logistics at Mars Using 150kW SEP Tug Cyclers

The ultimate goal of the Evolvable Mars Campaign is to build up a sustainable outpost at Mars that would be continually staffed with rotating crews. During this stage of human Mars exploration, it would be necessary to provision the crews with equipment and supplies both before and during their missions. In this paper, we study the use of 150 kW reusable SEP tugs as a means to deliver elements both to orbit and to the surface. The SEP tugs would make use of technology currently being developed for the proposed Asteroid Robotic Redirect Mission (ARRM). They would also be used to deliver food and supplies to sustain the crews similar to resupply missions for the International Space Station. The SEP tugs envisioned would be staged at a quasi-stable Lunar Near Rectilinear Orbit (NRO). The tugs would then mate with cargo vessels and xenon propellant being delivered by an SLS launch vehicle and continue on to Mars orbit where the cargo is delivered and the SEP tug returns to NRO to repeat the process. It was found that it is more efficient to deliver surface cargo via direct launch and entry versus using the tug cycler. Thousands of optimized low-thrust trajectories were simulated in order to create “Bacon plots” (like porkchop plots, but for low thrust transfers) in order to map out potential trajectories for dates from 2039 to 2052. This study maps out the buildup of a surface outpost as well as the necessary orbital and surface resupply launches in order to maintain it. In the steady state, a cadence of 9 cargo launches is required every 4 years to sustain the human outpost.

Post, Kevin↗

Milestone 1.2.15: Feasibility of In Situ Accelerated Aluminum Coupon Radiolysis by Synchrotron X-Rays

Extended dry storage of aluminum-clad spent nuclear fuel (ASNF) requires an assessment of the radiolytic generation of molecular hydrogen gas (H2) from the ASNF’s corrosion layers. H2 accumulation can potentially lead to storage canister embrittlement/rupture and the accumulation of flammable gas mixtures in accident scenarios. To date, computational models have been developed for the prediction of radiolytic H2 generation from samples exposed to lower absorbed dose regimes (up to ~3 MGy). However, greater accuracy is needed for higher absorbed doses (> 25 MGy) where the concentration of H2 ultimately reaches a steady-state. Data in this area is limited due to the amount of time taken (months) to accumulate such high gamma doses. Furthermore, a recent study observed radiation-induced damage to the surface of pre-corroded aluminum alloy coupons at high absorbed gamma doses. Despite this observation and its implications, there is currently no computational connection between the radiolytic formation of H2 and changes in the composition and morphology of the cladding’s corrosion layers with absorbed dose. To develop a better understanding of the processes and effects described above, the present study aimed to evaluate the feasibility of leveraging synchrotron x-ray capabilities for coupled accelerated irradiation and in situ surface characterization of ASNF alloys. To that end, National Synchrotron Light Source II (NSLS-II) beamtime was secured and a series of aluminum alloy 1100 (AA1100) wires, prepared under a variety of conditions, were interrogated using ex situ x-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques at Idaho National Laboratory (INL) and ex situ and in situ XRD capabilities at the NSLS-II x-ray powder diffraction (XPD) beamline. Results indicated that performing synchrotron XRD using the available XPD beamline configuration can provide sufficient radiation dose but cannot discern changes in the composition and morphology of sample corrosion layers. Additionally, ex situ XRD and SEM data from Idaho National Laboratory (INL) indicated that the pre-corrosion of AA1100 wires did not generate an appreciably thick corrosion layer, as compared with previous aluminum coupon studies. These thinner surface corrosion layers were not sufficient for detection by both ex situ and in situ synchrotron XRD, specifically under consideration of the NSLS-II XPD beamline’s configuration primarily capturing information from the sample’s bulk material, i.e., aluminum metal. In summary, the use of the NSLS-II XPD beamline for promoting accelerated radiolysis and in situ surface characterization is not appropriate for this program’s goals.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗

Direct simulations of wall-bounded compressible turbulence

A study has been initiated on the effects of compressibility (Mach number) on turbulent boundary layers. An understanding of both qualitative (turbulence structures, physical processes) and quantitative (turbulence statistics) effects are desired. This understanding should lead to better turbulence models for applications involving supersonic wall-bounded flows. Direct numerical simulations of an idealized problem will be used to accomplish these objectives. Among several possibilities we chose plane Couette flow with constant-temperature walls as the first problem to be studied. The lack of mean streamwise pressure gradient plus isothermal walls implies that both horizontal directions can be assumed to be homogeneous and that the flow can reach a statistically steady state. Together, these features greatly simplify the calculations and analyses of the results. To date, an algorithm has been developed and implemented (but not yet fully tested) for the accurate solution of the Navier-Stokes equations with the assumptions noted above.

Buell, J. C.↗

Initial development of a generic fluoride salt-cooled reactor model

Fluoride high-temperature reactors (FHRs) are high-temperature, low-pressure reactor concepts that use tri-structural isotropic (TRISO) fuel and molten fluoride salt coolant. These reactors have the potential to provide both electrical power and high-temperature process heat. We used generic FHR parameters for a pebble-bed FHR to develop an initial model with fresh fuel for a generic FHR (gFHR) in MELCOR and SCALE (NEWT and KENO). In this paper, we present the development of our gFHR models, which will serve as the baseline for a sensitivity and uncertainty analysis to quantify the range of possible source terms for FHRs in severe accidents. We present MELCOR results for fuel and coolant temperatures through the core, a nodalization study for the steady-state thermal hydraulic model, and development of reactor physics models in SCALE. As this work progresses, these models will be used to calculate source terms for a loss-of-forced-flow accident and to conduct a sensitivity study on this accident to establish a range of possible source terms. SCALE will provide reactor physics parameters like isotopic inventory, decay heat generation, and temperature coefficients of reactivity. Using the uncertainty quantification tools within SCALE, we will generate distributions for those parameters and will use the uncertainty quantification code RAVEN or DAKOTA to sample those distributions in MELCOR to quantify the impact of reactor physics and thermal hydraulic uncertainties on FHR source terms. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Preliminary proof-of-concept of real-time divertor heat flux control from infrared cameras with nitrogen injection in the DIII-D tokamak

In future tokamak reactors like ITER and the Fusion Pilot Plant (FPP), real-time feedback control of heat flux to the plasma-facing components (PFC) will be critical for steady-state operation. This work presents the first experimental demonstration of real-time divertor heat flux estimation with infrared thermography and feedback control with impurity seeding on the DIII-D tokamak. The flexible infrastructure of the Plasma Control System (PCS) on DIII-D makes this new capability possible. The PCS software runs on a gateway computer system, and five real-time compute nodes. An array of low latency streaming digitizers from D-TACQ Solutions connects to these real-time computers to collect and process data, and send commands to actuators during plasma discharges. This system handles the signal IO from the tokamak and allows the PCS to utilize the diagnostic data necessary to perform control in real-time. Feedback control on heat flux was accomplished by feeding infrared camera data from the “Infrared TV” (IRTV) camera to a custom-developed User Datagram Protocol (UDP) server. This server transmits infrared data to a newly developed PCS algorithm that estimates the heat flux to PFC. Here, a proportional integral derivative (PID) controller minimizes the error between a heat flux reference and the real-time estimate by injecting nitrogen gas into the divertor.

Algorithm development↗

Multi-input, multi-output regulator design for constant disturbances and non-zero set points with application to automatic landing in a crosswind

Undesirable steady offsets result when a stationary, linear regulator using state feedback is subjected to constant disturbances and/or non-zero setspoints. To eliminate these offsets, the disturbances and non-zero setpoints can be fed forward to the control. Only when the number of outputs is less than or equal to the number of control inputs can the outputs be maintained at arbitrary non-zero setpoints. The state and the disturbance may be estimated using a constant gain Kalman filter or by modeling the constant disturbances as exponentially correlated processes with long correlation times.

Holley, W. E.↗

Intermediate-phase engineering via dimethylammonium cation additive for stable perovskite solar cells

Achieving the long-term stability of perovskite solar cells is arguably the most important challenge required to enable widespread commercialization. Understanding the perovskite crystallization process and its direct impact on device stability is critical to achieving this goal. The commonly employed dimethyl-formamide/dimethyl-sulfoxide solvent preparation method results in a poor crystal quality and microstructure of the polycrystalline perovskite films. In this work, we introduce a high-temperature dimethyl-sulfoxide-free processing method that utilizes dimethylammonium chloride as an additive to control the perovskite intermediate precursor phases. By controlling the crystallization sequence, we tune the grain size, texturing, orientation (corner-up versus face-up) and crystallinity of the formamidinium (FA)/caesium (FA) y Cs 1-y Pb(I x Br 1-x )3 perovskite system. A population of encapsulated devices showed improved operational stability, with a median T80 lifetime (the time over which the device power conversion efficiency decreases to 80% of its initial value) for the steady-state power conversion efficiency of 1,190 hours, and a champion device showed a T80 of 1,410 hours, under simulated sunlight at 65 degrees C in air, under open-circuit conditions. In conclusion, this work highlights the importance of material quality in achieving the long-term operational stability of perovskite optoelectronic devices.

14 SOLAR ENERGY↗

PACT Center: Perovskite PV Accelerator for Commercializing Technologies (Final Technical Report)

The Perovskite PV Accelerator for Commercializing Technologies (PACT) center was established in July 2021 as a national resource to accelerate the commercialization of perovskite photovoltaic (PV) technology in the United States. Since its inception, PACT has been led by Sandia National Laboratories (Sandia) in partnership with the National Laboratory of the Rockies (NLR), formerly known as NREL. From FY20-FY23, Los Alamos National Laboratory (LANL), CFV Labs, Black & Veatch (B&V), and the Electric Power Research Institute (EPRI) were part of the project team. LANL brought expertise in perovskite PV device designs and processing, CFV Labs (now GroundWork Renewables) provided initial indoor and outdoor measurement hardware technology, B&V led the initial effort on perovskite PV bankability, and EPRI worked on reviewing testing standards, identifying commercialization gaps, and helping to run PACT’s Industry Advisory Board, a group including representatives from commercial testing labs, independent engineering firms, insurance companies, state regulators, and electric utilities. To source perovskite PV module samples, PACT contracted with the University of North Carolina (UNC), the University of Toledo, the University of Washington, and SLAC/Stanford University to provide a steady stream of research-grade perovskite mini modules, enabling protocol development in advance of commercial module availability. The project period ran from July 1, 2021, through December 31, 2025, including a No Cost Extension. Starting in FY25, the project was continued as a Core Capability in the Lab Call portfolio and continues at a reduced budget with only Sandia and NLR as funded recipients. Notably, starting in FY25 PACT expanded its scope beyond MHP modules to accept all emerging PV mini module technologies for testing, including organic PV (OPV) and all-thin-film tandems, with the aim of supporting commercialization across the broader emerging PV ecosystem. With this change in scope the program was renamed the PV Accelerator for Commercializing Technologies, dropping perovskite from the name.

14 SOLAR ENERGY↗

Modulation of terrestrial ion escape flux composition /by low-altitude acceleration and charge exchange chemistry/

Motivated by recent observations of highly variable hot plasma composition in the magnetosphere, control of the ionospheric escape flux composition by low-altitude particle dynamics and ion chemistry has been investigated for an e(-), H(+), O(+) ionosphere. It is found that the fraction of the steady state escape flux which is O(+) can be controlled very sensitively by the occurrence of parallel or transverse ion acceleration at altitudes below the altitude where the neutral oxygen density falls rapidly below the neutral hydrogen density and the ionospheric source of O(+) tends to be rapidly converted by charge exchange to H(+). The acceleration is required both to overcome the gravitational confinement of O(+) and to violate charge exchange equilibrium so that the neutral hydrogen atmosphere appears 'optically' thin to escaping O(+). Constraints are placed on the acceleration processes, and it is shown that O(+) escape is facilitated by observed ionospheric responses to magnetic activity.

Moore, T. E.↗

Plastic flow of plasma sprayed ceramics

The plastic flow of plasma-sprayed ZrO2-8Y2O3 ceramic has been measured at temperatures up to 1250 C and compared to the plastic flow of pressed and sintered ZrO2-8Y2O3. Plasma spraying of binary oxide ceramics is found to result in a metastable state which is inelastic at high temperature but can also be stabilized or devitrified through heat treating so as to decrease plastic properties. The mechanical properties of the as-plasma sprayed and devitrified ceramic sheet material was measured. An improved algorithm that incorporates the inherently nonlinear thermomechanical field equations was used to determine the influence of inelastic material behavior on the thermomechanical response of ceramic coated seal components. Significant creep was found during the thermal shock and steady heating periods with insufficient time during thermal quench to reverse the process, thereby inducing significant residual stresses into the components.

Padovan, J.↗

A Summary of the Quasi-Steady Acceleration Environment on-Board STS-94 (MSL-1)

The continuous free-fall state of a low Earth orbit experienced by NASA's Orbiters results in a unique reduced gravity environment. While microgravity science experiments are conducted in this reduced gravity environment, various accelerometer systems measure and record the microgravity acceleration environment for real-time and post-flight correlation with microgravity science data. This overall microgravity acceleration environment is comprised of quasi-steady, oscillatory, and transient contributions. The First Microgravity Science Laboratory (MSL-1) payload was dedicated to experiments studying various microgravity science disciplines, including combustion, fluid physics, and materials processing. In support of the MSL-1 payload, two systems capable of measuring the quasi-steady acceleration environment were flown: the Orbital Acceleration Research Experiment (OARE) and the Microgravity Measurement Assembly (MMA) system's Accelerometre Spatiale Triaxiale most evident in the quasi-steady acceleration regime. Utilizing such quasi-steady events, a comparison and summary of the quasi-steady acceleration environment for STS-94 will be presented

McPherson, Kevin M.↗