Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing methods”

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,783 records · Page 99

A History of Jet Noise Research at the National Aeronautics and Space Administration

This paper reviews jet noise research conducted at the National Aeronautics and Space Administration (NASA) from the early 1950s to the present day. Research conducted by NASA’s predecessor, the National Advisory Committee for Aeronautics (NACA), and early years of NASA focused on turbojet noise, where a common approach for reducing jet noise was to limit the jet exit velocity to speeds that provided acceptable noise levels. Suppressors tested during this time resulted in thrust losses that were too severe to be implemented. With the introduction of turbofan engines in the 1960s, NASA shifted research to programs for both subsonic and supersonic aircraft applications with specific noise reduction goals. Subsonic research focused on increasing the bypass ratio of the engine to reduce the jet exit velocity of the core exhaust and adding mixers to the dual exhaust streams. Advances in computational methods improved aerodynamic designs and jet noise prediction tools. Supersonic applications proved to be more troublesome as programs aimed at large commercial transports required higher specific thrust engines. Changing the engine cycle to reduce jet noise was not compatible with mission range and speed requirements. Research for supersonic commercial aircraft remains an area of interest today at NASA.

Jet Noise↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

Path-Adaptive Guidance Algorithm Trades for a Two-Stage Lunar Descent Vehicle

For the next generation of NASA’s missions, explicit, path-adaptive descent guidance algorithms must provide the stability and customizability required for a safe and efficient descent to the lunar surface, while also meeting program and vehicle constraints. Several descent algorithms have been flown and tested for single-stage landers through the Apollo and Altair programs, but thus far little analysis has been conducted involving the application of these algorithms to a two-stage descent vehicle. Due to payload mass and fairing constraints of the existing fleet of launch vehicles, multi-stage descent architectures are a unique option for achieving the greatest possible mass to lunar surface. This paper seeks to compare the performance of guidance configurations of a lunar lander system consisting of two stages, one of which separates partway through descent. Through development of this paper, an optimization suite has been written that is specifically designed for optimizing planetary non-atmospheric two-stage descent trajectories, and is used as a baseline to compare the guidance algorithms tested. Time-to-go computational methods and ignition logic routines that may be employed in a lunar environment are also discussed. Preliminary results are presented that show relative performance metrics for a range of different guidance algorithm configurations.

Jason M Everett↗

The Evolution of Deep Space Navigation: 2012–2014

The exploration of the planets of the solar system using robotic vehicles has been underway since the early 1960s. During this time the navigational capabili-ties employed have increased greatly in accuracy, as required by the scientific objectives of the missions and as enabled by improvements in technology. This paper is the seventh in a chronological sequence dealing with the evolution of deep space navigation. The time interval covered extends from 2012 to 2014. The paper focuses on the observational techniques that have been used to obtain navigational information, propellant-efficient means for modifying spacecraft trajectories, and the computational methods that have been employed, tracing their evolution through 13 planetary missions.

Wood, Lincoln J↗

The Evolution of Deep Space Navigation: 2014–2016

The exploration of the planets of the solar system using robotic vehicles has been underway since the early 1960s. During this time the navigational capabilities employed have increased greatly in accuracy, as required by the scientific objectives of the missions and as enabled by improvements in technology. This paper is the eighth in a chronological sequence dealing with the evolution of deep space navigation. The time interval covered extends from 2014 to 2016. The paper focuses on the observational techniques that have been used to obtain navigational information, propellant-efficient means for modifying spacecraft trajectories, and the computational methods that have been employed, tracing their evolution through 11 planetary missions.

Wood, Lincoln J↗

Packing and flow particle simulations

Granular material is present across natural and industrial processes on Earth and other planets. Granular particles show up in space exploration (lunar regolith), avalanches (boulders), food (coffee), construction (concrete powder) and manufacturing (additive manufacturing powder, battery slurries). Important phenomena emerge from large collections of granular particles. The size scales of the particles in granular material, 10 μm diameter or larger, makes particle-based simulations a tractable computational method. This seminar will present the packing and flow of granular matter using particle-based discrete element modeling simulations. Particles modeled with rotational friction only require as few as 2.6 contacts for mechanical stability, as opposed to 6 contacts for frictionless particles. Optimal parameters for in-space manufacturing particulate material are identified. Specifically, the maximum density and contacts occurs for a large-to-small particle volume ratio of 0.265. Stress and contact fabric fluctuations of flowing dry granular matter have power-law scaling with strain rate, and a kink. The pressure-dependent slope change kink could identify the transition between slower, quasistatic and faster, inertial flows.

granular↗

Packing and flow particle simulations

