Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “closure model”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An Atomic Mix Closure for Stochastic Media Transport Problems

The accurate modeling of stochastic media in radiation transport calculations remains on ongoing research challenge. The Levermore-Pomraning (LP) closure for the statistical transport equation is simple and inexpensive, but it is often inaccurate. More accurate closures have been proposed, but they can be prohibitively expensive. In this work we present a relatively inexpensive closure that we have observed to be more accurate than LP in almost all cases studied. Obtaining solutions to stochastic transport problems can be quite difficult. Transport calculations on an ensemble of explicit realizations generated by Monte Carlo sampling can be prohibitively expensive. Atomic mix or LP calculations are relatively inexpensive, but they can be inaccurate. The present work illustrates an approach that is less expensive than Monte Carlo sampling yet more accurate than atomic mix or LP approximations. We still need to analyze the reasons why the atomic mix closures derived for thicker problems lead to source iteration instabilities. It is not clear if the closure itself induces an ill-posed problem or if it is merely the iterative process that is problematic. We hope to examine alternative source terms to drive the subsidiary calculations. We also want to extend the work to multigroup and multidimensional problems. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1997

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗
DOE OSTI2015

Using ASCEM Modeling and Visualization to Inform Stakeholders of Contaminant Plume Evolution and Remediation Efficacy at F-Basin Savannah River, SC – 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along with its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various toolsets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS’s monitoring activities and operational decisions.

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

Using ASCEM modeling and visualization to inform stakeholders of contaminant plume evolution and remediation efficacy at F-Basin Savannah river, SC - 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along the its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various tool-sets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS's monitoring activities and operational decisions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
DOE OSTI2017

Variability in GRMHD Simulations of Sgr A{sup ∗}: Implications for EHT Closure Phase Observations

Closure phases along different baseline triangles carry a large amount of information regarding the structures of the images of black holes in interferometric observations with the Event Horizon Telescope. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A* to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of the small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the funnel-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the u–v space that show a high level of variability will need to be identified and excised.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

RESOLVING THE INNER JET STRUCTURE OF 1924-292 WITH THE EVENT HORIZON TELESCOPE

We present the first 1.3 mm (230 GHz) very long baseline interferometry model image of an active galactic nucleus (AGN) jet using closure phase techniques with a four-element array. The model image of the quasar 1924-292 was obtained with four telescopes at three observatories: the James Clerk Maxwell Telescope on Mauna Kea in Hawaii, the Arizona Radio Observatory's Submillimeter Telescope in Arizona, and two telescopes of the Combined Array for Research in Millimeter-wave Astronomy in California in 2009 April. With the greatly improved resolution compared with previous observations and robust closure phase measurement, the inner jet structure of 1924-292 was spatially resolved. The inner jet extends to the northwest along a position angle of -53 Degree-Sign at a distance of 0.38 mas from the tentatively identified core, in agreement with the inner jet structure inferred from lower frequencies, and making a position angle difference of {approx}80 Degree-Sign with respect to the centimeter jet. The size of the compact core is 0.15 pc with a brightness temperature of 1.2 Multiplication-Sign 10{sup 11} K. Compared with those measured at lower frequencies, the low brightness temperature may argue in favor of the decelerating jet model or particle-cascade models. The successful measurement of closure phase paves the way for imaging and time resolving Sgr A* and nearby AGNs with the Event Horizon Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Linearly exact parallel closures for slab geometry

Parallel closures are obtained by solving a linearized kinetic equation with a model collision operator using the Fourier transform method. The closures expressed in wave number space are exact for time-dependent linear problems to within the limits of the model collision operator. In the adiabatic, collisionless limit, an inverse Fourier transform is performed to obtain integral (nonlocal) parallel closures in real space; parallel heat flow and viscosity closures for density, temperature, and flow velocity equations replace Braginskii's parallel closure relations, and parallel flow velocity and heat flow closures for density and temperature equations replace Spitzer's parallel transport relations. It is verified that the closures reproduce the exact linear response function of Hammett and Perkins [Phys. Rev. Lett. 64, 3019 (1990)] for Landau damping given a temperature gradient. In contrast to their approximate closures where the vanishing viscosity coefficient numerically gives an exact response, our closures relate the heat flow and nonvanishing viscosity to temperature and flow velocity (gradients)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Turbulent transport and mixing in transitional Rayleigh-Taylor unstable flow: A priori assessment of gradient-diffusion and similarity modeling

