Search NASASearch

SEARCH · Search NASA

Results for “Meshed Terrain Thermal Model”

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.

730 records · Page 8

Controls for Electrified Aircraft Propulsion

Advanced electrified aircraft propulsion (EAP) concepts with integrated power, propulsion, and thermal systems require the development of equally sophisticated controllers to fly safe and efficient missions. Hybrid electric aircraft utilize electric machines mechanically coupled to the engine shafts to extract and insert power for a variety of purposes. Additionally, electric machines may drive propulsive fans pulling from a combination of on-board energy storage devices and engine extracted power. NASA’s investment in hybrid electric aircraft hardware-in-the-loop testing enables controls research on representative engine models using a novel emulation and scaling methodology. High-voltage, high-power electrical powertrain presents several technical risks at altitude. High-power requires high efficiency to minimize losses. Superconducting electric machines and power distribution research at NASA has identified key challenges with potential controls solutions. These risks necessitate the development of controllers robust to model uncertainty, disturbances, and fault conditions. System health management and fault detection schemes play an important role in a multi-layer controls approach that utilizes a supervisor which oversees inner loops for the highly coupled power, propulsion, and thermal systems. Such a supervisory controller enables coordinates energy transfer between engines, energy storage devices, electric machines, propulsors, and heat exchangers. Existing methods have shown an improvement in engine operability using the hybrid electric powertrain.

Halle E Buescher

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics

Prediction of prompt NO(x) in hydrocarbon air flames

The gas turbine industry is directing particular attention to very low NOx combustors, whether for aircraft or land based CCGT systems. These low NOx combustors frequently use liquid fuels or natural gas burning under very lean premixed conditions with air or under rich-lean conditions, although only the first case is studied here. In land based systems, diluted steam or nitrogen are sometimes injected into the combustion chambers to reduce flame temperature. The NOx emissions from such systems are the product of three chemical mechanisms which are interrelated: the hydrocarbon prompt NO, the thermal NO (extended Zeldovich mechanism), and the nitrous oxide route to NO. Formation of NO2 from NO also occurs, as well as emission of carbon monoxide and unburnt hydrocarbons. When the fuel-oxidant proportion decreases towards leaner conditions, flame temperatures are lowered, resulting in the total NO being reduced and the thermal-NO contribution greatly diminished to the benefit of the remaining two mechanisms of NO formation. While knowledge of the elementary reactions and their chemical kinetics concerning methane and simple hydrocarbons combustion has existed for a number of years, its use for computer modeling is limited to simple flow dynamics configurations. Nevertheless, understanding of such combusting flows under a wide range of experimental conditions allows for analogies or speculations with more complex actual systems. Such understanding can be achieved by means of one dimensional laminar premixed flame modeling, with a full chemical mechanism which incorporates the three routes of NO formation. Complementary to this understanding is the modeling of the actual combustion system using a full description of the fluid dynamics coupled with a reduced chemical scheme, which is then compared against the first model. The objective of this investigation is to evaluate the relative importance of the three mechanisms of NO formation in lean premixed methane-air combustion with increasing pressure using the one dimensional plug flow package PREMIX, and to test the validity of a three dimensional model with a global chemical mechanism against the one dimensional model in the atmospheric pressure case. Methane is chosen because it is the only mechanism which is reasonably well known and is a good guide to the behavior of other hydrocarbons. The mixture ratio chosen is richer than that in lean gas turbines, but the combustion of this mixture with the low preheat gives realistic gas turbine final flame temperatures. Conditions of NO2 formation are also analyzed in the one dimensional model and results extrapolated to the case of gas turbines.

Valerie Dupont

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity

Use of Assistive Technology to Augment API Capabilities

Application Programming Interfaces (APIs) allow for access to data and capabilities of computer applications by developers or users with experience in computer programming. Recent development with both Thermal Desktop and ESATAN-TMS have provided APIs to allow users to develop their own capabilities that interface with the Graphical User Interfaces (GUI) or manipulate the thermal model data. However, these APIs are only as good as the breadth of features in the native code accessible through the API; if a particular code’s feature is not accessible through the API, then users have very limited options besides waiting for updates to the API that expose the necessary functionality, particularly if model data access or user action, such as a button click, is required. However, Assistive Technology features that allow for differently-abled users to more fully experience a software’s capabilities may be creatively utilized to gain further access to data and capabilities not yet exposed by the API. This paper describes the process to augment the features of the OpenTD API via assistive technology and describes how to identify the application instance, navigate GUI elements, updates values on forms, and execute actions such as selecting a listbox item or clicking a button. It concludes with identifying some of the pitfalls to avoid and describes methods to best implement this approach.

