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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Boundary-element shape sensitivity analysis for thermal problems with nonlinear boundary conditions

Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. 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 are presented that have closed-form solutions.

Kane, James H.↗
NASA NTRS1988

Structural analyses of engine wall cooling concepts and materials

The severe thermal environments under which hypersonic aircraft such as the National Aerospace Plane (NASP) will operate require cooling of the engine walls, especially in the combustor. A preliminary assessment is made of some candidate materials based on structural analyses for a number of convective cooling configurations. Three materials were studied: graphite/copper and tungsten/copper composite alloys with 50 percent fiber volume fractions and a wrought cobalt-base superalloy, Haynes 188. Anisotropic mechanical and thermal properties for the composites were obtained from a computer code, ICAN, which determines the composite material properties from the individual properties of the fiber and matrix materials. The structural analyses were performed by using the MARC nonlinear finite element code. Heat transfer analyses were conducted to calculate the metal temperature distributions.

Kaufman, Albert↗
NASA NTRS1990

Cellular structure of lean hydrogen flames in microgravity

Detailed, time-dependent, two-dimensional numerical simulations of premixed laminar flames have been used to study the initiation and subsequent development of cellular structures in lean hydrogen-air flames. The model includes detailed hydrogen-oxygen combustion with 24 elementary reactions of eight reactive species and a nitrogen diluent, molecular diffusion of all species, thermal conduction, viscosity, and convection. This model has been used to study the nonlinear evolution of cellular flame structure and shows that cell splitting, as observed in experiments, can be predicted numerically for sufficiently reactive mixtures. The structures that evolved also resembled the cellular structures observed in experiments. The present study shows that the 'cell-split limit' postulated from experimental observations is an intrinsic property of the mixture and that external factors such as heat losses are not necessary to cause this limit.

Patnaik, G.↗
NASA NTRS1994

Multiscale low-frequency circulation modes in the global atmosphere

In this paper, fundamental multiscale circulation modes in the global atmosphere are identified with the objective of providing better understanding of atmospheric low-frequency variabilities over a wide range of spatial and temporal scales. With the use of a combination of rotated principal component technique, singular spectrum analysis, and phase space portraits, three categories of basic multiscale modes in the atmosphere are found. The first is the interannual-mode (IAM), which is dominated by time scales longer than a year and can be attributed to heating and circulation anomalies associated with the coupled tropical ocean-atmosphere, in particular the El Nino-Southern Oscillation. The second is a set of tropical intraseasonal modes consisting of three separate multiscale patterns (ISO-1, -2, -3) related to tropical heating that can be identified with the different phases of the Madden-Julian Oscillation (MJO), including its teleconnection to the extratropics. The ISO spatial and temporal patterns suggest that the extratropical wave train in the North Pacific and North America is related to heating over the Maritime Continent and that the evolution of the MJO around the equator may require forcing from the extratropics spawning convection over the Indian Ocean. The third category represents extratropical intraseasonal oscillations arising from internal dynamics of the basic-state circulation. In the Northern Hemisphere, there are two distinct circulation modes with multiple frequencies in this category: the Pacific/North America (PNA) and the North Atlantic/Eurasia (NAE). In the Southern Hemisphere, two phase-locked modes (PSA-1 and PSA-2) are found depicting an eastward propagating wave train from eastern Australia, via the Pacific South America to the South Atlantic. The extratropical modes exhibit temporal characteristics such as phase locking and harmonic oscillations possibly associated with quadratically nonlinear dynamical systems. Additionally, the observed monthly and seasonal anomalies arise from a complex interplay of the various multiscale low-frequency modes. The relative dominance of the different modes varies widely from month to month and from year to year. On the monthly time scale, while one or two mechanisms may dominate in one year, no single mechanism seems to dominate for all years. There are indications that when the IAM, that is, ENSO heating patterns are strong, the extratropical modes may be suppressed and vice versa. For the seasonal mean, the interannual mode tends to dominate and the contribution from the PNA remains quite significant.

Lau, K.-M.↗
NASA NTRS1990

Magnetohydrodynamic modeling of the solar corona

The ideal and resistive magnetohydrodynamic (MHD) model is used to examine the dynamics and structure of the solar corona. When the coronal magnetic field is deformed by photospheric flow it can evolve to states that become unstable to ideal MHD modes. The nonlinear evolution of these instabilities can lead to the generation of current sheets, field line reconnection, and energy release. The disruption of an arcade field and the kinking of coronal loops is described. The braiding of the large-scale coronal field by convective photospheric motions develops fine-scale structure in the magnetic field and leads to the development of intense current filaments. The resistive dissipation of these currents can provide an efficient coronal heating mechanism.

