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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗
NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibrium states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, C. G.↗
DOE OSTI2021

PORFLOW Modeling of Vadose Zone Flow and Transport for E-Area Intermediate Level Vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOW TM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater is accessible to members of the public following site closure.

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

PORFLOW Modeling of Vadose Zone Flow and Transport for the E- Area Low Activity Waste Vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Low Activity Waste Vault (LAWV) and local vadose zone has been developed using the PORFLOW TM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the LAWV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be accessed by members of the public following site closure.

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

Turbulence energy and diffusion transport of third-moments in a separating and reattaching flow

For accurate prediction of turbulent flow in separated and reattaching regions, it is necessary to incorporate second- and third-moments of turbulent fluctuations. The turbulence energy and the energy dissipation rate equations are modified by incorporating second-order closure. Moreover, a transport equation model for the third-order closure with a near-wall correction is developed for the evaluation of the diffusive action of the second-moments. After comparison of the results with experimental data, it is concluded that the models developed here improve the prediction of triple-velocity correlations in both recirculating and redeveloping flow regions.

Amano, R. S.↗
NASA NTRS2000

Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing

A novel mapping closure approximation (MCA) technique is developed to construct a model for the conditional dissipation rate (CDR) of a scalar in homogeneous turbulence. It is shown that the CDR model from amplitude mapping closure is incorrect in asymptotic behavior for unsymmetric binary mixing. The correct asymptotic behavior can be described by the CDR model formulated by the MCA technique. The MCA approach is outlined for constructing successive approximation to probability density function (PDF) and conditional moment.

He, Gouwei↗
NASA NTRS2006

Jet Noise Diagnostics Supporting Statistical Noise Prediction Methods

The primary focus of my presentation is the development of the jet noise prediction code JeNo with most examples coming from the experimental work that drove the theoretical development and validation. JeNo is a statistical jet noise prediction code, based upon the Lilley acoustic analogy. Our approach uses time-average 2-D or 3-D mean and turbulent statistics of the flow as input. The output is source distributions and spectral directivity. NASA has been investing in development of statistical jet noise prediction tools because these seem to fit the middle ground that allows enough flexibility and fidelity for jet noise source diagnostics while having reasonable computational requirements. These tools rely on Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solutions as input for computing far-field spectral directivity using an acoustic analogy. There are many ways acoustic analogies can be created, each with a series of assumptions and models, many often taken unknowingly. And the resulting prediction can be easily reverse-engineered by altering the models contained within. However, only an approach which is mathematically sound, with assumptions validated and modeled quantities checked against direct measurement will give consistently correct answers. Many quantities are modeled in acoustic analogies precisely because they have been impossible to measure or calculate, making this requirement a difficult task. The NASA team has spent considerable effort identifying all the assumptions and models used to take the Navier-Stokes equations to the point of a statistical calculation via an acoustic analogy very similar to that proposed by Lilley. Assumptions have been identified and experiments have been developed to test these assumptions. In some cases this has resulted in assumptions being changed. Beginning with the CFD used as input to the acoustic analogy, models for turbulence closure used in RANS CFD codes have been explored and compared against measurements of mean and rms velocity statistics over a range of jet speeds and temperatures. Models for flow parameters used in the acoustic analogy, most notably the space-time correlations of velocity, have been compared against direct measurements, and modified to better fit the observed data. These measurements have been extremely challenging for hot, high speed jets, and represent a sizeable investment in instrumentation development. As an intermediate check that the analysis is predicting the physics intended, phased arrays have been employed to measure source distributions for a wide range of jet cases. And finally, careful far-field spectral directivity measurements have been taken for final validation of the prediction code. Examples of each of these experimental efforts will be presented. The main result of these efforts is a noise prediction code, named JeNo, which is in middevelopment. JeNo is able to consistently predict spectral directivity, including aft angle directivity, for subsonic cold jets of most geometries. Current development on JeNo is focused on extending its capability to hot jets, requiring inclusion of a previously neglected second source associated with thermal fluctuations. A secondary result of the intensive experimentation is the archiving of various flow statistics applicable to other acoustic analogies and to development of time-resolved prediction methods. These will be of lasting value as we look ahead at future challenges to the aeroacoustic experimentalist.

Bridges, James E.↗
NASA NTRS2004

Apparent I-Xe Cooling Rates of Chondrules Compared with Silicates from the Colomera Iron Meteorite

In I-Xe dating, a regular pattern of increasing Xe-129/128Xe-128 ratio with increasing extraction temperature is often observed. If one makes the crude assumption that the temperatures at which the Xe is extracted in the laboratory is approximately the same as the temperature at which those sites closed 4.6 Ga ago, a (zeroth order) model cooling rate can be found. In order to test and refine this model we can apply the cooling theory of Dodson to those extraction steps approaching the I-Xe isochrons. Using an Arrhenius plot for these temperature fractions, and assuming that an only single phase is involved, the effective diffusion parameters can be estimated (frequency factor and activation energy). From the apparent (zeroth order) cooling rate, the closure temperature can be estimated from the Dodson equation. This model closure temperature can then be compared with the actual laboratory temperature at which the isochron begins. The ratio of the closure temperature and the temperature corresponding to the start of the isochron provides the ratio of the two temperature scales, incorporation and extraction. The actual cooling rate is then given by the apparent (zeroth order) cooling rate times the temperature scale factor. Figure 1 shows Arrhenius plots for I-Xe data

Hohenberg, C. M.↗
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 OSTI2020

