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 55 records · Page 3

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, C. G.↗

The application of Newman crack-closure model to predicting fatigue crack growth

Newman crack-closure model and the relevant crack growth program were applied to the analysis of crack growth under constant amplitude and aircraft spectrum loading on a number of aluminum alloy materials. The analysis was performed for available test data of 2219-T851, 2024-T3, 2024-T351, 7075-T651, 2324-T39, and 7150-T651 aluminum materials. The results showed that the constraint factor is a significant factor in the method. The determination of the constraint factor is discussed. For constant amplitude loading, satisfactory crack growth lives could be predicted. For the above aluminum specimens, the ratio of predicted to experimental lives, Np/Nt, ranged from 0.74 to 1.36. The mean value of Np/Nt was 0.97. For a specified complex spectrum loading, predicted crack growth lives are not in very good agreement with the test data. Further effort is needed to correctly simulate the transition between plane strain and plane stress conditions, existing near the crack tip.

Si, Erjian↗

A dynamic closure modeling framework for model order reduction of geophysical flows

Here in this paper, a dynamic closure modeling approach has been derived to stabilize the projection-based reduced order models in the long-term evolution of forced-dissipative dynamical systems. To simplify our derivation without losing generalizability, the proposed reduced order modeling (ROM) framework is first constructed by Galerkin projection of the single-layer quasigeostrophic equation, a standard prototype of large-scale general circulation models, onto a set of dominant proper orthogonal decomposition modes. We then propose an eddy viscosity closure approach to stabilize the resulting surrogate model considering the analogy between large eddy simulation (LES) and truncated modal projection. Our efforts, in particular, include the translation of the dynamic subgrid-scale model into our ROM setting by defining a test truncation similar to the test filtering in LES. Finally, the a posteriori analysis shows that our approach is remarkably accurate, allowing us to integrate simulations over long time intervals at a nominally small computational overhead.

42 ENGINEERING↗

Problems in simulating the stratocumulus-topped boundary layer with a third-order closure model

The Andre et al. (1976, 1978) third-order closure model, in which the time rate of change terms, the relaxation and rapid effects for pressure-related terms, and the clipping approximation are used along with the quasi-normal closure, is invoked in the study of turbulence in a cloudy layer that is radiatively cooled from above. A spurious oscillation whose greatest amplitude lies near the inversion is shown by analysis to arise from the mean gradient and buoyancy terms of the triple-moment equations. An attempt is made to damp the oscillation through the introduction of diffusion terms into the triple-moment equations. The results obtained are noted to be sensitive to the ad hoc eddy coefficient applied in the third-moment equations.

Moeng, C.-H.↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This report summarizes the work done to validate the droplet entrainment and de-entrainment models as well as two-phase closure models in the CTF code by comparison with experimental data obtained at Riso National Laboratory. The Riso data included a series of over 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions were set so that the majority of cases were in the annular/mist ow regime. Measurements included liquid lm ow rate, droplet ow rate, lm thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet ow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor and entrained liquid ow fractions, pressure drop per unit length, and liquid lm thickness. The entire process of generating CTF input decks, running cases, extracting data, and generating comparison plots was scripted using Python and Matplotlib for a completely automated validation process. All test cases and scripting tools have been committed to the COBRA-TF master repository and selected cases have been added to the continuous testing system to serve as regression tests. The dierences between the CTF- and experimentally-calculated ow fraction values were con- sistent with previous calculations by Wurtz, who applied the same entrainment correlation to the same data. It has been found that CTF's entrainment/de-entrainment predictive capability in the annular/mist ow regime for this particular facility is comparable to the licensed industry code, COBRAG. While lm and droplet predictions are generally good, it has been found that accuracy is diminished at lower ow qualities. This nding is consistent with the noted deciencies in the Wurtz entrainment model employed by CTF. The CTF predicted two-phase pressure drop in the annular/mist ow regime has been found to be highly inaccurate, exhibiting a clear bias with respect to the experimental data. This inaccuracy led to an investigation that revealed deciencies in the implementation of the annular/mist interfacial friction model, which should be investigated further in the future. Looking to published COBRAG results for this same facility reveal it exhibits no bias with regard to experimental pressure drop results. In addition to the problems with pressure drop prediction, the lm thickness was also signicantly under-predicted by CTF compared to both experimental data and Wurtz's analytical calculations. Film thickness is calculated using a simple geometric relationship and lm void fraction in CTF, which is dependent on slip ratio and interfacial friction. It is possible that the issues aecting the pressure drop and lm void prediction are related.

