SEARCH · Search NASA

Results for “nonlinear convective heating”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Hybrid perturbation/Bubnov-Galerkin technique for nonlinear thermal analysis

A two step hybrid analysis technique to predict the nonlinear steady state temperature distribution in structures and solids is presented. The technique is based on the regular perturbation expansion and the classical Bubnov-Galerkin approximation. The functions are obtained by using the regular perturbation method. These functions are selected as coordinate functions and the classical Bubnov-Galerkin technique is used to compute their amplitudes. The potential of the proposed hybrid technique for the solution of nonlinear thermal problems is discussed. The effectiveness of this technique is demonstrated by the effects of conduction, convection, and radiation modes of heat transfer. It is indicated that the hybrid technique overcomes the two major drawbacks of the classical techniques: (1) the requirement of using a small parameter in the regular perturbation method; and (2) the arbitrariness in the choice of the coordinate functions in the Bubnov-Galerkin technique. The proposed technique extends the range of applicability of the regular perturbation method and enhances the effectiveness of the Bubnov-Galerkin technique.

Noor, A. K.↗

NASA NTRS1980

Ion heating by strong electrostatic ion cyclotron turbulence

A theory of the ion heating due to electrostatic ion cyclotron (EIC) waves in the auroral zone is presented. Due to the slowly convecting nature of the EIC mode, quasi-linear plateau formation cannot stabilize the waves, and growth occurs until the nonlinear mechanisms of ion resonance broadening and electron trapping provide saturation. The large amplitude and coherent nature of the resulting wave imply that quasi-linear theory provides only a lower limit to the ion heating. An upper bound on the heating rate is derived using a time-average model of ion dynamics in the coherent waves. The effects of ion heating in the presence of the magnetic gradient force and parallel electric fields are considered, with the result that perpendicular energies over 100 eV are easily attainable from a 1 eV source plasma. Perpendicular heating in the absence of a parallel electric field yields conical ion distributions, which in the presence of an electric field become field-aligned beams.

Lysak, R. L.↗

NASA NTRS1980

Thermal Pollution Mathematical Model. Volume 3: User's Manual for One-Dimensional Numerical Model for the Seasonal Thermocline

A user's manual for a one dimensional thermal model to predict the temperature profiles of a deep body of water for any number of annual cycles is presented. The model is essentially a set of partial differential equations which are solved by finite difference methods using a high speed digital computer. The model features the effects of area change with depth, nonlinear interaction of wind generated turbulence and buoyancy, adsorption of radiative heat flux below the surface, thermal discharges, and the effects of vertical convection caused by discharge. The main assumption in the formulation is horizontal homogeneity. The environmental impact of thermal discharges from power plants is emphasized. Although the model is applicable to most lakes, a specific site (Lake Keowee, S.C.) application is described in detail. The programs are written in FORTRAN 5.

Lee, S. S.↗

DOE OSTI2024

Development of Modified Perturbation Solutions to the One-Phase Stefan Problems With a Convective Boundary

The classical Stefan problem is used to track the moving solid-liquid interface during the freezing process. Perturbation theory has often been applied to find an approximate analytical solution due to the nonlinearity of the moving interface. However, the Stefan number (i.e., the sensible over latent heat) must be small and usually less than 0.01 to assume the perturbation expansion, which in turn limits the thermal engineering applications. In this study, a modified perturbation solution is developed by adding a correction term after the leading-order solution to be valid for a much wider range of Stefan numbers (i.e., 0.01 less than or equal to Ste less than or equal to 1). Specifically, a one-phase Stefan problem is first formulated subjected to a convective boundary in the Cartesian, cylindrical, and spherical coordinate systems for inward solidification. The leading-order solution is calculated based on the regular perturbation theory, while the correction term is obtained using the Monte-Carlo method and a multi-variant regression. Results show that the correction term has a linear relationship with the Stefan number and is not significantly influenced by the Biot number. The proposed modified perturbation solution can accurately and rapidly predict the nonlinear moving interface motion for the freezing process.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Computer simulation of the Farley-Buneman instability and anomalous electron heating in the auroral ionosphere

The nonlinear saturation of the Farley-Buneman instability in a collisional plasma is studied by a 2 1/2 dimensional electrostatic particle simulation which includes inelastic and elastic collisions of electrons and elastic collision of ions with neutrals. In the simulation, a uniform convection electric field is applied externally so that the relative velocity between the electrons and ions is greater than the ion sound speed and destabilizes the instability. A nonlinear frequency shift from higher to lower frequencies and diffusion of the wave spectrum in two dimensional wave number space are found. It is found that the dominant mechanism for electron heating is due to an enhanced effective electron collision frequency and hence enhanced resistive heating as suggested by Primdahl (1986) and not due to the heating of electrons by the electric field of the waves parallel to the magnetic field. For the ionospheric conditions discussed by Schlegel and St.-Maurice (1981)an anomalous heating rate of about 4 x 10 to the 7th W/cu cm is found.

