Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal system modeling”

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 865 records · Page 48

Overview of Glenn Mechanical Components Branch Research

Mr. James Zakrajsek, chief of the Mechanical Components Branch, gave an overview of research conducted by the branch. Branch members perform basic research on mechanical components and systems, including gears and bearings, turbine seals, structural and thermal barrier seals, and space mechanisms. The research is focused on propulsion systems for present and advanced aerospace vehicles. For rotorcraft and conventional aircraft, we conduct research to develop technology needed to enable the design of low noise, ultra safe geared drive systems. We develop and validate analytical models for gear crack propagation, gear dynamics and noise, gear diagnostics, bearing dynamics, and thermal analyses of gear systems using experimental data from various component test rigs. In seal research we develop and test advanced turbine seal concepts to increase efficiency and durability of turbine engines. We perform experimental and analytical research to develop advanced thermal barrier seals and structural seals for current and next generation space vehicles. Our space mechanisms research involves fundamental investigation of lubricants, materials, components and mechanisms for deep space and planetary environments.

Zakrajsek, James↗

Model Validation for Bigelow Expandable Activity Module (BEAM) with Stowage

The Bigelow Expandable Activity Module (BEAM) was a collaboration between NASA and Bigelow Aerospace to develop an expandable habitat technology that could be used for future space exploration missions such as the Artemis program. On April 26, 2016, the BEAM was launched and berthed to the International Space Station (ISS) on the aft port of Node 3. After the originally planned two-year mission life, BEAM was approved for use as a stowage module for a 109-cargotransfer bag equivalent (CTBE) to alleviate stowage limitations onboard the ISS. The present work builds on the previous 2017 model validation to flight data that demonstrated that BEAM could generally accommodate cargo stowage without adversely affecting BEAM thermal control, primarily for condensation management. Once the stowage was configured in early 2019, model validation efforts resumed. The thermal model revisions included the CTBE stowage and updated CFD-derived heat transfer coefficients for various intra-module ventilation (IMV) flow rates. The coefficients were provided by the Crew and Thermal Systems Division (EC) at NASA JSC and Jacobs Technology. The model is now revalidated and able to predict minimum temperatures with good agreement to flight temperatures (within 1 °C) for both historic and recurring minimum values. Additionally, the model is partially validated for reduced flow rates. With this model, the team was able to provide new flight rule recommendations for condensation and IMV operations management.

Sydney Taylor↗

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗

System Modeling of a Lunar Molten Regolith Electrolysis Plant

