Search NASASearch

SEARCH · Search NASA

Results for “Distributed Autonomy”

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.

131 records · Page 6

Design and Optimization of an Open Rotor Stage

Open rotor propulsors enable greatly increased effective bypass ratios over current aircraft engines, resulting in significantly improved efficiencies. One of the most promising of these geometries is the “open fan”, which consists of a single rotor followed by a single stator. This configuration will be referred to as an “open rotor stage” for design work at NASA. In this work we consider the multidisciplinary design and optimization of an open rotor stage using a variety of 2D and 3D techniques. The 2D aerodynamic techniques used for design include a 2D Trefftz plane method, axisymmetric throughflow, and quasi-3D blade-to-blade CFD. Aerodynamic performance is prioritized, while structural requirements are met and some acoustic considerations are included. Between the rotor, stator, and hub profiles, a total of 101 design variables are used. Differences in the process of designing an open rotor stage vs. designing other high-speed propellers are also discussed. The current design, the NASA Optimized Propulsor for Efficiency and Noise (N-OPEN-1), achieves 84.7% net efficiency as predicted by 3D CFD, with further improvements expected to be possible. This design makes effective use of a slower rotor and a less tip-loaded work distribution than typical single-rotation propellers, allowing the stator to produce 14.8% of the thrust. Laminar flow on the blades was responsible for a 1.7% efficiency benefit and warrants further study.

Open rotors

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

Design and Optimization of an Open Rotor Stage

Open rotor propulsors enable greatly increased effective bypass ratios over current aircraft engines, resulting in significantly improved efficiencies. One of the most promising of these geometries is the “open fan”, which consists of a single rotor followed by a single stator. This configuration will be referred to as an “open rotor stage” for design work at NASA. In this work we consider the multidisciplinary design and optimization of an open rotor stage using a variety of 2D and 3D techniques. The 2D aerodynamic techniques used for design include a 2D Trefftz plane method, axisymmetric throughflow, and quasi-3D blade-to-blade CFD. Aerodynamic performance is prioritized, while structural requirements are met and some acoustic considerations are included. Between the rotor, stator, and hub profiles, a total of 101 design variables are used. Differences in the process of designing an open rotor stage vs. designing other high-speed propellers are also discussed. The current design, the NASA Optimized Propulsor for Efficiency and Noise (N-OPEN-1), achieves 84.7% net efficiency as predicted by 3D CFD, with further improvements expected to be possible. This design makes effective use of a slower rotor and a less tip-loaded work distribution than typical single-rotation propellers, allowing the stator to produce 14.8% of the thrust. Laminar flow on the blades was responsible for a 1.7% efficiency benefit and warrants further study.

turbomachinery

Physics-Based Modeling and Simulation of Emerging Battery Technologies for Aerospace

Recently there is a growing interest in the aviation sector to reduce air and noise pollution. Electrochemical energy storage devices such as batteries coupled with a distributed electric propulsion system can reduce noise concerns as well as emissions and allow the concept of Urban Air Mobility to come to fruition. The battery performance needs to improve considerably compared from current state-of-art Li-ion battery (about 200Wh/Kg) to realize all-electric passenger jet for short flights (up to 690 miles). NASA is exploring lithium-oxygen battery chemistry to power hybrid (battery powered electrical system) and all-electric aircraft for short distance and long-distance flights. Li-O2 is one of the advanced Li-ion technologies that promise to provide specific energy of more than 750Wh/Kg. For this presentation, we present our work on improving power density of Li-O2 batteries through the use of multiphysics simulations. Next, a path is outlined to port these model to simulate performance for a new battery chemistry for space application, Li-CO2. Li-CO2 uses carbon dioxide as the active material instead of oxygen. Although this technology is in its early development, the offers two benefits: it can be used as a CO2 scrubber, as oxygen is one of the by-products on charging, and as a backup or a standalone battery for various Mars or Venus missions, where the carbon dioxide content in the atmosphere is high and need battery to operate at higher temperatures.

