Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gas Dynamics”

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.

267 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

CFD-Assisted Nodal Modeling of Sloshing in a Cryogenic Propellant Tank

During autogenous pressurization, tank sloshing causes a significant increase in pressurant consumption to maintain constant ullage pressure during draining of the tank. This increase in pressurant consumption is caused by a significant increase in condensation at the liquid-vapor interface. Sloshing strongly affects the liquid side heat transfer coefficient and thereby the condensation rate. Traditionally, sloshing is modeled by CFD code using the VOF (Volume of Fluid) method to track the liquid vapor interface during sloshing. CFD calculations require a very fine grid to accurately compute the heat and mass transfer at the interface. Therefore, computations are time consuming and prohibit performing many parametric studies often needed during the design of a new system. This paper describes an alternative approach by developing a CFD-assisted nodal model to predict system parameters more economically with reasonable accuracy. In this approach, a multi-node model of tank pressurization was developed using GFSSP. A multi-node model was needed to account for stratification. The model computes heat and mass transfer at the interface to calculate the condensation rate. The liquid side heat transfer is computed using the parameters of sloshing dynamics such as frequency, wave amplitude, and interface area. The parameters of sloshing dynamics are computed by the CFD code LOCI-Stream. The model predictions were compared with test data for several cases.

Cryogenic Tank Sloshing↗
NASA NTRS2026

Evidence of an Extended Period of Aqueous Alteration on the Bennu Parent Body Revealed By I-Xe Analyses

NASA’s OSIRIS-REx mission delivered to Earth 121.6 g of regolith material from asteroid Bennu on 24th September 2023. We report multi-step laser pyrolysis xenon isotopic analysis to investigate the I-Xe system in individual particles separated from aggregate material (a mixture of small particles of undefined lithology). The high sensitivity of the RELAX mass spectrometer allows us to observe small excesses of individual isotopes and separate out different components, which may not be possible using a conventional noble gas mass spectrometer. The temperature resolution of RELAX is sufficient to seek an I-Xe isochron from step heating analysis of a sample mass <100 μg, which can contribute to our understanding of the nature and timing of aqueous alteration on Bennu’s parent body. We also aim to better estimate and understand the iodine concentration of CI chondrite-like material.

S A Crowther↗
NASA NTRS1997

Gravitational Effects in Smoldering and SHS

Smolder waves and SHS (self-propagating high-temperature synthesis) waves are both examples of filtration combustion waves propagating in porous media. Smoldering combustion is important for the study of fire safety. Smoldering itself can cause damage, its products are toxic and it can also lead to the more dangerous gas phase combustion which corresponds to faster propagation at higher temperatures. In SHS, a porous solid sample, consisting of a finely ground powder mixture of reactants, is ignited at one end. A high temperature thermal wave, having a frontal structure, then propagates through the sample converting reactants to products. The SHS technology appears to enjoy a number of advantages over the conventional technology, in which the sample is placed in a furnace and 'baked' until it is 'well done'. The advantages include shorter synthesis times, greater economy, in that the internal energy of the reactions is employed rather than the costly external energy of the furnace, purer products, simpler equipment and no intrinsic limitation on the size of the sample to be synthesized, as exists in the conventional technology.

Bernard J Matkowsky↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗
NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗
NASA NTRSDate not supplied

Solvent-Free Preparation of High Energy, Binder-Free Electrodes Enabled by Dry Compressible Holey Graphene

Graphene is an atomically thick sheet consisting of a graphitic carbon network with excellent mechanical strength, electrical and thermal conductivity, and chemical stability. Holey graphene, a structural derivative of graphene, has an array of through-the-thickness holes across the lateral surface of the nanosheet. The presence of these holes has minimal detrimental effect on the graphene properties and leads to enhanced performance in applications such as electronics, sensors, and energy storage. For example, these holes allow more facile cross-plane ion and gas transport than intact graphene, making holey graphene an ideal electrode material for electrochemical energy storage. This presentation will focus on the ability of holey graphene to be compression molded into robust articles or architectures under solvent-free conditions without the need for potentially parasitic binders. The unique dry compressibility of holey graphene has enabled facile fabrication of high mass loading electrodes with both high density and high porosity, which have found use in supercapacitors and various high-energy battery systems such as lithium-oxygen, lithium-sulfur, and lithium-selenium batteries.

Yi Lin↗
NASA NTRS1985

Structural Integrity and Durability of Reusable Space Propulsion Systems

