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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Modelling and measurement of crack closure and crack growth following overloads and underloads

Ignoring crack growth retardation following overloads can result in overly conservative life predictions in structures subjected to variable amplitude fatigue loading. Crack closure is believed to contribute to the crack growth retardation, although the specific closure mechanism is dabatable. The delay period and corresponding crack growth rate transients following overload and overload/underload cycles were systematically measured as a function of load ratio and overload magnitude. These responses are correlated in terms of the local 'driving force' for crack growth, i.e. the effective stress intensity factor range. Experimental results are compared with the predictions of a Dugdale-type (1960) crack closure model, and improvements in the model are suggested.

Dexter, R. J.↗
DOE OSTI2016

Newton-Krylov Method in Applications of Solving Two-phase Problems using Drift Flux Model

In this paper, Jacobian-free Newton-Krylov (JFNK) method is investigated in an application to implicitly solve two-phase flow problems using four-equation drift flux model. The closure models include the EPRI drift flux correlations and additional constitutive models to determine flow regimes, wall boiling and interfacial heat/mass transfer, two-phase flow wall friction, etc. Different from many traditional computer codes, fully implicit methods are used for the time integration. The resulted nonlinear discretized equation system is solved using the JFNK method. Expensive and error prone derivation and implementation of the analytical Jacobian matrix are avoided. Numerical results are successfully validated using existing experimental data on flow boiling under forced convection conditions in both pipe and rod bundle geometries.

42 ENGINEERING↗
DOE OSTI2018

Device-scale CFD modeling of gas-liquid multiphase flow and amine absorption for CO2 capture

In this paper we aim to develop a validated device-scale CFD model that can predict quantitatively both hydrodynamics and CO2 capture efficiency for an amine-based solvent absorber column with random Pall ring packing. A Eulerian porous-media approach and a two-fluid model were employed, in which the momentum and mass transfer equations were closed by literature-based empirical closure models. We proposed a hierarchical approach for calibrating the parameters in the closure models to make them accurate for the packed column. Specifically, a parameter for momentum transfer in the closure was first calibrated based on data from a single experiment. With this calibrated parameter, a parameter in the closure for mass transfer was next calibrated under a single operating condition. Last, the closure of the wetting area was calibrated for each gas velocity at three different liquid flow rates. For each calibration, cross validations were pursued using the experimental data under operating conditions different from those used for calibrations. This hierarchical approach can be generally applied to develop validated device-scale CFD models for different absorption columns.

Pan, Wenxiao↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
NASA NTRS1993

Development of a recursion RNG-based turbulence model

Reynolds stress closure models based on the recursion renormalization group theory are developed for the prediction of turbulent separated flows. The proposed model uses a finite wavenumber truncation scheme to account for the spectral distribution of energy. In particular, the model incorporates effects of both local and nonlocal interactions. The nonlocal interactions are shown to yield a contribution identical to that from the epsilon-renormalization group (RNG), while the local interactions introduce higher order dispersive effects. A formal analysis of the model is presented and its ability to accurately predict separated flows is analyzed from a combined theoretical and computational stand point. Turbulent flow past a backward facing step is chosen as a test case and the results obtained based on detailed computations demonstrate that the proposed recursion -RNG model with finite cut-off wavenumber can yield very good predictions for the backstep problem.

Zhou, YE↗
DOE OSTI2018

Gen-0 BWR Modeling Application and Validation (L3 Milestone THM.P18.01)

