Search NASASearch

SEARCH · Search NASA

Results for “Polymer Matrix Composites”

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.

198 records · Page 10

Evaluation of Particle Degradation of a Ytterbium Disilicate Gas Turbine Coating in a Combustion Environment

Ytterbium disilicate (Yb2Si2O7) is a state-of-the-art topcoat material used in environmental barrier coatings (EBCs) to protect SiC-SiC composites from water vapor corrosion in gas turbine engines. However, it’s degradation by thermochemical and thermomechanical particle interactions may limit its use in next generation engine systems. Several static furnace tests have been performed on Yb2Si2O7 to elucidate methods of degradation, although dynamic testing in an engine-relevant environment has not been carried out. Ultimately, testing and evaluation of materials that incorporates all environmental damage modes will allow for an accurate assessment of long-term durability and operating lifetime. This work details the investigation of particle-induced erosion and corrosion of Yb2Si2O7 in NASA Glenn’s Erosion Burner Rig Facility. Coating durability was investigated as a function of test temperature, particle size, and deposition rate. Thermal cycling was also incorporated into testing to determine the interplay of heating/cooling cycles with particle degradation. Overall, these analyses will be used to determine coating robustness in an engine relevant environment.

Jamesa L Stokes

Evaluation of Particle Degradation of a Ytterbium Disilicate Gas Turbine Coating in a Combustion Environment

Ytterbium disilicate (Yb2Si2O7) is a state-of-the-art topcoat material used in environmental barrier coatings (EBCs) to protect SiC-SiC composites from water vapor corrosion in gas turbine engines. However, it’s degradation by thermochemical and thermomechanical particle interactions may limit its use in next generation engine systems. Several static furnace tests have been performed on Yb2Si2O7 to elucidate methods of degradation, although dynamic testing in an engine-relevant environment has not been carried out. Ultimately, testing and evaluation of materials that incorporates all environmental damage modes will allow for an accurate assessment of long-term durability and operating lifetime. This work details the investigation of particle-induced erosion and corrosion of Yb2Si2O7 in NASA Glenn’s Erosion Burner Rig Facility. Coating durability was investigated as a function of test temperature, particle size, and deposition rate. Thermal cycling was also incorporated into testing to determine the interplay of heating/cooling cycles with particle degradation. Overall, these analyses will be used to determine coating robustness in an engine relevant environment.

Jamesa L Stokes

List of Commercial and Advanced Developmental Niobium-Based Alloys of the Space Age

This report compiles a list of commercial and developmental niobium-based alloys developed during the Space Age (late 1950s through mid-1970s) for extreme-temperature applications, including rocket engine thrust chambers, hypersonic re-entry thermal protection systems, and space fission reactor loops. Niobium (Nb) was widely pursued because it provided the lowest density (~8.6 g/cc) among the primary refractory metals, a high melting temperature (~2470°C), exceptional low-temperature ductility, good formability, and compatibility with liquid alkali metals. An evaluation of physical metallurgy mechanisms, focusing on solid-solution strengthening via heavy refractory solutes (W, Mo, Ta), dual-purpose reactive solutes (Hf, Zr, Ti), and dispersion strengthening using carbides, nitrides, and oxides is presented. Additionally, the report compares Western and Soviet Union metallurgical approaches, explaining how supply chain factors and manufacturing infrastructure influenced element selection, interstitial chemistry, and alloy identification/naming conventions. Cataloging these historical alloy chemical compositions serves as a foundational reference for modern alloy additive manufacturing, thermodynamic CALPHAD modeling, and machine-learning discovery pipelines for next-generation extreme-environment niobium-based alloys.

Physical Metallurgy

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

NASA JSC’s Simulant Development Lab Capabilities and Artemis Testing

The Simulant Development Lab (SDL) is a multifunctional collaborative workspace that supports the development, curation, analysis, testing, and distribution of planetary regolith simulants – including lunar, Martian, asteroidal, and other granular materials. The lab provides a multidisciplinary setting for scientific characterization of simulant physical properties and for engineering evaluations conducted with simulant test beds. To enable this work, the SDL curates and maintains a stock of more than 35 metric tons of simulant material. To evaluate these materials and support testing goals, the lab is equipped with a comprehensive suite of processing tools and analytical instruments. These capabilities enable the SDL’s mission at NASA’s Johnson Space Center to distribute, develop, process, characterize, and test regolith simulants for mission relevant applications. Through controlled and repeatable testing environments that replicate the physical and compositional properties of lunar regolith, the SDL supports Artemis hardware maturation, providing safe, Earth‑based analogs for evaluating systems that must withstand regolith dust interactions, physical wear and abrasion, and operational loads. The facility’s extensive simulant inventory and integrated geological and engineering test infrastructure accelerate technology readiness for Artemis and future exploration campaigns (e.g., future crewed or robotic missions to Mars).