Introduction: In-Situ Resource Utilization (ISRU) is the process of extracting local resources to produce commodities for propulsion, life support systems, and off-planet construction rather than transporting consumables from Earth. Molten Regolith Electrolysis (MRE) is a novel ISRU method of extracting oxygen gas and metal alloy from lunar regolith. The MRE process involves placing lunar regolith between two electrodes, through which current is passed, to melt the regolith and reduce the metal oxide constituents by direct electrolysis (e.g. FeO, SiO2, MgO, Al2O3) into oxygen gas and metal alloys. The oxygen is liquefied and used as propellant for landers, while the metals (e.g. Ferro-alloys) are further processed and used in structural building materials and parts manufacturing. A system model was developed that accounted for the major processes of an MRE plant (from excavation of raw materials to storage of products) to assess the feasibility of a lunar MRE plant. The System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team utilized its previously documented system sizing model, the Mission Analysis and Integration Tool (MAIT) [1] as framework of the system model. MAIT uses MATLAB/Simulink to integrate subsystem models into a complete system model of the MRE plant. Total mass, volume, and power requirements were computed for numerous iterations of a MRE plant. System Model: Figure 1: MRE Plant Block Diagram The regolith excavation model determines the mass and power needed to excavate sufficient regolith. The preheating auger initiates the regolith heating process before regolith enters the MRE re-actor to reduce the energy required to turn the solid into a molten liquid. The MRE reactor is modeled in COMSOL Multiphysics and based on the research by Dominguez, Sibille, and Schreiner [2, 3, 4]. This preliminary reactor model provides an accurate calculation of thermal equilibrium during electrochemical operation of the reactor system to assess the optimal mass and power required to process the inlet flow of regolith. The model also computes the outlet flowrates of oxygen and molten products. For this analysis, the primary components of the metal alloy considered were iron and silicon. The oxygen is then purified using an Yttrium Stabilized Zirconia (YSZ) electrode, followed by liquefaction using a 90K cryocooler to be stored as liquid oxygen in insulated cylindrical tanks. In future iterations of the system model, the molten metal tapped from the MRE reactor will undergo additional processing or refinement. However, downstream handling of metals is currently a technology gap that is missing a high TRL subsystem model. Therefore, for this analysis, the accumulated metal alloy stream terminates after leaving the MRE reactor. Study Goals: This analysis investigates multiple input variables to the system to determine the sensitivity of a (near) complete plant at full-scale. This preliminary investigation ran parametric sweeps on the MRE reactor geometry, electrical current supply, layers of multi-layer insulation (MLI) on the reactor, size of the electrodes in the oxygen purification model, and regolith composition (based on landing site location). Three production targets of oxygen (1,000, 10,000, and 50,000 kg/yr) were investigated for this analysis. The parametric sweeps conducted in this analysis provide valuable insight into the expected impact of the various model inputs on plant size. This information can be used to identify the most critical components of the plant and guide future decisions on allocating funding for research and development, providing subsystem developers with appropriate interfaces with downstream and upstream processes, and assessing the overall feasibility of MRE when compared to other ISRU plants. References: [1] Carlson, A. et al. (2024) ICES, ICES-2024-53. [2] Dominguez, D.A., and Sibille, L. (2011) AIAA, AIAA-2011-700. [3] Schreiner, S.S. (2015) MIT, Dissertation. [4] Schreiner, S.S. et al. (2016) ASR, 57(7), pp.1585-1603.

ISRU↗

Characterization of Titan Entry Radiative Heating in the Low Density Electric Arc Shock Tube

Prediction of aftbody radiative heating for Dragonfly’s entry at Titan requires accurate models to inform the design of its thermal protection systems. Previously reported investigations of shock layer aftbody radiation in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center focused on the lower altitude/higher density part of the trajectory. This paper presents the first optical emission spectroscopy measurements of an ongoing test series in the new aluminum Low Density Shock Tube (LDST), recently installed at EAST, focusing on the radiative heating characterization through the higher altitude portion of Dragonfly’s entry path. Experimental data collected at a fill pressure of 0.1 Torr and shock velocities between 5.2 and 6.2 km/s show favorable agreement to similar conditions previously tested in the 4" High-Velocity Shock Tube, as well as improvements in the signal to noise ratio and available test time. Simulations employing DPLR and NEQAIR underpredict the experimental radiance, with increasing discrepancy from the NIR towards the UV wavelength range, thus motivating further efforts in improving the current non-Boltzmann model for the CN molecule.

A Fagnani↗

National Modeling of Geothermal District Energy Systems with Ambient-Temperature Loops Using dGeo: Preprint

Geothermal district energy systems (DES) with ambient-temperature loops, also known as thermal energy networks, are one option for decarbonizing space heating and cooling loads. Geothermal fifth-generation DES include an "ambient" temperature thermal loop that connects heat pumps at each building with thermal balancing sources such as geothermal borehole fields. Heating and cooling are provided via a water-source heat pump at each end-user. This project seeks to analyze the nationwide potential for ambient-temperature loop districts by creating a new module within the Distributed Geothermal Market Demand Model (dGeo). dGeo is an agent-based modeling tool for distributed geothermal resources; it can investigate potential on a nationwide or statewide scale using geospatial data for all 50 states and thermal demands for existing buildings. This process allows for high-level estimates of technical and economic potential for ambient-temperature loop districts across the United States. Using GHEDesigner, a lookup table was created to size borehole fields for different thermal loads and ground conditions experienced across the country. A cost and financing structure, along with incentives, were applied. Cost estimates include costs for the distribution network, borehole field installation and operation, and circulation pump operation, while savings are calculated based on agent energy bills. This newly developed module can be used for assessing which areas of the country have the highest potential for agent benefits from ambient-temperature loop installation and assess the impact of different costing and pricing future scenarios. While the code is still under development and nationwide simulations are ongoing, initial results for two states are provided. Future work includes expanding the module to consider mixed residential and commercial districts and considering multiple costing scenarios.