This milestone report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the thermal hydraulics Methods (THM) Level 3 milestone L3:THM.P18.01 Gen-0 BWR Modeling Application and Validation. The extension of existing boiling closure designed for Pressurized Water Reactors (PWRs) to the high void fraction condition of Boiling Water Reactors (BWRs) is one of the challenges for Multiphase Computational Fluid Dynamics (M-CFD) applications. In the efforts aimed at developing the advanced boiling closure for BWRs, the baseline closure model, so-called ‘Zero’ closure, has been developed by leveraging the experience gained through the Numerical Nuclear Reactor project. This baseline closure model is physics-based, simplest and robust, but it does not contain any explicit descriptions of the complex high void fraction mechanisms. This is a continuation of the CASL thermal-hydraulic study to develop and assess the CFD-BWR models. In FY-17, the applicability has been evaluated against the international OECD/NRC BFBT benchmark data. In the present work, we assessed the capability of CFD-BWR model for void drift phenomena in the modern BWR fuel designs which make use of both full-length and part-length rods.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

A Reynolds-Averaged Navier-Stokes Perspective for the High Lift Common Research Model Using the LAVA Framework

An in-depth assessment of the High Lift-Common Research Model (HL-CRM) aerodynamic performance is made using a Reynolds-Averaged Navier-Stokes (RANS) methodology. Results presented here were produced in support of the 4th AIAA High Lift Prediction Workshop(HLPW-4), with additional simulation results and analysis used to support findings from the workshop. Mesh resolution, RANS closure models, and tunnel effects were among the sources of simulation sensitivity addressed, all in the pursuit of evaluating best-practice simulation methods for complex high-lift configurations. In general, free-air simulations using RANS yield accurate predictions of integrated loads in the pre-stall regime for the HL-CRM, while strong sensitivity to simulation method is observed near CLmax. Additionally, a study assessing changes in aerodynamic performance due to flap deflection angle indicated strong sensitivity to turbulence model closure and mesh resolution, particularly at the highest deflection setting. Several available corrections for the Spalart-Allmaras (SA) turbulence model were tested and generally introduced more error relative to the experiment in comparison to the baseline SA closure model. The Standard Menter Baseline Two-Equation (k-ωBSL) turbulence model was also explored, which arguably yielded results highly consistent with experiment, but at a substantially higher computational cost. Therefore, as a compromise between computational cost and accuracy, the baseline SA closure model was selected for additional HL-CRM studies that assess alternative solution methods. These strategies include “warm start” approaches, which initialize the simulation from a previously converged solution at a nearby angle of attack, unsteady calculations, and simulations that model the experimental wind tunnel. Results from these studies would reinforce the shortcomings of RANS in predicting aircraft performance at high-lift conditions, namely at high-αand highly loaded flap conditions, while also answering questions about appropriate modeling techniques to minimize the errors associated with RANS where possible

TTT↗
NASA NTRS1976

Remarks on the construction of a second-order closure description of turbulent reacting flows

In this paper an attempt is made to illustrate how second-order closure techniques might be used to make detailed calculations of turbulent diffusion flames. Two of the problems of closure that are peculiar to turbulent chemically reacting flows are discussed in some detail. These are the problem of computing in a consistent way all the correlations of fluctuations in scalar quantities needed for closure and the problem of the generation of turbulence by chemical reactions. Tentative closure models are put forward for both of these problems. Although the closure model for fluctuations in scalar quantities can, to a certain extent, be checked against idealized numerical examples for which exact solutions are known, the validity of both the closure models proposed here must, in the end, be determined by detailed comparison of computed results with experimental data.

DONALDSON C. DUP.↗
DOE OSTI2017

BWR Full Fuel Assembly Testing and Validation

This report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the Thermal Hydraulics Methods (THM) Level 3 milestone THM.CLS.P15.10: BWR Full Fuel Assembly Testing and Validation. The extension of existing boiling closure designed for Pressurized Water Reactors (PWRs) to the high void fraction condition of Boiling Water Reactors (BWRs) is one of the challenges for Multiphase Computational Fluid Dynamics (M-CFD) applications. In the efforts aimed at developing the advanced boiling closure for BWRs, the baseline closure model, so-called ‘Zero’ closure, has been developed based on the GEN-I closure for PWRs which was used as a starting point. Since this simplest closure model is robust but does not contain any explicit description of the complex high void fraction mechanisms, the applicability of M-CFD needs to be evaluated before developing advanced closures for BWRs. To support this need, the high fidelity M-CFD simulation with Zero closure has been conducted for the BWR full fuel assembly. The capabilities of M-CFD with Zero closure were evaluated against the international OECD/NRC BFBT benchmark data. The scalability of the M-CFD application to the nuclear fuel assembly was demonstrated with the INL High Performance Computing (HPC) system. This report presents the recent view of the best practices on M-CFD for BWRs to provide a useful reference to identify and overcome common challenges in the M-CFD applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2003

