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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗
NASA NTRS2015

Revisiting Turbulence Model Validation for High-Mach Number Axisymmetric Compression Corner Flows

Two axisymmetric shock-wave/boundary-layer interaction (SWBLI) cases are used to benchmark one- and two-equation Reynolds-averaged Navier-Stokes (RANS) turbulence models. This validation exercise was executed in the philosophy of the NASA Turbulence Modeling Resource and the AIAA Turbulence Model Benchmarking Working Group. Both SWBLI cases are from the experiments of Kussoy and Horstman for axisymmetric compression corner geometries with SWBLI inducing flares of 20 and 30 degrees, respectively. The freestream Mach number was approximately 7. The RANS closures examined are the Spalart-Allmaras one-equation model and the Menter family of kappa − omega two equation models including the Baseline and Shear Stress Transport formulations. The Wind-US and CFL3D RANS solvers are employed to simulate the SWBLI cases. Comparisons of RANS solutions to experimental data are made for a boundary layer survey plane just upstream of the SWBLI region. In the SWBLI region, comparisons of surface pressure and heat transfer are made. The effects of inflow modeling strategy, grid resolution, grid orthogonality, turbulent Prandtl number, and code-to-code variations are also addressed.

Shock Wave↗
NASA NTRSDate not supplied

Exploring Radiation, Water and Energy Closure in Key Deep Tropical Regions Over the Last 3 Decades

The closure of global energy and water budgets at various time scales depends heavily upon the on the processes of water and energy exchange in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance in key tropical regions is observed of water budget terms between observations and atmospheric reanalysis systems, implying that processes are not simulated well. For instance, periods of imbalance in the water budget may be at least partly attributed to precipitation estimates due to the transition between cloud and precipitation systems in which ice cloud affects microwave radar estimates (Lietmann-Niimi et al. 2024, submitted). A recent paper by Hsaio et al (2023, pre-published version) found that the transition process is characterized by wind shear and longwave cloud radiative feedback signatures as well. This is just one example how these processes also affect the radiative fluxes and thus the energy budget at the TOA. In this presentation, we assess the net radiative budget and components in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES EBAF/SYN1deg Ed4.2, ISSCP FH, etc.) and describing the temporal variability of energy balance terms across the deep tropics for the period from 1988 to near present in six key oceanic regions. These six key tropical oceanic regions are West Indian Ocean, the south Indian Ocean east, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean. In these regions we assess the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear to assess the uncertainties of energy closure in these regions. These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. These differences also affect closure terms and inferences on systematic uncertainties within the observational estimates. We conclude with recommendations for continuing work

Radiation budget↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2026

Large-scale calculations of 𝛽-decay rates and implications for 𝑟-process nucleosynthesis

Nuclear 𝛽 decay is a key element of the astrophysical rapid neutron capture process (𝑟 process). In this work, we present state-of-the-art global 𝛽-decay calculations based on the quantified relativistic nuclear energy density functional theory and the deformed proton-neutron quasiparticle random-phase approximation. Our analysis considers contributions from allowed and first-forbidden transitions. We used two point-coupling functionals with carefully calibrated time-odd terms and isoscalar pairing strength. The new calculations display consistent results for both employed functionals, especially near the neutron drip line, suggesting slower 𝛽 decays past the 𝑁=126 neutron shell closure than in commonly used 𝛽-decay models. The new rates, along with the existing rates based on the recent nonrelativistic global calculations, are found to slow down the synthesis of heavy elements in the 𝑟 process and significantly reduce the contribution of neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Characterization of an electrically conductive proppant for fracture diagnostics

Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.

Geochemistry & Geophysics↗
DOE OSTI2020

Solid Waste Release Calculations for the Composite Analysis Baseline Assessment

The purpose of the Hanford Site Composite Analysis (CA) is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The need for a waste form release model addressing multiple solid waste forms is documented in CP-60410, Hanford Site Composite Analysis Technical Approach Description: Waste Form Release.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1988

Theoretical study of turbulent channel flow: Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr--Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O--S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend stongly on the structure of the turbulence spectrun at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗
NASA NTRS1990

Theoretical study of turbulent channel flow - Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O-S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend strongly on the structure of the turbulence spectrum at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗
NASA NTRS1990

On near-wall turbulent flow modelling

