SEARCH · Search NASA
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.
A user's guide is presented for the transient version of the thermoregulatory model. The model is designed to simulate the transient response of the human thermoregulatory system to thermal inputs. The model consists of 41 compartments over which the terms of the heat balance are computed. The control mechanisms which are identified are sweating, vaso-constriction and vasodilation.
A large set of experiments performed with the NASA Ames Mars General Circulation Model is analyzed to determine the properties, structure, and dynamics of the simulated transient baroclinic eddies. There is strong transient baroclinic eddy activity in the extratropics of the Northern Hemisphere during the northern autumn, winter, and spring seasons. The eddy activity remains strong for very large dust loadings, though it shifts northward. The eastward propagating eddies are characterized by zonal wavenumbers of 1-4 and periods of about 2-10 days. The properties of the GCM baroclinic eddies in the northern extratropics are compared in detail with analogous properties inferred from Viking Lander meteorology observations.
Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.
This paper describes the modeling approach used to simulate the transient behavior of a salinity gradient solar pond. A system of finite difference equations are used to generate the time dependent temperature and salinity profiles within the pond. The stability of the pond, as determined by the capacity of the resulting salinity profile to suppress thermal convection within the primary gradient region of the pond, is continually monitored and when necessary adjustments are made to the thickness of the gradient zone. Results of the model are then compared to measurements taken during two representative seasonal periods at the University of Texas at El Paso's (UTEP's) research solar pond.
A spherically symmetric adiabatic single-fluid model is outlined for simulating the nonlinear time-dependent response of the corona to solar events that are simulated by perturbations in the appropriate physical variables from their steady-state values at the coronal base. Several observed features of the coronal transient that occurred on June 10, 1973, are simulated by using a particular steady-state solar wind, a specific combination of density and temperature perturbations, and a particular time dependence of the perturbations. A different steady-state solar wind, a perturbation of shorter duration, and other perturbation combinations are also employed to simulate the same transient so that the effect of each quantity can be determined. It is found that the model cannot adequately simulate all the observational results for the investigated transient, that the steady-state solar wind is relatively unimportant in such numerical simulations, and that studies which attempt only to reproduce observed shock trajectories may lead to erroneous conclusions regarding the physics of the solar event that produced the transient.
The equations of motion necessary to calculate the transient response of a multimass flexible rotor supported by nonlinear, damped bearings are derived from energy principles. Rotor excitation may be the result of imbalance, internal friction, rotor acceleration, nonlinear forces due to any number of bearing or seal stations, and gyroscopic couples developed from skewed disk effects. The method of solution for transient response simulation is discussed in detail and is based on extensive evaluation of numerical methods available for transient analysis. Examples of the application of transient response for the analysis of rotor bearing systems are presented and compared to actual machine performance. Recommendations for the use and extension of the present system simulation model are discussed.
The direct Monte Carlo method is used to simulate the transient motion of a confined rarefied gas as a container wall is rapidly heated or cooled. Results show that, for the rapid wall temperature changes, the gas response is better measured by an acoustic, rather than conduction, timescale, with steady state conditions being achieved after about ten wave crossing times.
Numerical simulations of the formation and propagation of mass ejection, loop transients in coronal streamers are discussed. The simulations are accomplished with numerical solutions of the single fluid, ideal MHD equations of motion in the meridional plane. The streamer is produced by simulating the relaxation of an initially radial hydrodynamic flow coupled with a dipole magnetic field. The simulated transient then results from an energy release at the base of the streamer. The legs of the loop transient produced remain essentially stationary while the loop expands mainly in the radial direction with velocities of 400 to 750 km s-1. Once the leading edge of the transient has passed out of the lower corona, the initial streamer configuration is restored after 15 to 24 hours. A second energy release 2 hours later than, and with an energy release identical to, the first does not produce a significant coronal disturbance.
Computerized FORTRAN equations simulation of linear and nonlinear circuits for symbolic analysis and transient response
A computer program has been developed to simulate the transient wall pressure field produced by a low speed fully turbulent boundary layer. The theoretical basis for the simulation has been discussed and preliminary results from a pressure simulation are presented.
The effect of thermal diffusively on the defect detectability in Carbon/Epoxy composite panels by transient thermography is presented in this paper. A series of Finite Element Models were constructed and analyzed to simulate the transient heat transfer phenomenon during Thermographic Non-destructive Evaluation (TNDE) of composite panels with square defects. Six common carbon fibers were considered. The models were built for composites with various combinations of fibers and volumetric ratios. Finite Element Analysis of these models showed the trends of the detectable range and the maximum thermal contrast versus the thermal diffusivity of various composites. Additionally, the trends of defect size to depth ratio and the thermal contrast has been investigated.
Time-dependent motions in the solar corona are subject to modification through conductive transport of thermal energy away from compression regions and into rarefaction regions. The evolution depends on how the energy is introduced - by an essentially thermal perturbation or by a density perturbation. The net effect of thermal conduction is to convert local kinetic energy in the transient into widely distributed thermal energy, which then may in turn be partially converted into kinetic energy of the bulk flow. Local motion in the transient is thereby damped. These effects are modeled here using a numerical solution of the time-dependent solar wind equations for single-fluid thermally conductive flow. The simulations of transient phenomena are made for ordinary collisional thermal conduction and compared with similar simulations done using polytropic flow. Between one and five solar radii, the gross character in the two cases is very similar. However, there are effects with thermal conduction that cannot be modeled with polytropic flow, illustrated here by temperature and velocity forerunners and a delayed and large-amplitude trailing velocity low. If thermal conduction is artificially inhibited by simple reduction of the thermal conduction coefficient, relative velocity amplitudes remain larger in approximate proportion to the amount the coefficient is reduced.
Combustion instability in ramjets is a complex phenomenon that involve nonlinear interaction between acoustic waves, vortex motion and unsteady heat release in the combustor. To numerically simulate this 3-D, transient phenomenon, enormous computer resources (time, memory and disk storage) are required. Although current generation vector supercomputers are capable of providing adequate resources for simulations of this nature, their high cost and limited availability, makes such machines less than satisfactory for routine use. The primary focus of this study is to assess the feasibility of using highly parallel computer systems as a cost-effective alternative for conducting such unsteady flow simulations. Towards this end, a large-eddy simulation model for combustion instability was implemented on the Intel iPSC/860 and a careful study was conducted to determine the benefits and the problems associated with the use of such machines for transient simulations. Details of this study along with the results obtained from the unsteady combustion simulations carried out on the iPSC/860 are discussed in this paper.
Transient heat transfer analysis of Apollo Service Module
The DYNGEN, a digital computer program for analyzing the steady state and transient performance of turbojet and turbofan engines, is described. The DYNGEN is based on earlier computer codes (SMOTE, GENENG, and GENENG 2) which are capable of calculating the steady state performance of turbojet and turbofan engines at design and off-design operating conditions. The DYNGEN has the combined capabilities of GENENG and GENENG 2 for calculating steady state performance; to these the further capability for calculating transient performance was added. The DYNGEN can be used to analyze one- and two-spool turbojet engines or two- and three-spool turbofan engines without modification to the basic program. A modified Euler method is used by DYNGEN to solve the differential equations which model the dynamics of the engine. This new method frees the programmer from having to minimize the number of equations which require iterative solution. As a result, some of the approximations normally used in transient engine simulations can be eliminated. This tends to produce better agreement when answers are compared with those from purely steady state simulations. The modified Euler method also permits the user to specify large time steps (about 0.10 sec) to be used in the solution of the differential equations. This saves computer execution time when long transients are run. Examples of the use of the program are included, and program results are compared with those from an existing hybrid-computer simulation of a two-spool turbofan.
A finite element model of an ATR42-300 commuter-class aircraft was developed and a crash simulation was executed. Analytical predictions were correlated with data obtained from a 30-feet per second (9.14-meters per second) vertical drop test of the aircraft. The purpose of the test was to evaluate the structural response of the aircraft when subjected to a severe, but survivable, impact. The aircraft was configured with seats, dummies, luggage, and other ballast. The wings were filled with 8,700 lb. (3,946 kilograms) of water to represent the fuel. The finite element model, which consisted of 57,643 nodes and 62,979 elements, was developed from direct measurements of the airframe geometry. The seats, dummies, luggage, simulated engines and fuel, and other ballast were represented using concentrated masses. The model was executed in LS-DYNA, a commercial finite element code for performing explicit transient dynamic simulations. Analytical predictions of structural deformation and selected time-history responses were correlated with experimental data from the drop test to validate the simulation.
Recently, heat transfer correlations based on liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching data were developed to improve the predictive accuracy of lumped node codes like SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP). After implementing these correlations into both programs, updated model runs showed strong improvement in LN2 pipe chilldown modeling but only modest improvement in LH2 modeling. Due to large differences in thermal and fluid properties between the two fluids, results indicated a need to develop a separate set of LH2-only correlations to improve the accuracy of the simulations. This paper presents a new set of two-phase convection heat transfer correlations based on LH2 pipe quenching data. A correlation to predict the bulk vapor temperature was developed after analysis showed that high amounts of thermal nonequilibrium of the liquid and vapor phases occurred during film boiling of LH2. Implemented in a numerical model, the new correlations achieve a mean absolute error of 19.5 K in the predicted wall temperature when compared to recent LH2 pipe chilldown data, an improvement of 40% over recent GFSSP predictions. This correlation set can be implemented in simulations of the transient LH2 chilldown process. Such simulations are useful for predicting the chilldown time and boil-off mass of LH2 for applications such as the transfer of LH2 from a ground storage tank to the rocket vehicle propellant tank, or through a rocket engine feedline during engine startup.
Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC/yr for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC/yr). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).