Data from a 1152×760×1280 direct numerical simulation of a Rayleigh-Taylor mixing layer modeled after a small-Atwood-number water-channel experiment is used to investigate the validity of gradient diffusion and similarity closures a priori. The budgets of the mean flow, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate transport equations across the mixing layer were previously analyzed at different evolution times to identify the most important transport and mixing mechanisms. Here a methodology is introduced to systematically estimate model coefficients as a function of time in the closures of the dynamically significant terms in the transport equations by minimizing the L2 norm of the difference between the model and correlations constructed using the simulation data. It is shown that gradient-diffusion and similarity closures used for the turbulent kinetic energy K, turbulent kinetic energy dissipation rate ε, heavy-fluid mass fraction variance S, and heavy-fluid mass fraction variance dissipation rate χ equations capture the shape of the exact, unclosed profiles well over the nonlinear and turbulent evolution regimes. Using order-of-magnitude estimates for the terms in the exact transport equations and their closure models, it is shown that several of the standard closures for the turbulent production and dissipation (destruction) must be modified to include Reynolds-number scalings appropriate for Rayleigh-Taylor flow at small to intermediate Reynolds numbers. The late-time, large Reynolds number coefficients are determined to be different from those used in shear flow applications and largely adopted in two-equation Reynolds-averaged Navier-Stokes (RANS) models of Rayleigh-Taylor turbulent mixing. In addition, it is shown that the predictions of the Boussinesq model for the Reynolds stress agree better with the data when additional buoyancy-related terms are included. It is shown that an unsteady RANS paradigm is needed to predict the transitional flow dynamics from early evolution times, analogous to the small Reynolds number modifications in RANS models of wall-bounded flows in which the production-to-dissipation ratio is far from equilibrium. Although the present study is specific to one particular flow and one set of initial conditions, the methodology could be applied to calibrations of other Rayleigh-Taylor flows with different initial conditions (which may give different results during the early-time, transitional flow stages, and perhaps asymptotic stage). Finally, the implications of these findings for developing high-fidelity eddy viscosity-based turbulent transport and mixing models of Rayleigh-Taylor turbulence are discussed.

42 ENGINEERING↗
DOE OSTI2015

Applying understanding of likelihood and consequences to demonstrate post-closure criticality safety of a geological disposal facility - 14377

Radioactive Waste Management (RWM) is responsible for implementing geological disposal of the UK's higher activity wastes. These wastes include significant quantities of fissile nuclides and, therefore, demonstration of criticality safety under disposal conditions forms an important component of RWM's disposal system safety case. At the time of disposal, controls (e.g. fissile material limits and/or water barriers) will ensure that waste packages are sub-critical. These controls, plus the containment provided by the engineered barrier system of a geological disposal facility (GDF) will also ensure that criticality is prevented for a long time after disposal. However, over post-closure timescales conditions in the GDF will evolve, and therefore underpinning knowledge is required to demonstrate: under what conditions criticality could occur; what the likelihood of such conditions developing is; and what the consequences would be if a post-closure criticality did occur. These requirements are set out in the UK environment agencies' Guidance on Requirements for Authorisation of a GDF. This paper discusses the disposal of Pressurised Water Reactor fuel, as an example of a waste type that may require disposal. It briefly selects relevant research to combine and summarise results from three recent RWM research projects. The first 'Likelihood of Criticality' considered under what conditions post-closure criticality could occur, whilst a second 'Modelling of Consequences of Hypothetical Criticality' focused on understanding the local physical consequences (temperature rise, power and rock response) of hypothetical 'what-if' post-closure criticality. The paper then explains how this understanding is applied to make safety arguments as to why post-closure criticality is not a significant concern. This is achieved via a third, stream of work that has considered the effects of a potential criticality on pathways that might give rise to a post-closure risk in the future, such as transport of radionuclides in groundwater. Collectively RWMs post-closure criticality safety work allows us to conclude in our generic environmental safety case that a post closure criticality is a low likelihood event, and even if one of these low likelihood events was to occur its consequences would be insignificant to GDF performance. This paper also serves to give context to two companion UK post-closure criticality safety papers which communicate further information on our modelling development and more detailed results underpinning our demonstration of post-closure criticality safety. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1987

Computation of turbulent channel flow using Large-Eddy Interaction Model

The objective of the paper is to investigate the nature and values of closure parameters appearing in the proposed Large-Eddy Interaction Model for prediction of turbulent flow field. Effects of two closure parameters on predicted Reynolds stresses and other turbulence structural quantities are examined for channel flows at two Reynolds numbers.

Hong, S. K.↗
NASA NTRS1988

A multiple-time-scale turbulence model based on variable partitioning of turbulent kinetic energy spectrum

The paper presents a multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method. Consideration is given to a class of turbulent boundary layer flows and of separated and/or swirling elliptic turbulent flows. For the separated and/or swirling turbulent flows, the present turbulence model yielded significantly improved computational results over those obtained with the standard k-epsilon turbulence model.

Kim, S.-W.↗
DOE OSTI2017

