Search NASA⌕ Search

SEARCH · Search NASA

Results for “point kinetics equations”

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.

206 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗
DOE OSTI2021

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Theoretical modeling of hepatitis C acute infection in liver-humanized mice support pre-clinical assessment of candidate viruses for controlled-human-infection studies

Designing and carrying out a controlled human infection (CHI) model for hepatitis C virus (HCV) is critical for vaccine development. However, key considerations for a CHI model protocol include understanding of the earliest viral-host kinetic events during the acute phase and susceptibility of the viral isolate under consideration for use in the CHI model to antiviral treatment before any infections in human volunteers can take place. Humanized mouse models lack adaptive immune responses but provide a unique opportunity to obtain quantitative understanding of early HCV kinetics and develop mathematical models to further understand viral and innate immune response dynamics during acute HCV infection. We show that the models reproduce the measured HCV kinetics in humanized mice, which are consistent with early acute HCV-host dynamics in immunocompetent chimpanzees. Our findings suggest that humanized mice are well-suited to support development of a CHI model. In-silico and in-vivo modeling estimates provide a starting point to characterize candidate viruses for testing in CHI model studies.

Agent-based modeling↗
DOE OSTI2022

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