Risk-Informed MultiPphysics BestEestimate Plus Uncertainties (BEPU) Application Development of RELAP5-3D Perturbation Model

The United States nuclear industry is facing a strong challenge to ensure maximum safety while enhancing economic benefit. Safety is a key parameter to all aspects related to light water reactor (LWR) nuclear power plants (NPPs), especially cost savings. Since the goal is to extend the lifetimes of these NPPs, the traditional deterministic safety concept may not guarantee a current economic asset. The Light Water Reactor Sustainability (LWRS) Program has been promoting a wide range of research and development (R&D) in this field to maximize the safety, economics, and performance of these NPPs through improved scientific understanding. One of the best practices to achieve this goal is to identify and optimize safety margins, which can lead to cost reduction. To do this, under the LWRS framework, the Risk-Informed Systems Analysis (RISA) Pathway will focus on the optimization of safety margin and minimization of uncertainties to ensure both safety and economics at the highest level. The RISA Pathway will provide enhanced capabilities for analyzing and characterizing LWR systems performance by developing and demonstrating methods, tools, and data to enable risk-informed margins management (RIMM). The goals of the RISA Pathway are twofold: (1) deploy the risk-informed tools and methods that enable better representation of safety margins and factors that contribute to cost and safety; and (2) conduct advanced risk assessment applications with industry to support margin management strategies that enable more cost-effective plant operation. The tools and methods provided by the RISA Pathway will support effective margin management for both active and passive safety systems, structures, and components (SSC) of an NPP. The tools and methods used in the RISA Pathway should have high confidence and highest technical maturity for and implementation to industry at its current setting. They should also have a capability to support risk-informed decision making for both probabilistic and deterministic elements of safety. The RISA Pathway will, therefore, perform a comprehensive assessment of verification and validation (V&V) status of RISA Toolkit to enhance credibility RISA Toolkit which be used by industry. This report summarizes RELAP5-3D development activities for best estimated and uncertainty (BEPU) capability to support ongoing RISA Pathway pilot projects for risk informed uncertainty quantification applications. Work scope includes (1) development of initial perturbation model for selected closure laws and (2) testing of suitable probability distribution functions (PDF) for the selected correlations. Developed BEPU modules are tested with selected separate effect cases and compared with conventional thermal-hydraulics model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1992

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

The presently obtained closed-form analytical expressions, which predict the limiting rate of mean reactant conversion in homogeneous turbulent flows under the influence of a binary reaction, are derived via the single-point pdf method based on amplitude mapping closure. With this model, the maximum rate of the mean reactant's decay can be conveniently expressed in terms of definite integrals of the parabolic cylinder functions. The results obtained are shown to be in good agreement with data generated by direct numerical simulations.

Madnia, C. K.↗
DOE OSTI2025

Resolving discrepancies in bang-time predictions for indirect-drive ICF experiments on the NIF: Insights from the Build-A-Hohlraum campaign

This study investigated discrepancies between measured and simulated x-ray drive in Indirect-Drive Inertial Confinement Fusion (ID-ICF) hohlraums at the National Ignition Facility. Despite advances in radiation-hydrodynamic simulations, a consistent “drive deficit” remains. Experimentally measured ID-ICF capsule bang-times are systematically 400–700 ps later than simulations predict. The Build-A-Hohlraum (BAH) campaign explored potential causes for this discrepancy by systematically varying hohlraum features, including laser entrance hole (LEH) windows, capsules, and gas fills. Overall, the agreement between simulated and experimental x-ray drive was found to be largely unaffected by these changes. The data allow us to exclude some hypotheses put forward to potentially explain the discrepancy. Errors in the local thermodynamic equilibrium (LTE) atomic modeling, errors in the modeling of LEH closure, and errors due to a lack of plasma species mix physics in simulations are shown to be inconsistent with our measurements. Instead, the data support the hypothesis that errors in NLTE emission modeling are a significant contributor to the discrepancy. X-ray emission in the 2–4 keV range is found to be approximately 30% lower than in simulations. This is accompanied by higher than predicted electron temperatures in the gold bubble region, pointing to errors in non-LTE modeling. Introducing an opacity multiplier of 0.87 on energy groups above 1.8 keV improves agreement with experimental data, reducing the bang-time discrepancy from 300 to 100 ps. These results underscore the need for refined NLTE opacity models to enhance the predictive power of hohlraum simulations.

Band emission↗
NASA NTRS1996

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗
DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

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

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1973

Effects of fineness and closure ratios on boattail drag of circular-arc-afterbody models with jet exhaust at Mach numbers up to 1.3

The effects of variations in fineness ratio and closure ratio on the boattail drag of circular-arc afterbody models were investigated at static conditions and at Mach numbers from 0.40 to 1.30. Jet total-pressure ratio was varied from jet off to about 6, depending on Mach number. Reynolds number based on model maximum diameter varied from about 1.2 million to 2.18 million. Angle of attack was varied from -4 deg to 8 deg. The results indicate that, at subsonic speeds before the drag rise and at a scheduled pressure ratio, the total drag (pressure plus calculated skin friction) is approximately equal for configurations with the same closure ratio. Also, for configurations with the same closure ratio, increasing fineness ratio results in an increased drag-rise Mach number; and for configurations with the same fineness ratio, increasing closure ratio results in an increased drag-rise Mach number.

Reubush, D. E.↗
NASA NTRS2007

Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures

The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.

Xu, Kuan-Man↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 1: Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be diferent from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites.

Wang, S. S.↗