Granular material is present across natural and industrial processes on Earth and other planets. Granular particles show up in space exploration (lunar regolith), avalanches (boulders), food (coffee), construction (concrete powder) and manufacturing (additive manufacturing powder, battery slurries). Important phenomena emerge from large collections of granular particles. The size scales of the particles in granular material, 10 μm diameter or larger, makes particle-based simulations a tractable computational method. This seminar will present the packing and flow of granular matter using particle-based discrete element modeling simulations. Particles modeled with rotational friction only require as few as 2.6 contacts for mechanical stability, as opposed to 6 contacts for frictionless particles. Optimal parameters for in-space manufacturing particulate material are identified. Specifically, the maximum density and contacts occurs for a large-to-small particle volume ratio of 0.265. Stress and contact fabric fluctuations of flowing dry granular matter have power-law scaling with strain rate, and a kink. The pressure-dependent slope change kink could identify the transition between slower, quasistatic and faster, inertial flows.

granular↗

Mixed-element USM3D Contributions to the 4th AIAA High-Lift Prediction Workshop

This paper discusses results of the mixed-element USM3D (USM3D-ME) simulations performed for the 4th AIAA High-Lift Prediction Workshop. The workshop was separated into six Technical Focus Groups to investigate the impact of geometry modeling, grid, and computational methods for predicting high lift flows. This work was performed under the Fixed Grid RANS Technical Focus Group. The primary geometry selected for the workshop was the High-Lift Common Research Model. The performed simulations included a flap deflection study and both a grid refinement study and pitch sweep for the nominal flap deflection configuration. The results show that USM3D-ME RANS solutions, generally, tends to underpredict the lift coefficient and to predict a less negative pitching moment relative to the experimental data. The predicted drag coefficient values agree better with experiment for smaller angles of attack but were observed to be larger than experiment for the largest angle of attack simulated. The results of the grid refinement study demonstrated a lack of grid convergence for the provided grid family. The results of the grid refinement study are consistent with the submissions to the Fixed Grid and Mesh Adaptation Technical Focus Groups. Grid convergence for the provided grid family remains elusive for the international community. A 2D Multielement Airfoil configuration was included to enable a turbulence model verification study, which illustrated favorable agreement between USM3D-ME and the solutions provided by other flow solvers.

CFD↗

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗

Mixed-Element USM3D Contributions to the 4th AIAA High-Lift Prediction Workshop

This paper discusses results of the mixed-element USM3D (USM3D-ME) simulations performed for the 4th AIAA High-Lift Prediction Workshop. The workshop was separated into six Technical Focus Groups to investigate the impact of geometry modeling, grid, and computational methods for predicting high lift flows. This work was performed under the Fixed Grid RANS Technical Focus Group. The primary geometry selected for the workshop was the High-Lift Common Research Model. The performed simulations included a flap deflection study and both a grid refinement study and pitch sweep for the nominal flap deflection configuration. The results show that USM3D-ME RANS solutions, generally, tends to underpredict the lift coefficient and to predict a less negative pitching moment relative to the experimental data. The predicted drag coefficient values agree better with experiment for smaller angles of attack but were observed to be larger than experiment for the largest angle of attack simulated. The results of the grid refinement study demonstrated a lack of grid convergence for the provided grid family. The results of the grid refinement study are consistent with the submissions to the Fixed Grid and Mesh Adaptation Technical Focus Groups. Grid convergence for the provided grid family remains elusive for the international community. A 2D Multielement Airfoil configuration was included to enable a turbulence model verification study, which illustrated favorable agreement between USM3D-ME and the solutions provided by other flow solvers.

CFD↗

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

This paper presents preliminary analysis using a finite volume based computational method for network flow analysis for the simulation of chill and fill of a cryogenic tank. An axisymmetric two-dimensional model of a spherical tank was developed using a nodal approach. No Vent Filling (NVF) was accomplished using a Thermodynamic Vent System assisted injector where the injector was cooled by the cold gas generated by Joule-Thompson expansion of the liquid propellant. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around spray droplets, and c) condensation of vapor at the interface of the cooled injector and ullage. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data and one-dimensional (1D) model predictions. It was demonstrated that the axisymmetric 2D model provided a more realistic prediction than the 1D model.

Cryogenics↗

Overview of NASA Research Activities in Shock Layer Kinetics and Radiation

Shock layer radiation is an important heating mechanism for entry probes to most planetary destinations and re-entry to Earth from beyond low Earth orbit. Our understanding of shock layer radiation phenomena has improved tremendously over the last decade thanks to NASA investment in fundamental radiation research through Entry Systems Modeling’s (ESM) Shock Layer Kinetics and Radiation (SLKR) task. This talk will overview the recent activities within SLKR, including validation of models through ground testing, flight instrumentation and remote observations; development of first principles ab initio calculations for reaction mechanisms and spectroscopic databases and advancing numerical and computational methods for prediction of radiation in reacting hypersonic flows. 1 Sr. Research Scientist, Aerothermodynamics Branch, and AIAA Associate Fellow.

Brett A. Cruden↗

Decoding the effects of synonymous variants

