Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closure problem”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Proposed Semi-Analytic Benchmark for Coupled Neutronics/Thermal-Hydraulics Feedback in Multiphysics Simulations

The past decade has seen considerable advancements in multiscale and multiphysics simulations that account for nonlinear effects of feedback between different physical processes within a single calculation. One area of particular interest has been multiphysics simulations that include both neutronics and thermal-hydraulics effects to accurately model the thermal feedback effects present in nuclear reactors. During this time, many groups have successfully modeled thermal feedback effects by coupling a variety of different types of neutron transport and thermal-hydraulic solvers together. However, formal verification and validation of these coupled code systems has proven challenging due to a lack of benchmark problems with well-characterized reference solutions. One approach is to confirm gross energy balance using bulk coolant properties (outlet temperature/enthalpy) for code verification purposes. Alternatively, code-to-code comparisons of results for a predefined reference/challenge problem have been used to establish reasonability of local results (neutron flux, heating, enthalpy rise) for various implementations, for example the VERA Core Physics Progression Problems. While community solutions are valuable for assessing the relative performance and fidelity of different implementations, they are not guaranteed to provide a measure of the absolute accuracy of a given implementation. In fact, these code-to-code comparisons can give a false sense of confidence by masking the effects of errors or approximations that are common to all of the solution methods tested. This is a particular concern when many methods derive from a common approach, such as Picard iteration. To this end, an analytical benchmark would be very valuable for verifying individual coupled physics implementations, as well as for providing a formal assessment of algorithm accuracy and sensitivity to model approximations such as mesh size and number of convergence iterations. However, creation of an analytical benchmark is difficult due to the complex nature of the eigenvalue form of the neutron transport equation (which does not permit a simple analytical solution for non-trivial cases) as well as the variety of different empirical closure relationships used in thermal-hydraulic modeling. In this paper we propose a benchmark for coupled neutronics/thermal-hydraulics simulations. The proposed benchmark overcomes the challenges associated with a multiphysics benchmark problem by using a fixed-source (rather than eigenvalue) formulation of the neutron transport equation and limiting fluid properties to the single-phase regime where the density and temperature of the coolant are related to enthalpy by a known relationship, while still preserving significant nuclear/thermal feedback effects within the simplified problem. The new benchmark is easy to simulate using existing coupled physics implementations, and allows for a direct comparison with a reference solution, as provided in this summary.

42 ENGINEERING↗
DOE OSTI2016

Final Technical Report

We have constructed and analyzed a series of idealized models of tropical convection interacting with large-scale circulations, with 25-50km resolution and with 1-2km cloud resolving resolution to set the stage for rigorous tests of convection closure schemes in high resolution global climate models. Much of the focus has been on the climatology of tropical cyclogenesis in rotating systems and the related problem of the spontaneous aggregation of convection in non-rotating systems. The PI (Held) will be delivering the honorary Bjerknes lecture at the Fall 2016 AGU meeting in December on this work. We have also provided new analyses of long-standing issues related to the interaction between convection and the large-scale circulation: Kelvin waves in the upper troposphere and lower stratosphere, water vapor transport into the stratosphere, and upper tropospheric temperature trends. The results of these analyses help to improve our understanding of processes, and provide tests for future high resolution global modeling. Our final goal of testing new convections schemes in next-generation global atmospheric models at GFDL has been left for future work due to the complexity of the idealized model results meant as tests for these models uncovered in this work and to computational resource limitations. 11 papers have been published with support from this grant, 2 are in review, and another major summary paper is in preparation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Internal Nozzle Flow Simulations of Gasoline-Like Fuels under Diesel Operating Conditions