ambient temperature loop↗

Scalability analysis of heavy-duty gas turbines using data-driven machine learning

With the increasing integration of variable renewable energy sources into power systems, the role of flexible power generation technologies like gas turbines (GT) in rapid grid balancing remains crucial. This sustained importance underscores the need for scaled and precise modeling of GT to ensure effective integration within evolving energy frameworks. While physics-driven GT models integrate thermodynamics, fluid dynamics, and combustion principles, they often rely on approximate mathematical representations to accommodate scaling that may not capture the actual complex dynamics for GTs and inertial effects associated to GTs with different ratings. In this study, a data-driven model is proposed using machine learning (ML) techniques to conduct GT scalability analysis and performance evaluation with high accuracy. The ML model, trained on data from various operating conditions and performance parameters, aims to uncover intricate relationships and patterns, resembling GT characteristics at different scales (ratings). The model is developed to capture complex system interaction and to adapt to changing operational scenarios at different capacities, providing valuable insights of power system dynamics. In this study, the real-time digital simulator platform was employed to generate training data for the ML model and assess its dynamic characteristics. The ultimate objective was to develop a detailed modeling framework based on governing equations and data-driven ML capable of predicting key performance indicators, in thermal systems such as GTs, including power output, speed, fuel consumption, and exhaust temperature under diverse operating conditions at different scales. The developed ML framework demonstrated high accuracy, with mean relative errors for GT power prediction, reference speed, exhaust temperature, and compressor pressure ratio (CPR) parameters consistently below 0.1% across typical load fluctuation scenarios. Maximum deviations were limited to approximately 0.5 K for exhaust temperature and 0.009 for CPR, underscoring the model’s ability to replicating dynamic GT behavior with high precision. The adaptability of the ML model enables its application across diverse operational conditions and its extension to other thermal systems. By leveraging advanced ML techniques, this study presents a robust and scalable modeling framework that enhances GT simulation precision, facilitating improved integration into evolving power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Performance Capability of Single-Cavity Vortex Gaseous Nuclear Rockets

An analysis was made to determine the maximum powerplant thrust-to-weight ratio possible with a single-cavity vortex gaseous reactor in which all the hydrogen propellant must diffuse through a fuel-rich region. An assumed radial temperature profile was used to represent conduction, convection, and radiation heat-transfer effects. The effect of hydrogen property changes due to dissociation and ionization was taken into account in a hydrodynamic computer program. It is shown that, even for extremely optimistic assumptions of reactor criticality and operating conditions, such a system is limited to reactor thrust-to-weight ratios of about 1.2 x 10(exp -3) for laminar flow. For turbulent flow, the maximum thrust-to-weight ratio is less than 10(exp -3). These low thrusts result from the fact that the hydrogen flow rate is limited by the diffusion process. The performance of a gas-core system with a specific impulse of 3000 seconds and a powerplant thrust-to-weight ratio of 10(exp -2) is shown to be equivalent to that of a 1000-second advanced solid-core system. It is therefore concluded that a single-cavity vortex gaseous reactor in which all the hydrogen must diffuse through the nuclear fuel is a low-thrust device and offers no improvement over a solid-core nuclear-rocket engine. To achieve higher thrust, additional hydrogen flow must be introduced in such a manner that it will by-pass the nuclear fuel. Obviously, such flow must be heated by thermal radiation. An illustrative model of a single-cavity vortex system employing supplementary flow of hydrogen through the core region is briefly examined. Such a system appears capable of thrust-to-weight ratios of approximately 1 to 10. For a high-impulse engine, this capability would be a considerable improvement over solid-core performance. Limits imposed by thermal radiation heat transfer to cavity walls are acknowledged but not evaluated. Alternate vortex concepts that employ many parallel vortices to achieve higher hydrogen flow rates offer the possibility of sufficiently high thrust-to-weight ratios, if they are not limited by short thermal-radiation path lengths.

