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 1,207 records · Page 67

Gravitational Effects on Distortion in Sintering

During sintering a powder compact gains strength through low-temperature interparticle bonding, usually induced by solid-state surface diffusion, followed by further strength contributions from high-temperature densification. In cases where a liquid phase forms, sintering densification is accelerated and shape retention is sustained while open pores remain and contribute capillary forces. Unfortunately, sintering densification requires the compact become thermally softened to a point where creep strain rates reach levels near 10(exp -2)/s when the liquid forms. On the other hand, thermal softening of the powder compact substantially reduces the strength at high temperatures. Therefore, the in situ strength evolution during sintering is a primary focus to separate compact densification (as required for high performance) with minimized distortion (as required for net-shaping). With respect to gravitation effects on distortion during sintering there are two points of substantial weakness - prior to significant interparticle bonding and during final pore closure. This research is focused on understanding the competition among interparticle neck growth, densification, thermal softening, grain boundary wetting, capillary effects associated with liquid wetting and residual porosity, and gravity. Most surprising is the apparent role of gravity, where the deviatoric stress acting on the powder structure induces skeletal formation that reduces distortion. In contrast with theory, microgravity samples exhibit more distortion yet fail to fully densify. Results are presented on the experimental concepts supporting an emerging model of sintering strength evolution that enables understanding of both distortion and densification. The experiments have relied on tungsten heavy alloys, various combinations of dihedral angle, pore size, initial porosity, liquid:solid ratio, and heating rates. On Earth, the dominant factor with respect to distortion is the starting body heterogeneity. Current modeling efforts are seeking some means to uniformly predict the distortion based on a starting pore structure heterogeneity parameter. Densification is largely unaffected by the initial pore structure, but distortion is highly variable, suggesting that nonuniform pore closure might be a significant parameter during sintering. With respect to flight experiments, plans include removal of the solid body forces acting on the solid grains, allowing stabilization of the pore structure and examination of the buoyancy effects with regard to distortion. In microgravity there is the surprising result that compacts do not densify, yet distort a factor that is contrary to all current sintering models. Densification without distortion during liquid phase sintering was achieved by manipulating microstructure and its evolution during sintering. Microstructure parameters such as the solid volume fraction, dihedral angle, initial porosity, and pore size were varied to measure densification and distortion behavior during LPS using W-Ni-Cu alloys. Green compacts were formed using ethylene-bis-stearamide as a pore-forming agent with the amount of polymer controlling the initial porosity. Different initial pore sizes were generated by varying the polymer particle size. Dihedral angle was varied by changing the Ni:Cu ratio in the alloys. Finally, the solid volume fraction was adjusted via the tungsten content. Distortion was quantified using profiles determined with a coordinate measuring machine to calculate a distortion parameter. Sintering results showed that solid volume fraction and dihedral angle are the dominant factors on densification and distortion during liquid phase sintering. Distortion decreases with increasing solid volume fraction and dihedral angle, while initial porosity and pore size have no observable effect on distortion at nearly full densification. Various strategies emerge to improve distortion control in liquid phase sintering.

German, Randall M.↗

Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces

This paper develops an interpretable, non-intrusive reduced-order modeling technique using regularized kernel interpolation. Existing non-intrusive approaches approximate the dynamics of a reduced-order model (ROM) by solving a data-driven least-squares regression problem for low-dimensional matrix operators. Our approach instead leverages regularized kernel interpolation, which yields an optimal approximation of the ROM dynamics from a user-defined reproducing kernel Hilbert space. We show that our kernel-based approach can produce interpretable ROMs whose structure mirrors full-order model structure by embedding judiciously chosen feature maps into the kernel. The approach is flexible and allows a combination of informed structure through feature maps and closure terms via more general nonlinear terms in the kernel. We also derive a computable a posteriori error bound that combines standard error estimates for intrusive projection-based ROMs and kernel interpolants. In conclusion, the approach is demonstrated in several numerical experiments that include comparisons to operator inference using both proper orthogonal decomposition and quadratic manifold dimension reduction.

Data-driven model reduction↗

Effect of specific non-pharmaceutical intervention policies on SARS-CoV-2 transmission in the counties of the United States

