Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Study”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Recycled Paper as a Source of Renewable Jet Fuel in the United States

Converting biomass into jet fuel involves more than the core chemical process. The overall process includes the logistics of harvesting and transporting the biomass, handling and preparing the material for processing, and processing and disposal of waste. All of these activities contribute to cost. Controlling cost involves more than developing efficient process chemistry. Choice of feedstock also has a significant impact on process economics. We consider chemical conversion of paper from municipal solid waste as a feedstock for the production of jet fuel and diesel. Paper has a significantly higher cellulose content than raw lignocellulosic biomass such as corn stover, so it requires less pretreatment to convert it into hydrocarbons than lignocellulosic biomass. Our techno-economic analysis showed that the cost of converting paper waste into jet fuel is about $1.00/gal less than jet fuel produced from corn stover. Although the cost of recycling paper into jet fuel is less than producing it from corn stover, the process is not competitive with petroleum. We estimated a minimum selling price of $3.97/gal for paper-derived jet fuel. Our sensitivity studies indicated that the biggest economic obstacle is the cost of cellulose hydrolysis. Direct hydrogenation of paper to sugar alcohols combined with increased economy of scale could make recycling paper jet fuel competitive.

09 BIOMASS FUELS↗
DOE OSTI2023

Nuclear data for space exploration

Understanding the harmful effects of galactic cosmic rays (GCRs) on space exploration requires a substantial amount of nuclear data. Specifically, the interaction of energetic GCR charged particles with spacecraft materials generates secondary radiations that, through energy deposition, can harm astronauts and electronic systems. By identifying the gaps in our knowledge of the relevant nuclear data—such as interaction cross sections—and identifying ways to fill those gaps—with measurements, compilations, evaluations, dissemination, reaction modeling, sensitivity studies, and uncertainty quantification—the safety and viability of space exploration can be improved. This work surveys the state of the art in this interdisciplinary field and identifies promising collaborative research topics that have significant potential to advance our understanding of the effects of the space radiation environment on space exploration.

42 ENGINEERING↗
DOE OSTI2022

Impacts of the Thermal Gradient on Inland Advecting Sea Breezes in the Southeastern United States

Sea breezes are frequently observed in the South Carolina/Georgia region of the Southeastern United States (SEUS) and can reach upwards of 150km inland. This region is unique among the places frequently affected by sea breeze due to it being a continental location with relatively flat topography. The thermal gradient between land and water environments is a factor in introducing the sea breeze, but its role in the inland extent of sea breeze propagation isn’t as well known. We investigate the role of the thermal gradient in previously catalogued sea breeze events observed at the Savannah River Site (SRS) by taking differences of temperature measurements at inland and coastal weather stations for the days that the events occurred. We saw that the temperature differences for those days were much higher than in the non-sea breeze days during the mornings and afternoon. Numerical models were also used to conduct a sensitivity study on a sea breeze case, using simple modifications of the temperature gradient. We found that while the modifications did not stop the generation of a sea breeze circulation, the extent of the inland propagation was dependent on the magnitude of the thermal gradient.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessing the Impact of Lightning NOx Emissions in CMAQ Using Lightning Flash Data from WWLLN over the Contiguous United States

Comparison of lightning flash data from the National Lightning Detection Network (NLDN) and from the World Wide Lightning Location Network (WWLLN) over the contiguous United States (CONUS) for the 2016–2018 period reveals temporally and spatially varying flash rates that would influence lightning NO x (LNO x ) production due to known detection efficiency differences especially during summer months over land (versus over ocean). However, the lightning flash density differences between the two networks show persistent seasonal patterns over geographical regions (e.g., land versus ocean). Since the NLDN data are considered to have higher accuracy (lightning detection with >95% efficiency), we developed scaling factors for the WWLLN flash data based on the ratios of WWLLN to NLDN flash data over time (months of year) and space. In this study, sensitivity simulations using the Community Multiscale Air Quality (CMAQ) model are performed utilizing the original data sets (both NLDN and WWLLN) and the scaled WWLLN flash data for LNO x production over the CONUS. The model performance of using the different lightning flash datasets for ambient O 3 and NO x mixing ratios that are directly impacted by LNO x emissions and the wet and dry deposition of oxidized nitrogen species that are indirectly impacted by LNO x emissions is assessed based on comparisons with ground-based observations, vertical profile measurements, and satellite products. During summer months, the original WWLLN data produced less LNO x emissions (due to its lower lightning detection efficiency) compared to the NLDN data, which resulted in less improvement in model performance than the simulation using NLDN data as compared to the simulation without any LNO x emissions. However, the scaled WWLLN data produced LNO x estimates and model performance comparable with the NLDN data, suggesting that scaled WWLLN may be used as a substitute for the NLDN data to provide LNO x estimates in air quality models when the NLDN data are not available (e.g., due to prohibitive cost or lack of spatial coverage).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Techno-Economic Analysis and Optimization of a Compressed-Air Energy Storage System Integrated with a Natural Gas Combined-Cycle Plant