Mehta, Mohit

Radio Noise Problems in Arctic Regions

Three main types of radio noise should be considered when establishing noise levels for communication purposes or in experiments utilizing radio techniques in Arctic regions. These are atmospheric and man-made noise, and precipitation static. National Bureau of Standards radio noise data obtained hourly at Arctic and Antarctic stations on eight .fixed frequencies between 51 kc and 20 Mc/s are analyzed to show the characteristics of the three types. The diurnal and seasonal variations of atmospheric radio noise at high latitudes are explainable in terms of changes in propagation factors and in the distribution of world thunderstorm activity. In both northern and southern Arctic regions atmospheric noise decreases during polar cap absorption (PCA), most markedly in the h.f. band. The probable magnitude of man-made noise in the h.f. band is estimated from data taken during PCA, when atmospheric noise is absent. At Thule man-made noise on 2.5 and 5 Mc/s appears to be 57 and 49 dB above kTB (where kTB is antenna thermal noise power), respectively, while at Byrd Station the values are about 20 and 12 dB respectively. Precipitation static is generated on exposed antennas during periods of blowing snow (blizzards). The noise power magnitude can be at least 50 dB above man-made and atmospheric levels, but the maximum enhancement cannot be ascertained because of missing data. In winter and spring months at Byrd Station radio data dependent upon exposed antennas may be lost 10-30 per cent of the time due to precipitation static alone.

Radio Noise

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Julia A. Shates

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)

Influence of Processing Parameters on the Mechanical Properties of 3D Printed Borosilicate Particulate Reinforced Polymer Composites

Emerging composite materials are expanding the potential of additive manufacturing and enabling applications previously restricted by traditional manufacturing methods. The multi-phase nature of these composite materials combined with the complex in-ternal geometry of additively manufactured parts have enabled unique behavior, and potentially new applications. Additionally, these materials can be pyrolyzed to create dense metal, ceramic, and glass parts with geometries typically not achievable by tra-ditional processes. Additive manufacturing of borosilicate glass-based systems can open new applications in nuclear engineering, astronomy, and bone regrowth therapy. To elucidate the process-parameter relationship of borosilicate-polylactic acid (PLA) composites, mechanical testing was conducted and compared with a pure polylactic acid polymer baseline. Test specimens were fabricated by fused-filament fabrication with minimal post-processing. Yield strength, ultimate strength, and elastic modulus were calculated from stress-strain curves. Optical and scanning electron microscopy were conducted to observe the specimen microstructure before and after testing. The highest compressive yield strength for the composite was 28.22 MPa, and the highest compressive yield strength for PLA was 49.30 MPa. Print orientation was found to benefit the composite material but have a detrimental effect on the pure matrix material. An elastic modulus of 2.66 GPa was recorded for the borosilicate-PLA composite at 100% infill, 1 shell wall, and layer lines parallel to compression axis. Microscopy revealed that lower modulus composite specimens had the particulates re-distributed within the matrix. Tensile testing was done according to a polymer testing standard, which caused difficulties obtaining consistent fracture within the gauge length.

mechanical testing

Evidence for Multiple Diagenetic Episodes in Ancient Fluvial‐Lacustrine Sedimentary Rocks in Gale Crater, Mars

The Curiosity rover's exploration of rocks and soils in Gale crater has provided diverse geochemical and mineralogical data sets, underscoring the complex geological history of the region. We report the crystalline, clay mineral, and amorphous phase distributions of four Gale crater rocks from an 80‐m stratigraphic interval. The mineralogy of the four samples is strongly influenced by aqueous alteration processes, including variations in water chemistries, redox, pH, and temperature. Localized hydrothermal events are evidenced by gray hematite and maturation of amorphous SiO 2 to opal‐CT. Low‐temperature diagenetic events are associated with fluctuating lake levels, evaporative events, and ground water infiltration. Among all mudstones analyzed in Gale crater, the diversity in diagenetic processes is primarily captured by the mineralogy and X‐ray amorphous chemistry of the drilled rocks. Variations indicate a transition from magnetite to hematite and an increase in matrix‐associated sulfates suggesting intensifying influence from oxic, diagenetic fluids up section. Furthermore, diagenetic fluid pathways are shown to be strongly affected by unconformities and sedimentary transitions, as evidenced by the intensity of alteration inferred from the mineralogy of sediments sampled adjacent to stratigraphic contacts.

