SEARCH · Engineering Papers
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.
Heat capacity expression due to chemical reaction, relating composition variations and binary diffusion coefficient with Lewis number, presenting integral equation for system enthalpy
A numerical simulation was developed to investigate the one dimensional heat transfer occurring in a system composed of a layered aircraft blade having an ice deposit on its surface. The finite difference representation of the heat conduction equations was done using the Crank-Nicolson implicit finite difference formulation. The simulation considers uniform or time dependent heat sources, from heaters which can be either point sources or of finite thickness. For the ice water phase change, a numerical method which approximates the latent heat effect by a large heat capacity over a small temperature interval was applied. The simulation describes the temperature profiles within the various layers of the de-icer pad, as well as the movement of the ice water interface. The simulation could also be used to predict the one dimensional temperature profiles in any composite slab having different boundary conditions.
Two computational methods are presented that enable a characterization of corrosion damage to be performed from thermal response measurements derived from a standard flash thermographic inspection. The first is based upon a one dimensional analytical solution to the heat diffusion equation and presumes the lateral extent of damage is large compared to the residual structural thickness, such that lateral heat diffusion effects can be considered insignificant. The second proposed method, based on a finite element optimization scheme, addresses the more general case where these conditions are not met. Results from an experimental application are given to illustrate the precision, robustness and practical efficacy of both methods.
This paper describes a finite volume procedure for network flow analysis in a thermofluid system. A flow network is defined as a group of inter-connected control volumes called “nodes” that are connected by “branches”. The mass and energy conservation equations are solved at the nodes and momentum conservation equations are solved at the branches. The flow network also includes solid nodes to account for fluid to solid heat transfer. The heat conduction equation is solved at the solid nodes in conjunction with the flow equations. The properties of a real fluid are calculated using a thermodynamic property program and used in the conservation equations. The system of equations describing the fluid-solid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution method. This procedure has been incorporated into a general-purpose computer program, the Generalized Fluid System Simulation Program, GFSSP. This paper also presents the application and verification of the method by comparison with test data for several applications that include (a) internal flow in a rocket engine turbo-pump, (b) pressurization and loading of a cryogenic propellant tank, (c) fluid transient during sudden opening of valve for priming of an evacuated feed line, and (d) chill down of a cryogenic transfer line with phase change and two phase flows.
A thermal-stress analysis of a wind turbine blade made of wood composite material is reported. First, the governing partial differential equation on heat conduction is derived, then, a finite element procedure using variational approach is developed for the solution of the governing equation. Thus, the temperature distribution throughout the blade is determined. Next, based on the temperature distribution, a finite element procedure using potential energy approach is applied to determine the thermal-stress distribution. A set of results is obtained through the use of a computer, which is considered to be satisfactory. All computer programs are contained in the report.
This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.
This study is to investigate the thermal stress of a wind turbine blade made of wood composite material. First, the governing partial differential equation on heat conduction is stated, then, a finite element procedure using a variational approach is employed for the solution of the governing equation. Thus, the temperature field throughout the blade is determined. Next, based on the temperature field, a finite element procedure using potential energy approach is applied to determine the thermal stress field. A set of results is obtained through the use of a computer, which is considered to be satisfactory.
A research effort was underway to study the use of two equation low Reynolds number turbulence models in predicting gas side heat transfer on turbine blades. The major objectives of this work are basicly threefold: study the predictive capabilities of two equation low Reynolds number turbulence models under the conditions characteristic of modern gas turbine blades; explore potential improvements to the models themselves as well as to the specification of initial conditions; and provide a comparison of the predictions of these models with the experimental data from a broad range of recently available turbine cascade experiments. The problems associated with predicting the boundary layer transition from laminar to turbulent flow are emphasized, as this may be the most serious deficiency of current modeling techniques. The results and conclusions of the first two phases are briefly described.
Equations derived to improve analysis of experimental data in flash method for measurement of thermal diffusivity and heat capacity. In flash method, pulse of radiant energy from flashlamp or other source deposited on front face of small, thin specimen, and temperature of rear face monitored as function of time. New equations account for both losses and exponentially decaying pulses.
Heat flux equations for magnetospheric electron temperature distribution
Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.
A solution is presented for the differential equation relating the combined effects of radiative and forced convective cooling for a heated sphere. The equation has the form where T and t are the variables temperature and time, respectively, and K sub o, T sub o, and H are constants. The solution can be used as a guideline for the design and understanding of space processing phenomena.
Two versions of the two-equation k-omega model and a shear stress transport (SST) model are used in a three-dimensional, multi-block, Navier-Stokes code to compare the detailed heat transfer measurements on a transonic turbine blade. It is found that the SST model resolves the passage vortex better on the suction side of the blade, thus yielding a better comparison with the experimental data than either of the k-w models. However, the comparison is still deficient on the suction side of the blade. Use of the SST model does require the computation of distance from a wall, which for a multiblock grid, such as in the present case, can be complicated. However, a relatively easy fix for this problem was devised. Also addressed are issues such as (1) computation of the production term in the turbulence equations for aerodynamic applications, and (2) the relation between the computational and experimental values for the turbulence length scale, and its influence on the passage vortex on the suction side of the turbine blade.
Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.
The current work is based on observations of boiling heat transfer over a continuous range of gravity levels between 0g to 1.8g and varying heater sizes with a fluorinert as the test liquid (FC-72/n-perfluorohexane). Variable gravity pool boiling heat transfer measurements over a wide range of gravity levels were made during parabolic flight campaigns as well as onboard the International Space Station. For large heaters and-or higher gravity conditions, buoyancy dominated boiling and heat transfer results were heater size independent. The power law coefficient for gravity in the heat transfer equation was found to be a function of wall temperature under these conditions. Under low gravity conditions and-or for smaller heaters, surface tension forces dominated and heat transfer results were heater size dependent. A pool boiling regime map differentiating buoyancy and surface tension dominated regimes was developed along with a unified framework that allowed for scaling of pool boiling over a wide range of gravity levels and heater sizes. The scaling laws developed in this study are expected to allow performance quantification of phase change based technologies under variable gravity environments eventually leading to their implementation in space based applications.
Observations of solar-wind He(++) can be useful for studying dynamical processes in the interplanetary medium. Several processes which may influence the kinetic temperature are considered, and an attempt is made to inquire whether they can account for the observed fact that the kinetic temperature exceeds the proton temperature. Resonant heating is discussed qualitatively, and equations of resonant heating are developed. The proton and He(++) heating rates are determined for the least damped magnetoacoustic wave under various conditions.
A closed-form solution for cylindrical thermal waves in metals is given based on the nonlinear hyperbolic system of energy-conservation and heat-flux relaxation equations. It is shown that heat released from a line source propagates radially outward with finite speed in the form of a thermal wave which exhibits a discontinuous wave front. Unique nonlinear thermal-wave solutions exist up to a critical amount of driving energy, i.e., for larger energy releases, the thermal flow becomes multivalued (occurrence of shock waves). By comparison, it is demonstrated that the parabolic thermal-wave theory gives, in general, a misleading picture of the profile and propagation of thermal waves and leads to physical (infinite speed of heat propagation) and mathematical (divergent energy integrals) difficulties. Attention is drawn to the importance of temporal heat-flux relaxation for the physical understanding of fast transient processes such as thermal waves and more general explosions and implosions.
A study has been made of the heat requirement for the cyclic de-icing of hollow steel propellers fitted with external blade heating shoes. Solutions to the equations for the heat flow in cyclic heating of propellers were obtained, using an electrical analogy. The study showed how the energy requirement for propeller de-icing with existing blade shoes could be decreased, and illustrated the effect of blade-shoe design on the energy requirement. It was demonstrated, for example, that by increasing the heating intensity and decreasing the heating period from those currently used the energy requirement could be decreased in the order of 60 percent. ' In addition, ft was shown that heating requirements could be decreased further, by as much as 60 percent, through proper design of the shoes. The' investigation also showed the energy requirement to increase with decreasing liquid-water content and air temperature. Uncertainties as to the exact values of convective heat-transfer coefficient prevailing over the surface of the blade and ice layer resulted in uncertainties of approximately proportional magnitude in the values of required heating intensity.