Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
NASA NTRS1977

A study of fatigue crack closure using electric potential and compliance techniques

The objective of this study was to compare closure data obtained on the same specimen with crack tip compliance gage and electric potential techniques. Equivalent results on closure using the two techniques were obtained on 7075-T651 aluminum center cracked panels. The results also indicated that closure was a function of stress ratio, specimen thickness and maximum applied stress intensity. Maximum stress intensity had a particularly strong effect on closure in the range of applied stresses used. This observed dependence of closure on specimen thickness and maximum stress intensity should account for many of the discrepancies in closure behavior found in the literature. Observations on potential system behavior during closure measurements are also reported.

Clarke, C. K.↗
NASA NTRS1995

Working toward a three-dimensional fatigue closure model for surface cracks

The first reliable elastic fracture mechanics solutions for a surface crack in a plate were obtained by Newman and Raju. The authors, both from the Mechanics of Materials Branch at NASA-Langley, used a highly detailed finite element solution requiring substantial computational resources. Computers have since become more powerful and available; however, many important related problems remain computationally expensive. The problem of three-dimensional fatigue crack growth taking into account plasticity-induced crack closure is one such problem. It is the goal of this research to provide an efficient method to account for three-dimensional crack closure in fatigue. Newman developed a two-dimensional plasticity-induced crack closure model for center cracked specimens. This model requires iterations to determine both the contact solution at each growth step and the extent of the plastic zone at the crack tip. A three-dimensional version of this model would require obtaining these nonlinear variables all along the crack front. This model must be efficient so that repeated calculations can be performed during crack growth simulations. The highly versatile line spring model (LSM) with contact, fatigue, and plasticity will form the basis of the closure model. There are several required additions to past work to address the three-dimensional crack closure problem. Initially, these additions will include (1) an improved LSM to more accurately obtain the crack opening displacement, stress intensity factors, and elastic T-stress near the ends of the surface crack; (2) a method to determine the extent of the plastic zone all along the crack front; (3) a method to determine the contact zone given a perfectly plastic layer of material on the crack surfaces; (4) a method to determine the magnitude of the compressive contact stress; and (5) a way to implement the degree of constraint along the curved crack front. During the summer ASEE program an enhanced LSM was developed. A method similar to that of 'strip synthesis' first introduced by Fujimoto was used. Briefly, the crack opening displacements of 'slices' of the surface crack in a direction parallel to the plate surface are considered in addition to the standard LSM approach that makes use of springs obtained from slices perpendicular to the plate surface. This enhancement is necessary so that an accurate three-dimensional representation of quantities such as contact zone size, plastic zone size, stress intensity factors, T-stress, and crack opening displacement can be determined. By combining results of previous investigations with the LSM, the problem of three-dimensional crack closure will be addressed. In addition to a closure model, the enhanced LSM can be used for many other problems including interacting surface cracks and fatigue crack growth of a through crack with a curved crack front.

Joseph, Paul F.↗
NASA NTRS2021

Aerosol–Ice Formation Closure: A Southern Great Plains Field Campaign

Prediction of ice formation in clouds presents one of the grand challenges in the atmospheric sciences. Immersion freezing initiated by ice-nucleating particles (INPs) is the dominant pathway of primary ice crystal formation in mixed-phase clouds, where supercooled water droplets and ice crystals coexist, with important implications for the hydrological cycle and climate. However, derivation of INP number concentrations from an ambient aerosol population in cloud-resolving and climate models remains highly uncertain. We conducted an aerosol-ice formation closure pilot study using a field-observational approach to evaluate the predictive capability of immersion freezing INPs. The closure study relies on co-located measurements of the ambient size-resolved and single-particle composition and INP number concentrations. The acquired particle data serve as input in several immersion freezing parameterizations, that are employed in cloud-resolving and climate models, for prediction of INP number concentrations. We discuss in detail one closure case study in which a front passed through the measurement site, resulting in a change of ambient particle and INP populations. We achieved closure in some circumstances within uncertainties, but we emphasize the need for freezing parameterization of potentially missing INP types and evaluation of the choice of parameterization to be employed. Overall, this closure pilot study aims to assess the level of parameter details and measurement strategies needed to achieve aerosol-ice formation closure. The closure approach is designed to accurately guide immersion freezing schemes in models, and ultimately identify the leading causes for climate model bias in INP predictions.

ice formation↗
DOE OSTI2026

Improbability of Post-Closure Criticality in Array of Transuranic Waste Containers after Compaction by Salt Creep at Waste Isolation Pilot Plant

