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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Flow interaction experiment. Volume 2: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗
NASA NTRS1991

Second-order near-wall turbulence closures - A review

Advances in second-order near-wall turbulence closures are summarized. All closures under consideration are based on high-Reynolds-number models. Most near-wall closures proposed to date attempt to modify the high-Reynolds-number models for the dissipation function and the pressure redistribution term so that the resultant models are applicable all the way to the wall. The asymptotic behavior of the near-wall closures is examined and compared with the proper near-wall behavior of the exact Reynolds-stress equations. It is found that three second-order near-wall closures give the best correlations with simulated turbulence statistics. However, their predictions of near-wall Reynolds-stress budgets are considered to be incorrect. A proposed modification to the dissipitation-rate equation remedies part of those predictions. It is concluded that further improvements are required if a complete replication of all the turbulence properties and Reynolds-stress budgets by a statistical model of turbulence is desirable.

So, R. M. C.↗
DOE OSTI2021

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2024

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Sub-grid scale model classification and blending through deep learning

In this article we detail the use of machine learning for spatio-temporally dynamic turbulence model classification and hybridization for large eddy simulations (LES) of turbulence. Here, our predictive framework is devised around the determination of local conditional probabilities for turbulence models that have varying underlying hypotheses. As a first deployment of this learning, we classify a point on our computational grid as that which requires the functional hypothesis, the structural hypothesis or no modelling at all. This ensures that the appropriate model is specified from a priori knowledge and an efficient balance of model characteristics is obtained in a particular flow computation. In addition, we also utilize the conditional-probability predictions of the same machine learning to blend turbulence models for another hybrid closure. Our test case for the demonstration of this concept is given by Kraichnan turbulence, which exhibits a strong interplay of enstrophy and energy cascades in the wavenumber domain. Our results indicate that the proposed methods lead to robust and stable closure and may potentially be used to combine the strengths of various models for complex flow phenomena prediction.

42 ENGINEERING↗
DOE OSTI2026

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗
DOE OSTI2014

Tracing the cosmic metal evolution in the low-redshift intergalactic medium

Using the Cosmic Origins Spectrograph aboard the Hubble Space Telescope, we measured the abundances of six ions (C III, C IV, Si III, Si IV, N V, and O VI) in the low-redshift (z ≤ 0.4) intergalactic medium (IGM). Both C IV and Si IV have increased in abundance by a factor of ∼10 from z ≈ 5.5 to the present. We derive ion mass densities, ρ{sub ion} ≡ Ω{sub ion}ρ{sub cr}, with Ω{sub ion} expressed relative to the closure density. Our models of mass-abundance ratios, (Si III/Si IV) =0.67{sub −0.19}{sup +0.35}, (C III/C IV) =0.70{sub −0.20}{sup +0.43}, and (Ω{sub C} {sub III}+Ω{sub C} {sub IV})/(Ω{sub Si} {sub III}+Ω{sub Si} {sub IV})=4.9{sub −1.1}{sup +2.2}, are consistent with the photoionization parameter log U = –1.5 ± 0.4, hydrogen photoionization rate Γ{sub H} = (8 ± 2) × 10{sup –14} s{sup –1} at z < 0.4, and specific intensity I {sub 0} = (3 ± 1) × 10{sup –23} erg cm{sup –2} s{sup –1} Hz{sup –1} sr{sup –1} at the Lyman limit. Consistent ionization corrections for C and Si are scaled to an ionizing photon flux Φ{sub 0} = 10{sup 4} cm{sup –2} s{sup –1}, baryon overdensity Δ {sub b} ≈ 200 ± 50, and ''alpha-enhancement'' (Si/C enhanced to three times its solar ratio). We compare these metal abundances to the expected IGM enrichment and abundances in higher photoionized states of carbon (C V) and silicon (Si V, Si VI, and Si VII). Our ionization modeling infers IGM metal densities of (5.4 ± 0.5) × 10{sup 5} M {sub ☉} Mpc{sup –3} in the photoionized Lyα forest traced by the C and Si ions and (9.1 ± 0.6) × 10{sup 5} M {sub ☉} Mpc{sup –3} in hotter gas traced by O VI. Combining both phases, the heavy elements in the IGM have mass density ρ {sub Z} = (1.5 ± 0.8) × 10{sup 6} M {sub ☉} Mpc{sup –3} or Ω {sub Z} ≈ 10{sup –5}. This represents 10% ± 5% of the metals produced by (6 ± 2) × 10{sup 8} M {sub ☉} Mpc{sup –3} of integrated star formation with yield y{sub m} = 0.025 ± 0.010. The missing metals at low redshift may reside within galaxies and in undetected ionized gas in galactic halos and circumgalactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Final report on DOE ASCR Grant DESC0001862 (ER25894)

