Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal coupling”

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.

329 records
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations↗
NASA NTRS1999

Coupled, Transient Simulations of the Interaction Between Power and Secondary Flowpaths in Gas Turbines

This paper presents a coupled analysis of the interaction between mainpath and secondary flowpaths in gas turbines using transient simulations. Some of the topics include: 1) Need for Coupled Analysis; 2) Primary-Secondary Coupling Schematic; 3) Secondary Flow Requirement; 4) Objectives of Present Methodology; 5) Current Methodologies Recap; 6) Proposed Coupled Code Methodology; 7) Description of SCISEAL Code; 8) Description of Turbo Code; 9) Code Coupling/Interface Issues; and 10) Current Interface Strategy. This paper is presented in viewgraph form.

M M Athavale↗
NASA NTRS1990

Discharge Transient Coupling in Large Space Power Systems

Experiments have shown that plasma environments can induce discharges in solar arrays. These plasmas simulate the environments found in low earth orbits where current plans call for operation of very large power systems. The discharges could be large enough to couple into the power system and possibly disrupt operations. Here, the general concepts of the discharge mechanism and the techniques of coupling are discussed. Data from both ground and flight experiments are reviewed to obtain an expected basis for the interactions. These concepts were applied to the Space Station solar array and distribution system as an example of the large space power system. The effect of discharges was found to be a function of the discharge site. For most sites in the array discharges would not seriously impact performance. One location at the negative end of the array was identified as a position where discharges could couple to charge stored in system capacitors. This latter case could impact performance.

N John Stevens↗
NASA NTRS1970

Prevention of Coupled Structure-Propulsion Instability (POGO) On The Space Shuttle

Many liquid-propellant rocket vehicles have experienced longitudinal vibration, nicknamed pogo after the jumping stick, due to an instability arising from interaction of the vehicle structure with the propulsion system. These vibrations can produce an intolerable environment for astronauts and equipment, can overload vehicle structure, and can lead to loss of propulsion performance. A NASA monograph concerned with the prevention of pogo is well on its way toward publication and should serve as a useful guide to design for the Space Shuttle. The paper will briefly describe the nature of the pogo phenomenon, present a summary of past instabilities, relate the possible deleterious effects of pogo, and discuss the contents of the design criteria monograph. Finally, research tasks are recommended for support of pogo development activities on the Space Shuttle.

S. Rubin↗
NASA NTRSDate not supplied

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↗
NASA NTRS1996

Minimum Heating Re-Entry Trajectories for Advanced Hypersonic Launch Vehicles

Optimal re-entry trajectories are generated for reusable launch vehicles which minimize: (1) the heat absorbed at the vehicle surface, (2) the lower surface temperature, and (3) the heat absorbed by the internal structure. The approach uses the energy state approximation technique and a finite control volume heat transfer code coupled to a flight path integration code. These trajectories are compared to the optimal re-entry trajectory minimizing the integrated convective heat rate to determine which trajectory produces the minimum internal structural temperatures for a given thermal protection system. Three different thermal protection systems are considered: tile, blanket, and metallic.

Robert Windhorst↗
NASA NTRS2012

Analysis of Screen Channel LAD Bubble Point Tests in Liquid Methane at Elevated Temperature

This paper examines the effect of varying the liquid temperature and pressure on the bubble point pressure for screen channel Liquid Acquisition Devices in cryogenic liquid methane using gaseous helium across a wide range of elevated pressures and temperatures. Testing of a 325 x 2300 Dutch Twill screen sample was conducted in the Cryogenic Components Lab 7 facility at the NASA Glenn Research Center in Cleveland, Ohio. Test conditions ranged from 105 to 160K and 0.0965 – 1.78 MPa. Bubble point is shown to be a strong function of the liquid temperature and a weak function of the amount of subcooling at the LAD screen. The model predicts well for saturated liquid but under predicts the subcooled data.

Reaction Control System↗
NASA NTRS2026

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↗
NASA NTRS2026

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↗
NASA NTRS2026

Development and Validation of a High-Vacuum Thermal Conductivity Testbed for Aerospace Interface Materials

Thermal Interface Materials (TIMs) are critical components in spacecraft thermal management systems, where thermal performance is strongly influenced by vacuum conditions, interface contact resistance, and layered metallic joint behavior. However, manufacturer-reported thermal conductivity values are often derived under idealized conditions and may not accurately represent performance within operational aerospace applications. To address this limitation, the Testbed for Advanced Interface Materials in Vacuum (TAIMV) was developed as a modular vacuum-compatible thermal conductivity characterization platform capable of evaluating aerospace-relevant TIM configurations under both ambient and high-vacuum environments. The testbed was derived from the ASTM C1044-16 guarded hot plate methodology and incorporates interchangeable layers of stainless steel coupon geometries, independently controlled main and guard heaters, embedded resistance temperature detectors (RTDs), thermocouples, multi-layer insulation (MLI), and a temperature-controlled cold plate to characterize through-thickness thermal gradients across layered interfaces. In the current configuration, interface compression is limited to the nominal contact pressure generated by the experimental stack assembly. Initial experimental campaigns were conducted at ambient pressure and below 1×10-5 torr for vacuum cases using multiple interface materials including Braycote 601EF and Krytox-based greases across a range of thermal operating conditions. In parallel, a coupled numerical Python thermal model was developed to predict temperature distribution throughout the stack while accounting for conduction, radiation, and parasitic heat transfer pathways and effective interface resistance effects. Experimental measurements and numerical predictions showed consistent thermal trends across multiple operating conditions and environmental states. Results also revealed measurable differences between ambient and vacuum thermal behavior, demonstrating the importance of interface resistance, parasitic heat transfer mechanisms, and stack geometry in determining effective thermal performance within layered thermal interfaces. The presented work establishes a foundation for future thermal model correlation efforts and expanded characterization of aerospace thermal interface materials under representative environmental conditions. Future work will focus on the integration of a load cell system to enable controlled pressure-dependent characterization of thermal interface materials under compressive loading. This capability will allow investigation of the influence of contact pressure on effective thermal conductivity, interface resistance, and thermal performance within layered aerospace thermal interfaces under representative operational conditions.

