Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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,423 records · Page 79

Numerical Simulation of AFP Nip Point Temperature Prediction for Complex Geometries

Material placement at the ideal nip point temperature over complex surfaces with uniformity across the width of the compaction rollers results in optimized part properties for Automated Fiber Placement (AFP) processes. However, current AFP systems utilize heat control models and methodologies, based on multiple process parameters such as feed-rate and orientation, that are mostly open-loop. Here, infrared (IR) heater input is calibrated as a function of process parameters during machine qualification. This work presents a numerical simulation to predict arrayed-infrared (AIR) emitter radiation onto a substrate that includes view factor implementation, IR radiative heat flow calculation, energy rate balance, and a transient heat transfer model. The purpose of this numerical model is to predict nip point temperature on complex surfaces, serving as a baseline for a new arrayed-infrared (AIR) thermoset heater to improve AFP process control. It is anticipated that this simulation will accurately control the temperature for high-speed AFP layup of complex geometries. An anticipated result of an AIR heater system is that material calibration and testing will be reduced as temperature is instantaneously monitored and controlled. Therefore, temperature across the roller width will be uniform during placement of complex parts, independent of their geometry.

Xia, Kaishu↗

Conceptualizing the Impact of Dust Contaminated Infrared Radiances on Data Assimilation for Numerical Weather Prediction

Numerical weather prediction systems depend on Hyperspectral Infrared Sounder (HIS) data, yet the impacts of dust-contaminated HIS radiances on weather forecasts has not been quantified. To determine the impact of dust aerosol on HIS radiance assimilation, we use a modified radiance assimilation system employing a one-dimensional variational assimilation system (1DVAR) developed under the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Numerical Weather Prediction–Satellite Application Facility (NWP-SAF) project, which uses the Radiative Transfer for TOVS (RTTOV). Dust aerosol impacts on analyzed temperature and moisture fields are quantified using synthetic HIS observations from rawinsonde, Micropulse Lidar Network (MPLNET), and Aerosol Robotic Network (AERONET). Specifically, a unit dust aerosol optical depth (AOD) contamination at 550 nm can introduce larger than 2.4 and 8.6 K peak biases in analyzed temperature and dewpoint, respectively, over our test domain. We hypothesize that aerosol observations, or even possibly forecasts from aerosol predication models, may be used operationally to mitigate dust induced temperature and moisture analysis biases through forward radiative transfer modeling.

Jared W. Marquis↗

Numerical Predictor-Corrector Based Guidance Scheme For Aero-Gravity Assist at Titan For Enceladus Missions

Aero-gravity assist is a spacecraft maneuver than can enable a vehicle to insert into an orbit of one planetary body using the atmosphere and gravity field of a secondary body, reducing the required propellant mass relative to a fully-propulsive orbit insertion. The Fully Numerical Predictor-Corrector Aerocapture Guidance algorithm was modified to work with a direct force control blunt body vehicle, executing an aero-gravity assist at Titan to enter a Saturnian orbit and conduct fly-by at Enceladus. Additionally, a proportional-integral-derivative controller was implemented to command sideslip angles for control of the orbital inclination. Numerical simulation showed the developed guidance scheme was capable of minimizing the energy and inclination error at atmospheric exit, allowing Saturnian moon tour trajectories and Enceladus exploration at a small delta-V cost, on the order of 100 m/s, in the nominal case.

Daniel Louis Engel↗

Numerical Predictor-Corrector Based Guidance Scheme For Aero-Gravity Assist at Titan For Enceladus Missions

Aero-gravity assist is a spacecraft maneuver than can enable a vehicle to insert into an orbit of one planetary body using the atmosphere and gravity field of a secondary body, reducing the required propellant mass relative to a fully-propulsive orbit insertion. The Fully Numerical Predictor-Corrector Aerocapture Guidance algorithm was modified to work with a direct force control blunt body vehicle, executing an aero-gravity assist at Titan to enter a Saturnian orbit and conduct fly-by at Enceladus. Additionally, a proportional-integral-derivative controller was implemented to command sideslip angles for control of the orbital inclination. Numerical simulation showed the developed guidance scheme was capable of minimizing the energy and inclination error at atmospheric exit, allowing Saturnian moon tour trajectories and Enceladus exploration at a small delta-V cost, on the order of 100 m/s, in the nominal case.

Daniel L. Engel↗

Turbulence, Transition, and Numerical Method Technologies

