Search NASA⌕ Search

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.

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

At least 37 records · Page 2

Self-adaptive closed constrained solution algorithms for nonlinear conduction

Self-adaptive solution algorithms are developed for nonlinear heat conduction problems encountered in analyzing materials for use in high temperature or cryogenic conditions. The nonlinear effects are noted to occur due to convection and radiation effects, as well as temperature-dependent properties of the materials. Incremental successive substitution (ISS) and Newton-Raphson (NR) procedures are treated as extrapolation schemes which have solution projections bounded by a hyperline with an externally applied thermal load vector arising from internal heat generation and boundary conditions. Closed constraints are formulated which improve the efficiency and stability of the procedures by employing closed ellipsoidal surfaces to control the size of successive iterations. Governing equations are defined for nonlinear finite element models, and comparisons are made of results using the the new method and the ISS and NR schemes for epoxy, PVC, and CuGe.

Padovan, J.↗

Numerical Simulations of Convectively Excited Gravity Waves

Magneto-convection and gravity waves are numerically simulated with a nonlinear, three-dimensional, time-dependent model of a stratified, rotating, spherical fluid shell heated from below. A Solar-like reference state is specified while global velocity, magnetic field, and thermodynamic perturbations are computed from the anelastic magnetohydrodynamic equations. Convective overshooting from the upper (superadiabatic) part of the shell excites gravity waves in the lower (subadiabatic) part. Due to differential rotation and Coriolis forces, convective cell patterns propagate eastward with a latitudinally dependent phase velocity. The structure of the excited wave motions in the stable region is more time-dependent than of the convective motions above. The magnetic field tends to be concentrated over giant-cell downdrafts in the convective zone but is affected very little by the wave motion in the stable region.

Glatzmaier, G. A.↗

The effect of ground hydrology on climate sensitivity to solar constant variations

A zonally averaged climate model is used to investigate the effects of different evaporation parameterizations on climate sensitivity to solar constant variations. The model is based on the annual mean zonally averaged hemispheric climate model of Ohring and Adler (1978), with some modifications on the heating parameterizations. The dynamic model utilizes the two-level quasi-geostrophic potential vorticity equations. The heating processes are related to solar radiation, long-wave radiation, convection, evaporation, latent heat release, and ocean transport. Evaporation parameterization has been improved by using a nonlinear formulation in which the Bowen ratio is determined by the vertical gradients of water vapor and temperature near the earth's surface.

Chow, S.-H.↗

Infinite Prandtl number thermal convection in a spherical shell

A Galerkin technique is used to calculate the steady-state axisymmetric nonlinear convective motions in an infinite-Prandtl-number Boussinesq fluid in a relatively thick spherical shell heated from below. A reasonably complete study of the properties of the even and general axisymmetric steady states is carried out for a range of moderately supercritical Rayleigh numbers. In addition, stability analyses are conducted to determine which form of axisymmetric steady convection is the preferred one and whether the axisymmetric steady flows are unstable to azimuthal perturbations.

Zebib, A.↗

Transient convection over the Amazon-Bolivia region and the dynamics of drought over Northeast Brazil

A two layer, nonlinear, equatorial beta-plane model, in p-coordinates is used to study the atmospheric response to a large scale prescribed heat source varying in time. The heat source is meant to represent a convective burst with total duration of approximately 48 hours over the Amazon/Bolivia region. The boundary conditions used are meridional velocity zero at 60 deg S, omega w = 0 at the top and zero geometric velocity at the lower boundary. Sensitivity study was done which includes initial state at rest, compared with realistic initial flow. The scale of the heat source is 1500 km in latitude and longitude and it is centered at 10 deg S. Special attention is paid to the distribution and intensity of the induced vertical motion. The model is integrated for two days and the preliminary results show agreement with the observed 200 mb flow. Of interest is the establishment of a trough and descending motion to the northeast of the heat source. A conjucture is thus made that the Amazon heat source and its fluctuations bear some relationship with the drought problem over Northeast Brazil.

Buchmann, J.↗

Two-dimensional hydrodynamic calculations of the nonlinear development of the Goldreich-Schubert-Fricke instability in a rotating annulus

