SEARCH · Engineering Papers

Results for “Fast Fluid Dynamics”

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.

176 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Fast particles in drift wave turbulence

This study aims to incorporate the effects of fast particles into our present fluid model for tokamak transport. The parameter ε f = ω / ω f, where ω is the mode frequency and ω f is the typical frequency of the fast particles, which enters as a factor in front of the fast particle response. Thus, for trapped fast particles, where ω f = ω pres the precession frequency of the fast particles, this parameter is of order 10 – 2 for drift waves, and thus, the fast particle response can be neglected. However, ε f will be of order 1 for fast particle modes such as in the fishbone instability. An important turbulence property, affecting both these limits, is resonance broadening. Effects of resonance broadening have recently been considered for fast particle instabilities, often coupled directly to the linear growth rate, while we here consider the original Dupree formulation where the turbulence directly drives a nonlinear frequency shift. Resonance broadening has a general tendency to counteract dissipative wave particle resonances. This has been observed for fast particle instabilities. Here, there is a resonant external source for the fast particles, so the instability survives if this source is dominant over the resonance broadening. For drift waves, however, external sources are not resonant since ε f << 1. Furthermore, the resonance broadening is able to remove the dissipative wave particle resonance completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Observation of shock-front separation in multi-ion-species collisional plasma shocks

We observe shock-front separation and species-dependent shock widths in multi-ion-species collisional plasma shocks, which are produced by obliquely merging plasma jets of a He/Ar mixture (97% He and 3% Ar by initial number density) on the plasma liner experiment. Visible plasma emission near the He-I 587.6-nm and Ar-II 476.5–514.5-nm lines is simultaneously recorded by splitting a single visible image of the shock into two different fast-framing cameras with different narrow bandpass filters (589 ± 5 nm for observing the He-I line and 500 ± 25 nm for the Ar-II lines). For conditions in these experiments (pre-shock ion and electron densities ≈5×10 14 cm –3 , ion and electron temperatures of ≈2.2 eV, and relative plasma-merging speed of 22 km/s), the observationally inferred magnitude of He/Ar shock-front separation and the shock widths themselves are <1 cm, which correspond to ~50 post-shock thermal ion–ion mean free paths. These experimental length scales are in reasonable qualitative and quantitative agreement with results from 1D multi-fluid simulations using the CHICAGO code. However, we note there are differences between the experimentally inferred and simulation-predicted ionization states and line emission intensities, particularly in the post-shock region. Overall, the experimental and simulation results are consistent with theoretical predictions that the lighter He ions diffuse farther ahead within the overall shock front than the heavier Ar ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗

DOE OSTI2024

