Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape analysis”

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.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Size distribution and scattering phase function of aerosol particles retrieved from sky brightness measurements

Ground-based measurements of the solar transmission and sky radiance in a horizontal plane through the Sun are taken in several geographical regions and aerosol types: dust in a desert transition zone in Israel, sulfate particles in Eastern and Western Europe, tropical aerosol in Brazil, and mixed continental/maritime aerosol in California. Stratospheric aerosol was introduced after the eruption of Mount Pinatubo in June 1991. Therefore measurements taken before the eruption are used to analyze the properties of tropospheric aerosol; measurements from 1992 are also used to detect the particle size and concentration of stratospheric aerosol. The measurements are used to retrieve the size distribution and the scattering phase function at large scattering angles of the undisturbed aerosol particles. The retrieved properties represent an average on the entire atmospheric column. A comparison between the retrieved phase function for a scattering angle of 120 deg, with phase function predicted from the retrieved size distribution, is used to test the assumption of particle homogeneity and sphericity in radiative transfer models (Mie theory). The effect was found to be small (20% +/- 15%). For the stratospheric aerosol (sulfates), as expected, the phase function was very well predicted using the Mie theory. A model with a power law distribution, based on the spectral dependence of the optical thickness, alpha, cannot estimate accurately the phase function (up to 50% error for lambda = 0.87 microns). Before the Pinatubo eruption the ratio between the volumes of sulfate and coarse particles was very well correlated with alpha. The Pinatubo stratospheric aerosol destroyed this correlation. The aerosol optical properties are compared with analysis of the size, shape, and composition of the individual particles by electron microscopy of in situ samples. The measured volume size distribution before the injection of stratospheric aerosol consistently show two modes, sulfate particles with r(sub m) less than 0.2 microns and coarse paritcles with r(sub m) greater than 0.7 microns. The 'window' in the tropospheric aerosol in this radius range was used to observe a stable stratospheric aerosol in 1992, with r(sub m) approximately 0.5 microns. A combination of such optical thickness and sky measurements can be used to assess the direct forcing and the climatic impact of aerosol. Systematic inversion for the key aerosol types (sulfates, smoke, dust, and maritime aerosol) of the size distribution and phase function can give the relationship between the aerosol physical and optical properties that can be used to compute the radiative forcing. This forcing can be validated in dedicated field experiments.

Kaufman, Y. J.↗
NASA NTRS1995

FTIR instrumentation to monitor vapors from Shuttle tile waterproofing materials

The Space Shuttle Thermal Protection System (TPS) tiles and blankets are waterproofed using DimethylEthoxySilane (DMEX) in the Orbiter Processing Facilities (OPF). DMES has a Threshold Limit Value (TLV) for exposure of personnel to vapor concentration in air of 0.5 ppm. The OPF high bay cannot be opened for normal work after a waterproofing operation until the DMES concentration is verified by measurement to be below the TLV. On several occasions the high bay has been kept closed for up to 8 hours following waterproofing operations due to high DMES measurements. In addition, the Miran 203 and Miran 1 BX infrared analyzers calibrated at different wavelengths gave different readings under the same conditions. There was reason to believe that some of the high DMES concentration readings were caused by interference form water and ethanol vapors. The Toxic Vapor Detection Laboratory (TVDL) was asked to test the existing DMES instruments and identify the best qualified instrument. In addition the TVDL was requested to develop instrumentation to ensure the OPF high bay could be opened safely as soon as possible after a waterproofing operation. A Fourier Transform Infrared (FTIR) spectrophotometer instrument developed for an earlier project was reprogrammed to measure DMES vapor along with ethanol, water, and several common solvent vapors. The FTIR was then used to perform a series of laboratory and field tests to evaluate the performance of the single wavelength IR instruments in use. The results demonstrated that the single wavelength IR instruments did respond to ethanol and water vapors, more or less depending on the analytical IR wavelength selected. The FTIR was able to separate the responses to DMES, water and ethanol, and give consistent readings for the DMES vapor concentration. The FTIR was then deployed to the OPF to monitor real waterproofing operations. The FTIR was also used to measure the time for DMES to evaporate from TPS tile under a range of humidity conditions in controlled laboratory tests. The combination of laboratory and field tests with the FTIR instrument demonstrated superior sensitivity, ability to reject interference from water and ethanol vapors, ruggedness to be transported from the lab to the OPF and set up without special procedures or degradation of performance. The multiple component vapor analysis algorithm was developed at KSC and incorporates automatic baseline correction and shape fitting of the spectra. The analysis for DMES, TetraMethylDiSiloxane (TMDS), ethanol, methanol, isopropanol, and baseline parameters uses 161 points per sample at 4 cm(exp -1) resolution, and processes an eight scan sample every ten seconds. The standard deviation of the measurements is 0.013 ppm and the upper linear limit is 125 ppm DMES. Based on successful demonstration of capabilities we produced three mobile instrument carts to be used in each OPF to support future waterproofing operations. The design and building of the 'DMES Carts' were accomplished in Fiscal year 1995.

