SEARCH · Engineering Papers

Results for “Finite-element”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Structural Configuration Systems Analysis for Advanced Aircraft Fuselage Concepts

Structural configuration analysis of an advanced aircraft fuselage concept is investigated. This concept is characterized by a double-bubble section fuselage with rear mounted engines. Based on lessons learned from structural systems analysis of unconventional aircraft, high-fidelity finite-element models (FEM) are developed for evaluating structural performance of three double-bubble section configurations. Structural sizing and stress analysis are applied for design improvement and weight reduction. Among the three double-bubble configurations, the double-D cross-section fuselage design was found to have a relatively lower structural weight. The structural FEM weights of these three double-bubble fuselage section concepts are also compared with several cylindrical fuselage models. Since these fuselage concepts are different in size, shape and material, the fuselage structural FEM weights are normalized by the corresponding passenger floor area for a relative comparison. This structural systems analysis indicates that an advanced composite double-D section fuselage may have a relative structural weight ratio advantage over a conventional aluminum fuselage. Ten commercial and conceptual aircraft fuselage structural weight estimates, which are empirically derived from the corresponding maximum takeoff gross weight, are also presented and compared with the FEM- based estimates for possible correlation. A conceptual full vehicle FEM model with a double-D fuselage is also developed for preliminary structural analysis and weight estimation.

Mukhopadhyay, Vivek↗

NASA NTRS2016

Applications of Displacement Transfer Functions to Deformed Shape Predictions of the G-III Swept-Wing Structure

In support of the Adaptive Compliant Trailing Edge [ACTE] project at the NASA Armstrong Flight Research Center, displacement transfer functions were applied to the swept wing of a Gulfstream G-III airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) to obtain deformed shape predictions. Four strainsensing lines (two on the lower surface, two on the upper surface) were used to calculate the deformed shape of the G III wing under bending and torsion. There being an insufficient number of surface strain sensors, the existing G III wing box finite element model was used to generate simulated surface strains for input to the displacement transfer functions. The resulting predicted deflections have good correlation with the finite-element generated deflections as well as the measured deflections from the ground load calibration test. The convergence study showed that the displacement prediction error at the G III wing tip can be reduced by increasing the number of strain stations (for each strain-sensing line) down to a minimum error of l.6 percent at 17 strain stations; using more than 17 strain stations yielded no benefit because the error slightly increased to 1.9% when 32 strain stations were used.

displacement theory↗

NASA NTRS2016

Buckling Design and Imperfection Sensitivity of Sandwich Composite Launch-Vehicle Shell Structures

Composite materials are increasingly being considered and used for launch-vehicle structures. For shell structures, such as interstages, skirts, and shrouds, honeycomb-core sandwich composites are often selected for their structural efficiency. Therefore, it is becoming increasingly important to understand the structural response, including buckling, of sandwich composite shell structures. Additionally, small geometric imperfections can significantly influence the buckling response, including considerably reducing the buckling load, of shell structures. Thus, both the response of the theoretically perfect structure and the buckling imperfection sensitivity must be considered during the design of such structures. To address the latter, empirically derived design factors, called buckling knockdown factors (KDFs), were developed by NASA in the 1960s to account for this buckling imperfection sensitivity during design. However, most of the test-article designs used in the development of these recommendations are not relevant to modern launch-vehicle constructions and material systems, and in particular, no composite test articles were considered. Herein, a two-part study on composite sandwich shells to (1) examine the relationship between the buckling knockdown factor and the areal mass of optimized designs, and (2) to interrogate the imperfection sensitivity of those optimized designs is presented. Four structures from recent NASA launch-vehicle development activities are considered. First, designs optimized for both strength and stability were generated for each of these structures using design optimization software and a range of buckling knockdown factors; it was found that the designed areal masses varied by between 6.1% and 19.6% over knockdown factors ranging from 0.6 to 0.9. Next, the buckling imperfection sensitivity of the optimized designs is explored using nonlinear finite-element analysis and the as-measured shape of a large-scale composite cylindrical shell. When compared with the current buckling design recommendations, the results suggest that the current recommendations are overly conservative and that the development of new recommendations could reduce the acreage areal mass of many composite sandwich shell designs by between 4% and 19%, depending on the structure.

Schultz, Marc R.↗

NASA NTRS2016

Benchmarking of Computational Models for NDE and SHM of Composites