Based on the rationale presented, post-closure nuclear criticality is improbable when room closure compacts containers disposing transuranic (TRU) waste emplaced at the Waste Isolation Pilot Plant (WIPP), an operating repository in bedded salt in southeastern New Mexico. As described in the original WIPP certification, a qualitative estimate of the probability of post-closure criticality in TRU waste produced from atomic energy defense activities has been low either because remote-handled TRU waste canisters are neutronically isolated by the bedded salt or because the low fissile mass in an array of contact-handled TRU waste drums cannot be compacted sufficiently by room closure from salt creep. These situations are still valid for the majority of TRU waste that is disposed at WIPP without any disposal constraints, as updated herein. This report also qualitatively evaluates the probability of criticality after disposal of TRU waste in pipe overpack containers (POC) where every POC in a shipment may have the maximum 200 fissile gram equivalent of 239Pu content. The probability of criticality for a disposal room filled with POCs is estimated during four representative phases of repository evolution: (1) a large salt block falls onto POCs in the first 20 years, (2) salt creep closes a disposal room in the first 1000 years without brine seepage and subsequent gas generation, which permits maximum compaction, (3) some brine seepage occurs into the closed room, which initiates consumption of the fiberboard (cellulose) impact absorber in the POC in the second 1000 years, and (4) full brine inundation of a room and consumption of all fiberboard thereafter. Salt-block fall in the first phase does not greatly disrupt three tiers of POCs. The compacted spacing of POCs in the later three repository conditions is calculated through high-fidelity, geomechanical modeling. Criticality evaluation of compacted 200-g 239 Pu spheres at the compacted spacing shows that neither 12-inch nor 6-inch POCs are critical after the first 1000 years, the second 1000 years, or thereafter as the sea of reflector material changes to represent the three repository conditions. Specifically, fiberboard and iron isolates 239 Pu while dry, and brine reduces the reactivity when a room is partially and fully inundated. Because POC behavior bounds behavior of other standard TRU waste containers, post-closure criticality caused by room closure compacting containers is omitted in the performance assessments for the 2019 and 2026 WIPP compliance re-certification applications to the US Environmental Protection Agency.

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

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

A study of fatigue crack closure using electric potential and compliance techniques

Compared are closure data produced on the same specimen by the crack tip compliance gage and electric potential techniques. Experiments on 7075-T651 aluminum center cracked panels produced equivalent results on closure using the two techniques. The results also indicated that closure is a function of stress ratio, specimen thickness and maximum applied stress intensity. Maximum stress intensity had a strong effect on closure in the range of applied stresses used. This dependence of closure on specimen thickness and maximum stress intensity accounts for many of the discrepencies in closure behavior reported in the literature.

Clarke, C. K.↗
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.↗
NASA NTRS1985

A comparison of three algebraic stress closures for combustor flow calculations

A comparison is made of the performance of two locally nonequilibrium and one equilibrium algebraic stress closures in calculating combustor flows. Effects of four different pressure-strain models on these closure models are also analyzed. The results show that the pressure-strain models have a much greater influence on the calculated mean velocity and turbulence field than the algebraic stress closures, and that the best mean strain model for the pressure-strain terms is that proposed by Launder, Reece and Rodi (1975). However, the equilibrium algebraic stress closure with the Rotta return-to-isotropy model (1951) for the pressure-strain terms gives as good a correlation with measurements as when the Launder et al. mean strain model is included in the pressure-strain model. Finally, comparison of the calculations with the standard k-epsilon closure results show that the algebraic stress closures are better suited for simple turbulent flow calculations.

Nikjooy, M.↗
NASA NTRS1991

Near-wall turbine closure for curved flows

At present, turbulence closures for curved flows only account for curvature effects in the fully turbulent region where the Reynolds number is large. The justification is that, near a wall, viscous effects dominate and curvature effects are only of secondary importance. Recent direct simulation data show that this assumption is not valid, even for simple two-dimensional fully developed turbulent curved channel flows. This paper presents an approach to develop a near-wall turbulence closure for wall-bounded curved flows. The curved channel direct simulation data is used as a guide to held develop such a closure. The proposed closure has the unique property of approaching conventional high-Reynolds-number Reynolds-stress closures far away from the wall. Hence, curvature effects in both the near-wall and the fully turbulent parts of the flow are accounted for properly. Validations of the closure are carried out with a set of low-Reynolds-number simulation data and with experimental measurements at high Reynolds number. Good agreement is obtained in both cases; in particular, the anisotropic behavior of the normal stresses and the shear stress behavior near the convex and concave walls.

So, R. M. C.↗
NASA NTRS1999

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗
NASA NTRS2003

Modeling the Interactions Between Multiple Crack Closure Mechanisms at Threshold

A fatigue crack closure model is developed that includes interactions between the three closure mechanisms most likely to occur at threshold; plasticity, roughness, and oxide. This model, herein referred to as the CROP model (for Closure, Roughness, Oxide, and Plasticity), also includes the effects of out-of plane cracking and multi-axial loading. These features make the CROP closure model uniquely suited for, but not limited to, threshold applications. Rough cracks are idealized here as two-dimensional sawtooths, whose geometry induces mixed-mode crack- tip stresses. Continuum mechanics and crack-tip dislocation concepts are combined to relate crack face displacements to crack-tip loads. Geometric criteria are used to determine closure loads from crack-face displacements. Finite element results, used to verify model predictions, provide critical information about the locations where crack closure occurs.

Newman, John A.↗
NASA NTRS1989

Quick actuating closure

A quick actuating closure for a pressure vessel 80 in which a wedge ring 30 with a conical outer surface 31 is moved forward to force shear blocks 40, with conical inner surfaces 41, radially outward to lock an end closure plug 70 within an opening 81 in the pressure vessel 80. A seal ring 60 and a preload ramp 50 sit between the shear blocks 40 and the end closure plug 70 to provide a backup sealing capability. Conical surfaces 44 and 55 of the preload ramp 50 and the shear blocks 40 interact to force the seal ring 60 into shoulders 73 and 85 in the end closure plug 70 and opening 81 to form a tight seal. The end closure plug 70 is unlocked by moving the wedge ring 30 rearward, which causes T-bars 32 of the wedge ring 30 riding within T -slots 42 of the shear blocks 40 to force them radially inward. The end closure plug 70 is then removed, allowing access to the interior of the pressure vessel 80.

White, III, Dorsey E.↗
DOE OSTI2026

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗
DOE OSTI2025

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