Spray formation in internal combustion engines with direct injection is strictly correlated with internal nozzle flow characteristics, which are in turn influenced by fuel physical properties and injector needle motion. This paper pre-sents a series of 3D simulations that model the in-nozzle flow in a 5-hole mini-sac diesel injector. Two gasoline-like naphtha fuels, namely full-range and light naphtha, were tested under operating conditions typical of diesel applica-tions and were compared with n-dodecane, selected from a palette used as diesel surrogates. Validated methodolo-gies from our previous work were employed to account for realistic needle motion. The multi-phase nature of the problem was described by the mixture model assumption with the Volume of Fluid method. Cavitation effects were included by means of the Homogeneous Relaxation Model and turbulence closure was achieved with the Standard k-ε model in an Unsteady Reynolds-Averaged Navier-Stokes formulation. The results revealed that injector perfor-mance and propensity to cavitation are influenced by the fuel properties. Analyses of several physical quantities were carried out to highlight the fuel-to-fuel differences in terms of mass flow rate, discharge coefficients, and fuel vapor volume fraction inside the orifices. A series of parametric investigations was also performed to assess the fuel response to varied fuel injection temperature, injection pressure, and cross-sectional orifice area. For all cases, the strict correlation between cavitation magnitude and saturation pressure was confirmed. Owing to their higher volatil-ity, the two gasoline-like fuels were characterized by higher cavitation across all the simulated conditions. Occur-rence of cavitation was mostly found at the needle seat and at the orifice inlets during the injection event’s transient, when very small gaps exist between the needle and its seat. This behavior tended to disappear at maximum needle lift, where cavitation was absent for all fuels. Differences in mass flow rate between the naphtha fuels and n-dodecane were measured and ascribed to the different densities of the three fuels. Nevertheless, they were found to be smaller than expected, owing to the lower viscosity of the gasoline-like fuels. This beneficial influence of the lower viscosity was shown to be less effective at higher temperature, where the relative viscosity differences de-creased.

Torelli, R.↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗
DOE OSTI2017

Turbulent Mixing Simulation via a Quantum Algorithm

Probability density function (PDF) methods have been very useful in describing many physical aspects of turbulent mixing. In applications of these methods, modeled PDF transport equations are commonly simulated via classical Monte Carlo techniques, which provide estimates of moments of the PDF at arbitrary accuracy. In this paper, recently developed techniques in quantum computing and quantum enhanced measurements (quantum metrology) are used to construct a quantum algorithm that accelerates the computation of such estimates. This quantum algorithm provides a quadratic speedup over classical Monte Carlo methods in terms of the number of repetitions needed to achieve the desired precision. This paper illustrates the power of this algorithm by considering a binary scalar mixing process modeled by means of the coalescence/dispersion (C/D) closure. The equation is first simulated using classical Monte Carlo methods, where error estimates for the computation of central moments are provided. Then the quantum algorithm for this problem is simulated by sampling from the same probability distribution as that of the output of a quantum computer, and it is shown that significantly fewer resources are required to achieve the same precision. Finally, the results demonstrate potential applications of future quantum computers for simulation of turbulent mixing, and large classes of related problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Multi-Material Closure Model for High-Order Finite Element Lagrangian Hydrodynamics

We present a new closure model for single fluid, multi-material Lagrangian hydrodynamics and its application to high-order finite element discretizations of these equations [1]. The model is general with respect to the number of materials, dimension and space and time discretizations. Knowledge about exact material interfaces is not required. Material indicator functions are evolved by a closure computation at each quadrature point of mixed cells, which can be viewed as a high-order variational generalization of the method of Tipton [2]. This computation is defined by the notion of partial non-instantaneous pressure equilibration, while the full pressure equilibration is achieved by both the closure model and the hydrodynamic motion. Exchange of internal energy between materials is derived through entropy considerations, that is, every material produces positive entropy, and the total entropy production is maximized in compression and minimized in expansion. Results are presented for standard one-dimensional two-material problems, followed by two-dimensional and three-dimensional multi-material high-velocity impact arbitrary Lagrangian–Eulerian calculations. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Panel Session 92: Reducing Risk to Enable Endstates

This panel focuses on creative and first-of-a-kind initiatives and approaches to improve DoE's options for achieving success in managing its most difficult wastes. Looking across the Complex, panelists will provide updates to some of the solutions being applied to DoE's challenging waste streams on sites including Hanford, Idaho, Savannah River, and other sites. Presentations proposed include: 1 - Tank Closure Cesium Removal (Mark Keefer, Savannah River Remediation (USA)); 2 - Evolution of Saltstone Disposal Units (Kent Rosenberger, Savannah River Remediation (USA)); 3 - CRESP's Landfill Partnership: Applied Research to Solve Challenging Problems for On-Site Disposal Facilities (Craig Benson, University of Virginia - CRESP (USA)); 4 - Reducing Risk to Enable Endstates (Frederick Hughes, Fluor Idaho, LLC (USA))

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Simulation of Wingtip Vortex Flows with Reynolds-Averaged Navier–Stokes and Scale-Resolving Simulation Methods