A two-day conference on the structural integrity and durability of reusable space propulsion systems was held on June 4 and 5, 1985, at the NASA Lewis Research Center. Presentations were made by industry, university, and government researchers organized into four sessions: aerothermodynamic loads; structural dynamics, fatigue, fracture, and constitutive modeling; and instrumentation. The principal objectives of the conference were to disseminate research results to date and future plans in each of the four areas. This publication contains the extended abstracts and the visual material presented during the conference.

Stanley J Marsik↗
NASA NTRS2003

Control of Thermo-Acoustics Instabilities: The Multi-Scale Extended Kalman Approach

"Multi-Scale Extended Kalman" (MSEK) is a novel model-based control approach recently found to be effective for suppressing combustion instabilities in gas turbines. A control law formulated in this approach for fuel modulation demonstrated steady suppression of a high-frequency combustion instability (less than 500Hz) in a liquid-fuel combustion test rig under engine-realistic conditions. To make-up for severe transport-delays on control effect, the MSEK controller combines a wavelet -like Multi-Scale analysis and an Extended Kalman Observer to predict the thermo-acoustic states of combustion pressure perturbations. The commanded fuel modulation is composed of a damper action based on the predicted states, and a tones suppression action based on the Multi-Scale estimation of thermal excitations and other transient disturbances. The controller performs automatic adjustments of the gain and phase of these actions to minimize the Time-Scale Averaged Variances of the pressures inside the combustion zone and upstream of the injector. The successful demonstration of Active Combustion Control with this MSEK controller completed an important NASA milestone for the current research in advanced combustion technologies.

Dzu K Le↗
NASA NTRSDate not supplied

Enabling Interoperability in Earth System Digital Twins (ESDT): Integrating Observations, Models, and AI for Actionable Insights Through NASA'S Intelligent Systems Technology Program

NASA’s Intelligent Systems Technology Program (IST) is driving a paradigm shift in Earth science through the development of Earth System Digital Twins (ESDT). These integrated information systems create a dynamic "digital replica" of the Earth by harmonizing continuous, multi-source observations with high-fidelity models and state-of-the-art artificial intelligence (AI) that enable “What now?”, “What next?”, and “What if?” scenario building. These scenarios are reflected in NASA IST’s series of ESDTs, from the Coastal Zone Digital Twin that integrates complex data on the current state of the Chesapeake Bay to the Terrestrial Environmental Rapid-Replication and Assimilation Hydrometeorological (TerraHydro) AI-based ESDT that forecasts water movement across Earth’s surface, to the Agriculture Land Information System (AgLIS) which can be used to assess optimal planting dates and crop yield estimates. By bridging the gap between vast data archives and actionable insights, these projects enable a system-of-systems approach to understanding complex, interacting Earth processes. This poster will highlight recent innovations and future directions from NASA’s ESDT initiatives: Continuous Data Assimilation & Multi-Source Fusion. A core requirement of the ESDT work is the transition from static models to dynamic "living" replicas. This involves creating frameworks for the continual assimilation of near-real-time data from uncoordinated, heterogeneous sources, including satellite observations and airborne assets, and ground-based Internet of Things (IoT) sensors. These systems link design, operational status, and environmental data, ensuring the digital twin accurately reflects the current state of the physical Earth system. High-Fidelity Hybrid Modeling & Computational Acceleration to enable interactive "what-if" explorations, programs are moving beyond traditional, slow physical solvers by developing fast surrogate machine learning models and Deep Generative Models (DGMs). These hybrid approaches use neural networks to emulate complex physics, such as cloud feedback or ocean dynamics, at a fraction of the original computing cost, often leveraging advanced hardware like Graphics Processing Units (GPUs) to achieve the necessary scale. Federated Ecosystems & Interoperable Frameworks rather than building isolated tools, NASA IST is moving toward federated ESDTs and reusable analytic collaborative frameworks. This theme focuses on interoperability standards and common ontologies that allow specialized digital twins to interact and share data. This system-of-systems architecture supports multi-discipline investigations, such as analyzing how upstream watershed changes impact downstream urban flooding or how wildfire emissions affect regional air quality. By leveraging these advancements, ESDTs empower researchers and decision-makers to conduct real-time analysis and run complex hypothetical scenarios, ultimately improving our understanding of Earth’s evolving systems and informing critical real-world applications.

Earth System↗
NASA NTRS1977

Research on battery-operated electric road vehicles

Mathematical analysis of battery-operated electric vehicles is presented. Attention is focused on assessing the influence of the battery on the mechanical and dynamical characteristics of dc electric motors with series and parallel excitation, as well as on evaluating the influence of the excitation mode and speed control system on the performance of the battery. The superiority of series excitation over parallel excitation with respect to vehicle performance is demonstrated. It is also shown that pulsed control of the electric motor, as compared to potentiometric control, provides a more effective use of the battery and decreases the cost of recharging.