97 MATHEMATICS AND COMPUTING↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This document summarizes the work to validate the droplet entrainment and de-entrainment models, as well as two-phase closure models in the CTF code [1] by comparison with experimental data obtained at Risø National Laboratory [2]. CTF is the modernized version of the subchannel code, COBRA-TF [3], being jointly developed by Oak Ridge National Laboratory and North Carolina State University for applications in the Consortium for Advanced Simulation of Light Water Reactors (CASL). The Risø data included a series of more than 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions wereset so that the majority of cases were in the annular/mist flow regime. Measurements included liquid film flow rate, droplet flow rate, film thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet flow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor flow fractions, and entrained liquid flow fractions, and pressure drop per unit length.

Wysocki, Aaron↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This document summarizes the work to validate the droplet entrainment and de-entrainment models, as well as two-phase closure models in the CTF (Coolant Boiling in Rod Arrays - Two Fluids) code by comparison with experimental data obtained at Risoe National Laboratory. CTF is the modernized version of the subchannel code, COBRA-TF, being jointly developed by Oak Ridge National Laboratory and North Carolina State University for applications in the Consortium for Advanced Simulation of Light Water Reactors (CASL). The Risoe data included a series of more than 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions were set so that the majority of cases was in the annular/mist flow regime. Measurements included liquid film flow rate, droplet flow rate, film thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet flow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor flow fractions, and entrained liquid flow fractions, and pressure drop per unit length.

42 ENGINEERING↗

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗

A study of Reynolds-stress closure model

A hybrid model of the Reynolds stress closure was developed. This model was tested for various sizes of step flow, and the computed Reynolds stress behavior was compared with experimental data. The third order closure model was reviewed. Transport equations for the triple velocity correlation were developed and implemented in a numerical code to evaluate the behavior of the triple velocity products in various regions of the flow field including recirculating, reattaching, and redeveloping flow domains.

Amano, R. S.↗

Differentiable physics-enabled closure modeling for Burgers’ turbulence

Abstract Data-driven turbulence modeling is experiencing a surge in interest following algorithmic and hardware developments in the data sciences. We discuss an approach using the differentiable physics paradigm that combines known physics with machine learning to develop closure models for Burgers’ turbulence. We consider the one-dimensional Burgers system as a prototypical test problem for modeling the unresolved terms in advection-dominated turbulence problems. We train a series of models that incorporate varying degrees of physical assumptions on an a posteriori loss function to test the efficacy of models across a range of system parameters, including viscosity, time, and grid resolution. We find that constraining models with inductive biases in the form of partial differential equations that contain known physics or existing closure approaches produces highly data-efficient, accurate, and generalizable models, outperforming state-of-the-art baselines. Addition of structure in the form of physics information also brings a level of interpretability to the models, potentially offering a stepping stone to the future of closure modeling.

97 MATHEMATICS AND COMPUTING↗

Predicting the growth of small and large cracks using a crack-closure model

An evaluation is conducted of the ability of a crack-closure model to predict the fatigue life of notched specimens under spectrum loading, for the cases of crack initiation and growth along the bore of a semicircular notch and the growth of large cracks in 2024-T3 Al alloy sheets under constant-amplitude loading. These cases are used to establish an initial defect size and relationship between the stress-intensity factor range and crack-growth rate. Experimental and calculation results for small-crack growth rates are in good agreement and exhibited a stress level effect at equal values of stress-intensity factor range.

Newman, J. C., Jr.↗

3D Cell-centered hydrodynamics with subscale closure model and multi-material remap

We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiscale modeling framework model (superparameterized CAM5) with a higher‐order turbulence closure: Model description and low‐cloud simulations

