Engineering PapersSearch

SEARCH · Engineering Papers

Results for “formulation generation”

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

A Diffusion‐Based Uncertainty Quantification Method to Advance E3SM Land Model Calibration

Abstract Calibrating land surface models and accurately quantifying their uncertainty are crucial for improving the reliability of simulations of complex environmental processes. This, in turn, advances our predictive understanding of ecosystems and supports climate‐resilient decision‐making. Traditional calibration methods, however, face challenges of high computational costs and difficulties in accurately quantifying parameter uncertainties. To address these issues, we develop a diffusion‐based uncertainty quantification (DBUQ) method. Unlike conventional generative diffusion methods, which are computationally expensive and memory‐intensive, DBUQ innovates by formulating a parameterized generative model and approximates this model through supervised learning, which enables quick generation of parameter posterior samples to quantify its uncertainty. DBUQ is effective, efficient, and general‐purpose, making it suitable for site‐specific ecosystem model calibration and broadly applicable for parameter uncertainty quantification across various earth system models. In this study, we applied DBUQ to calibrate the Energy Exascale Earth System Model land model at the Missouri Ozark AmeriFlux forest site. Results indicated that DBUQ produced accurate parameter posterior distributions similar to those from Markov Chain Monte Carlo sampling but with 30 times less computing time. This significant improvement in efficiency suggests that DBUQ can enable rapid, site‐level model calibration at a global scale, enhancing our predictive understanding of climate impacts on terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES

Orbiter/carrier separation for the ALT free flight no. 1 reference trajectories. Mission planning, mission analysis and software formulation

Details of the generation of the separation trajectories are discussed. The analysis culminated in definition of separation trajectories between physical separation and orbiter/carrier vortex clearance. Specifications, assumptions and analytical approach used to generate the separation trajectories are presented. Results of the analytical approach are evaluated. Conclusions and recommendations are summarized. Supporting references are listed.

Glenn, G. M.

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS

Modelling the Shuttle Remote Manipulator System: Another flexible model

High fidelity elastic system modeling algorithms are discussed. The particular system studied is the Space Shuttle Remote Manipulator System (RMS) undergoing full articulated motion. The model incorporates flexibility via a methodology the author has been developing. The technique is based in variational principles, so rigorous boundary condition generation and weak formulations for the associated partial differential equations are realized, yet the analyst need not integrate by parts. The methodology is formulated using vector-dyad notation with minimal use of tensor notation, therefore the technique is believed to be affable to practicing engineers. The objectives of this work are as follows: (1) determine the efficacy of the modeling method; and (2) determine if the method affords an analyst advantages in the overall modeling and simulation task. Generated out of necessity were Mathematica algorithms that quasi-automate the modeling procedure and simulation development. The project was divided into sections as follows: (1) model development of a simplified manipulator; (2) model development of the full-freedom RMS including a flexible movable base on a six degree of freedom orbiter (a rigid-body is attached to the manipulator end-effector); (3) simulation development for item 2; and (4) comparison to the currently used model of the flexible RMS in the Structures and Mechanics Division of NASA JSC. At the time of the writing of this report, items 3 and 4 above were not complete.

Barhorst, Alan A.

An Optimization-Based Coupling of Reduced Order Models with an Efficient Reduced Adjoint Basis Generation Approach

Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.

coupled problems

A Preliminary Study on the Vapor/Mist Phase Lubrication of a Spur Gearbox

Organophosphates have been the primary compounds used in vapor/mist phase lubrication studies involving ferrous bearing material. Experimental results have indicated that the initial formation of an iron phosphate film on a rubbing ferrous surface, followed by the growth (by cationic diffusion) of a lubricious pyrophosphate-type coating over the iron phosphate, is the reason organophosphates work well as vapor/mist phase lubricants. Recent work, however, has shown that this mechanism leads to the depletion of surface iron atoms and to eventual lubrication failure. A new organophosphate formulation was developed which circumvents surface iron depletion. This formulation was tested by generating an iron phosphate coating on an aluminum surface. The new formulation was then used to vapor/mist phase lubricate a spur gearbox in a preliminary study.

Morales, Wilfredo

Overview of Proposed ISRU Technology Development

ISRU involves any hardware or operation that harnesses and utilizes in-situ resources (natural and discarded) to create products and services for robotic and human exploration: Assessment of physical, mineral chemical, and volatile water resources, terrain, geology, and environment (orbital and local). Production of replacement parts, complex products, machines, and integrated systems from feedstock derived from one or more processed resources. Civil engineering, infrastructure emplacement, and structure construction using materials produced from in situ resources. Radiation shields, landing pads, roads, berms, habitats, etc. Generation and storage of electrical, thermal, and chemical energy with in situ derived materials. Solar arrays, thermal wadis, chemical batteries, etc. ISRU is a disruptive capability: Enables more affordable exploration than todays paradigm. Allows more sustainable architectures to be developed. Understand the ripple effect in the other Exploration Elements: MAV: propellant selection, higher rendezvous altitude (higher DV capable with ISRU propellants). EDL: significantly reduces required landed mass. Life Support: reduce amount of ECLSS closure, reduce trash mass carried through propulsive maneuvers. Power: ISRU drives electrical requirements, reactant and regeneration for fuel cells for landers, rovers, and habitat backup. Every Exploration Element except ISRU has some flight heritage (power, propulsion, habitats, landers, life support, etc.) ISRU will require a flight demonstration mission on Mars before it will be included in the critical path. Mission needs to be concluded at least 10 years before first human landed mission to ensure lessons learned can be incorporated into final design. ISRU Formulation team has generated a (still incomplete) list of over 75 technical questions on more than 40 components and subsystems that need to be answered before the right ISRU system will be ready for this flight demo.