Varpetian, V. S.↗
NASA NTRS1981

The 50 AMP-hour nickel cadmium battery manual

The battery is designed with a minimum battery to cell weight ratio consistent with adequate containment for operating conditions and dynamic environments and minimized weight. The battery is fully qualified and the environments to which it was successfully subjected were selected by NASA Goddard to cover a wide range of probable uses. The battery is suitable for either near-Earth geosynchronous missions, is compatible with passive or active thermal control systems and may be electrically controlled by a variety of changing routines. The initial application of the 50 A.H. Battery is a near-Earth mission aboard the LANDSAT D Satellite.

Webb, D. A.↗
NASA NTRS1990

Interactions Between the Space Station and the Environment: A Preliminary Assessment of EMI

A review of the interactions between proposed Space Station systems/payloads and the environment that contribute to electromagnetic interference was performed. Seven prime sources of interference have been identified. These are: The Space Station power system; active experiments such as beam injection; ASTROMAG; ram and wake density gradients; pick up ions produced by vented or offgassed clouds; waves produced by current loops that include the plasma and structure; arcing from high voltage solar arrays (or possible ESD in polar orbit). This review indicates that: minimizing leakage current from the 20 kHz power system to the structure; keeping the surfaces of the Space Station structure, arrays, and radiators nonconducting; minimizing venting of payloads or systems to non-operational periods; careful placement of payloads sensitive to magnetic field perturbations or wake noise; and designing an operational timeline compatible with experiment requirement are the most effective means of minimizing the effects of this interference. High degrees of uncertainty exist in the estimates of magnitudes of gas emission induced EMI, radiation of 20 kHz and harmonics, ASTROMAG induced interference, and arc threshold/frequency of the solar array. These processes demand further attention so that mitigation efforts are properly calibrated.

G B Murphy↗
NASA NTRSDate not supplied

Exploration Extravehicular Mobility Unit 11 Foot Vacuum Chamber Upgrades Test Results

The Exploration Extravehicular Mobility Unit (xEMU) uncrewed 11 foot vacuum chamber testing evaluated the capabilities of the 11 foot vacuum chamber facility to support advanced spacesuit testing. The government reference design xEMU spacesuit provided a high-fidelity test article to demonstrate 11 foot vacuum chamber capabilities which included: gas loading of the chamber at varying simulated metabolic rates and open loop suit abort operations, Intravehicular Activity (IVA) vacuum access, consumables recharge, IVA vehicle-provided thermal loop cooling, and IVA vehicle-provided power. To demonstrate the xEMU airlock operations transitioning from IVA to EVA conditions without a test subject in the suit, test support equipment was developed to remotely actuate both the Exploration, Servicing, and Cooling Umbilical (ESCU) and the vacuum access umbilical. This test also evaluated the performance of the Exploration Portable Life Support System (xPLSS) at vacuum conditions. Data was collected and analyzed for carbon dioxide (CO2) scrubbing performance of the Rapid Cycle Amine (RCA) swingbed, for thermal regulation performance of the Suit Water Membrane Evaporator (SWME), and for sensor performance across the xPLSS. This paper will detail the findings of the testing performed with these upgrades which discussed previously laid out in ICES-2025-342.

Robert Marsch↗
NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗
NASA NTRS1984

Nickel-iron battery system safety

The generated flow rates of gaseous hydrogen and gaseous oxygen from an electrical vehicle nickel-iron battery system were determined and used to evaluate the flame quenching capabilities of several candidate devices to prevent flame propagation within batteries having central watering/venting systems. The battery generated hydrogen and oxygen gases were measured for a complete charge and discharge cycle. The data correlates well with accepted theory during strong overcharge conditions indicating that the measurements are valid for other portions of the cycle. Tests confirm that the gas mixture in the cells is always flammable regardless of the battery status. The literature indicated that a conventional flame arrestor would not be effective over the broad spectrum of gassing conditions presented by a nickel-iron battery. Four different types of protective devices were evaluated. A foam-metal arrestor design was successful in quenching gaseous hydrogen and gaseous oxygen flames, however; the application of this flame arrestor to individual cell or module protection in a battery is problematic. A possible rearrangement of the watering/venting system to accept the partial protection of simple one-way valves is presented which, in combination with the successful foam-metal arrestor as main vent protection, could result in a significant improvement in battery protection.

Saltat, R. C.↗
NASA NTRS1995

Validation of the SINDA/FLUINT Code Using Several Analytical Solutions

The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

John R Keller↗
NASA NTRSDate not supplied

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