Application Programming Interface

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya

Use of Assistive Technology to Augment API Capabilities

Application Programming Interfaces (APIs) allow for access to data and capabilities of computer applications by developers or users with experience in computer programming. Recent development with both Thermal Desktop and ESATAN-TMS have provided APIs to allow users to develop their own capabilities that interface with the Graphical User Interfaces (GUI) or manipulate the thermal model data. However, these APIs are only as good as the breadth of features in the native code accessible through the API. If a particular code’s feature is not accessible through the API, then users have very limited options besides waiting for updates to the API that expose the necessary functionality, particularly if model data access or user action, such as a button click, is required. However, Assistive Technology features that allow for users with a disability to more fully experience a software’s capabilities may be creatively utilized to gain further access to data and capabilities not yet exposed by the API. This paper describes the process to augment the features of the OpenTD API via assistive technology and describes how to identify the application instance, navigate GUI elements, updates values on forms, and execute actions such as selecting a listbox item or clicking a button. It concludes with identifying some of the pitfalls to avoid and describes methods to best implement this approach.

Application Programming Interface

Preliminary Findings of the Experimental Development Unit Cold Flow Test for a Generation Zero Nuclear Propulsion Engine

Nuclear thermal propulsion (NTP) technology will greatly benefit human travel to Mars by significantly shortening transit times, improving crew safety, and providing more mission flexibility than traditional chemical rockets. As part of DRACO follow-on work to develop, build, and fly a generation zero NTP engine, a full scale flight-like experimental design unit (EDU) reactor was constructed to collect sufficient on-ground performance data to characterize flow induced vibrations (FIV) of critical reactor structures/components, inform development of the engine and reactor control algorithm, and collect pressure drop and flow distribution data across the reactor. The fluid conditions for the test program were designed to achieve system responses equivalent to that of an operational engine through all phases of engine operation including reactor startup, mainstage operation, reactor shutdown, and reactor cooldown. Over 100 tests were executed, flowing either GN2 or GHe through the EDU at varying flow rates and pressures. This experiment provided early validation of flow behavior and vibration risks before nuclear testing, boosted critical subsystem TRLs, informed design iterations, and reduced future test costs. The steady-state flow parameters for the experiment were modeled in Ansys Thermal Desktop, allowing rapid tuning and experiment-informed updates to a flight-like test matrix. The EDU dynamic environment was characterized with accelerometers, strain gauges, and high-frequency pressure transducers all sampled at 20 kHz. While many narrow-band oscillations were identified, no significant FIV occurred; the reactor structural responses tend to be enveloped by typical launch vehicle ascent vibration environments (defined up to 2 kHz), although significant energy is also present at higher frequencies.

Flow Induced Vibration

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

Under Pressure: The Artemis I Heatshield Char Loss Investigation and Artemis II’s Successful Entry

Following the successful Artemis I skip reentry on December 11, 2022, unexpected char liberation from the Orion heat- shield prompted formation of an Anomaly Response Team to determine root cause and establish flight rationale and corrective actions for subsequent missions. Through a comprehensive investigation including detailed hardware analysis, modeling & simulation, and ground testing, the team determined that the low permeability Avcoat experienced extreme internal gas pressure buildup during the skip entry that could not adequately outgas, leading to crack formation and char liberation. Key discoveries included significant performance differences between permeable and impermeable heatshield regions and successful replication of the anomaly through ground testing. Based on extensive ground testing and analysis, Artemis II flew a modified trajectory without the skip entry to minimize crack-inducing conditions, and successfully splashed down on April 10, 2026 with significantly reduced char loss. For Artemis III and beyond, a corrective action was implemented to produce a more permeable version of Avcoat that will substantially reduce internal pressure accumulation during entry.

Artemis II

Numerical Investigation of Heat Transfer and Fluid Flow within Electrochemical Hydrogen Peroxide Generation Unit