As was outlined in the proposal, the overall objective of this project is to investigate advanced optimization algorithms for solving entropy maximization problems. Entropy maximization is a general tool for reconstructing the least biased estimate of the state of a complex system, based on available information. In transport and kinetic theory, entropy maximization is used to derive closures for moment models in order to reduce the complexity of a kinetic description while maintaining many of its fundamental features. Application areas include gas dynamics, charged-particle transport, thermal radiative transfer, and neutron transport. Some such applications are of prime importance to the Department of Energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Nanoscale Morphology Evolution Under Ion Irradiation

We showed that the half-century-old paradigm of morphological instability under irradiation due to the curvature-dependence of the sputter yield, can account neither for the phase diagram nor the amplification or decay rates that we measure in the simplest possible experimental system -- an elemental semiconductor with an amorphous surface under noble-gas ion irradiation; We showed that a model of pattern formation based on the impact-induced redistribution of atoms that do not get sputtered away explains our experimental observations; We developed a first-principles, parameter-free approach for predicting morphology evolution, starting with molecular dynamics simulations of single ion impacts, lasting picoseconds, and upscaling through a rigorous crater-function formalism to develop a partial differential equation that predicts morphology evolution on time scales more than twelve orders of magnitude longer than can be covered by the molecular dynamics; We performed the first quantitative comparison of the contributions to morphological instability from sputter removal and from impact-induced redistribution of atoms that are removed, and showed that the former is negligible compared to the latter; We established a new paradigm for impact-induced morphology evolution based on crater functions that incorporate both redistribution and sputter effects; and We developed a model of nanopore closure by irradiation-induced stress and irradiationenhanced fluidity, for the near-surface irradiation regime in which nuclear stopping predominates, and showed that it explains many aspects of pore closure kinetics that we measure experimentally.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Diel cycle impacts on the chemical and light absorption properties of organic carbon aerosol from wildfires in the western United States

Abstract. Organic aerosol (OA) emissions from biomass burning have been the subject of intense research in recent years, involving a combination of field campaigns and laboratory studies. These efforts have aimed at improving our limited understanding of the diverse processes and pathways involved in the atmospheric processing and evolution of OA properties, culminating in their accurate parameterizations in climate and chemical transport models. To bring closure between laboratory and field studies, wildfire plumes in the western United States were sampled and characterized for their chemical and optical properties during the ground-based segment of the 2019 Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) field campaign. Using a custom-developed multiwavelength integrated photoacoustic-nephelometer spectrometer in conjunction with a suite of instruments, including an oxidation flow reactor equipped to generate hydroxyl (OH⚫) or nitrate (NO3⚫) radicals to mimic daytime or nighttime oxidative aging processes, we investigated the effects of multiple equivalent hours of OH⚫ or NO3⚫ exposure on the chemical composition and mass absorption cross-sections (MAC(λ)) at 488 and 561 nm of OA emitted from wildfires in Arizona and Oregon. We found that OH⚫ exposure induced a slight initial increase in absorption corresponding to short timescales; however, at longer timescales, the wavelength-dependent MAC(λ) decreased by a factor of 0.72 ± 0.08, consistent with previous laboratory studies and reports of photobleaching. On the other hand, NO3⚫ exposure increased MAC(λ) by a factor of up to 1.69 ± 0.38. We also noted some sensitivity of aerosol aging to different fire conditions between Arizona and Oregon. The MAC(λ) enhancement following NO3⚫ exposure was found to correlate with an enhancement in CHO1N and CHOgt1N ion families measured by an Aerodyne aerosol mass spectrometer.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Low-energy charged particle observations in the 5-20 Jupiter-radius region of the Jovian magnetosphere

Ion (greater than 0.5 MeV) and electron (greater than 30 keV) measurements made by the low-energy charged particle instrument during the Voyager 1 and 2 traversals of the 5-20 Jupiter-radius region of the Jovian magnetosphere are presented. The spatial morphology of particle intensities, energy spectra, and composition is emphasized. Diffusive radial transport is also discussed. The Jovian magnetosphere seemed to be much more disturbed during the Voyager 2 passage than during that of Voyager 1. Significant inbound-outbound asymmetries of the radial profiles of intensities are observed; an appropriate magnetic field model to provide closure has not been found. Low-energy electrons are not enhanced or depleted in the Io torus region except at the inner edge, about 5.2 Jupiter-radius, where they sharply decrease. This may be due to enhanced electron loss associated with a region of increased plasma density.

Armstrong, T. P.↗
NASA NTRS1988

The turbulent heat flux in low Mach number flows with large density variations