Mattson, C. B.↗
NASA NTRS1996

Dynamics and Statics of Nonaxisymmetric Liquid Bridges

Theoretical and experimental investigation of the stability of nonaxisymmetric and nonaxisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes, are conducted. Numerical analysis of the stability of nonaxisymmetric bridges between unequal disks for various orientations of the gravity vector is performed. Experimental and theoretical investigation of large (nonaxisymmetric) oscillations and breaking of liquid bridges are also conducted.

Alexander, J. Iwan D.↗
NASA NTRS1997

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS1999

Directional Solidification of Bi-Sn on USMP-4

The experiments used MEPHISTO hardware to study the solidification and melting behavior of bismuth alloyed with 1 at% tin. Three samples, each approximately 900 mm long and 6mm in diameter, were used. A portion of each sample also included a 2 mm diameter growth capillary, to assist in the formation of a single grain. One sample provided the Seebeck voltage generated during melting and freezing processes. Another provided temperature data and Peltier pulsed demarcation of the interface shape for post flight analysis. The third sample provided resistance and growth velocity measurements, as well as additional thermal data. The third sample was also quenched at the end of the mission to preserve the composition of the liquid near the interface for post flight determination. A total of 450mm of directionally solidified samples were preserved for post mission structural and compositional characterization. Substantial differences were observed in the Seebeck signal between the ground-based experiments and the space-based experiments. The temperature gradient in the liquid for the ground-based experiments was significantly lower than the temperature gradient in the liquid for the space-based experiments.

Abbaschian, Reza↗
NASA NTRS1993

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗
NASA NTRS1997

An Improved Treatment of External Boundary for Three-Dimensional Flow Computations

We present an innovative numerical approach for setting highly accurate nonlocal boundary conditions at the external computational boundaries when calculating three-dimensional compressible viscous flows over finite bodies. The approach is based on application of the difference potentials method by V. S. Ryaben'kii and extends our previous technique developed for the two-dimensional case. The new boundary conditions methodology has been successfully combined with the NASA-developed code TLNS3D and used for the analysis of wing-shaped configurations in subsonic and transonic flow regimes. As demonstrated by the computational experiments, the improved external boundary conditions allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable speedup of convergence of the multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS1996

Interface Configuration Experiments (ICE) Explore the Effects of Microgravity on Fluids