The close collaboration in the validation experimentation effort was excellent, the complementary research efforts related tofundamental numerical model development, and the activities related to effective HPC utilization on near term architectures that are coming down the pipeline were exactly what NASA should be doing. Even though the TQR panel cannot say that the success criteria were fully met, we do feel that the research funded by TTT/RCA was successful and useful in determining the next suite of CFD validation tests and the problematic flow conditions that continue to be difficult to predict using our current suite of computational tools. Given the TQR Panel’s review criteria: 1. The deliverables meet the technical requirements; 2.The appropriate technical approaches were followed in producing the deliverables; and 3. All technical caveats and concerns have been identified and the deliverables: a.Development of more accurate physics-based methods (e.g., higher moment closure); b. Large eddy simulation (LES); c. Advanced numerical methods; d. Transition prediction and modeling; e.Validation experiments; f. Multidisciplinary analysis and design (high fidelity). TQR panel feels that the technical processes were all completed and that the review criteria and all the deliverables were met. Thus, the Milestone TACP-2016-001 has been completed.

Mujeeb R Malik↗

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗

Numerical Modeling and Design of Gaseous Propellant Acquisition System for the Triton Hopper

Triton, the largest moon around Neptune, is a captured Kuiper Belt Object with a dynamic rarified atmosphere, active geysers, and unique “cantaloupe” like terrain. Little is known about Kuiper Belt Objects in the scientific community with most of the information coming from Pluto. Currently, NASA is designing an unmanned probe to explore the atmosphere and surface of Triton. This vehicle is unique in that it repeatedly takes off and lands (hops) using nitrogen propellant collected from its surroundings on Triton. This paper presents preliminary conceptual system design and numerical modeling of the atmospheric collection process. Four concepts are presented: a cryocooler-based system, a pump-based system, a sublimation and pump-based system, and a cryocooler and sublimation system. Numerical modeling takes into account the mass and energy transfer between all the relevant nodes in the system and is able to calculate the amount of nitrogen ice collected and ice accumulation rate as a function of time for each concept. The models are ultimately used to assess the feasibility of traversing from equator to pole in under two years using a gaseous propellant acquisition system for the Hopper.

Hopper↗

Conceptualizing the Impact of Dust-Contaminated Infrared Radiances on Data Assimilation for Numerical Weather Prediction

Numerical weather prediction systems depend on Hyperspectral Infrared Sounder (HIS) data, yet the impacts of dust-contaminated HIS radiances on weather forecasts has not been quantified. To determine the impact of dust aerosol on HIS radiance assimilation, we use a modified radiance assimilation system employing a one-dimensional variational assimilation system (1DVAR) developed under the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Numerical Weather Prediction–Satellite Application Facility (NWP-SAF) project, which uses the Radiative Transfer for TOVS (RTTOV). Dust aerosol impacts on analyzed temperature and moisture fields are quantified using synthetic HIS observations from rawinsonde, Micropulse Lidar Network (MPLNET), and Aerosol Robotic Network (AERONET). Specifically, a unit dust aerosol optical depth (AOD) contamination at 550 nm can introduce larger than 2.4 and 8.6 K peak biases in analyzed temperature and dewpoint, respectively, over our test domain. We hypothesize that aerosol observations, or even possibly forecasts from aerosol predication models, may be used operationally to mitigate dust induced temperature and moisture analysis biases through forward radiative transfer modeling.

Jared W Marquis↗

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗

Numerical Fluxes and Shocks

Accurately describing shocks within numerical partial differential equation solvers has been a significant computational hurdle for the field of applied mathematics. In our work we present and discuss the background and history of how shocks have been dealt with. We have analyzed numerical schemes and fluxes for the purpose of education and improved computational speed and accuracy. We demonstrate the need for specific methods to deal with shocks by showcasing some methods in a simple one dimensional case for burger’s equation. We aim to extend our knowledge to the euler equations.

Partial Differential Equations↗

Numerical Prediction of Heat Transfer Coefficients during Xenon Tank Fill in Microgravity

NASA’s Gateway will serve as an orbital outpost to enable sustained human lunar exploration and facilitate access to destinations beyond the Moon. To maintain sustainable deep space missions, on-orbit propellant transfer and refueling must be realized. One of Gateway’s critical components, the Power and Propulsion Element (PPE), will use solar electric propulsion to perform attitude control maneuvers and orbit transfers. Xenon will be used as the propellant for the PPE and stored in a composite overwrapped pressure vessel (COPV). During the propellant refueling process, xenon will compress causing the fluid and tank walls to warm. There is concern that the bond between the tank liner and composite may degrade at elevated temperatures, along with the possibility of the tank exceeding design pressure. While on-orbit testing is being proposed, numerical models are being used to predict the thermal response and reduce risk during a refueling operation in microgravity where natural convection is diminished. Accurate numerical models may be used to inform decisions for follow-on testing, or design and operation of refueling architectures. Challenges of modeling xenon tank fill in microgravity include non-ideal gas behavior, operating near the fluid’s critical point, and lack of experimental flight data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict averaged heat transfer coefficients during tank filling. At the anticipated injection flow rates, the CFD results showed that forced convection dominates in microgravity. Initial results from a multi-node Thermal Desktop model that assumed a purely conducting fluid (Nu=1) over-predicted the CFD fluid temperature and pressure, while a similar model allowing for natural convection under-predicted. The heat transfer coefficients predicted by CFD were implemented into the multi-node model, and the COPV thermal response from the two computational approaches were compared with good agreement.