Abstract In this study, a higher‐order turbulence closure scheme, called Cloud Layers Unified By Binormals (CLUBB), is implemented into a Multiscale Modeling Framework (MMF) model to improve low‐cloud simulations. The performance of CLUBB in MMF simulations with two different microphysics configurations (one‐moment cloud microphysics without aerosol treatment and two‐moment cloud microphysics coupled with aerosol treatment) is evaluated against observations and further compared with results from the Community Atmosphere Model, Version 5 (CAM5) with conventional cloud parameterizations. CLUBB is found to improve low‐cloud simulations in the MMF, and the improvement is particularly evident in the stratocumulus‐to‐cumulus transition regions. Compared to the single‐moment cloud microphysics, CLUBB with two‐moment microphysics produces clouds that are closer to the coast and agrees better with observations. In the stratocumulus‐to‐cumulus transition regions, CLUBB with two‐moment cloud microphysics produces short‐wave cloud forcing in better agreement with observations, while CLUBB with single‐moment cloud microphysics overestimates short‐wave cloud forcing. CLUBB is further found to produce quantitatively similar improvements in the MMF and CAM5, with slightly better performance in the MMF simulations (e.g., MMF with CLUBB generally produces low clouds that are closer to the coast than CAM5 with CLUBB). Improved low‐cloud simulations in MMF make it an even more attractive tool for studying aerosol‐cloud‐precipitation interactions.

54 ENVIRONMENTAL SCIENCES↗

A preliminary compressible second-order closure model for high speed flows

A preliminary version of a compressible second-order closure model that was developed in connection with the National Aero-Space Plane Project is presented. The model requires the solution of transport equations for the Favre-averaged Reynolds stress tensor and dissipation rate. Gradient transport hypotheses are used for the Reynolds heat flux, mass flux, and turbulent diffusion terms. Some brief remarks are made about the direction of future research to generalize the model.

Speziale, Charles G.↗

A critique of some recent second-order turbulence closure models for compressible boundary layers

Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.

Rubesin, M. W.↗

A new turbulence closure model for boundary layer flows with strong adverse pressure gradients and separation

A new turbulence closure model designed specifically to treat two-dimensional, turbulent boundary layers with strong adverse pressure gradients and attendant separation, is presented. The influence of history effects are modeled by using an ordinary differential equation (ODE) derived from the turbulence kinetic-energy equation, to describe the streamwise development of the maximum Reynolds shear stress in conjunction with an assumed eddy-viscosity distribution which has as its velocity scale the maximum Reynolds shear stress. In the outer part of the boundary layer, the eddy viscosity is treated as a free parameter which is adjusted in order to satisfy the ODE for the maximum shear stress. Because of this, the model s not simply an eddy-viscosity model, but contains features of a Reynolds-stress model. Comparisons with experiments are presented which clearly show the proposed model to be superior to the Cebeci-Smith model in treating strongly retarded and separated flows. In contrast to two-equation, eddy-viscosity models, it requires only slightly more computational effort than simple models like the Cebeci-Smith model.

Johnson, D. A.↗

Point-particle drag, lift, and torque closure models using machine learning: Hierarchical approach and interpretability

Developing deterministic neighborhood-informed point-particle closure models using machine learning has garnered interest recently from the dispersed multiphase flow community. The robustness of neural models for this complex multibody problem is hindered by the availability of particle-resolved data. Here, the present work addresses this unavoidable limitation of data paucity by implementing two strategies: (1) by using a rotation and reflection equivariant neural network and (2) by pursuing a physics-based hierarchical machine learning approach. The resulting machine-learned models are observed to achieve a maximum accuracy of 85% and 96% in the prediction of neighbor-induced force and torque fluctuations, respectively, for a wide range of Reynolds number and volume fraction conditions considered. Furthermore, we pursue force and torque network architectures that provide universal prediction spanning a wide range of Reynolds number (0.25 ≤ Re ≤250) and particle volume fraction (0 ≤ φ ≤0.4). The hierarchical nature of the approach enables improved prediction of quantities such as streamwise torque, by going beyond binary interactions to include trinary interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