The Interface Configuration Experiment (ICE) is actually a series of experiments that explore the striking behavior of liquid-vapor interfaces (i.e., fluid surfaces) in a low gravity environment under which major shifts in liquid position can arise from small changes in container shape or contact angle. Although these experiments are designed to test current mathematical theory, there are numerous practical applications that could result from these studies. When designing fluid management systems for space-based operations, it is important to be able to predict the locations and configurations that fluids will assume in containers under low-gravity conditions. The increased ability to predict, and hence control, fluid interfaces is vital to systems and/or processes where capillary forces play a significant role both in space and on the Earth. Some of these applications are in general coating processes (paints, pesticides, printing, etc.), fluid transport in porous media (ground water flows, oil recovery, etc.), liquid propellant systems in space (liquid fuel and oxygen), capillary-pumped loops and heat pipes, and space-based life-support systems. In space, almost every fluid system is affected, if not dominated, by capillarity. Knowledge of the liquid-vapor interface behavior, and in particular the interface shape from which any analysis must begin, is required as a foundation to predict how these fluids will react in microgravity and on Earth. With such knowledge, system designs can be optimized, thereby decreasing costs and complexity, while increasing performance and reliability. ICE has increased, and will continue to increase this knowledge, as it probes the specific peculiarities of current theory upon which our current understanding of these effects is based. Several versions of ICE were conducted in NASA Lewis Research Center's drop towers and on the space shuttle during the first and second United States Microgravity Laboratory missions (USML-1 and USML-2). Additional tests are planned for the space shuttle and for the Russian Mir space station. These studies will focus on interfacial problems concerning surface existence, uniqueness, configuration, stability, and flow characteristics.

Source record↗
NASA NTRS1998

Liquid-Vapor Interface Configurations Investigated in Low Gravity

The Interface Configuration Experiment (ICE) is part of a multifaceted study that is exploring the often striking behavior of liquid-vapor interfaces in low-gravity environments. Although the experiment was posed largely as a test of current mathematical theory, applications of the results should be manifold. In space almost every fluid system is affected, if not dominated, by capillarity (the effects of surface tension). As a result, knowledge of fluid interface behavior, in particular an equilibrium interface shape from which any analysis must begin, is fundamental--from the control of liquid fuels and oxygen in storage tanks to the design and development of inspace thermal systems, such as heat pipes and capillary pumped loops. ICE has increased, and should continue to increase, such knowledge as it probes the specific peculiarities of current theory upon which our present understanding rests. Several versions of ICE have been conducted in the drop towers at the NASA Lewis Research Center, on the space shuttles during the first and second United States Microgravity Laboratory missions (USML-1 and USML-2), and most recently aboard the Russian Mir space station. These studies focused on interfacial problems concerning the existence, uniqueness, configuration, stability, and flow characteristics of liquid-vapor interfaces. Results to date have clearly demonstrated the value of the present theory and the extent to which it can predict the behavior of capillary systems.

Concus, Paul↗
NASA NTRS2005

XMM-Newton Observations of Four Millisecond Pulsars

I present an analysis of the XMM-Newton observations of four millisecond pulsars, J0437-4715, J2124-3358, J1024-0719, and J0034-0534. The new data provide strong evidence of thermal emission in the X-ray flux detected from the first three objects. This thermal component is best interpreted as radiation from pulsar polar caps covered with a nonmagnetic hydrogen atmosphere. A nonthermal power-law component, dominating at energies E greater than or equal to 3 keV, can also be present in the detected X-ray emission. For PSR J0437-4715, the timing analysis reveals that the shape and pulsed fraction of the pulsar light curves are energy dependent. This, together with the results obtained from the phase-resolved spectroscopy, supports the two-component (thermal plus nonthermal) interpretation of the pulsar's X-ray radiation. Highly significant pulsations have been found in the X-ray flux of PSRs 52124-3358 and 51024-0719. For PSR 50034-0534, a possible X-ray counterpart of the radio pulsar has been suggested. The inferred properties of the detected thermal emission are compared with predictions of radio pulsar models.

Zavlin, Vyacheslav E.↗
NASA NTRS2008

High-precision Non-Contact Measurement of Creep of Ultra-High Temperature Materials for Aerospace