To address the rising electricity demand and greenhouse gas concentration in the environment, considerable effort is being carried out across the globe on installing and operating renewable energy sources. However, the renewable energy production is affected by diurnal and seasonal variability. To ensure that the electric grid remains reliable and resilient even for the high penetration of renewables into the grid, various types of energy storage systems are being investigated. In this paper, a compressed-air energy storage (CAES) system integrated with a natural gas combined-cycle (NGCC) power plant is investigated where air is extracted from the gas turbine compressor or injected back into the gas turbine combustor when it is optimal to do so. First-principles dynamic models of the NGCC plant and CAES are developed along with the development of an economic model. The dynamic optimization of the integrated system is undertaken in the Python/Pyomo platform for maximizing the net present value (NPV). NPV optimization is undertaken for 14 regions/cases considering year-long locational marginal price (LMP) data with a 1 h interval. Design variables such as the storage capacity and storage pressure, as well as the operating variables such as the power plant load, air injection rate, and air extraction rate, are optimized. Results show that the integrated CAES system has a higher NPV than the NGCC-only system for all 14 regions, thus indicating the potential deployment of the integrated system under the assumption of the availability of caverns in close proximity to the NGCC plant. The levelized cost of storage is found to be in the range of 136–145 $/MWh. Roundtrip efficiency is found to be between 74.6–82.5%. A sensitivity study with respect to LMP shows that the LMP profile has a significant impact on the extent of air injection/extraction while capital expenditure reduction has a negligible effect.

25 ENERGY STORAGE↗
DOE OSTI2022

Explanation and Application of the Evolving Contact Traction Fields in Shallow Foundation Systems

The present paper provides a qualitative discussion of the evolution of contact traction fields beneath rigid shallow foundations resting on granular materials. A phenomenological similarity is recognized in the measured contact traction fields of rigid footings and at the bases of sandpiles. This observation leads to the hypothesis that the stress distributions are brought about by the same physical phenomena, namely the development of arching effects through force chains and mobilized intergranular friction. A set of semi-empirical equations are suggested for the normal and tangential components of this contact traction based on past experimental measurements and phenomenological assumptions of frictional behaviors at the foundation system scale. These equations are then applied to the prescribed boundary conditions for the analysis of the settlement, resistance, and stress fields in supporting granular materials beneath the footing. A parametric sensitivity study is performed on the proposed modelling method, highlighting solutions to the boundary-value problems in an isotropic, homogeneous elastic half-space.

42 ENGINEERING↗
DOE OSTI2020

Multi-Size Proppant Pumping Schedule of Hydraulic Fracturing: Application to a MP-PIC Model of Unconventional Reservoir for Enhanced Gas Production

Slickwater hydraulic fracturing is becoming a prevalent approach to economically recovering shale hydrocarbon. It is very important to understand the proppant’s transport behavior during slickwater hydraulic fracturing treatment for effective creation of a desired propped fracture geometry. The currently available models are either oversimplified or have been performed at limited length scales to avoid high computational requirements. Another limitation is that the currently available hydraulic fracturing simulators are developed using only single-sized proppant particles. Motivated by this, in this work, a computationally efficient, three-dimensional, multiphase particle-in-cell (MP-PIC) model was employed to simulate the multi-size proppant transport in a field-scale geometry using the Eulerian–Lagrangian framework. Instead of tracking each particle, groups of particles (called parcels) are tracked, which allows one to simulate the proppant transport in field-scale geometries at an affordable computational cost. Then, we found from our sensitivity study that pumping schedules significantly affect propped fracture surface area and average fracture conductivity, thereby influencing shale gas production. Motivated by these results, we propose an optimization framework using the MP-PIC model to design the multi-size proppant pumping schedule that maximizes shale gas production from unconventional reservoirs for given fracturing resources.

42 ENGINEERING↗
DOE OSTI2021

Sensitivity of Neutron-Rich Nuclear Isomer Behavior to Uncertainties in Direct Transitions