K-ε Turbulence Model Parameter Estimates Using an Approximate Self-similar Jet-in-Crossflow Solution

The k-ε turbulence model has been described as perhaps “the most widely used complete turbulence model.” This family of heuristic Reynolds Averaged Navier-Stokes (RANS) turbulence closures is supported by a suite of model parameters that have been estimated by demanding the satisfaction of well-established canonical flows such as homogeneous shear flow, log-law behavior, etc. While this procedure does yield a set of so-called nominal parameters, it is abundantly clear that they do not provide a universally satisfactory turbulence model that is capable of simulating complex flows. Recent work on the Bayesian calibration of the k-ε model using jet-in-crossflow wind tunnel data has yielded parameter estimates that are far more predictive than nominal parameter values. In this paper, we develop a self-similar asymptotic solution for axisymmetric jet-in-crossflow interactions and derive analytical estimates of the parameters that were inferred using Bayesian calibration. The self-similar method utilizes a near field approach to estimate the turbulence model parameters while retaining the classical far-field scaling to model flow field quantities. Our parameter values are seen to be far more predictive than the nominal values, as checked using RANS simulations and experimental measurements. They are also closer to the Bayesian estimates than the nominal parameters. A traditional simplified jet trajectory model is explicitly related to the turbulence model parameters and is shown to yield good agreement with measurement when utilizing the analytical derived turbulence model coefficients. Finally, the close agreement between the turbulence model coefficients obtained via Bayesian calibration and the analytically estimated coefficients derived in this paper is consistent with the contention that the Bayesian calibration approach is firmly rooted in the underlying physical description.

42 ENGINEERING↗
NASA NTRS2015

Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes

The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.

Abdol-Hamid, Khaled S.↗
DOE OSTI2012

Implementation and Validation of the BHR Turbulence Model in the FLAG Hydrocode

The BHR-2 turbulence model, developed at Los Alamos National Laboratory for variable density and compressible flows, is implemented in an Arbitrary Lagrangian-Eulerian hydrocode, FLAG. The BHR-2 formulation is discussed, with emphasis on its connection to multi-component flow formulations that underlie FLAG's treatment of multi-species flow. One-dimensional and two-dimensional validation tests are performed and compared to experiment and Eulerian simulations. Turbulence is an often studied and ubiquitous phenomenon in nature, and modeling its effects is essential in many practical applications. Specifically the behavior of turbulence in the presence of strong density gradients and compressibility is of fundamental importance in applications ranging from Inertial Confinement Fusion (ICF) [1], supernovae [2], and atmospheric flows. The BHR closure approach [3] seeks to model the physical processes at work in variable density turbulence including Kelvin-Helmholtz (KH) [4], Rayleigh-Taylor (RT) [5], and Richtmyer-Meshkov (RM) [6], driven turbulence. The effectiveness of the BHR-2 implementation has been demonstrated for variable density mixing in the KH, RT, and RM cases in an Eulerian framework [7]. The primary motivation of the present work is to implement the BHR-2 turbulence model in the Arbitrary Lagrangian-Eulerian (ALE) hydrodynamics code FLAG. The goal is not only to demonstrate results in agreement with previous Eulerian calculations, but also document behavior that arises from the underlying differences in code philosophy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

A statistical analysis of the dependency of closure assumptions in cumulus parameterization on the horizontal resolution

Simulated data from the UCLA cumulus ensemble model are used to investigate the quasi-universal validity of closure assumptions used in existing cumulus parameterizations. A closure assumption is quasi-universally valid if it is sensitive neither to convective cloud regimes nor to horizontal resolutions of large-scale/mesoscale models. The dependency of three types of closure assumptions, as classified by Arakawa and Chen, on the horizontal resolution is addressed in this study. Type I is the constraint on the coupling of the time tendencies of large-scale temperature and water vapor mixing ratio. Type II is the constraint on the coupling of cumulus heating and cumulus drying. Type III is a direct constraint on the intensity of a cumulus ensemble. The macroscopic behavior of simulated cumulus convection is first compared with the observed behavior in view of Type I and Type II closure assumptions using 'quick-look' and canonical correlation analyses. It is found that they are statistically similar to each other. The three types of closure assumptions are further examined with simulated data averaged over selected subdomain sizes ranging from 64 to 512 km. It is found that the dependency of Type I and Type II closure assumptions on the horizontal resolution is very weak and that Type III closure assumption is somewhat dependent upon the horizontal resolution. The influences of convective and mesoscale processes on the closure assumptions are also addressed by comparing the structures of canonical components with the corresponding vertical profiles in the convective and stratiform regions of cumulus ensembles analyzed directly from simulated data. The implication of these results for cumulus parameterization is discussed.

Xu, Kuan-Man↗