For high-temperature applications (greater than 2,000 C) such as solid rocket motors, hypersonic aircraft, nuclear electric/thermal propulsion for spacecraft, and more efficient jet engines, creep becomes one of the most important design factors to be considered. Conventional creep-testing methods, where the specimen and test apparatus are in contact with each other, are limited to temperatures approximately 1,700 C. Development of alloys for higher-temperature applications is limited by the availability of testing methods at temperatures above 2000 C. Development of alloys for applications requiring a long service life at temperatures as low as 1500 C, such as the next generation of jet turbine superalloys, is limited by the difficulty of accelerated testing at temperatures above 1700 C. For these reasons, a new, non-contact creep-measurement technique is needed for higher temperature applications. A new non-contact method for creep measurements of ultra-high-temperature metals and ceramics has been developed and validated. Using the electrostatic levitation (ESL) facility at NASA Marshall Space Flight Center, a spherical sample is rotated quickly enough to cause creep deformation due to centrifugal acceleration. Very accurate measurement of the deformed shape through digital image analysis allows the stress exponent n to be determined very precisely from a single test, rather than from numerous conventional tests. Validation tests on single-crystal niobium spheres showed excellent agreement with conventional tests at 1985 C; however the non-contact method provides much greater precision while using only about 40 milligrams of material. This method is being applied to materials including metals and ceramics for non-eroding throats in solid rockets and next-generation superalloys for turbine engines. Recent advances in the method and the current state of these new measurements will be presented.

Rogers, Jan R.↗
NASA NTRS2010

Use of ILTV Control Laws for LaNCETS Flight Research

A report discusses the Lift and Nozzle Change Effects on Tail Shock (LaNCETS) test to investigate the effects of lift distribution and nozzle-area ratio changes on tail shock strength of an F-15 aircraft. Specific research objectives are to obtain inflight shock strength for multiple combinations of nozzle-area ratio and lift distribution; compare results with preflight prediction tools; and update predictive tools with flight results. The objectives from a stability and control perspective are to ensure adequate aircraft stability for the changes in lift distribution and plume shape, and ensure manageable transient from engaging and disengaging the ILTV research control laws. In order to change the lift distribution and plume shape of the F-15 aircraft, a decade-old Inner Loop Thrust Vectoring (ILTV) research control law was used. Flight envelope expansion was performed for the test configuration and flight conditions prior to the probing test points. The approach for achieving the research objectives was to utilize the unique capabilities of NASA's NF-15B-837 aircraft to allow the adjustment of the nozzle-area ratio and/or canard positions by engaging the ILTV research control laws. The ILTV control laws provide the ability to add trim command biases to canard positions, nozzle area ratios, and thrust vectoring through the use of datasets. Datasets consist of programmed test inputs (PTIs) that define trims to change the nozzle-area ratio and/or canard positions. The trims are applied as increments to the normally commanded positions. A LaNCETS non-linear, six-degrees-of-freedom simulation capable of realtime pilot-in-the-loop, hardware-in-the-loop, and non-real-time batch support was developed and validated. Prior to first flight, extensive simulation analyses were performed to show adequate stability margins with the changes in lift distribution and plume shape. Additionally, engagement/disengagement transient analysis was also performed to show manageable transients.

Moua, Cheng↗
NASA NTRS2016

Impact of Scattering Model on Disdrometer Derived Attenuation Scaling

NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP#5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 gigahertz attenuation from the disdrometer and the 20 gigahertz time-series as well as to directly measure the 40 gigahertz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data. In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer-derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.

Radio Transmission↗
NASA NTRS2016

Impact of Scattering Model on Disdrometer Derived Attenuation Scaling

NASA Glenn Research Center (GRC), the Air Force Research Laboratory (AFRL), and the Politecnico di Milano (POLIMI) are currently entering the third year of a joint propagation study in Milan, Italy utilizing the 20 and 40 GHz beacons of the Alphasat TDP5 Aldo Paraboni scientific payload. The Ka- and Q-band beacon receivers were installed at the POLIMI campus in June of 2014 and provide direct measurements of signal attenuation at each frequency. Collocated weather instrumentation provides concurrent measurement of atmospheric conditions at the receiver; included among these weather instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which records droplet size distributions (DSD) and droplet velocity distributions (DVD) during precipitation events. This information can be used to derive the specific attenuation at frequencies of interest and thereby scale measured attenuation data from one frequency to another. Given the ability to both predict the 40 GHz attenuation from the disdrometer and the 20 GHz timeseries as well as to directly measure the 40 GHz attenuation with the beacon receiver, the Milan terminal is uniquely able to assess these scaling techniques and refine the methods used to infer attenuation from disdrometer data.In order to derive specific attenuation from the DSD, the forward scattering coefficient must be computed. In previous work, this has been done using the Mie scattering model, however, this assumes a spherical droplet shape. The primary goal of this analysis is to assess the impact of the scattering model and droplet shape on disdrometer derived attenuation predictions by comparing the use of the Mie scattering model to the use of the T-matrix method, which does not assume a spherical droplet. In particular, this paper will investigate the impact of these two scattering approaches on the error of the resulting predictions as well as on the relationship between prediction error and rain rate.

