Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Level structure of light neutron-rich La isotopes beyond the 𝑁 = 82 shell closure

Here, the high spin excited states of Lanthanum isotopes 140–143 La, above the 𝑁 = 82 closed shell, were populated in fission reactions. The prompt 𝛾-ray transitions were measured using two complementary methods: (a) in coincidence with the isotopically identified fragments produced in the fission of the 238 U + 9 Be system using the Variable Mode Spectrometer (VAMOS++) and the Advanced Gamma Tracking Array (AGATA) spectrometer, and (b) high statistics threefold 𝛾−𝛾−𝛾 and fourfold 𝛾−𝛾−𝛾−𝛾 coincidence data from the spontaneous fission of 252 Cf using the Gammasphere. This work reports the first identification of a pair of parity doublet structures in 143 La and the new high spin level structure in 140–142 La from prompt 𝛾-ray spectroscopy. The level structures are interpreted in terms of the systematics of neighboring odd-𝑍 nuclei above the 𝑍 = 50 shell closure and large-scale shell model calculations. The present results indicate the presence of stable octupole deformation in 143 La. The excitation energy pattern and their comparison with neighboring isotones, moving away from the 𝑁 = 82 closed shell, point towards a transition from single-particle structures to an alternating parity rotational band structure in the La isotopic chain.

Navin, A. [Centre National de la Recherche Scienti↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Triple-velocity products in a channel with a backward-facing step

Attention is given to the evaluation of tripple velocity products in the reattaching and developing region behind a backward facing step. The change in tripple velocity products is significant in the wake region, yielding substantial variation in the diffusion rate of the Reynolds stresses. Four models of the third-order closure are examined, and the results are compared with the experimental data for Chandrsuda and Bradshaw (1980).

Amano, R. S.↗
NASA NTRS1987

Evaluation of a theory for pressure-strain rate

A theoretical expression for the slow part (the nonlinear fluctuation part) of the pressure-strain rate is compared with simulations of anisotropic homogeneous flows. The objective is to determine the quantitative accuracy of the theory and to test its prediction that the generalized Rotta coefficient, a non-dimensionalized ratio of slow term to the Reynolds stress anisotropy, varies with direction and can be negative. Comparisons are made between theoretical and simulation values of the slow term itself and of the generalized Rotta coefficients. The implications of the comparison for two-point closure theories and for Reynolds stress modeling are pointed out.

Weinstock, J.↗
NASA NTRS1989

Discussion of turbulence modelling: Past and future

The full text of a paper presented at the Whither Turbulence Workshop (Cornell University, March 22-24, 1989) on past and future trends in turbulence modeling is provided. It is argued that Reynolds stress models are likely to remain the preferred approach for technological applications for at least the next few decades. In general agreement with the Launder position paper, it is further argued that among the variety of Reynolds stress models in use, second-order closures constitute by far the most promising approach. However, some needed improvements in the specification of the turbulent length scale are emphasized. The central points of the paper are illustrated by examples from homogeneous turbulence.

Speziale, Charles G.↗
NASA NTRS1990

Fatigue crack growth in unidirectional metal matrix composite

The weight function method was used to determine the effective stress intensity factor and the crack opening profile for a fatigue tested composite which exhibited fiber bridging. The bridging mechanism was modeled using two approaches; the crack closure approach and the shear lag approach. The numerically determined stress intensity factor values from both methods were compared and correlated with the experimentally obtained crack growth rates for SiC/Ti-15-3 (0)(sub 8) oriented composites. The near crack tip opening profile was also determined for both methods and compared with the experimentally obtained measurements.

Ghosn, Louis J.↗
NASA NTRS1991

Computational study of the aerodynamics and control by blowing of asymmetric vortical flows over delta wings

Some of the work is described which was done in a study of the flow field produced by tangential leading edge blowing on a 60 deg. delta wing. The flow is studied computationally by solving the Thin Layer Navier-Stokes equations. Steady state flow fields are calculated for various angles of attack and yaw, with and without the presence of tangential leading edge blowing. The effectiveness of blowing as a rolling moment control mechanism to extend the envelope of controllability is illustrated at pre- and post-stall angles of attack. The numerical grid is generated using algebraic grid generation and various interpolation and blending techniques. The jet emanates from a slot with linearly varying thickness and is introduced into the flow field using the concept of an actuator plane, thereby not requiring resolution of the jet slot geometry. The Baldwin-Lomax algebraic turbulence model is used to provide turbulent closure. The computational results are compared with those of experiments.

Craig, Ken↗
NASA NTRS1990

Fatigue crack growth in unidirectional metal matrix composite

The weight function method was used to determine the effective stress intensity factor and the crack opening profile for a fatigue tested composite which exhibited fiber bridging. The bridging mechanism was modeled using two approaches; the crack closure approach and the shear lag approach. The numerically determined stress intensity factor values from both methods were compared and correlated with the experimentally obtained crack growth rates for SiC/Ti-15-3 (0) (sub 8) oriented composites. The near crack tip opening profile was also determined for both methods and compared with the experimentally obtained measurements.

Ghosn, L. J.↗
NASA NTRS2000

North Atlantic Aerosol Properties for Radiative Impact Assessments. Derived from Column Closure Analyses in TARFOX and ACE-2

Aerosol effects on atmospheric radiative fluxes provide a forcing function that can change the climate in potentially significant ways. This aerosol radiative forcing is a major source of uncertainty in understanding the climate change of the past century and predicting future climate. To help reduce this uncertainty, the 1996 Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX) and the 1997 Aerosol Characterization Experiment (ACE-2) measured the properties and radiative effects of aerosols over the Atlantic Ocean. Both experiments used remote and in situ measurements from aircraft and the surface, coordinated with overpasses by a variety of satellite radiometers. TARFOX focused on the urban-industrial haze plume flowing from the United States over the western Atlantic, whereas ACE-2 studied aerosols over the eastern Atlantic from both Europe and Africa. These aerosols often have a marked impact on satellite-measured radiances. However, accurate derivation of flux changes, or radiative forcing, from the satellite measured radiances or retrieved aerosol optical depths (AODs) remains a difficult challenge. Here we summarize key initial results from TARFOX and ACE-2, with a focus on closure analyses that yield aerosol microphysical models for use in improved assessments of flux changes. We show how one such model gives computed radiative flux sensitivities (dF/dAOD) that agree with values measured in TARFOX and preliminary values computed for the polluted marine boundary layer in ACE-2. A companion paper uses the model to compute aerosol-induced flux changes over the North Atlantic from AVHRR-derived AOD fields.