Mikic, Zoran↗
NASA NTRS1977

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗
NASA NTRS1989

Three-dimensional thermal convection in a spherical shell

Nonlinear convective solutions are presented for a shell with properties characteristic of the earth's whole mantle for Rayleigh numbers up to 70,000. The solutions are validated numerically, and two distinct convective patterns (cubic and tetrahedral) are identified which are closely related to the geometric planforms predicted by the analytical theories of slightly supercritical spherical convection. A quantitative analysis of the horizontal and vertical structures of the velocity and temperature fields of the solutions is presented, and their heat transport properties are examined, including the total heat flow and the spatial distribution of the heat flux at the shell boundaries.

Bercovici, D.↗
NASA NTRS1986

Efficient numerical techniques for complex fluid flows

The central feature in any flow prediction method is the treatment of the coupling between the momentum and continuity equations. In natural-convection flows, the energy equation also becomes strongly coupled with the momentum equations. Because of the nonlinear nature of the coupling, these equations are solved iteratively. Iterative methods are often prone to slow convergence, divergence, and extreme sensitivity to underrelaxation factors. The aim of the present research is to develop more efficient and reliable solution schemes for the coupled flow equations. Such schemes will significantly reduce the expense of computing complex flows encountered in combustion chambers, gas turbines, heat exchangers, and other practical equipment. In the work completed so far, a technique employing norm reduction in conjunction with the successive-substitution and Newton-Raphson techniques was developed. Also, a block-correction procedure for the flow equations is currently being formulated and tested.

Patankar, Suhas V.↗
NASA NTRS1992

Treating convection in sequential solvers

The treatment of the convection terms in the sequential solver, a standard procedure found in virtually all pressure based algorithms, to compute the flow problems with sharp gradients and source terms is investigated. Both scalar model problems and one-dimensional gas dynamics equations have been used to study the various issues involved. Different approaches including the use of nonlinear filtering techniques and adoption of TVD type schemes have been investigated. Special treatments of the source terms such as pressure gradients and heat release have also been devised, yielding insight and improved accuracy of the numerical procedure adopted.

Shyy, Wei↗
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
NASA NTRS1983

Multiple buoyancy driven flows in a vertical cylinder heated from below

The structure of axisymmetric buoyancy-driven convection in a vertical cylinder heated from below is probed by finite element solution of the Boussinesq equations coupled with computed-implemented perturbation techniques for detecting and tracking multiple flows and for determining flow stability. Results are reported for fluids with Prandtl number of one and for cylinders with aspect ratio (Lambda) (defined as the height to radius of the cylinder) between 0.5 and 2.25. Extensive calculations of the neutral stability curve for the static solution and of the nonlinear motions along the bifurcating flow families show a continuous evolution of the primary cellular motion from a single toroidal cell to two and three cells nested radially in the cylinder, instead of the sharp transitions found for a cylinder with shear-free sidewalls. The smooth transitions in flow structure with Rayleigh number and lambda are explained by nonlinear connectivity between the first two bifurcating flow families formed either by a secondary bifurcation point for Lambda or = Lambda * approximately 0.80 or by a limit point for Lambda Lambda *. The transition between these two modes may be described by the theory of multiple limit point bifurcation.

Yamaguchi, Y.↗
NASA NTRS1985

Convective motions in a spherical shell

We compute the axisymmetric convective motions that exist in a spherical shell heated from below with inner to outer radius ratio equal to 0.5. The boundaries are stress-free and gravity is directly proportional to radius. Accurate solutions at large Rayleigh numbers, O(100000), are made feasible by a spectral method that employs diagonal-mode truncation. By examining the stability of axisymmetric motions it is inferred that the preferred form of convection varies dramatically according to the value of the Rayleigh number. While axisymmetric motions with different patterns may exist for modestly nonlinear convection, only a single motion persists at sufficiently large values of the Rayleigh number. This circulation is symmetric about the equator and has two meridional cells with rising motion at the poles. Instability of this single axisymmetric motion determines that the preferred pattern of three-dimensional convection has one azimuthal wave.

Zebib, A.↗
NASA NTRS1983

Nonlinear anelastic modal theory for solar convection

