Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “significant quantity”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluation of DNA Extraction Efficiency in Diverse Algae Strains Using Commercial Kits and Lysis Approaches

Efficient DNA extraction is essential for accurately monitoring microalgae communities in large-scale cultivation systems such as raceway ponds and wastewater ponds. Traditional phenol chloroform extracts are a staple in microbiology but are obsolete for routine sampling due to its high toxicity reagents and time intensive setups. Commercial DNA extraction kits are more favorable for the microbes found in these ponds, but lack specific kits made for these communities. Little is known about which kits perform the best, leading researchers to use a variety of different kits with inconsistent results. This project compared one precipitation based commercial kit (Lucigen Masterpure) and five wash based kits (Monarch, Zymo Quick-DNA, and three Qiagen DNeasy kits) using four brackish algae strains to determine which methods yield the greatest quantity and quality of genomic DNA. Extractions were evaluated using the manufacturers protocol, and additional pretreatment options were administered before a single kit to compare its potential in being added routinely before extractions. Pretreatment options included both cryogenic freeze-thawing and heat incubation using enzymes. DNA was quantified using Qubit fluorometry and NanoDrop purity ratios. Overall, the Qiagen PowerWater kit provided the highest DNA yield and purity, but at a significantly higher cost then the precipitation-based kit (MasterPure). It was also noted that while the precipitation-based kit was significantly cheaper, provided similar results, it took significantly more time to complete a single run. Cryogenic pretreatment (6x cycles) increased average DNA yields by up to 80%, whereas enzymatic pretreatment most improved purity ratios without substantially improving quantity. The results suggest that it may be more cost and time efficient to use Qiagen kits with the addition of lysis pretreatments to procure better results. Future works includes developing a better system to efficiently collect multi variable data, and to upscale to artificial polycultures using similar methodologies alongside sequencing to confirm kit results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Implicit learning of convective organization explains precipitation stochasticity

Accurate prediction of precipitation intensity is crucial for both human and natural systems, especially in a warming climate more prone to extreme precipitation. Yet, climate models fail to accurately predict precipitation intensity, particularly extremes. One missing piece of information in traditional climate model parameterizations is subgrid-scale cloud structure and organization, which affects precipitation intensity and stochasticity at coarse resolution. Here, using global storm-resolving simulations and machine learning, we show that, by implicitly learning subgrid organization, we can accurately predict precipitation variability and stochasticity with a low-dimensional set of latent variables. Using a neural network to parameterize coarse-grained precipitation, we find that the overall behavior of precipitation is reasonably predictable using large-scale quantities only; however, the neural network cannot predict the variability of precipitation (R 2 ~ 0.45) and underestimates precipitation extremes. The performance is significantly improved when the network is informed by our organization metric, correctly predicting precipitation extremes and spatial variability (R 2 ~ 0.9). The organization metric is implicitly learned by training the algorithm on a high-resolution precipitable water field, encoding the degree of subgrid organization. The organization metric shows large hysteresis, emphasizing the role of memory created by subgrid-scale structures. We demonstrate that this organization metric can be predicted as a simple memory process from information available at the previous time steps. These findings stress the role of organization and memory in accurate prediction of precipitation intensity and extremes and the necessity of parameterizing subgrid-scale convective organization in climate models to better project future changes of water cycle and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Energy efficiency and renewable energy utilization in desalination systems

Desalination technologies and industry have advanced significantly in the past two decades to meet the growing freshwater demands stimulated by the compounding issues of both water quality and quantity in many regions of the world. As desalination processes are energy demanding, there have been many efforts dedicated to improve energy-efficiency of the process units, enhance energy conservation and recovery, and increase renewable energy integration in desalination plants. This review article highlights recent key advances and discusses possible venues for further development in desalination energy portfolio to reduce specific energy consumption and, via integration with solar energy, to minimize the environmental footprint associated with freshwater production. First, an overview of current desalination technologies and their energy requirements are presented followed by a discussion on opportunities for improving energy efficiency and energy recovery in both membrane and thermal desalination technologies. Then, various combinations of renewable energy driven desalination plants are discussed with some recent highlights in solar energy driven membrane, thermal and hybrid desalination processes. Technological readiness levels for novel desalination processes, their perceived impact and expected near-future developments in renewable energy integrated desalination technologies are presented. Lastly, the potential for solar driven desalination as a cost-competitive freshwater supply alternative is discussed

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1978