The simulation of wingtip vortex flows is reported on with Reynolds-averaged Navier–Stokes (RANS) equations and scale-resolving simulation (SRS) models. These range from turbulent viscosity and Reynolds-stress model (RSM) RANS closures, to hybrid and bridging SRS methods. The aim of the study was threefold: assess the relevance of replicating the experimental flow conditions, evaluate the numerical requisites of each mathematical model, and determine their modeling accuracy in prediction of a representative wingtip vortex flow. The selected problem is the flow around a NACA 0012 wing at 10 deg of angle of attack and Reynolds number of 4.60×10 6 . The results confirm the relevance of reproducing the experimental flow conditions, in particular at the inlet boundary where the flow is not uniform. For this reason, the evaluation of modeling errors using the available measurements requires the specification of the experimental inlet conditions. On the other hand, the quantification of the modeling error indicates that solely RANS–RSM can achieve an accurate representation of the flow dynamics. Whereas the well-recognized limitations of turbulent viscosity RANS closures to deal with solid-body rotation lead to the overprediction of turbulence and consequent rapid diffusion of the wingtip vortex, the predictive use of SRS methods may reveal excessively complex and demanding.

42 ENGINEERING↗
DOE OSTI2017

Problems Related to Temperature-Induced Processes in Long-Term Repositories for Nuclear Waste - 17514

Concepts for the final disposal of high level radioactive waste in deep geological formations must consider the heat pulse produced from radioactive decay during operations and in the post-closure phase as an important factor for the integrity of barrier systems. Several safety cases in different countries limit the maximum temperature at the boundary between the technical barrier/canister surface to the host rock/buffer material to 100 C. Increases in maximum temperature have very diverse effects with respect to the type, intensity and range of the heat pulse. These depend, among other factors, on i) the amount of heat emitted, ii) the type of host rock including site-specific characteristics, iii) the safety and repository concept and further site-specific boundary conditions. The German Commission on the Disposal of High-level Radioactive Waste (Commission) recommended, for precautionary reasons, that the temperature limit at the outer surface of waste containers should not exceed 100 C. Higher temperatures could only be considered if the maximum physically possible temperatures in the respective host rock have been reliably set as a result of research findings. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A fluid-kinetic framework for self-consistent runaway-electron simulations

The problem of self-consistently coupling kinetic runaway-electron physics to the macroscopic evolution of the plasma is addressed by dividing the electron population into a bulk and a tail. A probabilistic closure is adopted to determine the coupling between the bulk and the tail populations, preserving them both as genuine, non-negative distribution functions. Here, macroscopic one-fluid equations and the kinetic equation for the runaway-electron population are then derived, now displaying sink and source terms due to transfer of electrons between the bulk and the tail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Reverse Initiation of the Reynolds RP-87 Detonator Resulting from Off-Normal Functioning of the Ultrafast Closure Valve System

A series of tests were executed to ensure complete UCV high explosive material consumption in the case of a system anomaly occurring — specifically a firing event where one or more RP-87 detonators do not function as a result of an arming and firing system problem or some other circumstance where a detonator does not function even though a firing signal is received at the bridgewire. A simplified linear array was utilized with 5 detonators per test firing only one detonator and maintaining the critical as-designed UCV dimensions and HE mass-to-volume ratios which can affect sympathetic initiation behavior. Tests were performed using both KTech Corporation's stated detonator gap of 0.060" as well as larger gaps to determine safety margin. The detonator's explosive contents were completely consumed in all cases where the gap was maintained within stated specifications. Only extreme gap conditions (0.30") indicated remaining explosive material. The mechanism by which the material is reacted does change with increasing gap however, transferring to a deflagration at larger gaps. All tests were performed at SNL/NM site 9930 in April 2013.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