Thermal Development Testing↗
NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗
NASA NTRS2026

Thermochemical Stability of Ca2Yb8(SiO4)6O2 Apatite in Presence of Molten Calcium-Magnesium-Aluminosilicate (CMAS)

Thermochemical stability of ytterbium silicon oxyapatite Ca 2 Yb 8 (SiO 4 ) 6 O 2 (CYbS) in the presence of molten calcium-magnesium aluminosilicate (CMAS) has been investigated at elevated temperatures for consideration as a thermal and environmental barrier coating (T/EBC) material. CYbS apatite powder was synthesized from the constituent oxides via a solid-state reaction method. Hot-pressed apatite substrates were exposed to molten CMAS at 1200, 1300, and 1400 °C for 1, 10, and 50 h. Development of phases in the interaction region of the heat-treated specimens was monitored using scanning electron microscopy, transmission electron microscopy, high-angle annular dark-field imaging, selected area electron diffraction, and energy dispersive x-ray spectroscopy. Monoclinic cyclosilicate Ca 3 Yb 2 (Si 3 O 9 ) 2 formed from interaction of CYbS apatite with CaO in the CMAS melt at the apatite-CMAS reaction front and continued to nucleate and grow within the residual CMAS in diffusion couples annealed for 1 to 50 h at 1200 °C as well as in those heat treated at 1300 °C for 1 h. Residual CMAS was depleted of Ca when cyclosilicate was present. Dendritic wollastonite (CaSiO3) was observed within the residual CMAS in couples annealed at 1200 and 1300 °C. Ingress of molten CMAS, because of its exponential decrease in viscosity, occurred through open pores and along the grain boundaries of the apatite substrates without any detectable chemical reaction at 1300 and 1400 °C. Results of this study indicate that Ca 2 Yb 8 (SiO 4 ) 6 O 2 apatite has the potential to mitigate the CMAS corrosion up to about 1200 °C but not at higher temperatures.

Narottam P Bansal↗
NASA NTRS2024

Wave Rotor Enhanced Nuclear (WREN) Propulsion: NASA Innovative Advanced Concepts (NIAC) - Phase I Final Report

Nuclear Thermal Propulsion (NTP) is identified as one of the preferred propulsion technologies for manned missions throughout the solar system (NASA MSFC).[1, 2] The state-ofthe-art NTP cycle is based on a solid core Nuclear Engine for Rocket Vehicle Application (NERVA)[3] class technology (Fig. 1) that is envisioned to provide a specific impulse of 900 seconds doubling chemical rocket performance (450 seconds). Even with this impressive increase, the NTP NERVA designs still have issues providing adequate initial to final mass fractions for high ΔV missions.[4] Nuclear Electric Propulsion (NEP) can provide extremely high Isp (2,000 to over 10,000 seconds) but with only low thrust and limits on mass to power ratios. The need for an electric power source also adds the issue of heat rejection in space where thermal energy conversion is at best 30-40% under ideal conditions. NASA Space Technology Mission Directorate (STMD) has recently expressed interest in finding advanced nuclear propulsion technology through the NASA Go:Thrust RFI.[5, 6] A novel Wave Rotor (WR) topping cycle has been proposed for our NASA NIAC concept. It promises to deliver similar thrust as NERVA class NTP propulsion, but with Isp in the 1,200-2,000 second range. Coupled with an NEP cycle, the duty cycle Isp can further be increased (1,800-4,000 seconds) with minimal addition of dry mass. This bimodal design enables fast transit trajectories for manned missions to Mars and revolutionizes the deep space exploration of our solar system.

Nuclear Thermal Propulsion↗
NASA NTRS2026

Deep RL for Fast Long-Horizon Operations Scheduling on NASA's Carruthers Geocorona Observatory Mission

Spacecraft operations scheduling is a highly constrained, long-horizon combinatorial optimization problem that traditionally relies on heuristics, constraint programming, or manual planning. We present a scalable deep reinforcement learning framework developed and deployed for NASA’s Carruthers Geocorona Observatory mission. Our framework introduces a macro-action abstraction known as activity blocks coupled with dynamic action-masking to navigate the intractably large search space and strictly enforce complex power, thermal, and instrument constraints. The resulting architecture generates globally feasible schedules with overwhelming probability, establishes operational trust, and executes a full training cycle in under six hours, circumventing the need for policy robustness by enabling rapid, on-demand retraining. Further, resulting schedules outperform baseline heuristics in scheduled science quality. The deep reinforcement learning framework was deployed as the default operational scheduler for the Carruthers Geocorona Observatory mission from the outset of the mission, demonstrating that deep reinforcement learning can be trusted for real spacecraft operations under complex, evolving constraints.

Geocorona↗
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↗