In-place recalibration technique applied to a capacitance-type system for measuring rotor blade tip clearance

The rotor blade tip clearance measurement system consists of a capacitance sensing probe with self contained tuning elements, a connecting coaxial cable, and remotely located electronics. Tests show that the accuracy of the system suffers from a strong dependence on probe tip temperature and humidity. A novel inplace recalibration technique was presented which partly overcomes this problem through a simple modification of the electronics that permits a scale factor correction. This technique, when applied to a commercial system significantly reduced errors under varying conditions of humidity and temperature. Equations were also found that characterize the important cable and probe design quantities.

Barranger, J. P.↗
NASA NTRS1996

Precision and Radiosonde Validation of Satellite Gridpoint Temperature Anomalies. Part I; MSU Channel 2: MSU Channel 2 - Pt. 1

In Part 1 of this study, monthly 2.5 deg. gridpoint anomalies in the TIROS-N (Television and Infrared Operational Satellite-N) series Microwave Sounding Unit (MSU) channel 2 brightness temperatures during 1979-88 are evaluated with multiple satellites and radiosonde data for their climate temperature monitoring capability. The MSU anomalies we computed about a 10-year mean annual cycle at each grid point, with the MSUs intercalibrated to a common arbitrary level. The intercalibrations remove relative biases between instruments of up to several tenths of a degree celsius. The monthly gridpoint anomaly agreement between concurrently operating satellites reveals single-satellite precision on generally better than 0.07 C in the tropics and better than 0.15 C at higher latitudes. Monthly anomalies in radiosonde channel 2 brightness temperatures computed with the radiative transfer equation compare very closely to the MSU measured anomalies in all climate zones, with correlations generally from 0.94 to 0.98 and standard errors of 0.15 C in the tropics to 0.30 C at high latitudes. Simplification of these radiative transfer calculations to a static weighting profile applied to the radiosonde temperature profile leads to an average degradation of only 0.02 deg. in the monthly skill. In terms of a more traditionally measured quantity, the MSU channel 2 anomalies match best with either the radiosonde 100-20-kPa or 100-15-kPa layer anomalies. No significant spurious trends were found in the 10-yr satellite dataset compared to the radiosondes that would indicate a calibration drift in either system. Thus, sequentially launched, overlapping passive microwave radiometers provide a useful system for monitoring intraseasonal to interannual climate anomalies and offer hope for monitoring of interdecadal trends from space. The Appendix includes previously unpublished details of the MSU gridpoint anomaly dataset construction. Part II of this study addresses the removal from channel 2 of the temperature influence above the 30-kPa level, providing a sharper and thus potentially more useful weighting function for monitoring lower tropospheric temperatures.

Spencer, Roy W.↗
NASA NTRS2000

A Common Calibration Source Framework for Fully-Polarimetric and Interferometric Radiometers

Two types of microwave radiometry--synthetic thinned array radiometry (STAR) and fully-polarimetric (FP) radiometry--have received increasing attention during the last several years. STAR radiometers offer a technological solution to achieving high spatial resolution imaging from orbit without requiring a filled aperture or a moving antenna, and FP radiometers measure extra polarization state information upon which entirely new or more robust geophysical retrieval algorithms can be based. Radiometer configurations used for both STAR and FP instruments share one fundamental feature that distinguishes them from more 'standard' radiometers, namely, they measure correlations between pairs of microwave signals. The calibration requirements for correlation radiometers are broader than those for standard radiometers. Quantities of interest include total powers, complex correlation coefficients, various offsets, and possible nonlinearities. A candidate for an ideal calibration source would be one that injects test signals with precisely controllable correlation coefficients and absolute powers simultaneously into a pair of receivers, permitting all of these calibration quantities to be measured. The complex nature of correlation radiometer calibration, coupled with certain inherent similarities between STAR and FP instruments, suggests significant leverage in addressing both problems together. Recognizing this, a project was recently begun at NASA Goddard Space Flight Center to develop a compact low-power subsystem for spaceflight STAR or FP receiver calibration. We present a common theoretical framework for the design of signals for a controlled correlation calibration source. A statistical model is described, along with temporal and spectral constraints on such signals. Finally, a method for realizing these signals is demonstrated using a Matlab-based implementation.

