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,261 records · Page 70

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↗

Mineral Dust Optical Properties for Remote Sensing and Global Modeling: A Review

Dust plays a key role in many Earth system processes and is ubiquitous in the Martian atmosphere. Various intensive field campaigns, laboratory analyses, space-based remote sensing missions, and global modeling efforts aim to characterize dust optical properties. This is a bountiful time for dust scientists, and yet the interpretation of retrievals and comparison to models remains complicated by various conflicting assumptions that are part of each algorithm. For example, the conversion of satellite radiance measurements into products like aerosol optical depth for model evaluation depends upon aerosol properties like particle size and shape that are often prescribed and not part of the retrieval. Conversely, the model calculation of aerosol optical depth often uses different assumptions. The goal of this review is to first document algorithmic assumptions by various satellite retrieval products and models, and identify where there is consistency and where there are differences. In general, the differences documented in this paper reflect uncertainties resulting from incomplete observational characterization of dust aerosols and limitations in our understanding. Second, we note what observations might reduce uncertainties in our knowledge and bring greater consistency to retrievals and models, allowing for a more rigorous and harmonious comparison. The lack of comprehensive and realistic shape models for dust is an outstanding issue, such that closure between forward modeling from particle refractive index, shape, and size and observed optical properties cannot be achieved. Limitations in the computational methods that must be applied to model scattering from complex shapes also makes accurate optical modeling for dust challenging. Field observations indicate the persistence of coarse and giant dust particles at higher altitudes and farther downwind from their source than previously expected. Remote sensing retrieval algorithms based on observations at visible wavelengths have limited sensitivity to these particles and generally do not consider them, although a recent product based on longwave radiances is encouraging. Current measurements of the refractive index of bulk dust and fundamental dust minerology components such as hematite vary widely, inhibiting attempts to represent the variability in dust optical properties and forcing, as expected from different major dust source regions on Earth that have varying mineralogical composition. Some remote sensing retrieval algorithms allow for limited refractive index variability in their inversion solutions through mixing with other fine mode aerosol models, or optimizing the single scattering albedo, but Earth system models surveyed for this paper assume a globally uniform, size-invariant refractive index. Although no Martian dust samples have yet been returned to Earth, remote sensing observations indicate that Martian dust is globally homogenous in composition, and a single spectral refractive index assumption has been widely adopted to represent Martian dust. The lack of comprehensive, statistically representative measurements of dust particle microphysical properties (size distribution, morphology, complex index of refraction spectra, internal structure heterogeneity), and the resulting optical properties, limits our ability to verify the fidelity of these assumptions. A chain of measurements is needed, ranging from characterizing individual dust mineralogy components (e.g., pure hematite and goethite) to in situ sampling of complex atmospheric aerosol mixtures. Such results could be applied to both remote sensing retrievals that characterize the optical properties of the total aerosol burden in the atmosphere from total radiance measurements, and to global models that represent the total aerosol burden in the atmosphere by building it up from the balance of individual aerosol sources and sinks.

Aeolian dust↗

A Summary of the 2023 Joint High-Speed Boundary Layer Transition Workshop

A joint workshop between the Air Force Office of Scientific Research, the Office of Naval Research, and the NASA Hypersonic Technology Project was held in June 2023 on the topic of high-speed boundary-layer transition. The purpose of this workshop was to gather subject matter experts in the areas of supersonic and hypersonic boundary-layer transition in the areas of geometries of interest, computational methods, facilities and diagnostics, flight experiments, flow control, and new approaches and to discuss the current state of the art in the field. As gaps and accomplishments were identified, roadmaps were also created to guide future investments in these areas. The top recommendations are (1) to define relevant common geometries of interest that can be worked by multiple research groups; (2) improve modeling of high-enthalpy flows involving chemistry and ablation, studies of receptivity of freestream disturbances and atmospheric particulates into the boundary layer, surface feature characterization and modeling of changes to geometry; (3) design and determine viable flow control methods; and (4) move towards incorporating boundary layer transition prediction into a standard practice for design of high-speed vehicles. To execute this last recommendation successfully, investments beyond fundamental (basic and early applied) research are needed.

boundary layer transition↗

Satellite Optical Remote Sensing of Clouds and Aerosols: From Particle Single-Scattering and Gaseous Absorption Through Radiative Transfer to Retrieval Products

Clouds and aerosols are fundamental regulators of Earth’s radiation budget and climate system, influencing both solar and terrestrial radiation through scattering, absorption, and emission processes. Accurate characterization of their physical and radiative properties from space requires a rigorous understanding of particle single-scattering, gaseous absorption, and radiative transfer in the atmosphere, as well as reliable inversion methods. This review synthesizes the physical foundations and algorithmic implementations of satellite-based passive optical remote sensing of clouds and aerosols, spanning the ultraviolet to thermal infrared spectral range. Beginning with electromagnetic scattering theory and state-of-the-art methods for computing single-scattering by nonspherical particles and computationally efficient methods for accounting for atmospheric absorption, we discuss the radiative transfer framework underpinning cloud and aerosol retrievals. The connection between single-scattering and multiple-scattering is rigorously formulated. We then summarize operational and research-grade retrieval techniques, including cloud masking and thermodynamic phase determination, CO₂ slicing for cloud-top pressure, the Nakajima-King shortwave bi-spectral, and infrared split-window approaches for cloud optical thickness and effective particle size, inversion algorithms for determining aerosol properties from multi-spectral and/or multi-angle radiometric and polarimetric measurements, and active-passive sensing synergy. Examples of the global cloud and aerosol climatologies are illustrated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging SpectroRadiometer (MISR). Furthermore, the unique strengths of active remote sensing techniques based on spaceborne lidar observations are briefly elaborated in the context of studying ice clouds composed of randomly and horizontally oriented ice crystals, which is a significant challenge for conventional passive remote sensing techniques. By connecting physical theory to practical retrievals, this review highlights both the maturity of current methodologies and the remaining challenges in reducing uncertainties in particle morphology, vertical structure, absorption, and aerosol-cloud interactions. Furthermore, the impact of artificial intelligence (AI) on atmospheric remote sensing is briefly addressed.

Aerosols↗