Nuclear isomers are populated in the rapid neutron capture process (r process) of nucleosynthesis. The r process may cover a wide range of temperatures, potentially starting from several tens of GK (several MeV) and then cooling as material is ejected from the event. As the r-process environment cools, isomers can freeze out of thermal equilibrium or be directly populated as astrophysically metastable isomers (astromers). Astromers can undergo reactions and decays at rates very different from the ground state, so they may need to be treated independently in nucleosythesis simulations. Two key behaviors of astromers—ground state ↔ isomer transition rates and thermalization temperatures—are determined by direct transition rates between pairs of nuclear states. We perform a sensitivity study to constrain the effects of unknown transitions on astromer behavior. Detailed balance ensures that ground → isomer and isomer → ground transitions are symmetric, so unknown transitions are equally impactful in both directions. We also introduce a categorization of astromers that describes their potential effects in hot environments. We provide a table of neutron-rich isomers that includes the astromer type, thermalization temperature, and key unmeasured transition rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Jupiter in the Ultraviolet: Acetylene and Ethane Abundances in the Stratosphere of Jupiter from Cassini Observations between 0.15 and 0.19 μm

At wavelengths between 0.15 and 0.19 μm, the far-ultraviolet spectrum of Jupiter is dominated by the scattered solar spectrum, attenuated by molecular absorptions primarily by acetylene and ethane, and to a lesser extent ammonia and phosphine. We describe the development of our radiative transfer code that enables the retrieval of abundances of these molecular species from ultraviolet reflectance spectra. As a proof-of-concept we present an analysis of Cassini Ultraviolet Imaging Spectrograph (UVIS) observations of the disk of Jupiter during the 2000/2001 flyby. The ultraviolet-retrieved acetylene abundances in the upper stratosphere are lower than those predicted by models based solely on infrared thermal emission from the mid-stratosphere observed by the Composite Infrared Spectrometer (CIRS), requiring an adjustment to the vertical profiles above 1 mbar. We produce a vertical acetylene abundance profile that is compatible with both CIRS and UVIS, with reduced abundances at pressures <1 mbar: the 0.1 mbar abundances are 1.21 ± 0.07 ppm for acetylene and 20.8 ± 5.1 ppm for ethane. Finally, we perform a sensitivity study for the JUICE ultraviolet spectrograph, which has extended wavelength coverage out to 0.21 μm, enabling the retrieval of ammonia and phosphine abundances, in addition to acetylene and ethane.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Abstract The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Astromers: Nuclear Isomers in Astrophysics

In this work, we develop a method to compute thermally mediated transition rates between the ground state and long-lived isomers in nuclei. We also establish criteria delimiting a thermalization temperature above which a nucleus may be considered a single species and below which it must be treated as two separate species: a ground-state species and an astrophysical isomer ("astromer") species. Below the thermalization temperature, the destruction rates dominate the internal transition rates between the ground state and the isomer. If the destruction rates also differ greatly from one another, the nuclear levels fall out of or fail to reach thermal equilibrium. Without thermal equilibrium, there may not be a safe assumption about the distribution of occupation probability among the nuclear levels when computing nuclear reaction rates. In these conditions, the isomer has astrophysical consequences and should be treated as a separate astromer species which evolves separately from the ground state in a nucleosynthesis network. We apply our transition-rate methods and perform sensitivity studies on a few well-known astromers. We also study transitions in several other isomers of likely astrophysical interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evaluating cloud liquid detection against Cloudnet using cloud radar Doppler spectra in a pre-trained artificial neural network