Grid-Independent Large-Eddy Simulation in Turbulent Channel Flow using Three-Dimensional Explicit Filtering

In this paper, turbulence-closure models are evaluated using the 'true' LES approach in turbulent channel flow. The study is an extension of the work presented by Gullbrand (2001), where fourth-order commutative filter functions are applied in three dimensions in a fourth-order finite-difference code. The true LES solution is the grid-independent solution to the filtered governing equations. The solution is obtained by keeping the filter width constant while the computational grid is refined. As the grid is refined, the solution converges towards the true LES solution. The true LES solution will depend on the filter width used, but will be independent of the grid resolution. In traditional LES, because the filter is implicit and directly connected to the grid spacing, the solution converges towards a direct numerical simulation (DNS) as the grid is refined, and not towards the solution of the filtered Navier-Stokes equations. The effect of turbulence-closure models is therefore difficult to determine in traditional LES because, as the grid is refined, more turbulence length scales are resolved and less influence from the models is expected. In contrast, in the true LES formulation, the explicit filter eliminates all scales that are smaller than the filter cutoff, regardless of the grid resolution. This ensures that the resolved length-scales do not vary as the grid resolution is changed. In true LES, the cell size must be smaller than or equal to the cutoff length scale of the filter function. The turbulence-closure models investigated are the dynamic Smagorinsky model (DSM), the dynamic mixed model (DMM), and the dynamic reconstruction model (DRM). These turbulence models were previously studied using two-dimensional explicit filtering in turbulent channel flow by Gullbrand & Chow (2002). The DSM by Germano et al. (1991) is used as the USFS model in all the simulations. This enables evaluation of different reconstruction models for the RSFS stresses. The DMM consists of the scale-similarity model (SSM) by Bardina et al. (1983), which is an RSFS model, in linear combination with the DSM. In the DRM, the RSFS stresses are modeled by using an estimate of the unfiltered velocity in the unclosed term, while the USFS stresses are modeled by the DSM. The DSM and the DMM are two commonly used turbulence-closure models, while the DRM is a more recent model.

Gullbrand, Jessica↗
DOE OSTI2023

Multifidelity Approach to Sensitivity Estimation in Large-Eddy Simulation

A novel approach to compute affordable approximate sensitivities in a large-eddy simulation (LES) is proposed and assessed. The approach is based on solving a Reynolds-averaged Navier–Stokes (RANS) problem that has been linearized around the mean LES solution, with closure modeling required for the linearized changes in turbulent Reynolds stresses. In the present study, the closure modeling is based on a linearized algebraic turbulence model. The method is assessed for the flow over a NACA0012 airfoil at a fixed angle of attack, with the Reynolds number as the variable parameter. The results show that, provided an accurate linearized closure model, the method predicts the correct sensitivity of the skin friction coefficient and the mean velocity field at the cost of a linearized RANS, which provides an important proof-of-concept for this approach. Here, the linearized algebraic turbulence model with standard model coefficients produces a reasonably accurate sensitivity, but the results also suggest that accuracy could be gained from recalibrating the model coefficients for this new use of the model.

42 ENGINEERING↗
NASA NTRS1991

An analysis of RNG based turbulence models for homogeneous shear flow