C N Achilles

Incremental Learning for Passive Microwave Precipitation Retrievals using Advanced Technology Microwave Sounder

Spaceborne passive microwave (PMW) radiometry is central to global precipitation monitoring, yet retrieval uncertainties remain substantial, particularly for cross-track sounders whose variable footprints and channel configurations are optimized for atmospheric temperature and moisture profiling rather than precipitation. Consequently, existing operational products often exhibit angular-dependent biases, limited effective swath utilization, unrealistic rainfall probability distributions, and systematic misclassification of precipitation phase. These limitations are further compounded by the scarcity of globally accurate and representative precipitation observations, as training data from the Dual-frequency Precipitation Radar (DPR) and the Cloud Profiling Radar (CPR) are spatially sparse, lack uniform global coverage, and exhibit heterogeneous error characteristics across precipitation regimes. To address these challenges, this study presents a supervised retrieval algorithm that incrementally trains an ensemble of extreme gradient-boosted decision trees by augmenting base learners with pre-training on reanalysis data and post-training on coincident DPR and CPR observations matched with the Advanced Technology Microwave Sounder (ATMS). By transferring prior information from reanalysis to posterior constraints from radar observations and adopting a sequential detection–estimation strategy for precipitation phase and rate retrieval, the proposed approach yields retrievals across the full ATMS swath that are largely free from persistent deficiencies in current Global Precipitation Measurement (GPM) passive microwave operational products. In particular, the method resolves bimodal artifacts in rainfall retrievals and mitigates systematic high-latitude snowfall biases, including overestimation across the Arctic and underestimation across the Antarctic. Validation against independent Multi-Radar Multi-Sensor (MRMS) data over the Contiguous United States (CONUS) further demonstrates improved performance in precipitation phase detection and rate estimation relative to both reanalysis and current GPM PMW products.

Mahyar Garshasbi

Towards an Improved Understanding of the Antarctic Coastal Zone and Its Contribution to Future Global Sea Level

Understanding the coastal zone of the Antarctic Ice Sheet (AIS), where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate and sea level. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the AIS, the Southern Ocean and their global connections in the coming centuries. The AIS remains the largest source of uncertainty in future sea-level projections. Bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, and is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography data sets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation data set of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and improving predictions of future ocean and climate changes.

Kenichi Matsuoka

High Energy Hydroforming Friction Stir Welded 2050 Al-Li Blanks for Aerospace Domes and Cones

Friction stir welded (FSW) aluminum blanks offer larger preform sizes than available from the mill in single width plates. Large blanks made from aerospace-grade aluminum alloys are essential for making large domes and cones for launch vehicles and crew vehicles. Previous research efforts have demonstrated >5m FSW blanks can be formed into large domes or cones using spin forming, but the subsequent solution heat treatment consistently led to abnormal grain growth (AGG). The AGG phenomenon occurs because the FSW blanks must be spin formed in the -O annealed temper to withstand the high strain process without cracking. Subsequently, the solution heat treatment process leads to the onset of AGG, whereby a small number of grains grow much larger than the surrounding matrix. Ultimately, the large grains and bimodal grain-size distribution leads to unacceptable ductility in the short-transverse direction of the part. Alternatively, explosive hydroforming is a nascent process that can be applied to large and thick aerospace-grade aluminum alloy blanks in the -T3 temper and only require a post-forming artificial aging process to a -T8. Avoiding the solution heat treatment (SHT) process mitigates AGG, and also eliminates the need for very rare drop-furnace quench baths that can accommodate a part size larger than 5m in diameter. This presentation summarizes a Lockheed Martin-led effort that demonstrated explosive hydroforming of a 1.9m diameter, 39mm thick 2050 aluminum lithium FSW blank. Details will be provided on the process flow, mechanical properties, and metallographic analyses. This work was performed in partnership with NASA Langley Research Center (LaRC) under a NASA Announcement of Collaboration Opportunities program.