Machida, S.↗

NASA NTRS1976

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

NASA NTRS1980

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

NASA NTRS1992

Boundary formulations for shape sensitivity of temperature dependent conductivity problems

Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.

Kane, James H.↗

NASA NTRS1995

Low-frequency oscillations in radiative-convective systems. Part 2: An idealized model

A simple model is used to examine the hypothesis that nonlinear interactions among atmospheric radiation, cumulus convection, and the surface moisture flux can result in a stationary, low-frequency (30-60 day period) oscillating heat source in the tropical atmosphere. The model produces low-frequency oscillations of temperature, moisture, and precipitation. The mechanism that produces these oscillations is identified through analyses of the model and its results. The relevance of this mechanism to understand the observed Madden-Julian oscillation in the tropical atmosphere over the Indian and western Pacific Oceans is discussed.

Hu, QI↗

NASA NTRS1980

Waves in space plasmas - The mirror trapping of hot auroral electrons

A brief review is given of the problem of precipitation of auroral electrons by electrostatic Bernstein waves. Since the magnetospheric loss cone is small, only moderately small intense levels of wave turbulence are required to remove any large anisotropy sources of free energy and to maintain a weakly anisotropic electron distribution on strong diffusion precipitation. The electrostatic electron cyclotron harmonic waves are nonconvectively unstable for weak loss cone anisotropies and over a large range of parameters for both the hot and cold distributions. Since the instability is nonconvective, weak wave growth can be maintained independent of the flux level of the hot electrons, i.e., the instability does not have the stably trapped flux limit imposed by convective amplification. Recent plasma numerical simulations show that the nonlinear evolution of this instability involves both the pitch angle diffusion of the hot electrons and the heating of the cold electrons.

Ashour-Abdalla, M.↗

NASA NTRS1983

Reduction methods for nonlinear steady-state thermal analysis

Hybrid analysis techniques and a problem adaptive computational algorithm are presented for predicting the nonlinear steady state temperature distribution in structures and solids. In the proposed techniques, the structure is discretized by using the finite element method. The vector of nodal temperatures is then expressed as a linear combination of a small number of global temperature modes (or basis vectors), and the Bubnov Galerkin technique is used to compute the amplitudes of the global modes. The global temperature modes (or basis vectors) are chosen to include the various order derivatives of the nodal temperature vector with respect to preselected path parameter(s). The potential of the proposed reduction methods for solution of large scale, nonlinear thermal problems is discussed and the effectiveness of the methods is demonstrated by means of numerical examples, including steady state conduction, convection, and radiation modes of heat transfer.

Noor, A. K.↗

NASA NTRS1996

Magnetothermal Convection in Nonconducting Diamagnetic and Paramagnetic Fluids

Nonuniform magnetic fields exert a magnetic body force on electrically nonconducting classical fluids. These include paramagnetic fluids such as gaseous and liquid oxygen and diamagnetic fluids such as helium. Recent experiments show that this force can overwhelm the force of gravity even at the surface of the earth; it can levitate liquids and gases, quench candle flames, block gas flows, and suppress heat transport. Thermal gradients render the magnetic force nonuniform through the temperature-dependent magnetic susceptibility. These thermal gradients can therefore drive magnetic convection analogous to buoyancy-driven convection. This magnetothermal convection can overwhelm convection driven by gravitational buoyancy in terrestrial experiments. The objectives of the proposed ground-based theoretical study are (a) to supply the magnetothermohydrodynamic theory necessary to understand these recent experiments and (b) to explore the consequences of nonuniform magnetic fields in microgravity. Even the linear theory for the onset of magnetothermal convection is lacking in the literature. We intend to supply the linear and nonlinear theory based on the thermohydrodynamic equations supplemented by the magnetic body force. We intend to investigate the effect of magnetic fields on gas blockage and heat transport in microgravity. Since magnetic fields provide a means of creating arbitrary, controllable body force distributions, we intend to investigate the possibility of using magnetic fields to position and control fluids in microgravity. We also intend to investigate the possibility of creating stationary terrestrial microgravity environments by using the magnetic force to effectively cancel gravity. These investigations may aid in the design of space-based heat-transfer, combustion, and human-life-support equipment.

Edwards, Boyd F.↗

NASA NTRS1999

Studies in Thermocapillary Convection of the Marangoni-Benard Type

A nonlinear feedback control strategy for delaying the onset and eliminating the subcritical nature of long-wavelength Marangoni-Benard convection is investigated based on an evolution equation. A control temperature is applied to the lower wall in a gas-liquid layer otherwise heated uniformly from below. It is shown that, if the interface deflection is assumed to be known via sensing as a function of both horizontal coordinates and time, a control temperature with a cubic-order polynomial dependence on the deflection is capable of delaying the onset as well as eliminating the subcritical instability altogether, at least on the basis of a weakly nonlinear analysis. The control coefficients required for stabilization are O(1) for both delaying onset indefinitely and eliminating subcritical instability.