In a recent paper, the authors compared the performance of a variety of turbulence models including the k-epsilon model and the second-order closure model based on Renormalization Group (RNG) Methods. The performance of these RNG models in homogeneous turbulent shear flow was found to be quite poor, apparently due to the value of the constant C(sub epsilon1) in the modeled dissipation rate equation which was substantially lower than its traditional value. However, recently a correction has been made in the RNG based calculation of C(sub epsilon1). It is shown that with the new value of C(sub epsilon1), the performance of the RNG k-epsilon model is substantially improved. On the other hand, while the predictions of the revised RNG second-order closure model are better, some lingering problems still remain which can be easily remedied by the addition of higher order terms.

Speziale, Charles G.↗
NASA NTRS1991

An analysis of RNG-based turbulence models for homogeneous shear flow

In a recent paper, the authors compared the performance of a variety of turbulence models including the k-epsilon model and the second-order closure model based on Renormalization Group (RNG) Methods. The performance of these RNG models in homogeneous turbulent shear flow was found to be quite poor, apparently due to the value of the constant C(sub epsilon 1) in the modeled dissipation rate equation which was substantially lower than its traditional value. However, recently a correction has been made in the RNG-based calculation of C(sub epsilon 1). It is shown that with the new value of C(sub epsilon 1), the performance of the RNG k-epsilon model is substantially improved. On the other hand, while the predictions of the revised RNG second-order closure model are better, some lingering problems still remain which can be easily remedied by the addition of higher order terms.

Speziale, Charles G.↗
NASA NTRS1992

Workshop on Engineering Turbulence Modeling

Discussed here is the future direction of various levels of engineering turbulence modeling related to computational fluid dynamics (CFD) computations for propulsion. For each level of computation, there are a few turbulence models which represent the state-of-the-art for that level. However, it is important to know their capabilities as well as their deficiencies in order to help engineers select and implement the appropriate models in their real world engineering calculations. This will also help turbulence modelers perceive the future directions for improving turbulence models. The focus is on one-point closure models (i.e., from algebraic models to higher order moment closure schemes and partial differential equation methods) which can be applied to CFD computations. However, other schemes helpful in developing one-point closure models, are also discussed.

Povinelli, Louis A.↗
NASA NTRS1992

Effects of constraint on crack growth under aircraft spectrum loading

The objective is to study the effects of constraint on fatigue crack growth under aircraft spectrum loading. A plasticity-induced crack-closure model that accounts for constraint variations during the transition from flat-to-slant crack growth was used to correlate crack-growth rate data under constant-amplitude loading and to calculate crack growth under simulated aircraft spectrum loading. The model was applied to several thin-sheet aluminum alloy materials. Under laboratory air conditions, the transition was shown to be related to the size of the cyclic plastic zone based on the effective stress-intensity factor range for several sheet materials and thicknesses. Results from three-dimensional, elastic-plastic, finite-element analyses of a flat, straight-through crack in a thin-sheet aluminum alloy specimen showed a constraint loss similar to that assumed in the model. Using test data and the closure model, the location of the constraint-loss regime in terms of growth rate and the value of the constraint factor at these rates were determined by trial and error. The model was then used to calculate crack growth under the TWIST spectrum. The calculated results agreed reasonably well with test data. In general, the model predicted shorter crack-growth lives than tests under the TWIST spectrum by about 40 percent. For the TWIST spectrum clipped at Level 3, the calculated lives were within about 20 percent. The results demonstrated that constraint variations, especially for thin-sheet alloys, should be accounted for to predict crack growth under typical aircraft spectra.

Newman, J. C., Jr.↗
DOE OSTI2018

Regional evapotranspiration from an image-based implementation of the Surface Temperature Initiated Closure (STIC1.2) model and its validation across an aridity gradient in the conterminous US