Russell, Philip A.↗
NASA NTRS2004

Radiative Transfer Modeling of a Large Pool Fire by Discrete Ordinates, Discrete Transfer, Ray Tracing, Monte Carlo and Moment Methods

Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).

Jensen, K. A.↗
NASA NTRS1986

On the modeling of low-Reynolds-number turbulence

A full Reynolds-stress closure that is capable of describing the flow all the way to the wall was formulated for turbulent flow through circular pipe. Since viscosity does not appear explicitly in the pressure redistribution terms, conventional high-number models for these terms are found to be applicable. However, the models for turbulent diffusion and viscous dissipation have to be modified to account for viscous diffusion near a wall. Two redistribution and two diffusion models are investigated for their effects on the model calculations. Wall correction to pressure redistribution modeling is also examined. Diffusion effects on calculated turbulent properties are further investigated by simplifying the transport equations to algebraic equations for Reynolds stress. Two approximations are explored. These are the equilibrium and nonequilibrium turbulence assumptions. Finally, the two-equation closure is also used to calculate the flow in question and the results compared with all the other model calculations. Fully developed pipe flows at two moderate Reynolds numbers are used to validate these model calculations.

So, R. M. C.↗
NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
NASA NTRS1987

Study of one-dimensional spectral dynamic equations of the Reynolds stresses in homogeneous anisotropic turbulence: Application to split-spectrum modeling

The CTR numerical data base generated by direct simulation of homogeneous anisotropic turbulence was used to calculate all of the terms in the spectral balance equations for the turbulent Reynolds stresses. The aim in not only to test the main closure assumptions used in the split-spectrum models, but also to try to devise improved hypotheses deduced from the statistical information. Numerical simulations of turbulent flows provide a large amount of data, a thought provoking wealth of information. The main advantage of this type of comparison is that a great variety of flows can be considered, and this is necessary to test closure hypotheses. Moreover various initial conditions can be introduced in the calculation, even if they are not experimentally feasible. All the terms in the spectral equations can be calculated. The limited Reynolds numbers of the simulations and the statistical noise caused by a small sample, particularly at the large scales, causes some difficulty in the interpretation of the results, but the method of approach proved to be a powerful tool for testing and improving spectral closures.

Schiestel, R.↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
NASA NTRS2024

Assessment of Dust Size Retrievals Based on AERONET: A Case Study of Radiative Closure From Visible-Near-Infrared to Thermal Infrared

Super-coarse dust particles (diameters >10 μm) are evidenced to be more abundant in the atmosphere than model estimates and contribute significantly to the dust climate impacts. Since super-coarse dust accounts for less dust extinction in the visible-to-near-infrared (VIS-NIR) than in the thermal infrared (TIR) spectral regime, they are suspected to be underestimated by remote sensing instruments operates only in VIS-NIR, including Aerosol Robotic Networks (AERONET), a widely used data set for dust model validation. In this study, we perform a radiative closure assessment using the AERONET-retrieved size distribution in comparison with the collocated Atmospheric Infrared Sounder (AIRS) TIR observations with comprehensive uncertainty analysis. The consistently warm bias in the comparisons suggests a potential underestimation of super-coarse dust in the AERONET retrievals due to the limited VIS-NIR sensitivity. An extra super-coarse mode included in the AERONET-retrieved size distribution helps improve the TIR closure without deteriorating the retrieval accuracy in the VIS-NIR.

thermal infrared (TIR)↗
NASA NTRS1981

Comparison of closed loop model with flight test results

An analytic technique capable of predicting the landing characteristics of proposed aircraft configurations in the early stages of design was developed. In this analysis, a linear pilot-aircraft closed loop model was evaluated using experimental data generated with the NT-33 variable stability in-flight simulator. The pilot dynamics are modeled as inner and outer servo loop closures around aircraft pitch attitude, and altitude rate-of-change respectively. The landing flare maneuver is of particular interest as recent experience with military and other highly augmented vehicles shows this task to be relatively demanding, and potentially a critical design point. A unique feature of the pilot model is the incorporation of an internal model of the pilot's desired flight path for the flare maneuver.

George, F. L.↗