Non-pharmaceutical interventions (NPIs) remain the only widely available tool for controlling the ongoing SARS-CoV-2 pandemic. We estimated weekly values of the effective basic reproductive number (R eff ) using a mechanistic metapopulation model and associated these with county-level characteristics and NPIs in the United States (US). Interventions that included school and leisure activities closure and nursing home visiting bans were all associated with a median R eff below 1 when combined with either stay at home orders (median R eff 0.97, 95% confidence interval (CI) 0.58–1.39) or face masks (median R eff 0.97, 95% CI 0.58–1.39). While direct causal effects of interventions remain unclear, our results suggest that relaxation of some NPIs will need to be counterbalanced by continuation and/or implementation of others.

60 APPLIED LIFE SCIENCES↗

β-Delayed One and Two Neutron Emission Probabilities Southeast of 132 Sn and the Odd-Even Systematics in r-Process Nuclide Abundances

Here, the β-delayed one- and two-neutron emission probabilities (P 1n and P 2n ) of 20 neutron-rich nuclei with N≥82 have been measured at the RIBF facility of the RIKEN Nishina Center. P 1n of 130,131 Ag, 133,134 Cd, 135,136 In, and 138,139 Sn were determined for the first time, and stringent upper limits were placed on P2n for nearly all cases. β-delayed two-neutron emission (β2n) was unambiguously identified in 133 Cd and 135,136 In, and their P 2n were measured. Weak β2n was also detected from 137,138 Sn. Our results highlight the effect of the N=82 and Z=50 shell closures on β-delayed neutron emission probability and provide stringent benchmarks for newly developed macroscopic-microscopic and self-consistent global models with the inclusion of a statistical treatment of neutron and γ emission. The impact of our measurements on r-process nucleosynthesis was studied in a neutron star merger scenario. Our P 1n and P 2n have a direct impact on the odd-even staggering of the final abundance, improving the agreement between calculated and observed Solar System abundances. The odd isotope fraction of Ba in r-process-enhanced (r-II) stars is also better reproduced using our new data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MELCOR Deposition Models Applied to Noble Metals

Noble metals are chemically inert elements that form particles that can circulate around the fuel loop of a liquid-fueled molten salt reactor (LF-MSR). Due to their health-hazardous abilities, it is important to track these elements and to understand their contributions to source term during a reaction, particularly for Te as it decays into iodine. This report describes three solid deposition models that were implemented into the MELCOR reactor modeling code. The models that were implemented were Dittus-Boelter for turbulent fluid deposition, Talbot thermophoretic transport equation and Richardson-Zaki gravitational deposition. These closures were used to calculate the mass transfer coefficient and was compared to those reported for the Molten Salt Reactor Experiment (MSRE). The impacts of salt, operation and particle parameters on how they affect solid deposition are described, and found that solid deposition is strongly dependent on salt flow and viscosity, as well as particle diffusion rate and size. The conclusions form the basis for activities to be carried out in FY26 and for pursuit of future activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A ray tracing model of gravity wave propagation and breakdown in the middle atmosphere

Gravity wave ray tracing and wave packet theory is used to parameterize wave breaking in the mesosphere. Rays are tracked by solving the group velocity equations, and the interaction with the basic state is determined by considering the evolution of the packet wave action density. The ray tracing approach has a number of advantages over the steady state parameterization as the effects of gravity wave focussing and refraction, local dissipation, and wave response to rapid changes in the mean flow are more realistically considered; however, if steady state conditions prevail, the method gives identical results. The ray tracing algorithm is tested using both interactive and noninteractive models of the basic state. In the interactive model, gravity wave interaction with the polar night jet on a beta-plane is considered. The algorithm produces realistic polar night jet closure for weak topographic forcing of gravity waves. Planetary scale waves forced by local transfer of wave action into the basic flow in turn transfer their wave action into the zonal mean flow. Highly refracted rays are also found not to contribute greatly to the climatology of the mesosphere, as their wave action is severely reduced by dissipation during their lateral travel.

Schoeberl, M. R.↗

Transient aerodynamic forces on a fighter model during simulated approach and landing with thrust reversers

Previous wind tunnel tests of fighter configurations have shown that thrust reverser jets can induce large, unsteady aerodynamic forces and moments during operation in ground proximity. This is a concern for STOL configurations using partial reversing to spoil the thrust while keeping the engine output near military (MIL) power during landing approach. A novel test technique to simulate approach and landing was developed under a cooperative Northrop/NASA/USAF program. The NASA LaRC Vortex Research Facility was used for the experiments in which a 7-percent F-18 model was moved horizontally at speeds of up to 100 feet per second over a ramp simulating an aircraft to ground rate of closure similar to a no-flare STOL approach and landing. This paper presents an analysis of data showing the effect of reverser jet orientation and jet dynamic pressure ratio on the transient forces for different angles of attack, and flap and horizontal tail deflection. It was found, for reverser jets acting parallel to the plane of symmetry, that the jets interacted strongly with the ground, starting approximately half a span above the ground board. Unsteady rolling moment transients, large enough to cause the probable upset of an aircraft, and strong normal force and pitching moment transients were measured. For jets directed 40 degrees outboard, the transients were similar to the jet-off case, implying only minor interaction.