Simulant Development Lab

The Alpha-Proton-X-Ray Spectrometer Deployment Mechanism: an Anthropomorphic Approach to Sensor Placement on Martian Rocks and Soil

On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Richard S Blomquist

Microgravity Game Changers: Recent Accomplishments of Microgravity Assisted Materials Experiments on the ISS

This paper discusses important new materials, devices and technology produced in low-Earth orbit (LEO) that have been developed with funding from NASA’s In Space Production Applications (InSPA) Portfolio. In coordination with the International Space Station (ISS) National Laboratory, InSPA supports innovative small businesses and other research & development organizations to develop microgravity-assisted applications that make significant impacts in a wide variety of industries on Earth. These technology development efforts are using the microgravity of orbital space to overcome gravity-induced defects on Earth to create high-value products for terrestrial markets in the areas of semiconductors, pharmaceuticals, and communications. Through 24/7/365 operation on ISS for 25 years, coupled with commercial space transportation and services to, in, and from orbit, the microgravity environment has enabled pioneering science and significant commercial applications, including novel medical advances. World-leading researchers, innovators, and space payload developers provide compelling evidence through demonstrations on the ISS that materials processing in microgravity improves outcomes for a wide range of materials and products at all levels of the organization of matter, from quantum entities such as Bose-Einstein condensates, to crystals, alloys, composites, thin films, and photonics, to stem cells, medical diagnoses, medical treatments, and medical devices for humans. As of 2026, numerous examples proved by demonstration on ISS show the potential of these technologies in various fields. In this paper, a summary of four important experiments is presented.

microgravity medical devices

Aerothermal Shape Optimization of Actively-Cooled Battery Packs Using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

heat transfer

Aerothermal Shape Optimization of Actively-Cooled Battery Packs using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

thermal management

The Canadian Galactic Emission Mapper: A New 8-10 GHz Telescope To Map the Polarization of the Northern Sky

The Canadian Galactic Emission Mapper (CGEM) radio telescope is mapping polarized Galactic foregrounds from 8-10 GHz to aid in the search for B-modes in the Cosmic Microwave Background (CMB). Here we describe the design and early on-sky characterization of the CGEM optics. CGEM employs an on-axis, hat feed optical design that is a body of revolution (BOR) and exhibits excellent intrinsic polarization purity. We describe how we have optimized the optics with a novel framework that can directly minimize simulated intensity to polarized intensity (T → P) leakage in angular power spectrum space. The optimized optics exhibit simulated T → P leakage that is orders of magnitude below anticipated CℓBB for r = 10 −3 when scaled from 8-10 GHz to the CMB observing window near 95 GHz. We go on to describe the mechanical design of the optics, including a novel secondary mirror support made from Astroquartz composite that has low loss, low dielectric constant, and preserves the BOR symmetry of the optics. We then showcase the early on-sky performance of the optics with observations of the Sun and satellites, which probe the beam to 40dB down from the peak. A beam model for the deployed CGEM optics based on electromagnetic simulations is in excellent agreement with the on-sky data.

cosmic microwave background

Isotopic Signature of Organic Molecules from Beyond the Solar System: An Enriched Methane D/H Ratio in the Interstellar Object 3I/ATLAS

Interstellar objects are interlopers from other planetary systems, and their volatile compositions provide a glimpse into planet formation around their host star. We present near-infrared spectra of the coma of interstellar object 3I/ATLAS measured with the James Webb Space Telescope. Our results demonstrate an unexpectedly high D/H = (3.31 ± 0.34)% for methane and represent an exceedingly rare detection of deuterated organic molecules in an interstellar object. This D/H ratio is a factor of 14 ± 2 higher than that measured in comet 67P/Churyumov-Gerasimenko by the Rosetta spacecraft, the only other comet for which CH3D has been detected, yet the ratio of deuteration in methane compared with water is consistent for both comets within 1.2σ. The D/H ratio in methane is observationally unconstrained in extrasolar sources to date, but the enriched ratio in 3I/ATLAS is most similar to those measured in other organic molecules toward primitive environments. The high D/H ratios of water and methane in 3I/ATLAS are a natural consequence of formation in a high D/H elemental ratio environment as a result of locally cold conditions in the protoplanetary disk and prior interstellar cloud. Thus, 3I/ATLAS formed in an environment very different from that in which our Sun and planets originated.

Interstellar Object 3I/ATLAS

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

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

computational materials

The Effect of Gravity on the Combustion of Bulk Metals

