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 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 NTRS2022

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow↗

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.↗

NASA NTRS1995

Dynamics of Vapour Bubbles in Nucleate Boiling: Basic Equations of Bubble Evolution - 1

We consider the behaviour of a vapour bubble formed at a nucleation site on a heated horizontal wall. There is no forced convection of an ambient liquid, and the bubble is presumably separated from the wall by a thin liquid microlayer. The energy conservation law results in a variational equation for the mechanical energy of the whole system consisting of the bubble and liquid. It leads to a set of two strongly nonlinear equations which govern bubble expansion and motion of its centre of mass. A supplementary equation to find out the vapour temperature follows from consideration of heat transfer to the bubble, both from the bulk of surrounding liquid and through the microlayer. The average thickness of the microlayer is shown to increase monotonously with time as the bubble meniscus spreads along the wall. Bubble expansion is driven by the pressure head between vapour inside and liquid far away from the bubble, with due allowance for surface tension and gravity effects. It is resisted by inertia of liquid being placed into motion as the bubble grows. The inertia originates also a force that presses the bubble to the wall. This force is counteracted by the buoyancy and an effective surface tension force that tends to transform the bubble into a sphere. The analysis brings about quite a new formulation of the familiar problem of bubble growth and detachment under conditions of nucleate pool boiling.

Buyevich, Yu A.↗

NASA NTRS1983

Ignition of a combustible half space

A half space of combustible material is subjected to an arbitrary energy flux at the boundary where convection heat loss is also allowed. An asymptotic analysis of the temperature growth reveals two conditions necessary for ignition to occur. Cases of both large and order unity Lewis number are shown to lead to a nonlinear integral equation governing the thermal runaway. Some global and asymptotic properties of the integral equation are obtained.

Olmstead, W. E.↗

NASA NTRS1992

Numerical models of mantle convection

An overview of numerical methods describing the structure and dynamics of the mantle is presented with attention given to novel 3D modeling techniques. The paper reviews 3D spherical and Cartesian models for constant viscosity emphasizing the assumptions regarding style of convection, time dependence, and implications for the mantle. Similarly treated are 3D Cartesian models with temperature-dependent viscosities, and briefly examined are models that are based on compressibility, nonlinear viscosity, or plates. Extensive illustrations are presented detailing: (1) temperature variations from models of 3D thermal convection in spherical shells; (2) thermal anomalies in equatorial cross sections; and (3) temperature variations in a spherical shell heated from within. The discussion relates the numerical results of the models with real mantle-convection events, and the simulations are shown to yield increasingly realistic representations of material behavior.

Schubert, G.↗

NASA NTRS1992

On the physics of waves in the solar atmosphere: Wave heating and wind acceleration

In the area of solar physics, new calculations of the acoustic wave energy fluxes generated in the solar convective zone was performed. The original theory developed was corrected by including a new frequency factor describing temporal variations of the turbulent energy spectrum. We have modified the original Stein code by including this new frequency factor, and tested the code extensively. Another possible source of the mechanical energy generated in the solar convective zone is the excitation of magnetic flux tube waves which can carry energy along the tubes far away from the region. The problem as to how efficiently those waves are generated in the Sun was recently solved. The propagation of nonlinear magnetic tube waves in the solar atmosphere was calculated, and mode coupling, shock formation, and heating of the local medium was studied. The wave trapping problems and evaluation of critical frequencies for wave reflection in the solar atmosphere was studied. It was shown that the role played by Alfven waves in the wind accelerations and the coronal hole heating is dominant. Presently, we are performing calculations of wave energy fluxes generated in late-type dwarf stars and studying physical processes responsible for the heating of stellar chromospheres and coronae. In the area of physics of waves, a new analytical approach for studying linear Alfven waves in smoothly nonuniform media was recently developed. This approach is presently being extended to study the propagation of linear and nonlinear magnetohydrodynamic (MHD) waves in stratified, nonisothermal and solar atmosphere. The Lighthill theory of sound generation to nonisothermal media (with a special temperature distribution) was extended. Energy cascade by nonlinear MHD waves and possible chaos driven by these waves are presently considered.

Musielak, Z. E.↗

NASA NTRS1994

Coupling of microprocesses and macroprocesses due to velocity shear: An application to the low-altitude ionosphere

Recent observations indicate that low-altitude (below 1500 km) ion energization and thermal ion upwelling are colocated in the convective flow reversal region. In this region the convective velocity V(sub perpendicular) is generally small but spatial gradients in V(sub perpendicular) can be large. As a result, Joule heating is small. The observed high level of ion heating (few electron volts or more) cannot be explained by classical Joule heating alone but requires additional heating sources such as plasma waves. At these lower altitudes, sources of free energy are not obvious and hence the nature of ion energization remains ill understood. The high degree of correlation of ion heating with shear in the convective velocity (Tsunoda et al., 1989) is suggestive of an important role of velocity shear in this phenomenon. We provide more recent evidence for this correlation and show that even a small amount of velocity shear in the transverse flow is sufficient to excite a large-scale Kelvin-Helmholtz mode, which can nonlinearly steepen and give rise to highly stressed regions of strongly sheared flows. Futhermore, these stressed regions of strongly sheared flows may seed plasma waves in the range of ion cyclotron to lower hybrid frequencies, which are potential sources for ion heating. This novle two-step mechanism for ion energization is applied to typical observations of low-altitude thermal ion upwelling events.

Ganguli, G.↗