Recent studies have highlighted the need for improved characterizations of aerodynamic conductance and temperature (g A and T 0 ) in thermal remote-sensing-based surface energy balance (SEB) models to reduce uncertainties in regional-scale evapotranspiration (ET) mapping. By integrating radiometric surface temperature (TR) into the Penman–Monteith (PM) equation and finding analytical solutions of gA and T0, this need was recently addressed by the Surface Temperature Initiated Closure (STIC) model. However, previous implementations of STIC were confined to the ecosystem-scale using flux tower observations of infrared temperature. This study demonstrates the first regional-scale implementation of the most recent version of the STIC model (STIC1.2) that integrates the Moderate Resolution Imaging Spectroradiometer (MODIS) derived TR and ancillary land surface variables in conjunction with NLDAS (North American Land Data Assimilation System) atmospheric variables into a combined structure of the PM and Shuttleworth–Wallace (SW) framework for estimating ET at 1 km × 1 km spatial resolution. Evaluation of STIC1.2 at 13 core AmeriFlux sites covering a broad spectrum of climates and biomes across an aridity gradient in the conterminous US suggests that STIC1.2 can provide spatially explicit ET maps with reliable accuracies from dry to wet extremes. When observed ET from one wet, one dry, and one normal precipitation year from all sites were combined, STIC1.2 explained 66 % of the variability in observed 8-day cumulative ET with a root mean square error (RMSE) of 7.4 mm/8-day, mean absolute error (MAE) of 5 mm/8-day, and percent bias (PBIAS) of -4 %. These error statistics showed relatively better accuracies than a widely used but previous version of the SEB-based Surface Energy Balance System (SEBS) model, which utilized a simple NDVI-based parameterization of surface roughness (z OM ), and the PM-based MOD16 ET. SEBS was found to overestimate (PBIAS = 28 %) and MOD16 was found to underestimate ET (PBIAS =-26 %). The performance of STIC1.2 was better in forest and grassland ecosystems as compared to cropland (20 % underestimation) and woody savanna (40 % overestimation). Model inter-comparison suggested that ET differences between the models are robustly correlated with g A and associated roughness length estimation uncertainties which are intrinsically connected to T R uncertainties, vapor pressure deficit (D A ), and vegetation cover. A consistent performance of STIC1.2 in a broad range of hydrological and biome categories, as well as the capacity to capture spatio-temporal ET signatures across an aridity gradient, points to the potential for this simplified analytical model for near-real-time ET mapping from regional to continental scales.

Bhattarai, Nishan (ORCID:0000000327493549)↗
NASA NTRS1993

Simulations of diffusion-reaction equations with implications to turbulent combustion modeling

An enhanced diffusion-reaction reaction system (DRS) is proposed as a statistical model for the evolution of multiple scalars undergoing mixing and reaction in an isotropic turbulence field. The DRS model is close enough to the scalar equations in a reacting flow that other statistical models of turbulent mixing that decouple the velocity field from scalar mixing and reaction (e.g. mapping closure model, assumed-pdf models) cannot distinguish the model equations from the original equations. Numerical simulations of DRS are performed for three scalars evolving from non-premixed initial conditions. A simple one-step reversible reaction is considered. The data from the simulations are used (1) to study the effect of chemical conversion on the evolution of scalar statistics, and (2) to evaluate other models (mapping-closure model, assumed multivariate beta-pdf model).

Girimaji, Sharath S.↗
NASA NTRS1992

Advances in fatigue life prediction methodology for metallic materials

The capabilities of a plasticity-induced crack-closure model to predict small- and large-crack growth rates, and in some cases total fatigue life, for four aluminum alloys and three titanium alloys under constant-amplitude, variable-amplitude, and spectrum loading are described. Equations to calculate a cyclic-plastic-zone corrected effective stress-intensity factor range from a cyclic J-integral and crack-closure analysis of large cracks were reviewed. The effective stress-intensity factor range against crack growth rate relations were used in the closure model to predict small- and large-crack growth under variable-amplitude and spectrum loading. Using the closure model and microstructural features, a total fatigue life prediction method is demonstrated for three aluminum alloys under various load histories.

Newman, J. C., Jr.↗