Search NASA⌕ Search

SEARCH · Search NASA

Results for “lumped parameter model”

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 163 records · Page 9

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗

Physics Informed Neural Nets for Systems Health Management

To facilitate and solve the prediction problem, awareness of the current health state of the system is key, since it is necessary to perform condition-based predictions. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditions. Development in data-driven algorithms in regression of complex nonlinear functions and classification tasks have generated a growing interest in artificial intelligence for industrial applications. Complex multi-physics models as well as digital twins, once purely built on physics and corresponding simplified lumped parameter iterations, can now benefit from machine learning algorithms to mitigate the lack of understanding of some complex behavior. The research work presents application of physics-informed neural nets application to a representative electric powertrain for unmanned aerial vehicles. The model is composed of physics-derived and empirical equations, integrated with connected networks that are strategically placed within the model to substitute equations that are subject to large uncertainty. Polynomial fit driven by heuristics or empirical observations can be substituted by more flexible networks that can minimize the error between model predictions and observations without being restricted to a predefined functional form. This modeling strategy allows training of networks deep inside the model and unknown parameters in a single learning stage.

Physics Informed↗

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

Building Thermal Models

This presentation is meant to be an overview of the model building process It is based on typical techniques (Monte Carlo Ray Tracing for radiation exchange, Lumped Parameter, Finite Difference for thermal solution) used by the aerospace industry This is not intended to be a "How to Use ThermalDesktop" course. It is intended to be a "How to Build Thermal Models" course and the techniques will be demonstrated using the capabilities of ThermalDesktop (TD). Other codes may or may not have similar capabilities. The General Model Building Process can be broken into four top level steps: 1. Build Model; 2. Check Model; 3. Execute Model; 4. Verify Results.

Modeling↗

On the performance of explicit and implicit algorithms for transient thermal analysis

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite element models of the configurations are discribed. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system and a model of the space shuttle orbiter wing. Calculations were carried out using the SPAR finite element program, the MITAS lumped parameter program and a special purpose finite element program incorporating the GEAR algorithms. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff. Careful attention to modeling detail such as avoiding thin or short high-conducting elements can sometimes reduce the stiffness to the extent that explicit methods become advantageous.

Adelman, H. M.↗

Recent advances in thermostructural finite element analysis

A comparison study is described which has indicated that the finite element and lumped parameter methods are comparable for steady-state conduction/convection analyses. The finite element model was advantageous because of the ease of model verification with computer graphics. Some problems encountered in constructing compatible thermal and structural models have indicated the need for further methodology developments to automate the transition between thermal and structural analyses. The need for more efficient finite element nonlinear transient algorithms has been cited. With further development in the thermal analysis area, the finite element method offers high potential for an integrated thermal/structural analysis capability.

Thornton, E. A.↗

Cascaded Perforates as One-Dimensional, Bulk Absorbers

Porous cell honeycomb liners for aircraft engine nacelles offer the possibility of exploiting extended reaction effects to improve liner attenuation bandwidth as generally attributed to the performance of bulk absorbers. This paper describes an analytical procedure, starting with an impedance prediction model for a single perforated plate, to estimate the bulk-absorber parameters for a cascade of such perforates - a first step to modeling a porous wall honeycomb structure. The objective is to build confidence in a lumped element impedance model, when applied to a uniformly-spaced set of porous plates to predict its .bulk. absorber properties. The model is based upon a modified version of the two-parameter flow resistance model of the form A + BV(sub inc), where A and B are physics-based, semi-empirical parameters that are adjusted to provide an optimum fit to a composite dataset from three plate porosities of 2.5, 5 and 10%. The composite dataset is achieved by reformulating the two-parameter flow resistance model into a .reduced pressure drop coefficient. dependency on perforate hole Reynolds number. The resulting impedance model is employed to calculate surface impedance spectra for N and 2N-layer perforate cascades. The well-known two-thickness method for experimental determination of bulk-absorber parameters is then applied to these .synthesized. data sets to predict the characteristic impedance and propagation constant for the perforate cascades. These results are then compared with experimental results reported in a companion paper.

Parrott, T. L.↗

Development status of SINDA/FLUINT and SINAPS