Report on Initial Sodium Testing on the Thermal Hydraulic Experimental Test Article (THETA) (Fiscal Year 2024 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that is used to develop sodium components and instrumentation as well as to acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility simulates nominal thermal hydraulic conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field may be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility was designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary coolant and secondary coolant system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. In fiscal year 2023, thermal stratification tests were completed with the primary system online, while the secondary system was being constructed [1]. These tests had shown that the core barrel and intermediate heat exchanger (IHX) outlet required increased thermal insulation. The THETA primary system was removed from METL, cleaned, thermal insulators installed, and then inserted into METL Test Vessel 4. At the time of this writing the THETA primary and secondary system are operational. During this fiscal year 100+ hours of testing was completed to characterize thermal hydraulic phenomena associated with steady state and transient conditions in a pool type liquid metal cooled reactor. A majority of the testing campaign was completed to satisfy the experimental data acquisition requirements for the GAIN Voucher with Oklo, CRADA 2021-21121. THETA is still operational at the time of this publication and future testing is planned for fiscal year 2025. Work is underway to publish existing and future data to an online database to facilitate collaboration with SFR engineers looking to validate their systems code or computational fluid dynamics models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Perspectives on Microfluidics for the Study of Asphaltenes in Upstream Hydrocarbon Production: A Minireview

The utilization of microfluidics has generated deep insights into asphaltene precipitation mechanisms and oil–water emulsion stabilization. Agglomeration and precipitation of asphaltenes can cause flow assurance problems during the extraction and transportation of crude oil. Change in temperature, pressure, reservoir conditions, and solvents can change the local environment, leading to asphaltene precipitation. Understanding asphaltene properties and precipitation pathways becomes critical in devising mitigation methods, demulsifiers, and suitable conditions during hydrocarbon processing. Microfluidics has helped in high throughput measurement studies, understanding critical processing conditions, fast demulsifier screening, and the effect of solvent concentration on deposition, generating useful information for utilization at the point of resource extraction facilitating improved resource management. It has become possible to capture the porous, complex nature of reservoir formations and the interaction of chemicals during precipitation through integrated analytics and visualization studies available only through microfluidics. The use of droplet microfluidics, with optical microscopy and high-speed imaging to study the oil–water interface, has resulted in greater understanding of the role of asphaltenes in interfacial properties and emulsion stabilization. Here, this minireview highlights the crucial aspects of microfluidics that have been used to understand physicochemical behavior and dynamics of asphaltene deposition. Some of the unique devices have been presented focusing on the key elements of microfluidics design, fabrication, and analysis, as the insight obtained from microfluidics strongly depends on the device design and the controllability of the experimental parameters. Successful implementation of microfluidics for efficient and controlled experiments, short analysis time scales and rapid screening, and generation of high-quality, reliable data that convey asphaltene deposition issues and interface behavior in emulsions shows the importance of microsystems for advancing knowledge in hydrocarbon production and processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗

DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗

DOE OSTI2023

Modeling DOE Standard Canister Configurations with Updated Surface Chemistry

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Conformally coated scaffold design using water-tolerant Pr 1.8 Ba 0.2 NiO 4.1 for protonic ceramic electrochemical cells with 5,000-h electrolysis stability

Protonic ceramic electrochemical cells (PCECs) have potential as long-duration energy storage systems. However, their operational stability is limited under industrially relevant conditions due to the intrinsic chemical instability of doped barium cerate-based electrolytes and oxygen electrodes against H 2 O, as well as the poor electrode–electrolyte interfacial contact. Here, in this study, we present a conformally coated scaffold (CCS) design to comprehensively address these issues. A porous proton-conducting scaffold is constructed and conformally coated with Pr 1.8 Ba 0.2 NiO 4.1 electrocatalyst, which has high chemical stability against H 2 O, triple conductivity and hydration capability, and protects vulnerable electrolytes from H 2 O. The CCS structure consolidates the electrode–electrolyte interfacial bonding to enable fast proton transfer in the percolated network. This design enables PCECs to reach electrolysis stability for 5,000 h at −1.5 A cm−2 and 600 °C in 40% H 2 O. This work provides a general strategy to stabilize PCECs and offers guidance for designing resilient and stable solid-state energy storage systems.

Protonic ceramic electrochemical cell↗

DOE OSTI2022

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Rapid Evaluation Framework for the CMIP7 Assessment Fast Track

As Earth system models (ESMs) grow in complexity and in volume of output data, there is an increasing need for rapid, comprehensive evaluation of their scientific performance. The upcoming Assessment Fast Track for the Seventh Phase of the Coupled Model Intercomparison Project (CMIP7) will require expeditious response for model analyses designed to inform and drive integrated Earth system assessments. To meet this challenge, the Rapid Evaluation Framework (REF), a community-driven platform for benchmarking and performance assessment of ESMs, was designed and developed. The initial implementation of the REF, constructed to meet the near-term needs of the CMIP7 Assessment Fast Track, builds upon four disparate community evaluation and benchmarking tools that are coupled together using the Coordinated Model Evaluation Capabilities (CMEC) framework. The REF runs within a containerized workflow for portability and reproducibility and is aimed at generating and organizing diagnostics covering a variety of model variables. The REF leverages well documented observational datasets to provide assessments of model fidelity across a collection of diagnostics. All diagnostics were identified and selected with community involvement and consultation. Operational integration with the Earth System Grid Federation (ESGF) will permit automated execution of the REF for selected diagnostics as soon as model output data are published on ESGF by the originating modeling centers. The REF is designed to be portable across a range of current computational platforms to facilitate use by modeling centers for assessing the evolution of model versions or gauging the relative performance of CMIP simulations before being published on ESGF. When integrated into production simulation workflows, results from the REF provide immediate quantitative feedback that allows model developers and scientists to quickly identify model biases and performance issues. After the REF is released to the community, its subsequent development and support will be prioritized by an international consortium of scientists and engineers, enabling a broader impact across Earth science disciplines. For instance, the REF will facilitate improvements to models and will enhance confidence in model projections through process-based selection of models based on their performance with respect to observations. Production of reproducible diagnostics and community-based assessments are key features of the REF. Furthermore, providing interoperability with existing evaluation packages assures that contributions from previous community efforts will be available for use in future model intercomparison projects.

Hoffman, Forrest [ORNL] (ORCID:0000000158024134)↗