Detection of liquid-containing cloud layers in thick mixed-phase clouds or multi-layer cloud situations from ground-based remote-sensing instruments still poses observational challenges, yet improvements are crucial since the existence of multi-layer liquid layers in mixed-phase cloud situations influences cloud radiative effects, cloud lifetime, and precipitation formation processes. Hydrometeor target classifications such as from Cloudnet that require a lidar signal for the classification of liquid are limited to the maximum height of lidar signal penetration and thus often lead to underestimations of liquid-containing cloud layers. Here we evaluate the Cloudnet liquid detection against the approach of Luke et al. (2010) which extracts morphological features in cloud-penetrating cloud radar Doppler spectra measurements in an artificial neural network (ANN) approach to classify liquid beyond full lidar signal attenuation based on the simulation of the two lidar parameters particle backscatter coefficient and particle depolarization ratio. We show that the ANN of Luke et al. (2010) which was trained under Arctic conditions can successfully be applied to observations at the mid-latitudes obtained during the 7-week-long ACCEPT field experiment in Cabauw, the Netherlands, in 2014. In a sensitivity study covering the whole duration of the ACCEPT campaign, different liquid-detection thresholds for ANN-predicted lidar variables are applied and evaluated against the Cloudnet target classification. Independent validation of the liquid mask from the standard Cloudnet target classification against the ANN-based technique is realized by comparisons to observations of microwave radiometer liquid-water path, ceilometer liquid-layer base altitude, and radiosonde relative humidity. In addition, a case-study comparison against the cloud feature mask detected by the space-borne lidar aboard the CALIPSO satellite is presented. Three conclusions were drawn from the investigation. First, it was found that the threshold selection criteria of liquid-related lidar backscatter and depolarization alone control the liquid detection considerably. Second, all threshold values used in the ANN framework were found to outperform the Cloudnet target classification for deep or multi-layer cloud situations where the lidar signal is fully attenuated within low liquid layers and the cloud radar is able to detect the microphysical fingerprint of liquid in higher cloud layers. Third, if lidar data are available, Cloudnet is at least as good as the ANN. The times when Cloudnet outperforms the ANN in liquid detections are often associated with situations where cloud dynamics smear the imprint of cloud microphysics on the radar Doppler spectra.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗
DOE OSTI2022

The sensitivity of the Fitch wind farm parameterization to a three-dimensional planetary boundary layer scheme

Abstract. Wind plant wake impacts can be estimated with a number of simulation methodologies, each with its own fidelity and sensitivity to model inputs. In turbine-free mesoscale simulations, hub-height wind speeds often significantly vary with the choice of a planetary boundary layer (PBL) scheme. However, the sensitivity of wind plant wakes to a PBL scheme has not been explored because, as of the Weather Research and Forecasting model v4.3.3, wake parameterizations were only compatible with one PBL scheme. We couple the Fitch wind farm parameterization with the new NCAR 3DPBL scheme and compare the resulting wakes to those simulated with a widely used PBL scheme. We simulate a wind plant in pseudo-steady states under idealized stable, neutral, and unstable conditions with matching hub-height wind speeds using two PBL schemes: MYNN and the NCAR 3DPBL. For these idealized scenarios, average hub-height wind speed losses within the plant differ between PBL schemes by between −0.20 and 0.22 m s−1, and correspondingly, capacity factors range between 39.5 %–53.8 %. These simulations suggest that PBL schemes represent a meaningful source of modeled wind resource uncertainty; therefore, we recommend incorporating PBL variability into future wind plant planning sensitivity studies as well as wind forecasting studies.

17 WIND ENERGY↗
DOE OSTI2023

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗
DOE OSTI2024

Modeling Dense Particle Flow in Multistage and Obstructed Flow Receivers Using High Fidelity Simulations

Particles are a leading contender for next-generation, concentrating solar power technologies, and the design of the particle receiver is critical to minimize the levelized cost of electricity. Falling particle receivers (FPRs) are a viable receiver concept, but many new designs feature complex particle obstructions that include dense discrete phase flows. This creates additional challenges for modeling as particle-to-particle interactions (i.e., collisions) and particle drag become more complex. To improve upon existing modeling strategies, a CFD-DEM simulation capability was created by coupling two independent codes: Sierra/Fuego and LAMMPS. A suitable receiver model was then defined using a traditional continuum-based model for the air and a granular model for the particle curtain. A sensitivity study was executed using this model to determine the relevance of different granular model inputs on important quantities of interest in obstructed flow FPRs: the particle velocity and curtain opacity. The study showed that the granular model inputs had little effect on the particle velocity magnitude and curtain opacity after an obstruction.

Mills, Brantley↗
NASA NTRS1971

Computational techniques for design optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

Computational techniques were developed and assimilated for the design optimization. The resulting computer program was then used to perform initial optimization and sensitivity studies on a typical thermal protection system (TPS) to demonstrate its application to the space shuttle TPS design. The program was developed in Fortran IV for the CDC 6400 but was subsequently converted to the Fortran V language to be used on the Univac 1108. The program allows for improvement and update of the performance prediction techniques. The program logic involves subroutines which handle the following basic functions: (1) a driver which calls for input, output, and communication between program and user and between the subroutines themselves; (2) thermodynamic analysis; (3) thermal stress analysis; (4) acoustic fatigue analysis; and (5) weights/cost analysis. In addition, a system total cost is predicted based on system weight and historical cost data of similar systems. Two basic types of input are provided, both of which are based on trajectory data. These are vehicle attitude (altitude, velocity, and angles of attack and sideslip), for external heat and pressure loads calculation, and heating rates and pressure loads as a function of time.

Source record↗