Ragsdale, Robert G.↗

World Record Demonstration of > 30% Thermophotovoltaic Conversion Efficiency

Thermophotovoltaic (TPV) devices are solid-state heat engines that convert thermal radiation into electricity using semiconductor diodes, and have applications in energy storage, primary power conversion, and waste heat recovery. Higher temperature emitters and wider band gap photovoltaics (PVs) generally afford higher TPV efficiencies. High temperature emitters are accessible in thermal energy storage applications. We have demonstrated a world record 31% +/- 2% TPV conversion efficiency with a 0.9 cm2 GaAs-based PV device under a 2430 °C thermal emitter, producing an electrical output power of 2.23 W. Critical to the result was the cell’s high reflectance of photon energies below the device band gap. Unlike solar PV, for which sub-band gap (SBG) light is lost, a TPV cell can reflect and recycle SBG light to the thermal emitter. The demonstration was made on a custom-built measurement platform in which a ~100 cm2 graphite thermal emitter was heated under vacuum up to 2430 °C. The TPV efficiency was evaluated by comparing the measured electrical output power of the TPV with the calorimetrically measured total net power incident on the cell. The measured TPV efficiency as a function of thermal emitter temperature was corroborated by our full system modeling predictions. As far as the authors are aware, this is the highest TPV conversion efficiency ever measured, and device improvements should yield > 40% efficiency in the near future.

27 ARPA - Advanced Research Projects Agency-Energy↗

Hypersonic three-dimensional nonequilibrium boundary-layer equations in generalized curvilinear coordinates

The basic governing equations for the second-order three-dimensional hypersonic thermal and chemical nonequilibrium boundary layer are derived by means of an order-of-magnitude analysis. A two-temperature concept is implemented into the system of boundary-layer equations by simplifying the rather complicated general three-temperature thermal gas model. The equations are written in a surface-oriented non-orthogonal curvilinear coordinate system, where two curvilinear coordinates are non-orthogonial and a third coordinate is normal to the surface. The equations are described with minimum use of tensor expressions arising from the coordinate transformation, to avoid unnecessary confusion for readers. The set of equations obtained will be suitable for the development of a three-dimensional nonequilibrium boundary-layer code. Such a code could be used to determine economically the aerodynamic/aerothermodynamic loads to the surfaces of hypersonic vehicles with general configurations. In addition, the basic equations for three-dimensional stagnation flow, of which solution is required as an initial value for space-marching integration of the boundary-layer equations, are given along with the boundary conditions, the boundary-layer parameters, and the inner-outer layer matching procedure. Expressions for the chemical reaction rates and the thermodynamic and transport properties in the thermal nonequilibrium environment are explicitly given.

Lee, Jong-Hun↗

Near Earth Asteroid Scout Solar Sail Thrust and Torque Model