Heat Transfer↗

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗

Development and Implementation of Large-Scale Numerical Models for Shape Memory Mars Spring Tires

The current work focuses on implementation of the user-defined shape memory alloy (SMA) model in the finite element analysis program ABAQUS for large-scale simulations of the Spring Tires made of SMA developed at the NASA Glenn Research Center. The main objective of this study was to improve and optimize the structural design of the SMA tires through in-depth numerical analysis and sensitivity studies. Various design variables (wire diam., coil diam., pitch, spring length, and bead angle) were varied to study their influence on the global load-displacement response of the tire construct. A detailed investigation of the three-dimensional stress states was also carried out to enhance understanding of the local changes as the tire goes through global deformation. It was concluded that a robust numerical model with a good predictive capability, together with a thoughtfully crafted sensitivity study, can result in reduced design iterations required to reach a desired tire performance. This will, in turn, lead to significantly reduced manufacturing time, required labor, and testing expense.

shape memory alloys↗

Comparing Planetary Helicopter Brownout Observations with Numerical Predictions for Ingenuity

Ingenuity is considered the Martian equivalent of the Wright Flyer. Its sustained flights on Mars proved that non-terrestrial powered flight is possible, even on planets with thin atmospheres. However, Ingenuity was small and carried no payload. Future planetary rotorcraft will be bigger and heavier to accomplish new scientific and exploratory goals. On dusty planets, the associated increase in rotor thrust poses a potential challenge well known to helicopter pilots who have flown in desert environments: dust clouds forming during take-off and landing due to rotor wash, referred to as brownout. This work evaluates the effectiveness of existing numerical simulation models in predicting the severity of this phenomenon for planetary rotorcraft. We compare the output of a CFD-based dust mobilization model with data recently gathered using image processing techniques of Ingenuity’s flights. The simulation results suggest that the sensitivity of the wall friction velocity to the rotorcraft height may be lower than previously assessed. We further evaluate the sensitivity of the total mass of mobilized dust to other model parameters such as the sandblasting efficiency and the saltation threshold velocity. Finally, we discuss the limitations of this comparison due to uncertainty in both the observation and the numerical model.

HEC↗

NUMERICAL NOISE STUDY IN EIC BEAM-BEAM SIMULATIONS

In the Electron-Ion Collider (EIC) design, a flat beam collision scheme is adopted to achieve 1034 cm -2 s -1 luminosity. We found that the vertical growth of the proton beam is much larger than of the round beam. In this article we present the numerical noise study about the number of macroparticles, the electron slice number, and the electron bunch length. Both weak-strong and strong-strong simulation methods are used. It turns out the proton emittance growth in the strong-strong simulation mainly comes from the numerical noise. This study helps us to perform beam-beam simulation correctly for the EIC.

43 PARTICLE ACCELERATORS↗

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗

A Numerical Modelling Study on Submerged Condensers for Heat Pump Water Heaters Using Low-GWP Refrigerants

Challenges to global warming require reduction of GHG emissions, and increased use of renewable energy. Heat pumps extract energy from renewable sources and convert them to more useable forms for space conditioning, space heating and for water heating. They can be powered by renewable sources of electricity and must be engineered to use low global warming potentials (GWP) refrigerants (e.g., propane, GWP = 3) for a sustainable economy. Therefore, an electric heat pump water heaters (HPWH) using propane as refrigerant are developed in Oak Ridge National Laboratory to replace water heaters using all electric resistance and natural gas. A submerged condenser directly contacting with water is employed in the HPWH, which eliminates heat transfer barrier of the tank wall, and also prevents heat loss from the condenser tubes to the surrounding air. In present work, a numerical model is developed to model the submerged condenser. A dynamic temperature distribution on the condenser tubes from experimental data is applied to the model to enhance the accuracy of the numerical model. The model is essential to guide the design of the submerged condenser to achieve the performance required by the HPWH.

Zhang, Mingkan↗

Porting numerical integration codes from CUDA to oneAPI: a case study

We present our experience in porting optimized CUDA implementations to oneAPI. We focus on the use case of numerical integration, particularly the CUDA implementations of PAGANI and $m$-Cubes. We faced several challenges that caused performance degradation in the oneAPI ports. These include differences in utilized registers per thread, compiler optimizations, and mappings of CUDA library calls to oneAPI equivalents. After addressing those challenges, we tested both the PAGANI and m-Cubes integrators on numerous integrands of various characteristics. To evaluate the quality of the ports, we collected performance metrics of the CUDA and oneAPI implementations on the Nvidia V100 GPU. We found that the oneAPI ports often achieve comparable performance to the CUDA versions, and that they are at most 10% slower.

97 MATHEMATICS AND COMPUTING↗