Humphreys, A. P.↗

Contact dynamics testing of automated three-point docking mechanism

An evaluation of an OMV docking mechanism, based on an adaptation of the Shuttle Flight Support System Pallet Berthing Mechanism was completed. The mechanism uses automatically actuated motorized latches to engage towel bars on the target satellite. LED sensors establish the towel bar position within the capture envelope and the latch capture commands are issued. Then, locking pawls engage the bar, locking and pre-loading the mechanism. Two series of tests were conducted to test nominal and failure mode captures and to evaluate design parameters such as LED sensor locations, automatic closure algorithms, latch closure velocity, position/velocity entry envelopes, and closure method. The first test series involved single latch testing on the Flat Floor Facility, the 6 DOF Facility, and an analytic simulation model. The intent was to compare results in order to validate the various facilities. Reasonably good agreement was achieved. The second test series repeated the single latch testing on the refurbished 6 DOF Facility to validate the facility modifications. The individual latches were tested under free-drift conditions for functionality and performance. Next, the three-latch configuration underwent parametric testing. Test results validated the improved fidelity of the 6 DOF Facility and verified successful docking at the required entry velocity. The tests determined the 'best' design parameter definitions and concluded that the locking pawls should not lock until all three latches completely close.

Rourke, Kenneth H.↗

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗

Afterbody Heating Characteristics of a Proposed Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 deg sideslip, and normal-shock density, ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75% of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three-dimensional separated wake flows.

Horvath, Thomas J.↗

Pressurization System Modeling for a Generic Bimese Two- Stage-to-Orbit Reusable Launch Vehicle

A pressurization system model was developed for a generic bimese Two-Stage-to-orbit Reusable Launch Vehicle using a cross-feed system and operating with densified propellants. The model was based on the pressurization system model for a crossfeed subscale water test article and was validated with test data obtained from the test article. The model consists of the liquid oxygen and liquid hydrogen pressurization models, each made up of two submodels, Booster and Orbiter tank pressurization models. The tanks are controlled within a 0.2-psi band and pressurized on the ground with ambient helium and autogenously in flight with gaseous oxygen and gaseous hydrogen. A 15-psi pressure difference is maintained between the Booster and Orbiter tanks to ensure crossfeed check valve closure before Booster separation. The analysis uses an ascent trajectory generated for a generic bimese vehicle and a tank configuration based on the Space Shuttle External Tank. It determines the flow rates required to pressurize the tanks on the ground and in flight, and demonstrates the model's capability to analyze the pressurization system performance of a full-scale bimese vehicle with densified propellants.

Mazurkivich, Pete↗

A Well-Clear Volume Based on Time to Entry Point

A well-clear volume is a key component of NASA's Separation Assurance concept for the integration of UAS in the NAS. This paper proposes a mathematical definition of the well-clear volume that uses, in addition to distance thresholds, a time threshold based on time to entry point (TEP). The mathematical model that results from this definition is more conservative than other candidate definitions of the wellclear volume that are based on range over closure rate and time to closest point of approach.

Narkawicz, Anthony J.↗

Downwelling Far-Infrared Emission Spectra Measured by FIRST at Cerro Toco, Chile and Table Mountain Facility, CA

The Far-Infrared Spectroscopy of the Troposphere (FIRST) instrument is a fourier transform spectrometer developed to measure the far-infrared portion of the Earth's longwave spectrum. Specifically, the purpose is measurement of the important, yet under-investigated 100 to 1000 cm-1 band. Presented here are measurements made by FIRST at two successful deployments in a ground-based configuration to measure downwelling longwave radiation. The initial deployment was to Cerro Toco, Chile. FIRST operated from August to October, 2009 as part of the Radiative Heating in Underexplored Bands Campaign (RHUBC-II) campaign. After recalibration, FIRST was deployed to the Table Mountain Facility from August to October, 2012. Dry days for both campaigns were chosen for analysis ? 09/24/2009 and 10/19/2012. Also available from both deployments are coincident radiosonde temperature and water vapor vertical profiles used as inputs to radiance calculations performed by two platforms, the AER Line-by-Line Radiative Transform Model (using the AER line database), and a research platform, MRTA utilizing HITRAN 2012. Thus, we have three streams of data for the closure studies, and results of the comparisons are presented showing good agreement in the optically thick regions of 200-400 cm-1 and 600-800 cm-1. Differences between measured and coincident calculated spectra between 400 and 600 cm-1 are less than 10 mW/m^2/sr/cm-1. Uncertainties accounted for in this study include uncertainty from multiple detectors, water vapor, temperature, line strength and width as presented in HITRAN 2012, and due to internal heating caused by black bodies. The residuals fall within the root-sum-square envelope of the uncertainties.