Long-term manned space missions require the onboard production of disinfectants essential for maintaining crew health and supporting life systems. Currently, disinfection aboard the International Space Station (ISS) relies on disposable wetted wipes, which are regularly resupplied from Earth. This approach imposes a significant burden on resupply logistics, storage, and waste management. To address these challenges and support future missions, efforts are underway to develop an in-situ solution that electrochemically generates hydrogen peroxide disinfectant using onboard resources. In collaboration with NASA, Faraday Technology, Inc. has advanced this concept through a series of Small Business Innovation Research (SBIR) projects, resulting in the development of a Peroxide Generation Unit (PGU). The PGU can produce up to 3 wt.% hydrogen peroxide on-demand at a rate of 1 liter per day, providing a sustainable alternative to Earth-dependent supplies. The resulting aqueous hydrogen peroxide (H₂O₂) is an effective disinfectant, safe for crew use, compatible with spacecraft systems, and free from volatiles, off-gassing, or residues. This innovation offers a reliable, efficient solution for onboard disinfection, reducing dependence on Earth-based resupply while ensuring the health and safety of space crews. Generating hydrogen peroxide at the required rate needs high voltages and currents, exceeding 20V and 2A respectively, which leads to significant heat generation from Joule heating. This temperature rise poses a risk to sensitive system components, especially critical and expensive membranes that can degrade under thermal stress. To mitigate this risk, the thermal, fluid, and electrical flows within the system are modeled computationally using the commercial software COMSOL. The numerical simulations are validated against experimental data from both sub-scale and alpha-scale systems. Once verified, the model is employed to identify thermal hotspots, investigate their underlying causes, and explore solutions to prevent them.

Life Support Systems (LSS)

Selecting an Encapsulant for an Aerospace Superconducting Machine

Achieving larger benefits from electrified propulsion in single aisle aircraft (the largest commercial aircraft market segment) necessitates that a significant fraction of their 20 MW or greater propulsion system be electrified. These aircraft would require up to 10-20 of today’s state-of-the art machines, which have been demonstrated up to 2.5 MW for aerospace use, and could benefit further from 5 MW+ designs. However, when conventionally cooled, high power density machines still produce a significant amount of waste heat, translating into heavy thermal management and detracting from system benefits. Highly efficient fully superconducting electric machines (or those with a cryogenically cooled, high purity copper or aluminum stator) are a possible solution but have multiple barriers to overcome. One such barrier is the durability of superconducting windings. Superconductors are typically brittle, and encapsulants used to hold windings in place are a significant source of stress due to mismatches in coefficient of thermally induced contraction between ambient and operational temperatures (20 to 77 K). Lack of compliance between the materials at operational temperatures also adds significant stress. These encapsulants also form a thermal barrier that can be detrimental to cooling the superconductors. This work describes an initial study of commercially available encapsulants as it relates to superconducting machine applications, along with initial modeling of thermal stress in superconducting windings for a superconducting or cryogenically cooled stator in a 5 MW electric machine.

encapsulant

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Battery Evaluation Profiles for X-57 and Future Urban Electric Aircraft

Battery energy density is one of the most critical design parameters for sizing all-electric aircraft, however it’s easily overestimated. Establishing the effective, usable energy density is confused by varying degrees of margin needed to account for structural and thermal management between different cell chemistry and pack designs. Therefore, a better methodology is needed to fairly compare emerging battery technologies for electric aircraft. Currently, there is a loss of critical information when vehicle trade studies are performed using “nominal” published cell-level performance metrics. Aircraft power demands rarely match these nominal power profiles, and aircraft designers lack the ability to accurately simulate the battery performance and temperature off-nominally unless the battery chemistry is well established. Conversely, battery suppliers have no generalized reference cases to publish more realistic performance metrics. This can lead to poor assumptions, such as aircraft studies assuming a fixed discharge efficiency of a battery, when in reality the usable energy in a pack is dependent on the power and thermal profile. Information needed to properly assess weight penalties for thermal management is also typically poorly characterized when assessing candidate batteries. This paper serves to better inform battery development, and similarly, provide aircraft designers with more realistic assumptions for applying knockdown margins in their designs. Detailed power and thermal performance estimates are provided, which provide a starting point for sizing power and thermal budgets using experimentally derived battery models. Results show that the X-57 battery-to-shaft efficiency is 77.3% for a particular optimized mission. Considering a 25% reserve on the battery capacity, this means that only roughly half of the original 55.3kWh ‘nominal’ pack energy can be converted to useful work during a mission. Further estimates on a clean-sheet VTOL optimization show an average 82.7% battery-to-shaft efficiency, using 98% peak efficiency inverters and 97.4% peak efficiency motors. Although higher battery efficiencies are possible, the resulting weight penalty negates improvement in vehicle performance. These trade-offs and resulting power profiles are provided as a starting point to better assess future battery designs.

Battery