Synonymous single nucleotide variants (sSNVs) are common in the human genome but are often overlooked. However, sSNVs can have significant biological impact and may lead to disease. Existing computational methods for evaluating the effect of sSNVs suffer from the lack of gold-standard training/evaluation data and exhibit over-reliance on sequence conservation signals. We developed synVep (synonymous Variant effect predictor), a machine learning-based method that overcomes both of these limitations. Our training data was a combination of variants reported by gnomAD (observed) and those unreported, but possible in the human genome (generated). We used positive-unlabeled learning to purify the generated variant set of any likely unobservable variants. We then trained two sequential extreme gradient boosting models to identify subsets of the remaining variants putatively enriched and depleted in effect. Our method attained 90% precision/recall on a previously unseen set of variants. Furthermore, although synVep does not explicitly use conservation, its scores correlated with evolutionary distances between orthologs in cross-species variation analysis. synVep was also able to differentiate pathogenic vs. benign variants, as well as splice-site disrupting variants (SDV) vs. non-SDVs. Thus, synVep provides an important improvement in annotation of sSNVs, allowing users to focus on variants that most likely harbor effects.

Zishuo Zeng↗

The Evolution of Deep Space Navigation: 2016-2018*

The exploration of the planets of the solar system using robotic vehicles has been underway since the early 1960s. During this time the navigational capabilities employed have increased greatly in accuracy, as required by the scientific objectives of the missions and as enabled by improvements in technology. This paper is the ninth in a chronological sequence dealing with the evolution of deep space navigation. The time interval covered extends from 2016 to 2018. The paper focuses on the observational techniques that have been used to obtain navigational information, propellant-efficient means for modifying spacecraft trajectories, and the computational methods that have been employed, tracing their evolution through 12 planetary missions.

Wood, Lincoln J.↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗

Path-Adaptive Guidance Algorithm Trades for a Two-Stage Lunar Descent Vehicle

For the next generation of NASA’s missions, the necessity for path-adaptive guidance algorithms has become clear in order to provide the stability and customizability required for a safe and efficient descent to the lunar surface, while meeting specified program and vehicle constraints. Several descent algorithms have been tested and flown for single-stage landers through the Apollo and Altair programs, but thus far little analysis has been conducted involving the application of these algorithms for a two-stage descent vehicle. Due to the limited payload mass constraints of the existing fleet of launch vehicles, multi-stage descent architectures have become a viable course of action. This paper seeks to compare the performance of guidance configurations for a lunar lander system consisting of two stages, one of which separates partway through descent. Through development of this paper, an optimization suite has been written that is specifically designed for optimizing planetary non-atmospheric two-stage descent trajectories, and is used as a comparison baseline for the guidance algorithms tested. Time-to -go computational methods and ignition logic routines that may be employed in a lunar environment are also discussed. Further work is to be completed on trajectory design trades as well as the effects of modifying guidance targets in simulation based on trajectories that are optimized for different performance indices.

Jason M Everett↗

In Pursuit of CFD-based Wind Tunnel Calibrations

Computational fluid dynamic simulations of models tested in wind tunnels require a high level of fidelity and accuracy, particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure a sufficient characterization of both the physical geometry of the wind tunnel, the thermodynamics of the tunnel, and flow conditions in the test section. The condition setting process among subsonic wind tunnels, closed or open throat, is generally the same, though they can differ in specific details of the thermodynamics. The derivation, application, and error estimation of condition setting and calibration is discussed. Computational flow solutions of the high-speed leg of the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel are used to evaluate the calibration process and discuss calibration curve sensitivities. Duplicating the characteristics of a particular wind tunnel is difficult at best. A calibration curve derived from the computational method is the most consistent method to use for tunnel condition setting.

wind tunnel↗

Ground Tests on the BOLT Geometry at Mach 7: Cross-Facility Comparison and Stability Analysis

Transition measurements have been obtained through two experimental campaigns conducted independently by the German Aerospace Center and the French Aerospace Lab with the French Alternative Energies and Atomic Energy Commission on subscale models of the BOLT-1 flight experiment geometry. This paper details a cross-facility comparison of measurements obtained at Mach 7, as well as subsequent computational analysis. Infrared (IR) thermography measurements obtained by both campaigns have facilitated a global comparison of the transition front across facilities at analagous conditions, which are found to be in good agreement. High-frequency surface pressure fluctuation data demonstrate significant amplification of instabilities with Mack-mode characteristics in the outboard regions of the acreage. These measurements are compared to stability analyses of varying fidelity. The computational methods employed to characterize the boundary-layer transition phenomena include the traditional line-marching implementation of the Parabolized Stability Equations (PSE), 2D eigenvalue analysis coupled with PSE, and Adaptive Mesh Refinement Wavepacket Tracking. While the 2D eigenvalue analysis is found to predict instabilities, which correlate in terms of frequency and acreage location to the experimental measurements, the predicted amplification for these instabilities is lower than would typically be expected for transition. Line-marching results for traveling crossflow produce the best match to the experimental transition front, with a consistent transition N factor of approximately 3-3.5.

Boundary Layer Transition↗