Kelly, R. E.↗

NASA NTRS1993

Analysis of Radiation-natural Convection Interactions in 1-g and low-g Environments using the Discrete Exchange Factor Method

A new numerical method is presented for the analysis of combined natural convection and radiation heat transfer with applications in many engineering situations such as materials processing, combustion and fire research. Because of the recent interest in the low gravity environment of space, attention is devoted to both 1-g and low-g applications. The two-dimensional mathematical model is represented by a set of coupled nonlinear integro-partial differential equations. Radiative exchange is formulated using the Discrete Exchange Factor method (DEF). This method considers point to point exchange and provides accurate results over a wide range of radiation parameters. Numerical results show that radiation significantly influences the flow and heat transfer in both low-g and 1-g applications. In the low-g environment, convection is weak, and radiation can easily become the dominant heat transfer mode. It is also shown that volumetric heating by radiation gives rise to an intricate cell pattern in the top heated enclosure.

Kassemi, M.↗

NASA NTRS1990

Cold plasma heating in the plasma sheet boundary layer - Theory and simulations

Satellite observations in recent years have confirmed that the plasma sheet boundary layer is a permanent feature of the earth's magnetotail located between the lobe and central plasma sheet during both quiet and active magnetic periods. Distinct features of the boundary layer include field aligned ion beams and intense electrostatic emissions known as broadband electrostatic noise. Since the plasma sheet boundary layer is a spatial feature of the magnetotail, within it will occur thermal mixing of the resident warm boundary layer plasma with inflowing (convecting) cold ionospheric plasma. A theoretical study involving linear theory and nonlinear numerical particle simulations is presented which examines ion beam instabilities in the presence of a thermally mixed hot and cold background plasma. It is found that the free energy in the ion beams can heat the cool ionospheric plasma to ambient plasma sheet boundary layer temperatures via broadband electrostatic noise. These results, along with recent observational reports that ionospheric outflow can account for measured plasma sheet densities, suggest that the ionospheric role in plasma sheet dynamics and content may be as large as the solar wind.

Schriver, David↗

NASA NTRS1973

Effects of interhemisphere transport on plasma temperatures at low latitudes.

The thermal balance of the equatorial plasma between 300 and 800 km is examined. Steady state nighttime calculations are made for O+, H+, and electrons. The following features are included: collisional heat transfer between ions, electrons, and neutrals; ion and electron thermal conduction along the field lines; curvature of the field lines; nonlinear advection due to field-aligned ion and electron motions; and convective compression or expansion due to field-aligned and E x B motions. The ion velocities necessary to calculate the effects of convection are obtained from the work of Moffett and Hanson, who include a meridional wind across the magnetic equator in their calculations. It is shown that field-aligned interhemisphere plasma flows appreciably affect the plasma temperatures.

Bailey, G. J.↗

NASA NTRS1975

Comparison of NASTRAN and MITAS nonlinear thermal analyses of a convectively cooled structure

Comparative steady state nonlinear thermal analyses of a scramjet fuel injection strut are presented. The analyses were performed using the NASTRAN finite element program and MITAS, a lumped-parameter thermal analyzer. The strut is subjected to aerodynamic heating on two sides and is internally cooled by hydrogen flowing from internal manifolds through heat exchangers bonded to the primary structure. Based on coolant temperatures determined by MITAS, NASTRAN predicted temperature distributions throughout the strut which were in close agreement with similar MITAS predictions.

Thornton, E. A.↗

NASA NTRS1991

A numerical modeling study of the interaction between the tides and the circulation forced by high-latitude plasma convection

A spectral, time-varying thermospheric general circulation model has been used to study the nonlinear interaction at high latitudes between the tides propagating into the thermosphere from below and the circulation induced by magnetospheric forcing and in situ solar heating. The model is discrete in the vertical with 27 layers spaced by half a scale height. In the horizontal, the fields are expanded in a series of spherical harmonics using a triangular truncation at wave number 31, equivalent to a homogeneous global resolution with a minimum wavelength of 1270 km. A hypothetical uniform grid point model would require a horizontal spacing of 417 km to describe the same minimum wavelength. In the high-latitude F region the tides affect the dusk vortex of the neutral flow very little, but the dawn vortex is either suppressed or amplified dependent upon the universal time and tidal phase. In the E region neutral flow, both the dusk and dawn vortices are shifted in local time by the tides, again as a function of universal time and tidal phase. At dusk a nonlinear amplification of the sunward winds occurs for certain combinations of parameters, and at dawn the winds may be completely suppressed. Below 120 km altitude the magnetospheric forcing creates a single cyclonic vortex which is also sensitive to the high-latitude tidal structure.

Mikkelsen, I. S.↗