Mars exploration

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING

Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations

Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.

42 ENGINEERING

Stochastic sea state for SRB studies

Ocean surface characteristics at two locations were studied for a space shuttle solid rocket booster ocean impact and recovery analysis. Probability distributions of wave heights, wave slopes, wave velocities, ocean currents, and 1-kilometer altitude winds were formulated. Procedures for generating ocean wave simulations are described.

Perlmutter, M.

Application of second-order turbulent modeling to the prediction of radiated aerodynamic sound

The Ribner formulation of the generation of aerodynamic sound is coupled with predictions of second-order velocity correlations and integral scale to estimate the sound radiated from several complicated jet flows. In particular, it is shown that the sound radiated from a cold swirling jet is greater than from its nonswirling equal thrust counterpart. The noise radiated from the flow field of a multitube suppressor was estimated and compared with an equal thrust diameter Gaussian jet. It is shown that the multitube concept is indeed quieter.

Bilanin, A. J.

The CCSDS Lossless Data Compression Algorithm for Space Applications

In the late 80's, when the author started working at the Goddard Space Flight Center (GSFC) for the National Aeronautics and Space Administration (NASA), several scientists there were in the process of formulating the next generation of Earth viewing science instruments, the Moderate Resolution Imaging Spectroradiometer (MODIS). The instrument would have over thirty spectral bands and would transmit enormous data through the communications channel. This was when the author was assigned the task of investigating lossless compression algorithms for space implementation to compress science data in order to reduce the requirement on bandwidth and storage.

Yeh, Pen-Shu

Ambipolar Electric Field, Photoelectrons, and Their Role in Atmospheric Escape From Hot Jupiters

Atmospheric mass loss from Hot Jupiters can be large due to the close proximity of these planets to their host star and the strong radiation the planetary atmosphere receives. On Earth, a major contribution to the acceleration of atmospheric ions comes from the vertical separation of ions and electrons, and the generation of the ambipolar electric field. This process, known as the "polar wind," is responsible for the transport of ionospheric constituents to Earth's magnetosphere, where they are well observed. The polar wind can also be enhanced by a relatively small fraction of super-thermal electrons (photoelectrons) generated by photoionization.We formulate a simplified calculation of the effect of the ambipolar electric field and the photoelectrons on the ion scale height in a generalized manner. We find that the ion scale height can be increased by a factor of 2-15 due to the polar wind effects. We also estimate a lower limit of an order of magnitude increase of the ion density and the atmospheric mass-loss rate when polar wind effects are included.

atmospheres

Second Generation Prototype Design and Testing for a High Altitude Venus Balloon

This paper describes the development of a second generation prototype balloon intended for flight in the upper atmosphere of Venus. The design of this new prototype incorporates lessons learned from the construction and testing of the first generation prototype, including finite element analyses of the balloon stresses and deformations, measured leak performance after handling and packaging, permeability and optical property measurements on material samples, and sulfuric acid testing. An improved design for the second generation prototype was formulated based on these results, although the spherical shape and 5.5 m diameter size were retained. The resulting balloon has a volume of 87 cubic meters and is capable of carrying a 45 kg payload at a 55 km altitude at Venus. The design and fabrication of the new prototype is described, along with test data for inflation and leakage performance.

mobility

The electron Boltzmann equation in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material show that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux but increases sharply in the presence of a sustainer electric field.

Hassan, H. A.

Electron distribution function in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material shows that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux.

Hassan, H. A.

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor

Generating Global Leaf Area Index from Landsat: Algorithm Formulation and Demonstration

This paper summarizes the implementation of a physically based algorithm for the retrieval of vegetation green Leaf Area Index (LAI) from Landsat surface reflectance data. The algorithm is based on the canopy spectral invariants theory and provides a computationally efficient way of parameterizing the Bidirectional Reflectance Factor (BRF) as a function of spatial resolution and wavelength. LAI retrievals from the application of this algorithm to aggregated Landsat surface reflectances are consistent with those of MODIS for homogeneous sites represented by different herbaceous and forest cover types. Example results illustrating the physics and performance of the algorithm suggest three key factors that influence the LAI retrieval process: 1) the atmospheric correction procedures to estimate surface reflectances; 2) the proximity of Landsatobserved surface reflectance and corresponding reflectances as characterized by the model simulation; and 3) the quality of the input land cover type in accurately delineating pure vegetated components as opposed to mixed pixels. Accounting for these factors, a pilot implementation of the LAI retrieval algorithm was demonstrated for the state of California utilizing the Global Land Survey (GLS) 2005 Landsat data archive. In a separate exercise, the performance of the LAI algorithm over California was evaluated by using the short-wave infrared band in addition to the red and near-infrared bands. Results show that the algorithm, while ingesting the short-wave infrared band, has the ability to delineate open canopies with understory effects and may provide useful information compared to a more traditional two-band retrieval. Future research will involve implementation of this algorithm at continental scales and a validation exercise will be performed in evaluating the accuracy of the 30-m LAI products at several field sites. ©

Landsat