statistics↗
NASA NTRS2015

Self- and Air-Broadened Line Shapes in the 2v3 P and R Branches of 12CH4

In this paper we report line shape parameters of 12CH4 for several hundred 2V(sub 3) transitions in the spectral regions 5891-5996 cm( exp -1) (P branch) and 6015-6115 cm(exp -1) (R branch). Air- and self-broadening coefficients were measured as a function of temperature; line mixing via off-diagonal relaxation matrix element coefficients was also obtained for 47 transition pairs. In total, nearly 1517 positions and intensities were retrieved, but many transitions were too weak for the line shape study. For this analysis, we used 25 high-resolution (0.0056 and 0.0067 cm(ex[ -1) and high signal-to-noise (S/N) spectra of high-purity 12CH4 and the same high-purity 12CH4 broadened by dry air recorded at different sample temperatures between 130 K and 295 K with the Bruker IFS 125HR Fourier transform spectrometer at JPL. Three different absorption cells were used (1) a White cell set to a path length of 13.09 m for room temperature data, (2) a single-pass 0.2038 m long coolable cell (for self-broadening) and (3) a multipass cell with 20.941 m total path coolable Herriott cell (for air-broadening). In total there were 13 spectra with pure 12CH4 (0.27-599 Torr) and 12 air-broadened spectra with total sample pressures of 80-805 Torr and volume mixing ratios (VMR) of methane between 0.18 and 1.0. An interactive multispectrum nonlinear least-squares technique was employed to fit the individual P10-P1 and R0-R10 manifolds in all the spectra simultaneously. Results obtained from the present analysis are compared to other recent measurements.

Devi, V. Malathy↗
NASA NTRS2022

Scattering By Relatively Small Oblate Spheroidal Drops of Water in the Rainbow Region: T-Matrix Results and Geometric Interpretation

Late in 1999 Michael Mishchenko carried out T-matrix calculations of light scattered by oblate spheroidal drops of water at Marston's request. The purpose of the investigation was to determine if the caustic features, known by that time to be visible in the far-field scattering near the primary rainbow region for horizontally illuminated relatively large drops, would also be present in T-matrix calculations of the scattering by small drops. The axis of rotational symmetry of the drop was taken to be vertical, and the illumination was taken to be horizontal. Scattering into the relevant angular region was evaluated for drop aspect ratios D/H = 1, 1.23, 1.311, and 1.37 where D denotes the equatorial diameter and H the vertical diameter. The results are limited to polarized monochromatic illumination with ka = 140 where a = D/2 and k = 2π/λ with λ being the wavelength of the illumination. As expected from geometrical analysis, a horizontal “V” shaped caustic is visible for D/H = 1.311. The T-matrix results, gray-tone intensity displays provided by Mishchenko, remain instructive.

Scattering by spheroids↗
NASA NTRS1975

Mean rates of the orbital elements of a satellite perturbed by a lens shaped mass concentration

Long arc gravity analysis of lunar orbiter tracking data in the past has been carried out with the help of averaged equations of motion, in which short period effects have been suppressed. This procedure has required that the harmonic terms in the gravity potential be averaged over an orbital period. In the present paper, this technique is extended to mass points and mass disks in the gravity field. This requires the evaluation of expressions for the mean rates of the orbit elements for a satellite perturbed by a lens shaped mass concentration. Corresponding expressions for the perturbations due to a mass point are obtained in the limit as the lens radius goes to zero. The derived equations have been programmed on the UNIVAC 1108 computer, and the results checked by numerical differencing.

Ananda, M.↗
NASA NTRSDate not supplied

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

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

Optimization↗