The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid to help prepare for human missions to Near Earth Asteroids. NEA Scout will launch as a secondary payload on the first SLS-Orion mission. NEA Scout will perform a small trim maneuver shortly after deploy from the spent SLS upper stage using a cold gas propulsion system, but from that point on will depend entirely on the solar sail for thrust. As such, it is important to accurately characterize the thrust of the sail in order to achieve mission success. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust; a flat plate model could potentially model thrust well enough to close mission design studies, but a three-dimensional solar sail is essential to control system design. The three-dimensional solar sail model revealed that thermal deformations of unshielded booms would create unacceptably large solar disturbance torques. The original large FEM model was used in control and mission simulations, but was resulted in simulations with prohibitive run times. This led us to adapt the Generalized Sail Model (GSM) of Rios-Reyes. A design reference sail model has been baselined for NEA Scout and has been used to design the mission and control system for the sailcraft. Additionally, since NEA Scout uses reaction wheels for attitude pointing and control, the solar torque model is essentially to successfully design the NEA Scout momentum management control system. We have also updated the estimate of diffusivity used for the aluminized sail material based on optical testing of wrinkled sail material. The model presented here represents the current state of the art of NASA's ability to model solar sail thrust and torque.

Heaton, Andy↗

Model predictive control for integrated control of air-conditioning and mechanical ventilation, lighting and shading systems

Modern buildings are increasingly automated and often equipped with multiple building services (e.g., air-conditioning and mechanical ventilation (ACMV), dynamic shading, dimmable lighting). These systems are conventionally controlled individually without considering their interactions, affecting the building’s overall energy inefficiency and occupant comfort. A model predictive control (MPC) system that features a multi-objective MPC scheme to enable coordinated control of multiple building services for overall optimized energy efficiency, indoor thermal and visual comfort, as well as a hybrid model for predicting indoor visual comfort and lighting power is proposed. The MPC system was implemented in a test facility having two identical, side-by-side experimental cells to facilitate comparison with a building management system (BMS) employing conventional reactive feedback control. The MPC system coordinated the control of the ACMV, dynamic façade and automated dimmable lighting systems in one cell while the BMS controlled the building services in the other cell in a conventional manner. The MPC side achieved 15.1–20.7% electricity consumption reduction, as compared to the BMS side. Simultaneously, the MPC system improved indoor thermal comfort by maintaining the room within the comfortable range (-0.5 < predicted mean vote < 0.5) for 98.3% of the time, up from 91.8% of the time on the BMS side. Visual comfort, measured by indoor daylight glare probability (DGP) and horizontal illuminance level at work plane height, was maintained for the entire test period on the MPC side, improving from having visual comfort for 94.5% and 85.7% of the time, respectively, on the BMS side.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling the TPF interferometer

The Terrestrial Planet Finder interferometer design concepts are large and complex systems that must operate in environments that are impractical to reproduce in preflight testing. The structurally-connected design is 36 meters long - longer than all but one thermal vacuum chamber in existence. The formation flying design will be comprised of up to five separate spacecraft, each with a sunshield over 15 meters on a side, and is designed to operate with formation sizes spanning over 100 meters to very close formations. System-level verification of the performance of the designs will need to rely on analytical modeling. The effort to model the many physical aspects of the designs under study is under way*. This paper describes the program of modeling for the TPF-I concepts. The program includes a number of types of models, such as the standard stand-alone optics, thermal, and structural models, as well as an end-to-end performance model of the project system called the Observatory Simulation. Aspects of each model are discussed including the purpose, methods of implementation (software applications), and approaches to validation. Program-level considerations (such as model-to-model integration and configuration management) are also discussed. Given that there are at least seven different organizations contributing to model developments and more than twenty separate models, these are special challenges.

Henry, Curt↗

Suppressing thermal transport in nonporous polymer hybrids by limiting thermally accessible vibrational modes

Achieving low thermal conductivity in nonporous polymer materials without compromising mechanical integrity remains a longstanding challenge. Conventional strategies, such as introducing porosity, are inherently limited in dense systems. Here, in this study, we demonstrate that thermal transport in polymer–organic filler hybrids is closely associated with the engineered availability of thermally accessible vibrational modes, rather than individual interfacial resistance or porosity. Using PU/organic filler hybrids as a model system, we show that incorporating rigid organic motifs shifts the vibrational density of states toward higher frequencies, where vibrational modes are weakly thermally populated at ambient conditions. This limits the number of heat-carrying channels, leading to suppressed thermal conductivity despite negligible individual interfacial resistance (∼10 −9 m 2 K W −1 ). These findings provide initial evidence that vibrational mode engineering may offer a promising molecular-level route for reducing thermal conductivity in dense polymer–organic hybrid systems.