Ultrasonic wave phenomena constitute the leading physical mechanism for nondestructive evaluation (NDE) and structural health monitoring (SHM) of solid composite materials such as carbon-fiber-reinforced polymer (CFRP) laminates. Computational models of ultrasonic guided-wave excitation, propagation, scattering, and detection in quasi-isotropic laminates can be extremely valuable in designing practically realizable NDE and SHM hardware and software with desired accuracy, reliability, efficiency, and coverage. This paper presents comparisons of guided-wave simulations for CFRP composites implemented using three different simulation codes: two commercial finite-element analysis packages, COMSOL and ABAQUS, and a custom code implementing the Elastodynamic Finite Integration Technique (EFIT). Comparisons are also made to experimental laser Doppler vibrometry data and theoretical dispersion curves.

composite materials↗

NASA NTRS2017

Curved Displacement Transfer Functions for Geometric Nonlinear Large Deformation Structure Shape Predictions

For shape predictions of structures under large geometrically nonlinear deformations, Curved Displacement Transfer Functions were formulated based on a curved displacement, traced by a material point from the undeformed position to deformed position. The embedded beam (depth-wise cross section of a structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution could be described with a piecewise linear or a piecewise nonlinear function. The discretization approach enabled piecewise integrations of the embedded-beam curvature equations to yield the Curved Displacement Transfer Functions, expressed in terms of embedded beam geometrical parameters and surface strains. By entering the surface strain data into the Displacement Transfer Functions, deflections along each embedded beam can be calculated at multiple points for mapping the overall structural deformed shapes. Finite-element linear and nonlinear analyses of a tapered cantilever tubular beam were performed to generate linear and nonlinear surface strains and the associated deflections to be used for validation. The shape prediction accuracies were then determined by comparing the theoretical deflections with the finiteelement- generated deflections. The results show that the newly developed Curved Displacement Transfer Functions are very accurate for shape predictions of structures under large geometrically nonlinear deformations.

displacement theory↗

NASA NTRS2016

Influence of Population Variation of Physiological Parameters in Computational Models of Space Physiology

The earliest manifestations of Visual Impairment and Intracranial Pressure (VIIP) syndrome become evident after months of spaceflight and include a variety of ophthalmic changes, including posterior globe flattening and distension of the optic nerve sheath. Prevailing evidence links the occurrence of VIIP to the cephalic fluid shift induced by microgravity and the subsequent pressure changes around the optic nerve and eye. Deducing the etiology of VIIP is challenging due to the wide range of physiological parameters that may be influenced by spaceflight and are required to address a realistic spectrum of physiological responses. Here, we report on the application of an efficient approach to interrogating physiological parameter space through computational modeling. Specifically, we assess the influence of uncertainty in input parameters for two models of VIIP syndrome: a lumped-parameter model (LPM) of the cardiovascular and central nervous systems, and a finite-element model (FEM) of the posterior eye, optic nerve head (ONH) and optic nerve sheath. Methods: To investigate the parameter space in each model, we employed Latin hypercube sampling partial rank correlation coefficient (LHSPRCC) strategies. LHS techniques outperform Monte Carlo approaches by enforcing efficient sampling across the entire range of all parameters. The PRCC method estimates the sensitivity of model outputs to these parameters while adjusting for the linear effects of all other inputs. The LPM analysis addressed uncertainties in 42 physiological parameters, such as initial compartmental volume and nominal compartment percentage of total cardiac output in the supine state, while the FEM evaluated the effects on biomechanical strain from uncertainties in 23 material and pressure parameters for the ocular anatomy. Results and Conclusion: The LPM analysis identified several key factors including high sensitivity to the initial fluid distribution. The FEM study found that intraocular pressure and intracranial pressure had dominant impact on the peak strains in the ONH and retro-laminar optic nerve, respectively; optic nerve and lamina cribrosa stiffness were also important. This investigation illustrates the ability of LHSPRCC to identify the most influential physiological parameters, which must therefore be well-characterized to produce the most accurate numerical results.

biomechanics↗

NASA NTRS2017

A Dynamic Calibration Method for Experimental and Analytical Hub Load Comparison

This paper presents the results from an ongoing effort to produce improved correlation between analytical hub force and moment prediction and those measured during wind-tunnel testing on the Aeroelastic Rotor Experimental System (ARES), a conventional rotor testbed commonly used at the Langley Transonic Dynamics Tunnel (TDT). A frequency-dependent transformation between loads at the rotor hub and outputs of the testbed balance is produced from frequency response functions measured during vibration testing of the system. The resulting transformation is used as a dynamic calibration of the balance to transform hub loads predicted by comprehensive analysis into predicted balance outputs. In addition to detailing the transformation process, this paper also presents a set of wind-tunnel test cases, with comparisons between the measured balance outputs and transformed predictions from the comprehensive analysis code CAMRAD II. The modal response of the testbed is discussed and compared to a detailed finite-element model. Results reveal that the modal response of the testbed exhibits a number of characteristics that make accurate dynamic balance predictions challenging, even with the use of the balance transformation.

Kreshock, Andrew R.↗

NASA NTRS2017

Microstructure Based Material-Sand Particulate Interactions and Assessment of Coatings for High Temperature Turbine Blades

Gas turbine engines for military/commercial fixed-wing and rotary wing aircraft use thermal barrier coatings in the high-temperature sections of the engine for improved efficiency and power. The desire to further make improvements in gas turbine engine efficiency and high power-density is driving the research and development of thermal barrier coatings and the effort of improving their tolerance to fine foreign particulates that may be contained in the intake air. Both commercial and military aircraft engines often are required to operate over sandy regions such as in the Middle-East nations, as well as over volcanic zones. For rotorcraft gas turbine engines, the sand ingestion is adverse during take-off, hovering near ground, and landing conditions. Although, most of the rotorcraft gas turbine engines are fitted with inlet particle separators, they are not 100 percent efficient in filtering fine sand particles of size 75 microns or below. The presence of these fine solid particles in the working fluid medium has an adverse effect on the durability of turbine blade thermal barrier coatings and overall performance of the engine. Typical turbine blade damages include blade coating wear, sand glazing, Calcia-Magnesia-Alumina-Silicate (CMAS) attack, oxidation, plugged cooling holes, all of which can cause rapid performance deterioration including loss of aircraft. The objective of this research is to understand the fine particle interactions with typical ceramic coatings of turbine blades at the microstructure level. A finite-element based microstructure modeling and analysis has been performed to investigate particle-surface interactions, and restitution characteristics. Experimentally, a set of tailored thermal barrier coatings and surface treatments were down-selected through hot burner rig tests and then applied to first stage nozzle vanes of the Gas Generator Turbine of a typical rotorcraft gas turbine engine. Laser Doppler velocity measurements were performed during hot burner rig testing to determine sand particle incoming velocities and their rebound characteristics upon impact on coated material targets. Further, engine sand ingestion tests were carried out to test the CMAS tolerance of the coated nozzle vanes. The findings from this on-going collaborative research to develop the next-gen sand tolerant coatings for turbine blades are presented in this paper.

Engine Tests↗

NASA NTRS2018

Buckling Response of a Large-Scale, Seamless, Orthogrid-Stiffened Metallic Cylinder

Results from the buckling test of a compression-loaded 8-ft-diameter seamless (i.e., without manufacturing joints), orthogrid-stiffened metallic cylinder are presented. This test was used to assess the buckling response and imperfection sensitivity characteristics of a seamless cylinder. In addition, the test article and test served as a technology demonstration to show the application of the flow forming manufacturing process to build more efficient buckling-critical structures by eliminating the welded joints that are traditionally used in the manufacturing of large metallic barrels. Pretest predictions of the cylinder buckling response were obtained using a finite-element model that included measured geometric imperfections. The buckling load predicted using this model was 697,000 lb, and the test article buckled at 743,000 lb (6% higher). After the test, the model was revised to account for measured variations in skin and stiffener geometry, nonuniform loading, and material properties. The revised model predicted a buckling load of 754,000 lb, which is within 1.5% of the tested buckling load. In addition, it was determined that the load carrying capability of the seamless cylinder is approximately 28% greater than a corresponding cylinder with welded joints.

Rudd, Michelle Tillotson↗

NASA NTRS2018

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

NASA NTRS2018

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗

NASA NTRS2018

Sensitivity Analysis for Multidisciplinary Systems (SAMS)

This report describes the research conducted under an interagency collaboration agreement between the Aerospace Systems Directorate of the Air Force Research Laboratory (AFRL/RQ) and the Computational AeroSciences Branch of NASA Langley (NASA LaRC). Both organizations have a long-term goal of developing a modular computational system for coupling fluids and structures to enable both analysis and optimization of aerospace vehicles. Ultimately, the system should support multiple solvers within the fluid and structure domains to allow the best combination for the task at hand, as well as to allow for institutional preferences of specific software components. Towards this goal, the current research was focused on enhancing the existing modal aeroelastic analysis in the NASA FUN3D (Fully-UNstructured three-dimensional CFD (Computational Fluid Dynamics) code) software (Biedron et al. 2018), as well as developing new aeroelastic analysis and optimization capabilities based on a non-linear finite-element method. The methods and enhancements described in this document pertain to FUN3D Version 13.4.

Biedron, Robert T.↗

NASA NTRS2018

Comet Hitchhiker: NIAC Phase 1 Final Report

Summary of Activities-Developed the Comet Hitchhiker concept, which is to hitch rides on small bodies (asteroids and comets) using a tethered spacecraft. (Section 2)-Identified five scientifically important missions that would be enabled or significantly benefited by the Comet Hitchhiker concept.The five mission concepts are: KBO rendezvous, Centaur rendezvous, Trojan rendezvous, Damocloid rendezvous, and Main asteroid belt tour to rendezvous with multiple (10) objects. (Section 3)-Derived the Space Hitchhike Equation, or "the rocket equation for hitchhiker", which relates the specific strength of tether, mass ratio, and V. (Section 4.1)-Performed in-depth feasibility analysis of the critical components of the concept through: Finite-element simulations of tether and spacecraft dynamics, as shown in Figure 1 (Section 4.4); Supercomputer simulations of the hypervelocity impact of harpoon on a small body, as shown in Figure 2. (Section 6)-Performed public outreach activities including the collaboration with a concept artist of the Museum of Science Fiction, exposure to media, and public presentations. (Section 8)

Deep Space↗

NASA NTRS2015

Aerodynamic Load Analysis of a Variable Camber Continuous Trailing Edge Flap System on a Flexible Wing Aircraft

This paper presents an investigative aerodynamic analysis conducted on the novel control surface known as a Variable Camber Continuous Trailing Edge Flap (VCCTEF). The VCCTEF is modeled as a control effector on the NASA Generic Transport Model (GTM) where wing flexibility is considered. Aerodynamic modeling of the aircraft is conducted using vortex-lattice method (VLM), and an aeroelastic model of the aircraft that utilizes a coupled finite-element analysis (FEA) vortex-lattice solution is employed. VLM solutions are used to determine quasi-steady aerodynamic loading over the aeroelastic wing structures with VCCTEF. The load data is used to calculate aerodynamic sensitivities to control surface deployment and is also integrated to determine overall hinge moments. This analysis is conducted for different flight conditions, where control sensitivities compare the VCCTEF effectiveness against conventional control surfaces. Hinge moment results provide insight into aeroelastic wing loads, and worst case hinge moments for the VCCTEF can be estimated. Results show that the VCCTEF offers greater control authority than conventional ailerons at cruise, and nominal flap settings are determined for a low-speed take-off condition where the VCCTEF maintains comparable control effectiveness. Worst case hinge moment values are presented, where for a possible VCCTEF configuration, the stiff wing model demonstrates up to a 5.90 percent increase in flap hinge moment relative to a rigid model, and a reduced stiffness model demonstrates up to an 11.42 percent increase.

Aeroelasticity↗

NASA NTRS2017

Magnetic Shield Design Modeling and Validation on SWOT Spacecraft Applied to Mars Flux Pinning Orbiting Sample Design

A modeling methodology was developed for use on the Surface Water and Ocean Topography (SWOT) mission for the computational modeling and design of a magnetic shield for a 63 A-sq m source. Shield options were modeled and tested across various design parameters and validated with measurement. Measurement results fell within 10% of simulation in most cases of concern with sources of error well understood. These methods results informed a subsequent modeling activity for a future Mars Sample Return (MSR) mission concept with stringent magnetic cleanliness requirements to retain geologic integrity of the samples. One option for capturing the sample from orbit requires rare earth magnets in close proximity to the Martian samples. Magnetic modeling with a finite-element method solver estimated the magnetic environment and established the need for shielding. Further modeling then determined the shielding necessary to meet magnetic requirements, followed by the design of a mass-optimized solution.

Gonzales, Edward↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

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↗