Solar envelope models are developed using single-mode anelastic equations as a description of turbulent convection which provide estimates for the variation with depth of the largest convective cellular flows, with horizontal sizes comparable to the total depth of the convection zone. These models can be used to describe compressible motions occurring over many density scale heights. Single-mode anelastic solutions are obtained for a solar envelope whose mean stratification is nearly adiabatic over most of its vertical extent because of the enthalpy flux explicitly carried by the big cell, while a subgrid scale representation of turbulent heat transport is incorporated into the treatment near the surface. It is shown that the single-mode equations allow two solutions for the same horizontal wavelength which are distinguished by the sense of the vertical velocity at the center of the three-dimensional cell. It is found that the upward directed flow experiences large pressure effects which can modify the density fluctuations so that the sense of the buoyancy force is changed, with buoyancy braking actually achieved near the top of the convection zone. It is suggested that such dynamical processes may explain why the amplitudes of flows related to the largest scales of convection are so weak in the solar atmosphere.

Latour, J.↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
NASA NTRS1985

Theoretical interpretation of thermal arrest data from differential thermal analysis measurements

A theory was developed to interpret the thermal arrest data of HgCdTe system from differential thermal analysis measurements. A nonlinear regression program is used to model fit the data. The purpose is to obtain a set of parameters which will give the information about average radiation coefficient, latent heat of fussion, and specific heat of the HgCdTe system at various compositions. The liquidus and solidus at various compositions are also obtained. The heat radiation from the furnace, heat radiation from the system, heat capacity from both the system and the quartz tube and the heat of fusion of the system are considered, however, it neglects heat transfer through conduction and convection.

Wang, J. C.↗
NASA NTRS1992

Observations of short large-amplitude magnetic structures at a quasi-parallel shock

Results of a detailed analysis of short large-amplitude magnetic structures (SLAMS) observed at an encounter of the quasi-parallel blow shock by the AMPTE UKS and IRM satellites are presented. Isolated SLAMS, surrounded by solar wind conditions, and embedded SLAMS, which lie within or form the boundary with regions of significant heating and deceleration, are identified. The duration, polarization, and other characteristics of SLAMS are all consistent with their growth directly out of the ULF wave field, including the common occurrence of an attached whistler as found in ULF shocklets. The plasma rest frame propagation speeds and two-spacecraft time delays for all cases show that the SLAMS attempt to propagate upstream against the oncoming flow, but are convected back downstream. The speeds and delays vary systematically with SLAMS amplitude in the way anticipated from nonlinear wave theory, as do their polarization features. Inter-SLAMS regions and boundary regions with solar wind contain hot deflected ions of lesser density than within the SLAMS.

Schwartz, Steven J.↗
NASA NTRS1991

Laboratory and theoretical models of planetary-scale instabilities and waves

Meteorologists and planetary astronomers interested in large-scale planetary and solar circulations recognize the importance of rotation and stratification in determining the character of these flows. The two outstanding problems of interest are: (1) the origins and nature of chaos in baroclinically unstable flows; and (2) the physical mechanisms responsible for high speed zonal winds and banding on the giant planets. The methods used to study these problems, and the insights gained, are useful in more general atmospheric and climate dynamic settings. Because the planetary curvature or beta-effect is crucial in the large scale nonlinear dynamics, the motions of rotating convecting liquids in spherical shells were studied using electrohydrodynamic polarization forces to generate radial gravity and centrally directed buoyancy forces in the laboratory. The Geophysical Fluid Flow Cell (GFFC) experiments performed on Spacelab 3 in 1985 were analyzed. The interpretation and extension of these results have led to the construction of efficient numerical models of rotating convection with an aim to understand the possible generation of zonal banding on Jupiter and the fate of banana cells in rapidly rotating convection as the heating is made strongly supercritical. Efforts to pose baroclinic wave experiments for future space missions using a modified version of the 1985 instrument have led us to develop theoretical and numerical models of baroclinic instability. Some surprising properties of both these models were discovered.

Hart, John E.↗
NASA NTRS1993

Marangoni instability in a liquid layer with two free surfaces

The onset of the Marangoni instability in a liquid layer with two free nearly insulating surfaces heated from below is studied. Linear stability analysis yields a condition for the emergence of a longwave or a finite wavelength instability from the quiescent equilibrium state. Using the method of asymptotic expansions, a weakly nonlinear evolution equation describing the spatiotemporal behavior of the velocity and temperature fields at the onset of the longwave instability is derived. The latter is given by delta(M) = 24, delta(M) being the difference between the upper and the lower Marangoni numbers. It is shown that in some parametric range one convective cell forms across the layer, while in other parametric domains two convective cells emerge between the two free surfaces.

Deissler, Robert J.↗