36 MATERIALS SCIENCE↗

Thermal stress effects in intermetallic matrix composites

Intermetallic matrix composites develop residual stresses from the large thermal expansion mismatch (delta-alpha) between the fibers and matrix. This work was undertaken to: establish improved techniques to measure these thermal stresses in IMC's; determine residual stresses in a variety of IMC systems by experiments and modeling; and, determine the effect of residual stresses on selected mechanical properties of an IMC. X ray diffraction (XRD), neutron diffraction (ND), synchrotron XRD (SXRD), and ultrasonics (US) techniques for measuring thermal stresses in IMC were examined and ND was selected as the most promising technique. ND was demonstrated on a variety of IMC systems encompassing Ti- and Ni-base matrices, SiC, W, and Al2O3 fibers, and different fiber fractions (Vf). Experimental results on these systems agreed with predictions of a concentric cylinder model. In SiC/Ti-base systems, little yielding was found and stresses were controlled primarily by delta-alpha and Vf. In Ni-base matrix systems, yield strength of the matrix and Vf controlled stress levels. The longitudinal residual stresses in SCS-6/Ti-24Al-llNb composite were modified by thermomechanical processing. Increasing residual stress decreased ultimate tensile strength in agreement with model predictions. Fiber pushout strength showed an unexpected inverse correlation with residual stress. In-plane shear yield strength showed no dependence on residual stress. Higher levels of residual tension led to higher fatigue crack growth rates, as suggested by matrix mean stress effects.

Wright, P. K.↗

Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. The quadrupole method has been shown as a viable technique for rapidedly simulating the thermal repsonse of layered systems. Often, the effort has focused on a one-dimentional models, in particular for improved analysis of thermal data. For composites, the in-plane heat diffusion often significantly impacts the thermal response of defects of interest, therefore three-dimensional simulations are desirable. This paper discusses the extension of the quadrupole methodology to perform simulations of thermographic responses in three-dimensional configuration. This enables the more realistic simulations of the thermal response of delaminiations in composites. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are known obtained from other sources.

Thermography↗

Simulation of Thermographic Responses of Delaminations in Composites with Quadrupole Method

Thermography has been shown to be a viable technique for inspection of composites. The quadrupole method is a valuable technique for rapidly simulating the thermal response of layered systems. Often, the effort has focused on a one-dimensional model, in particular for improved analysis of thermal data. For composites, the in-plane heat diffusion often significantly impacts the thermal response of defects of interest, therefore three-dimensional simulations are desirable. This paper discusses the extension of the quadrupole methodology to perform simulations of thermographic responses in three-dimensional configurations. This enables the more realistic simulations of the thermal response of delaminations in composites. The simulations are compared to finite element simulations of the same inspection configurations.

Thermography↗

A simulation model for wind energy storage systems. Volume 1: Technical report

A comprehensive computer program for the modeling of wind energy and storage systems utilizing any combination of five types of storage (pumped hydro, battery, thermal, flywheel and pneumatic) was developed. The level of detail of Simulation Model for Wind Energy Storage (SIMWEST) is consistent with a role of evaluating the economic feasibility as well as the general performance of wind energy systems. The software package consists of two basic programs and a library of system, environmental, and load components. The first program is a precompiler which generates computer models (in FORTRAN) of complex wind source storage application systems, from user specifications using the respective library components. The second program provides the techno-economic system analysis with the respective I/O, the integration of systems dynamics, and the iteration for conveyance of variables. SIMWEST program, as described, runs on the UNIVAC 1100 series computers.

Warren, A. W.↗