SINDA/FLUINT (Systems Improved Numerical Differencing Analyzer/Fluid Integrator, formerly SINDA '85) is a computer code used to analyze thermal/fluid systems that can be represented in lumped parameter form. In addition to conduction and radiation heat transfer, the code is capable of modeling both single- and two-phase flow networks, their associated hardware, and their heat transfer processes. In this paper, recent improvements to SINDA/FLUINT are described, as are those in progress that will be available in the fall of 1992 in Version 2.5. Also, a preview of planned enhancements is provided. This paper also introduces SINAPS (SINDA Application Programming System), a powerful graphical pre- and postprocessor that will also be available in the fall of 1992.

Cullimore, Brent A.↗

Fire Detection Tradeoffs as a Function of Vehicle Parameters

Fire survivability depends on the detection of and response to a fire before it has produced an unacceptable environment in the vehicle. This detection time is the result of interplay between the fire burning and growth rates; the vehicle size; the detection system design; the transport time to the detector (controlled by the level of mixing in the vehicle); and the rate at which the life support system filters the atmosphere, potentially removing the detected species or particles. Given the large differences in critical vehicle parameters (volume, mixing rate and filtration rate) the detection approach that works for a large vehicle (e.g. the ISS) may not be the best choice for a smaller crew capsule. This paper examines the impact of vehicle size and environmental control and life support system parameters on the detectability of fires in comparison to the hazard they present. A lumped element model was developed that considers smoke, heat, and toxic product release rates in comparison to mixing and filtration rates in the vehicle. Recent work has quantified the production rate of smoke and several hazardous species from overheated spacecraft polymers. These results are used as the input data set in the lumped element model in combination with the transport behavior of major toxic products released by overheating spacecraft materials to evaluate the necessary alarm thresholds to enable appropriate response to the fire hazard.

microgravity↗

Fire Detection Tradeoffs as a Function of Vehicle Parameters

Fire survivability depends on the detection of and response to a fire before it has produced an unacceptable environment in the vehicle. This detection time is the result of interplay between the fire burning and growth rates; the vehicle size; the detection system design; the transport time to the detector (controlled by the level of mixing in the vehicle); and the rate at which the life support system filters the atmosphere, potentially removing the detected species or particles. Given the large differences in critical vehicle parameters (volume, mixing rate and filtration rate) the detection approach that works for a large vehicle (e.g. the ISS) may not be the best choice for a smaller crew capsule. This paper examines the impact of vehicle size and environmental control and life support system parameters on the detectability of fires in comparison to the hazard they present. A lumped element model was developed that considers smoke, heat, and toxic product release rates in comparison to mixing and filtration rates in the vehicle. Recent work has quantified the production rate of smoke and several hazardous species from overheated spacecraft polymers. These results are used as the input data set in the lumped element model in combination with the transport behavior of major toxic products released by overheating spacecraft materials to evaluate the necessary alarm thresholds to enable appropriate response to the fire hazard.

microgravity↗

Improved test methods for determining lightning-induced voltages in aircraft

A lumped parameter transmission line with a surge impedance matching that of the aircraft and its return lines was evaluated as a replacement for earlier current generators. Various test circuit parameters were evaluated using a 1/10 scale relative geometric model. Induced voltage response was evaluated by taking measurements on the NASA-Dryden Digital Fly by Wire F-8 aircraft. Return conductor arrangements as well as other circuit changes were also evaluated, with all induced voltage measurements being made on the same circuit for comparison purposes. The lumped parameter transmission line generates a concave front current wave with the peak di/dt near the peak of the current wave which is more representative of lightning. However, the induced voltage measurements when scaled by appropriate scale factors (peak current or di/dt) resulting from both techniques yield comparable results.

Crouch, K. E.↗

Using Piezoelectric Devices to Transmit Power through Walls

A method denoted wireless acoustic-electric feed-through (WAEF) has been conceived for transmitting power and/or data signals through walls or other solid objects made of a variety of elastic materials that could be electrically conductive or nonconductive. WAEF would make it unnecessary to use wires, optical fibers, tubes, or other discrete wall-penetrating signal-transmitting components, thereby eliminating the potential for structural weakening or leakage at such penetrations. Avoidance of such penetrations could be essential in some applications in which maintenance of pressure, vacuum, or chemical or biological isolation is required. In a basic WAEF setup, a transmitting piezoelectric transducer on one side of a wall would be driven at resonance to excite ultrasonic vibrations in the wall. A receiving piezoelectric transducer on the opposite side of the wall would convert the vibrations back to an ultrasonic AC electric signal, which would then be detected and otherwise processed in a manner that would depend on the modulation (if any) applied to the signal and whether the signal was used to transmit power, data, or both. An electromechanical-network model has been derived as a computationally efficient means of analyzing and designing a WAEF system. This model is a variant of a prior model, known in the piezoelectric-transducer art as Mason's equivalent-circuit model, in which the electrical and mechanical dynamics, including electromechanical couplings, are expressed as electrical circuit elements that can include inductors, capacitors, and lumped-parameter complex impedances. The real parts of the complex impedances are used to account for dielectric, mechanical, and coupling losses in all components (including all piezoelectric-transducer, wall, and intermediate material layers). In an application to a three-layer piezoelectric structure, this model was shown to yield the same results as do solutions of the wave equations of piezoelectricity and acoustic propagation in their full complexity.

Sherrit, Stewart↗

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Kuo, K. K.↗

Combustion instability coupling with feed system acoustics

High frequency combustion instability has recently been observed by Rocketdyne in a 40K thrust methane/LOX combustion chamber. The oscillations had frequencies as high as 14,000 Hz with pressure amplitudes in the LOX dome of 500 psi at a chamber pressures of 2,000 psi. At this frequency the wave length associated with a period of oscillation is 2.3 inches in LOX and 1.4 inches in methane. These distances are comparable to the lengths of the injector elements which requires that acoustic waves be considered in the feed systems rather than using lumped parameters as is normally considered for feed system coupled oscillations. To expand the capability of existing models, the Feiler and Heidmann feed system coupled instability model was modified to include acoustic oscillations in the feed system. Similarly the vaporization controlled instability model of Heidmann and Wieber was modified to include flow oscillations that would be produced by feed system coupling. The major elements that control oscillations in a rocket combustion chamber are shown and discussed.

Priem, Richard J.↗

Experience in distributed parameter modeling of the Spacecraft Control Laboratory Experiment (SCOLE) structure

The Spacecraft Control Laboratory Experiment (SCOLE) configuration is used to compare exact and approximate solutions of the partial differential equations which define its structural dynamics. The need for a proof model for evaluating competing control laws demands that solutions be generated which not only exhibit accurate modal characteristics, but precise static deflections as well. Because precise pointing is required, the motion of the end bodies of the Shuttle-attached antenna must be known with great accuracy. Modal models are attractive because of their stable solutions but require hundreds of modes to obtain a static deflection accuracy of only one percent. Although proportional damping in bending agrees well with experimental results using the SCOLE experimental apparatus, modes which involve both torsion and bending differ significantly from proportional damping. A lumped mass model is used to generate exact static deflections, but only approximate modal characteristics. Asymptotic solutions to the distributed parameter system approximate very accurately the modal characteristics at high mode numbers. Ways are examined for refining the approximate solutions by applying a first-order variation and by employing singular perturbation techniques which are usually limited to ordinary differential equations. The most accurate solutions of the distributed parameter model of SCOLE are obtained by combining exact and asymptotic solutions.

Taylor, L. W.↗

Quantification of Dynamic Model Validation Metrics Using Uncertainty Propagation from Requirements

The Space Launch System, NASA's new large launch vehicle for long range space exploration, is presently in the final design and construction phases, with the first launch scheduled for 2019. A dynamic model of the system has been created and is critical for calculation of interface loads and natural frequencies and mode shapes for guidance, navigation, and control (GNC). Because of the program and schedule constraints, a single modal test of the SLS will be performed while bolted down to the Mobile Launch Pad just before the first launch. A Monte Carlo and optimization scheme will be performed to create thousands of possible models based on given dispersions in model properties and to determine which model best fits the natural frequencies and mode shapes from modal test. However, the question still remains as to whether this model is acceptable for the loads and GNC requirements. An uncertainty propagation and quantification (UP and UQ) technique to develop a quantitative set of validation metrics that is based on the flight requirements has therefore been developed and is discussed in this paper. There has been considerable research on UQ and UP and validation in the literature, but very little on propagating the uncertainties from requirements, so most validation metrics are "rules-of-thumb;" this research seeks to come up with more reason-based metrics. One of the main assumptions used to achieve this task is that the uncertainty in the modeling of the fixed boundary condition is accurate, so therefore that same uncertainty can be used in propagating the fixed-test configuration to the free-free actual configuration. The second main technique applied here is the usage of the limit-state formulation to quantify the final probabilistic parameters and to compare them with the requirements. These techniques are explored with a simple lumped spring-mass system and a simplified SLS model. When completed, it is anticipated that this requirements-based validation metric will provide a quantified confidence and probability of success for the final SLS dynamics model, which will be critical for a successful launch program, and can be applied in the many other industries where an accurate dynamic model is required.

Brown, Andrew M.↗

A Simple Model for Rotating Detonation Rocket Engine Sizing and Performance Estimates

A Rotating Detonation Rocket Engine (RDRE) model is described which characterizes the device as an essentially infinite number of circumferentially arranged, sequentially firing pulse detonation engine (PDE) tubes. The PDE tubes are in turn treated as lumped-parameter chambers, each of which executes a dynamic (i.e. time dependent) Atkinson cycle with an allowance for a finite Mach number during the refill portion of the cycle. The formulation results in several free parameters which are used to essentially tune the model such that it closely matches higher fidelity computational fluid dynamic RDRE simulations in terms of performance, flow rates, etc. The simplicity of the model allows for straightforward combination with component models (e.g. turbopumps, cooling jackets, etc.) as well as implementation within larger vehicle simulations. Such capability allows for realistic mission analyses and benefits studies which will help to determine the best application for the RDRE concept.

detonation↗

A Simple Model for Rotating Detonation Rocket Engine Sizing and Performance Estimates

A Rotating Detonation Rocket Engine (RDRE) model is described which characterizes the device as an essentially infinite number of circumferentially arranged, sequentially firing pulse detonation engine (PDE) tubes. The PDE tubes are in turn treated as lumped-parameter chambers, each of which executes a dynamic (i.e. time dependent) Atkinson cycle with an allowance for a finite Mach number during the refill portion of the cycle. The formulation results in several free parameters which are used to essentially tune the model such that it closely matches higher fidelity computational fluid dynamic RDRE simulations in terms of performance, flow rates, etc. The simplicity of the model allows for straightforward combination with component models (e.g. turbopumps, cooling jackets, etc.) as well as implementation within larger vehicle simulations. Such capability allows for realistic mission analyses and benefits studies which will help to determine the best application for the RDRE concept.

detonation↗