In recent years, metal combustion studies at the University of Colorado have focused on the effects of gravity (g) on the ignition and burning behavior of bulk metals. The impetus behind this effort is the understanding of the ignition conditions and flammability properties of structural metals found in oxygen (O,sub>2 ) systems for space applications. Since spacecraft are subjected to higher-than-1g loads during launch and reentry and to a zero-gravity environment while in orbit, the study of ignition and combustion of bulk metals at different gravitational accelerations is of great practical concern. From the scientific standpoint, studies conducted under low gravity conditions provide simplified boundary conditions, since buoyancy is removed, and make possible the identification of fundamental ignition and combustion mechanisms. This investigation is intended to provide experimental verification of the influence of natural convection on the burning behavior of metals. In addition, the study offers the first findings of the influence of gravity on ignition of bulk metals and on the combustion mechanism and structure of metal-oxygen, vapor-phase diffusion flames in a buoyancy-free environment. Titanium (Ti) and magnesium (Mg) metals were chosen because of their importance as elements of structural materials and their simple chemical composition-pure metals instead of multicomponent alloys to simplify chemical and spectroscopic analyses. In addition, these elements present the two different combustion modes observed in metals: heterogeneous or surface burning (for Ti) and homogeneous or gas-phase reaction (for Mg). Finally, Mg, Ti, and their oxides exhibit a wide range of thermophysical and chemical properties. Metal surface temperature profiles, critical and ignition temperatures, propagation rates, burning times, and spectroscopic measurements are obtained under normal and reduced gravity. Visual evidence of all phenomena is provided by high-speed photography.

Melvyn C Branch

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

Designing a Propylene-Glycol Coolant Servicer System for Gateway’s Internal Active Thermal Control System

A human spacecraft ATCS—especially one using single-phase coolant loops exposed to cabin atmospheric conditions—requires periodic degassing and refilling to support long-duration missions of 15 to 30 years. During initial fill operations, system maintenance, gas permeation, and quick-disconnect mating or de-mating, small amounts of gas can gradually enter the coolant system over time. This can lead to degraded heat transfer performance, pump cavitation, and potentially pump vapor lock if a significant gas volume accumulates over time. Additionally, system leaks and routine fluid sampling can gradually reduce accumulator volumes to unacceptable levels, requiring periodic refills. In more severe cases, catastrophic changes in fluid composition may necessitate emergency draining, refilling, and degassing to ensure continued system functionality. These risks were identified and mitigated on the ISS ITCS through the development of a dual-membrane degasser ORU and a Fluid Servicer System ORU for coolant refilling. These systems were developed for the fully water-based ISS ITCS Coolant[1]. The Gateway space station, the first permanent human habitat in lunar orbit, uses a propylene-glycol/water coolant mixture, which has significantly different fluid properties compared to pure water. Because microgravity degassing technologies are sensitive to fluid surface tension and viscosity, existing ISS hardware is not suitable for servicing a propylene-glycol-based TCS. Unlike the ISS, the Gateway will operate in a higher-radiation environment and must meet stricter mass constraints due to its location outside of low earth orbit. This Government Furnished Equipment (GFE) flight hardware project aims to develop a lightweight, radiation-resistant Coolant Servicer System (CSS) capable of degassing and refilling Gateway’s propylene-glycol water-based IATCS.

Propylene Glycol Water

Designing a Propylene-Glycol Coolant Servicer System for Gateway’s Internal Active Thermal Control System

A human spacecraft ATCS—especially one using single-phase coolant loops exposed to cabin atmospheric conditions—requires periodic degassing and refilling to support long-duration missions of 15 to 30 years. During initial fill operations, system maintenance, gas permeation, and quick-disconnect mating or de-mating, small amounts of gas can gradually enter the coolant system over time. This can lead to degraded heat transfer performance, pump cavitation, and potentially pump vapor lock if a significant gas volume accumulates over time. Additionally, system leaks and routine fluid sampling can gradually reduce accumulator volumes to unacceptable levels, requiring periodic refills. In more severe cases, catastrophic changes in fluid composition may necessitate emergency draining, refilling, and degassing to ensure continued system functionality. These risks were identified and mitigated on the ISS ITCS through the development of a dual-membrane degasser ORU and a Fluid Servicer System ORU for coolant refilling. These systems were developed for the fully water-based ISS ITCS Coolant[1]. The Gateway space station, the first permanent human habitat in lunar orbit, uses a propylene-glycol/water coolant mixture, which has significantly different fluid properties compared to pure water. Because microgravity degassing technologies are sensitive to fluid surface tension and viscosity, existing ISS hardware is not suitable for servicing a propylene-glycol-based TCS. Unlike the ISS, the Gateway will operate in a higher-radiation environment and must meet stricter mass constraints due to its location outside of low earth orbit. This Government Furnished Equipment (GFE) flight hardware project aims to develop a lightweight, radiation-resistant Coolant Servicer System (CSS) capable of degassing and refilling Gateway’s propylene-glycol water-based IATCS.

Propylene Glycol Water

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