Kim, Edward J.↗
NASA NTRS2001

Thresholds in Atmosphere-Soil Moisture Interactions: Results from Climate Model Studies

The potential predictability of the effects of warm season soil moisture anomalies over the central U.S. has been investigated using a series of GCM (Global Climate Model) experiments with the NCAR (National Center for Atmospheric Research) CCM3 (Community Climate Model version 3)/LSM (Land Surface Model). Three different types of experiments have been made, all starting in either March (representing precursor conditions) or June (conditions at the onset of the warm season): (1) 'anomaly' runs with large, exaggerated initial soil moisture reductions, aimed at evaluating the physical mechanisms by which soil moisture can affect the atmosphere; (2) 'predictability' runs aimed at evaluating whether typical soil moisture initial anomalies (indicative of year-to-year variability) can have a significant effect, and if so, for how long; (3) 'threshold' runs aimed at evaluating if a soil moisture anomaly must be of a specific size (i.e., a threshold crossed) before a significant impact on the atmosphere is seen. The 'anomaly' runs show a large, long-lasting response in soil moisture and also quantities such as surface temperature, sea level pressure, and precipitation; effects persist for at least a year. The 'predictability' runs, on the other hand, show very little impact of the initial soil moisture anomalies on the subsequent evolution of soil moisture and other atmospheric parameters; internal variability is most important, with the initial state of the atmosphere (representing remote effects such as SST anomalies) playing a more minor role. The 'threshold' runs, devised to help resolve the dichotomy in 'anomaly' and 'predictability' results, suggest that, at least in CCM3/LSM, the vertical profile of soil moisture is the most important factor, and that deep soil zone anomalies exert a more powerful, long-lasting effect than do anomalies in the near surface soil zone. We therefore suggest that soil moisture feedbacks may be more important in explaining prolonged decadal to century-long droughts evident in the historic and recent prehistoric records, but less important on a seasonal to interannual time-scale.

Oglesby, Robert J.↗
NASA NTRS2006

Micromechanics Modeling of Functionally Graded Interphase Regions in Carbon Nanotube-Polymer Composites

The effective elastic properties of a unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes are modeled herein through the use of a multi-scale method involving the molecular dynamics/equivalent continuum and micromechanics methods. The specific lamina representative volume element studied consists of a carbon fiber surrounded by a region of epoxy containing a radially varying concentration of carbon nanotubes which is then embedded in the pure epoxy matrix. The variable concentration of carbon nanotubes surrounding the carbon fiber results in a functionally graded interphase region as the properties of the interphase region vary according to the carbon nanotube volume fraction. Molecular dynamics and equivalent continuum methods are used to assess the local effective properties of the carbon nanotube/epoxy comprising the interphase region. Micromechanics in the form of the Mori-Tanaka method are then applied to obtain the global effective properties of the graded interphase region wherein the carbon nanotubes are randomly oriented. Finally, the multi-layer composite cylinders micromechanics approach is used to obtain the effective lamina properties from the lamina representative volume element. It was found that even very small quantities of carbon nanotubes (0.36% of lamina by volume) coating the surface of the carbon fibers in the lamina can have a significant effect (8% increase) on the transverse properties of the lamina (E22, k23, G23 and G12) with almost no affect on the lamina properties in the fiber direction (E11 and v12).