The behavior of various terms in the Reynolds-stress transport equations is evaluated in view results of near-wall turbulence characterizations which show that all components of the velocity-pressure gradient correlation vanish at the wall. A methodology for deriving an asymptotically correct model of the velocity-pressure gradient correlation is proposed, and a model capable of approaching the high-Reynolds number model for pressure distribution far from the wall is derived. A parametric study of the model constants introduced by the near-wall closure reveals that one constant in the dissipation-rate equation is Reynolds number-dependent; with this modification, excellent agreement is obtained with measured and simulated near-wall turbulence statistics.

Lai, Y. G.↗
NASA NTRS1992

Long-term orbital evolution of short-period comets found in Project Cosmo-DICE

Orbital evolutions of about 160 short-period (SP) comets are numerically integrated for 4400 years in the framework of a realistic dynamical model. By the round-trip error in closure test, a reliable time space of the integrated orbits is estimated for each comet. Majority of the SP comets with their Tisserand's constant(J) between 2.8 and 3.1 are found to evolve within the past 1000-2000 years from the orbits whose perihelia are near the Jovian orbit to the orbits with perihelia of 1-2 AU. This evolution is much more rapid than that expected from Monte Carlo simulations based on symmetric distribution of planetary perturbations, thus suggesting that asymmetry of perturbation distribution play an important role in cometary evolution. Several comets are shown to evolve from the near-Saturn orbits and then to be handed over under the control of Jupiter. We also find that a few comets were captured from long-period orbits (a = 75-125 AU) via only a few close encounters with Jupiter. It is confirmed that the captured SP comets of low-inclination with 2.7 less than J less than 3.1 show more or less strong chaotic behavior. On the other hand, comets with longer orbital period and/or of high inclination reveal slow or quasi-periodic orbital evolution.

Nakamura, Tsuko↗
NASA NTRS1994

Rotorcraft control system design for uncertain vehicle dynamics using quantitative feedback theory

Quantitative Feedback Theory describes a frequency-domain technique for the design of multi-input, multi-output control systems which must meet time or frequency domain performance criteria when specified uncertainty exists in the linear description of the vehicle dynamics. This theory is applied to the design of the longitudinal flight control system for a linear model of the BO-105C rotorcraft. Uncertainty in the vehicle model is due to the variation in the vehicle dynamics over a range of airspeeds from 0-100 kts. For purposes of exposition, the vehicle description contains no rotor or actuator dynamics. The design example indicates the manner in which significant uncertainty exists in the vehicle model. The advantage of using a sequential loop closure technique to reduce the cost of feedback is demonstrated by example.

Hess, R. A.↗
NASA NTRS2000

Long-Term Global Morphology of Gravity Wave Activity Using UARS Data

Gravity waves in satellite data from CRISTA and MLS are studied in depth this quarter. Results this quarter are somewhat limited due to the PI'S heavy involvement throughout this reporting period in on-site forecasting of mountain wave-induced turbulence for the NASA's ER-2 research aircraft at Kiruna, Sweden during the SAGE Ill Ozone Loss and Validation Experiment (SOLVE). Results reported concentrate on further mesoscale modeling studies of mountain waves over the southern Andes, evident in CRISTA and MLS data. Two-dimensional mesoscale model simulations are extended through generalization of model equations to include both rotation and a first-order turbulence closure scheme. Results of three experiments are analyzed in depth and submitted for publication. We also commence simulations with a three-dimensional mesoscale model (MM5) and present preliminary results for the CRISTA 1 period near southern South America. Combination of ground-based temperature data at 87 km from two sites with global HRDl data was continued this quarter, showing stationary planetary wave structures. This work was also submitted for publication.

Eckermann, Stephen D.↗
NASA NTRS2017

Elastomeric Thermal Insulation Design Considerations in Long, Aluminized Solid Rocket Motors

An all-new sounding rocket was designed at NASA's Marshall Space Flight Center that featured an aft finocyl, aluminized solid propellant grain and silica-filled ethylene-propylene-diene monomer (SFEPDM) internal insulation. Upon the initial static firing of the first of this new design, the solid rocket motor (SRM) case failed thermally just upstream of the aft closure early in the burn time. Subsequent fluid modeling indicated that the high-velocity combustion-product jets emanating from the fin-slots in the propellant grain were likely inducing a strongly swirling flow, thus substantially increasing the severity of the convective environment on the exposed portion of the SFEPDM insulation in this region. The aft portion of the fin-slots in another of the motors were filled with propellant to eliminate the possibility of both direct jet impingement on the exposed SFEPDM and the appearance of strongly swirling flow in the aft region of the motor. When static-fired, this motor's case still failed in the same axial location, and, though somewhat later than for the first static firing, still in less than 1/3rd of the desired burn duration. These results indicate that the extreme material decomposition rates of the SFEPDM in this application are not due to gas-phase convection or shear but rather to interactions with burning aluminum or alumina slag. Further comparisons with between SFEPDM performance in this design and that in other hot-fire tests provide insight into the mechanisms of SFEPDM decomposition in SRM aft domes that can guide the upcoming redesign effort, as well as other future SRM designs. These data also highlight the current limitations of modeling elastomeric insulators solely with diffusion-controlled, gas-phase thermochemistry in SRM regions with significant viscous shear and/or condense-phase impingement or flow.