A transport equation has been derived which is the difference between the volume- and mass-averaged velocities and is simply related to the turbulent heat flux phi sup h. Using this equation and an assumption analogous to the drift flux approximation of two-phase flow modeling, an algebraic closure relation for phi sup h that exibits fluxes due to directed transport proportional to -del anti p and due to gradient transport proportional to -del tau has been obtained.

Orourke, Peter J.↗
NASA NTRS1990

Numerical simulation of an F-16A at angle of attack

The transonic flowfield around an F-16A fighter configuration at a moderate incidence angle is simulated by solving the Navier-Stokes equations on a single-block grid. The numerical solution matches experimental freestream conditions with a mach number of 0.85, 16 degrees angle of attack, and a characteristic Reynolds number of 12.75 million. MacCormack's explicit algorithm is used in conjuction with a local time step and consecutive mesh refinement procedure to accelerate numerical convergence. The Baldwin-Lomax algebraic model provides turbulent closure. Computed surface pressure distributions and the aircraft lift coefficient compare favorably with wind tunnel data. The drag coefficient in the simulation overpredicts the experimental value by 8 percent.

Huband, G. W.↗
NASA NTRS1991

Computation of turbulent flow in a thin liquid layer of fluid involving a hydraulic jump

Numerically computed flow fields and free surface height distributions are presented for the flow of a thin layer of liquid adjacent to a solid horizontal surface that encounters a hydraulic jump. Two kinds of flow configurations are considered: two-dimensional plane flow and axisymmetric radial flow. The computations used a boundary-fitted moving grid method with a k-epsilon model for the closure of turbulence. The free surface height was determined by an optimization procedure which minimized the error in the pressure distribution on the free surface. It was also checked against an approximate procedure involving integration of the governing equations and use of the MacCormack predictor-corrector method. The computed film height also compared reasonably well with previous experiments. A region of recirculating flow was found to be present adjacent to the solid boundary near the location of the jump, which was caused by a rapid deceleration of the flow.

Rahman, M. M.↗
NASA NTRS1995

PDF methods for turbulent reactive flows

Viewgraphs are presented on computation of turbulent combustion, governing equations, closure problem, PDF modeling of turbulent reactive flows, validation cases, current projects, and collaboration with industry and technology transfer.

Hsu, Andrew T.↗
DOE OSTI2016

Improving performance of two-group interfacial area transport equation models using genetic algorithms and high resolution experimental data

In the two-fluid transport model, the coupling of mass, momentum and energy transfer between phases is highly dependent on interfacial area transfer terms. Several research efforts in the past have been focused on the development of an interfacial area transport equation model (IATE), in order to eliminate the drawbacks of static flow regime maps currently used in best-estimate thermal-hydraulic system codes. The IATE attempts to model the dynamics that are involved in two phase flows by accounting for the different interaction mechanisms affecting bubble transport in the flow. In general, current IATE models perform poorly in high void fraction flow regimes. As both models depend on closure relations that require empirically determined coefficients, a genetic algorithm is implemented to study their optimization. The preliminary study indicates that improvements are achievable by application of the genetic algorithm. Optimization of coefficients for an entire set of high-resolution wire mesh sensor (WMS) data provided a single set of optimal coefficients. The application of WMS-optimized coefficients to external tests improved IATE performance for similar hydraulic diameters and provided a neutral change in performance for significantly smaller/larger hydraulic diameters. This is a crucial outcome as an optimized set of coefficients that provide significant improvement for a subset of experimental data - is shown to benefit performance in external data for which the coefficients are not specifically optimized. (authors)

42 ENGINEERING↗
NASA NTRS1998

An Evaluation of the Plasticity-Induced Crack-Closure Concept and Measurement Methods

An assessment of the plasticity-induced crack-closure concept is made, in light of some of the questions that have been raised on the validity of the concept, and the assumptions that have been made concerning crack-dp damage below the crack-opening stress. The impact of using other crack-tip parameters, such as the cyclic crack-tip displacement, to model crack-growth rate behavior was studied. Crack-growth simulations, using a crack-closure model, showed a close relation between traditional Delta K eff, and the cyclic crack-tip displacement (Delta eff) for an aluminum alloy and a steel. Evaluations of the cyclic hysteresis energy demonstrated that the cyclic plastic damage below the crack-opening stress was negligible in the Paris crack-growth regime. Some of the standard and newly proposed remote measurement methods to determine the 'effective' crack-tip driving parameter were evaluated on middle-crack tension specimens. A potential source of the Kmax effect on crack-growth rates was studied on an aluminum alloy. Results showed that the ratio of Kmax to Kc had a strong effect on crack-growth rates at high stress ratios and at low stress ratios for very high stress levels. The crack-closure concept and the traditional crack-growth rate equations were able to correlate and predict crack-growth rates under these extreme conditions.

Newman, James C., Jr.↗