Seidel, Gary D.↗
NASA NTRS2007

ISRU at a Lunar Outpost: Implementation and Opportunities for Partnerships and Commercial Development

The NASA Lunar Architecture Team (LAT), which was commissioned to help answer the question "how" will humans return to the Moon, and the Synthesis Team and the recently released Global Exploration Strategy, which was commissioned to help answer the question "why" will humans return to the Moon and go on to Mars have identified the ability to extract and use in-situ resources as important to extending human frontiers, reduce dependence on Earth, and further economic and commercial expansion into space. The extraction and processing of space resources into useful products, known as In-Situ Resource Utilization (ISRU), can have a substantial impact on mission and architecture concepts. In particular, the ability to make propellants, life support consumables, and fuel cell reagents can significantly reduce the cost, mass, and risk of sustained human activities beyond Earth. Potential lunar resources include solar wind implanted volatiles, vast quantities of metal and mineral oxides, possible water/ice at the poles, abundant solar energy, regions of permanent light and darkness, the vacuum of space itself, and even scavenging leftover descent propellants and/or trash and waste from human crew activities. Suitable processing can transform these raw resources into useful materials and products. The establishment of a human lunar Outpost, as proposed by NASA at the 2nd Space Exploration Conference in Houston in December 2006, opens up the possibility for the first time of breaking our reliance on Earth supplied consumables and learn to "live off the land". The ISRU phasing and capability incorporation strategy developed during LAT Phase I & II is based on the premise that while ISRU is a critical capability and key to successful implementation of the US Vision for Space Exploration, it is also an unproven capability for human lunar exploration and can not be put in the critical path of architecture success until it has been proven. Therefore, ISRU needs to take incremental steps toward the desired end state. However, at the same time, the lunar architecture needs to be open enough to take advantage of ISRU when proven available.

Sanders, Gerald B.↗
NASA NTRSDate not supplied

RTG Radiator Efficiency in the Presence of Lunar Dust

The accumulation of Lunar dust on various surfaces is the subject of much recent study, yet is still poorly understood. Recent plans to return humans to the Moon as part of the Artemis program mean that the presence of robust infrastructure to support a permanent human presence is inevitable. Power supply systems, such as Radioisotope Thermoelectric Generators are a likely part of that infrastructure, but require radiative thermal dissipation through a series of radiator fins to function effectively. In this work, the results of a parametric study of the impact on radiator efficiency due to various quantities and types of lunar dust simulants are presented, determining that the impact from low levels of lunar dust accumulation is not significantly detrimental to the efficiency of Multi-Mission Radioisotope Thermoelectric Generator performance in a vacuum, but increased fidelity of data would improve these findings.

Dust↗
DOE OSTI2024

State of Innovation 2024: Paving the Way for Low-Carbon Cement and Concrete

Concrete production in the U.S. accounted for nearly 393 million cubic yards in 2023, with $38.8 billion in revenue. While cement, the key ingredient in concrete, is only 10%-15% of concrete's mixture by mass, 8% of total global emissions come from the production of cement. A significant challenge lies in the fact that roughly 51% of concrete emissions stem from the material calcination process of cement production. Stakeholders ranging from cement producers to government agencies are beginning to take measures to significantly reduce concrete emissions by 2050 and are seeking novel technologies from the startup community. Strategies to lower carbon emissions include reducing quantities of cement in concrete formulas; optimizing digital and automated production; lower temperature processing; carbon capture, utilization, and durable storage; and other cost-saving approaches to energy and material efficiencies. New materials based on carbon mineralization and novel cement chemistries hold the potential to reach net-zero or even carbon-negative concrete production; however, there is presently no readily available substitute that can replicate concrete's unique properties and versatility at the volume demanded by construction worldwide. The nature of concrete's raw materials, diverse applications, and scale of demand means that complete decarbonization will rely on a combination of innovative production methods and novel cement chemistries. Low-carbon solutions will need to compete economically with traditional concrete to become viable in a high-volume commodity market, although consumer demand and regulation will play important roles. Despite these challenges, this analysis reveals that venture capital (VC) investments in low-carbon concrete reflect a growing awareness of the decarbonized cement market opportunity. Between 2022 and 2023 emerging low-carbon technologies garnered more than $700 million in VC investments, representing a growing share of investment in the built environment. An investment gap appears after Series A for technology solutions, demonstrating the sector's potential, as well as the need for additional performance assurance and technology incubation.