Martin, Heath T.↗
NASA NTRS2018

TPSAS-NF1676L-30451-DND

The fundamental laws of the physics of fluid flow in the continuum regime are embedded in a somewhat "simple" set of partial differential equations describing the transport of mass, momentum, energy, and/orentropy. These basic equations can be manipulated to construct other important physical quantities such as vorticity that help further elucidate flow physics. Additionally, these equations can be nondimensionalized to identify important nondimensional groups governing compressible, reacting fluid flows, most prominently Reynolds, Mach, Froude, Prandtl, Lewis, Schmidt, Damkohler, and Zeldovich numbers. In this lecture, we present these governing equations and discuss the significance of the different terms and nondimensional groups. Furthermore, since these equations cannot be solved directly on current supercomputers for turbulent flows at Reynolds numbers of engineering interest, a scale reduction by averaging or filtering is introduced, which leads to a "closure" problem. The closure is then provided by introducing some common models. The current introduction is by no means complete, but it and the references herein may serve as a good starting point for learning about the most common approaches to turbulent reacting flows.

Tomasz (Tom) G Drozda↗
NASA NTRSDate not supplied

Validation of Delamination Growth Predictions under Fatigue Loading using the Single Leg Bending Specimen

The Single Leg Bending Specimen (SLB) was chosen for a study to compare experimentally determined delamination growth with analysis predictions. First, fatigue tests were performed on SLB specimens with several different starter delamination lengths. Specimens were cycled in load control and delamination length increase was continuously recorded with increasing number of cycles. Test results from different load levels were combined to capture the entire delamination growth law (Paris Law) over a wide range of energy release rates. Second, a recently developed interpolation scheme was used to predict the Paris Law associated with the SLB specimen. This interpolation scheme allows the calculation of the growth rates for a given mixed-mode ratio from known single-mode growth laws. Third, finite element analyses were performed using two-dimensional models of the SLB specimen. New Virtual Crack Closure Technique based analysis tools for automated delamination growth analysis under cyclic loading in Abaqus were used for growth predictions. Fourth, the predictions from analyses were compared to experimental fatigue data. Remaining analysis shortcomings were highlighted.

Ronald Krueger↗
NASA NTRS2020

Mechanical twinning of monazite expels radiogenic Pb

Mechanical twins form by the simple shear of the crystal lattice during deformation In order to test the potential of narrow twins in monazite to record the timing of their formation, a ~1,700 Ma monazite grain deformed at ~980 Ma was investigated by electron backscattered diffraction (EBSD), transmission electron microscopy (TEM) and atom probe tomography (APT). APT 208Pb/232Th ages indicate that the twin was entirely reset by radiogenic Pb-loss during its formation at conditions far below the monazite closure temperature. The results are consistent with a model where Pb is liberated during rupture of REE-O bonds in the large REEO9 polyhedra during twinning. Liberated Pb is likely to migrate along fast diffusion pathways such as crystal defects. The combination of a quantitative microstructural investigation and nanogeochronology provides a new approach for understanding the history of accessory phases.

monazite↗
NASA NTRSDate not supplied

Numerical Analysis of Cavitation-Induced Two-Phase Fluid Transients in Cryogenic Feedlines Using Finite Volume Procedure

This article presents a numerical investigation of the cavitation induced fluid transient in cryogenic feedline due to sudden closure of the valve A one dimensional numerical model based on finite volume procedure consisting of nodes and branches using Generalized Fluid System Simulation Program ( has been employed for the prediction of pressure variations The unsteady friction formulation accurately determines wall shear stress, and a homogeneous mixture model is used for two phase flow modeling Numerical results obtained with the present formulation are compared with experimental results and Method of Characteristic ( based results It is observed that low residual in flow variable calculation improves the solution providing numerically stable results GFSSP is able to reproduce the fluid transient peak due to vapor collapse Subsequently, it can retain the shape of pressure wave for two phase fluid transients, but it is dampening faster compared to experiments Further investigation is required with set of experiments to validate the present numerical approach for different experimental data.

Fluid Transient↗