Friction stir welded blank

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

Towards an improved understanding of the Antarctic coastal zone and its contribution to future global sea level

Understanding the coastal zone of the Antarctic Ice Sheet, where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the coming centuries. The Antarctic Ice Sheet remains the largest source of uncertainty in future sea-level projections. Insufficient knowledge of bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, but is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography datasets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation dataset of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and enhancing predictions of future ocean and climate changes.

Kenichi Matsuoka

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

Coma Physics of an Interstellar Object: JWST Spatial-Spectral Mapping of 3I/ATLAS

We report a survey of molecular emission from cometary volatiles using the James Webb Space Telescope (JWST) toward interstellar object 3I/ATLAS carried out on UT 2025 December 22 and 23 at a heliocentric distance (\rh{}) of $2.37-2.41$ au. These measurements of CO, \ce{CO2}, \ce{H2O}, \ce{CH3OH}, and \ce{CH4} sampled molecular chemistry in 3I/ATLAS as it receded from its encounter with our Sun and entered the vicinity of the \ce{H2O} ice line --- the region between \rh{} = $2-3$ au where the temperature becomes too low for H$_2$O to vigorously sublime and CO and \ce{CO2} begin to control the overall activity. CO was the most abundant molecule, followed by \ce{H2O} and \ce{CO2}, whose molecular abundances with respect to CO were $(40.5\pm3.1)\%$ and ($41.6\pm0.3)\%$, respectively. This work presents spatial-spectral maps of column density and rotational temperature as a function of distance from the nucleus for all detected species. The spatial distributions of both quantities were highly anisotropic for the apolar species in the coma of 3I/ATLAS, yet were more nearly symmetric for the polar molecules. These results demonstrate how volatiles were segregated in the nucleus ices of 3I/ATLAS and reveal heating and cooling mechanisms in its coma. Derived maps of the ortho-to-para ratio (OPR) for \ce{H2O} were flat with increasing distance from the nucleus and consistent with a coma-averaged value $\mathrm{OPR}=2.7\pm0.2$, slightly less than the expected equilibrium value of three.

Nathan X Roth

Understanding and Controlling Hybrid Electric Gas Turbine Engine Transient Dynamics

The electrification of the gas turbine engine is known to increase the flexibility of aircraft architectures by enabling the generation of electrical power to distribute to other electrically based, thrust producing subsystems. It also has potential for direct performance benefits in the gas turbine engine itself, both at steady state and dynamically. Although the design focus of the gas turbine engine performance is primarily at steady state, it is often the instabilities occurring during transients that cause disequilibrium and constrain performance improvements. Instabilities arise due to the disequilibrium of the energy storage mechanisms within the traditional engine system and likewise for the electrified engine system. The primary energy storage mechanisms in the traditional system are the rotational inertia, gas path volumes, the thermal masses that make up the mechanical structure, and now, the electrical power system will provide additional contributions. Understanding the effect each of these energy storage mechanisms has on the others and controlling them appropriately allows for the suppression of state changes within the turbomachinery components to the degree that the components remain near steady state, thus reducing the disequilibrium within the system. The ability to tightly regulate the state changes of turbomachinery components, such as the compressor, minimizes the excursion of the compressor operating point from the operating line (operability), allowing for higher performing, more efficient compressor designs by decreasing the amount of stall margin needed for safe engine power level changes. Preliminary studies with the electrification of the turbine engine have shown that this is possible, and this concept can lead to design trades benefitting engine performance, weight, and volume. This paper explores in detail the energy storage mechanisms and control approaches for coordinating their state changes, and ultimately proposes that a higher performing, more efficient compressor design can result.

transient dynamics

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