cement↗
DOE OSTI2023

Control optimization for parametric Hamiltonians by pulse reconstruction

Optimal control techniques provide a means to tailor the control pulses required to generate customized quantum gates, which helps to improve the resilience of quantum simulations to gate errors and device noise. However, the significant amount of (classical) computation required to generate customized gates can quickly undermine the effectiveness of this approach, especially when pulse optimization needs to be iterated. We propose a method to reduce the computational time required to generate the control pulse for a Hamiltonian that is parametrically dependent on a time-varying quantity. We use simple interpolation schemes to accurately reconstruct the control pulses from a set of pulses obtained in advance for a discrete set of predetermined parameter values. We obtain a reconstruction with very high fidelity and a significant reduction in computational effort. We report the results of the application of the proposed method to device-level quantum simulations of the unitary (real) time evolution of two interacting neutrons based on superconducting qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2020

Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics

Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. Here, we propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo (MLMF) estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models (zero- and one-dimensional) to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these two low-fidelity model alternatives and our high-fidelity three-dimensional models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. Our goal is to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to single fidelity Monte Carlo estimators and to a multilevel Monte Carlo approach based only on three-dimensional simulations, but leveraging multiple spatial resolutions. We demonstrate significant, on the order of 10 to 100 times, reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities such as outlet pressure and flow show larger reductions than local quantities, such as those relating to wall shear stress, as the latter rely more heavily on the highest fidelity model evaluations. Similarly, healthy models show larger reductions than diseased models. In all cases, our workflow coupling Dakota’s MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2002

Correlation of the Coronal Mass Ejection Productivity of Solar Active Regions with Measures of their Global Nonpotentiality from Vector Magnetograms: Baseline Results

Conventional magnetograms and chromospheric and coronal images show qualitatively that the fastest coronal mass ejections (CMEs) are magnetic explosions from sunspot active regions where the magnetic field is globally strongly sheared and twisted from its minimum-energy potential configuration. We present measurements from active region vector magnetograms that start to quantify the dependence of an active region's CME productivity on the global nonpotentiality of its magnetic field. From each of 17 magnetograms of 12 bipolar active regions, we measured the size of the active region (the magnetic flux content, phi) and three separate measures of the global nonpotentiality (L(sub SS), the length of strong-shear, strong-field main neutral line: I(sub N), the net electric current connecting one polarity to the other; and alpha = (mu)I(sub N)/phi), a flux normalized measure of the field twist). From these measurements and the observed CME productivity of the active regions, we find that: (1) All three measures of global nonpotentiality are statistically correlated with the active region flux content and with each other; (2) All three measures of global nonpotentiality are significantly correlated with CME productivity. The flux content correlates with CME productivity, but at a lower statistically significant confidence level (less than 95%); (3) The net current is less closely correlated with CME productivity than alpha and the correlation of CME productivity with flux content is even weaker. If these differences in correlation strength, and a significant correlation of alpha with flux content, persist to larger active regions, this would imply that the size of active regions does not affect CME productivity except through global nonpotentiality; and (4) For each of the four global magnetic quantities, the correlation with CME productivity is stronger for a two-day time window for the CME production than for windows half as wide or twice as wide. This plausibly is a result of the most counterproductive active regions producing less than one CME per day, and from the active region's evolution often significantly changing the global nonpotentiality over the course of several days. These results establish that measures of active region global nonpotentiality from vector magnetograms (such as L(sub SS), I(sub N), and alpha) should be useful for prediction a active region CMEs.

Falconer, D. A.↗
NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗
NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