Two-dimensional nonlinear hydrodynamic calculations are presented which may help assess the effectiveness of the instability in transporting angular momentum in the equatorial zones of stars and planets which are stably stratified with respect to convection. The calculations were made by numerically integrating the 2D axisymmetric Navier-Stokes equations, including viscosity and heat conduction. The instability was followed into the nonlinear regime. The maximum rms velocity amplitude was found to correlate well with the product of the linear growth rate and radial length scale of the instability, consistent with the idea that the instability grows to an amplitude such that an eddy turnover time becomes equal to the growth time defined by the inverse of the growth rate. The time scale for angular momentum to be redistributed to a state of marginal stability was consistent with this picture. The results suggest that in physical situations a state of marginal stability will be maintained, since departures from such a state will be rapidly corrected.

Korycansky, D. G.↗

Invisibility of Solar Active Region Umbra-to-Umbra Coronal Loops: New Evidence that Magnetoconvection Drives Solar-Stellar Coronal Heating

Coronal heating generally increases with increasing magnetic field strength: the EUV/X-ray corona in active regions is 10--100 times more luminous and 2--4 times hotter than that in quiet regions and coronal holes, which are heated to only about 1.5 MK, and have fields that are 10--100 times weaker than that in active regions. From a comparison of a nonlinear force-free model of the three-dimensional active region coronal field to observed extreme-ultraviolet loops, we find that (1) umbra-to-umbra coronal loops, despite being rooted in the strongest magnetic flux, are invisible, and (2) the brightest loops have one foot in an umbra or penumbra and the other foot in another sunspot's penumbra or in unipolar or mixed-polarity plage. The invisibility of umbra-to-umbra loops is new evidence that magnetoconvection drives solar-stellar coronal heating: evidently, the strong umbral field at both ends quenches the magnetoconvection and hence the heating. Our results from EUV observations and nonlinear force-free modeling of coronal magnetic field imply that, for any coronal loop on the Sun or on any other convective star, as long as the field can be braided by convection in at least one loop foot, the stronger the field in the loop, the stronger the coronal heating.

Tiwari, Sanjiv K.↗

Effects from equation of state and rheology in dissipative heating in compressible mantle convection

The effects of compressibility on mantle convection are considered, incorporating the effects of equations of state and rheology in the dissipative heating term of the energy equation. The ways in which compression may raise the interior mantle temperature are explicitly demonstrated, and it is shown how this effect can be used to constrain some of the intrinsic parameters associated with the equation of state in the mantle. It is concluded that the coupling between variable viscosity and equation of state in dissipative heating is potentially an important mechanism in mantle convection. These findings emphasize that rheology, equation of state, and radiogenic heating are all linked to each other by nonlinear thermomechanical couplings.

Yuen, David A.↗

Physical vapor transport of mercurous chloride under a nonlinear thermal profile

Our study investigates numerically the flow field characteristics during the growth of mercurous chloride (Hg2Cl2) crystals in a rectangular ampoule under terrestrial and microgravity conditions for a nonlinear thermal gradient. With a residual gas lighter than the nutrient, the solutal Grashof number is dominant. We observe that in tilted configurations, when solutal convection is dominant, the maximum transport rate occurs at approximately 40 percent. For the vertical configurations, we were able to obtain solutions only for the cases either below the critical Rayleigh numbers or the stabilized configurations. The total mass flux decreases exponentially with an increase of pressure of residual gas, but it increases following a power law with the temperature difference driving the transport. The nonlinear thermal gradient appears to destabilize the flow field when thermal convection is dominant for both vertical top-heated and bottom-heated configurations. However, when the solutal Grashof number is dominant, the density gradient resulting from the solutal gradient appears to stabilize the flow for the bottom-heated configuration. The flow field for the top-heated configuration is destabilized for high Grashof numbers. The microgravity environment provides a means for lowering convection. For gravity levels of 10(exp -3) g(0) or less, the Stefan wind drives the flow, and no recirculating cell is predicted.

Mennetrier, Christophe↗

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗

Adiabaticity and viscosity in deep mantle convection

A study has been conducted of steady convection with adiabatic and viscous heating for variable viscosity in the Boussinesq limit using the mean-field theory. A strong nonlinear coupling is found between the thermodynamic constants governing adiabatic heating and the rheological parameters. The range of rheological values for which adiabaticity would occur throughout the mantle has been established. Too large an activation volume, greater than 6 cu cm/mol for the cases examined, would produce unreasonably high temperature at the bottom of the mantle (greater than 6000 K) and superadiabatic gradients, especially in the lower mantle. Radiogenic heating plays a profound role in controlling dynamically mantle temperatures. Present values for the averaged mantle heat production would yield objectionably high temperatures in the lower mantle.

Quareni, F.↗

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗

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

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

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

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

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

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