J Mast↗

Hot Isostatic Pressing for Reduction of Defects in Additively Manufactured Inconel-718 Preliminary Simulation Approach and Results

Selective laser melting (SLM) of Inconel-718 (IN718) is a promising advanced manufacturing technology; however, challenges exist for fracture-critical parts where pores induced during the manufacturing process reduce fatigue life of the material. Hot isostatic pressing (HIP) is a common post processing technique used to reduce the part porosity by exposing parts to high temperature and pressures. Pore closure depends on the pore shape, trapped gas, and HIP parameters, and the process can only be optimized by understanding the effect of these variables. An analysis capability is described to better understand the conditions necessary for pore closure, and the effect of several variables on porosity are explored and compared to a baseline case. The importance of valid material response models is identified as a significant factor in model accuracy; especially the creep material response model which, in the baseline case, predicts that creep strain accounts for 87.7% of the total strain at the conclusion of the HIP cycle. Lessons learned are described, in addition to future work including model validation.

Joshua M Fody↗

Aerosol properties computed from aircraft-based observations during the ACE- Asia campaign: A case study of lidar ratio closure and aerosol radiative effects - 2

For a vertical profile with three distinct layers (marine boundary, pollution and dust), observed during the ACE-Asia campaign, we carried out a comparison between the modeled lidar ratio vertical profile and that obtained from collocated airborne NASA AATS-14 sunphotometer and shipborne Micro-Pulse Lidar (MPL) measurements. Vertically resolved lidar ratio was calculated from two size distribution vertical profiles - one obtained by inversion of sunphotometer-derived extinction spectra, and one measured in-situ - combined with the same refractive index model based on aerosol chemical composition. The aerosol model implies single scattering albedos of 0.78 - 0.81 and 0.93 - 0.96 at 0.523 microns (the wavelength of the lidar measurements), in the pollution and dust layers, respectively. The lidar ratios calculated from the two size distribution profiles have close values in the dust layer; they are however, significantly lower than the lidar ratios derived from combined lidar and sunphotometer measurements, most probably due to the use of a simple nonspherical model with a single particle shape in our calculations. In the pollution layer, the two size distribution profiles yield generally different lidar ratios. The retrieved size distributions yield a lidar ratio which is in better agreement with that derived from lidar/sunphotometer measurements in this layer, with still large differences at certain altitudes (the largest relative difference was 46%). We explain these differences by non-uniqueness of the result of the size distribution retrieval and lack of information on vertical variability of particle refractive index. Radiative transfer calculations for this profile showed significant atmospheric radiative forcing, which occurred mainly in the pollution layer. We demonstrate that if the extinction profile is known then information on the vertical structure of absorption and asymmetry parameter is not significant for estimating forcing at TOA and the surface, while it is of importance for estimating vertical profiles of radiative forcing and heating rates.

Kuzmanoski, Maja↗

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

Control and modeling of a CELSS (Controlled Ecological Life Support System)

Research topics that arise from the conceptualization of control for closed life support systems which are life support systems in which all or most of the mass is recycled are discussed. Modeling and control of uncertain and poorly defined systems, resource allocation in closed life support systems, and control structures or systems with delay and closure are emphasized.

Auslander, D. M.↗

The Harang discontinuity and magnetospheric forecasting

This paper describes the physics of the Harang discontinuity and its application to forecasting of the magnetospheric environment. Gradient/curvature depletion of energetic ions from earthward drifting flux tubes results in upward Birkeland currents. Ionospheric closure rotates the dawn-dusk convection electric field resulting in the Harang reversal. This physical mechanism has been incorporated into the Rice Convection Model (RCM) and can serve as a fundamental part of the forecasting of the magnetospheric system. The inclusion of this physics within the RCM results in a physical picture of the Harang discontinuity which did not previously exist. Model results are close to observations and represent a significant improvement in predicted electric field and Birkeland current patterns.

Erickson, G. M.↗