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

The Effect of Surface Topography on the Nonlinear Dynamics of Rossby Waves

Boussinesq convection in rotating systems attracts a sustained attention of the fluid dynamics community, because it has intricate non-linear dynamics (Cross & Hohenberg 1993) and plays an important role in geophysical and astrophysical applications, such as the motion of the liquid outer core of Earth, the Red Spot in Jupiter, the giant cells in the Sun etc. (Alridge et al. 1990). A fundamental distinction between the real geo- and astrophysical problems and the idealized laboratory studies is that natural systems are inhomogeneous (Alridge et al. 1990). Heterogeneities modulate the flow and influence significantly the dynamics of convective patterns (Alridge et al. 1990; Hide 1971). The effect of modulations on pattern formation and transition to turbulence in Boussinesq convection is far from being completely understood (Cross & Hohenberg 1993; Aranson & Kramer 2002). It is generally accepted that in the liquid outer core of the Earth the transport of the angular momentum and internal heat occurs via thermal Rossby waves (Zhang et al. 2001; Kuang & Bloxham 1999). These waves been visualized in laboratory experiments in rotating liquid-filled spheres and concentric spherical shells (Zhang et al. 2001; Kuang & Bloxham 1999). The basic dynamical features of Rossby waves have been reproduced in a cylindrical annulus, a system much simpler than the spherical ones (Busse & Or 1986; Or & Busse 1987). For convection in a cylindrical annulus, the fluid motion is two-dimensional, and gravity is replaced by a centrifugal force, (Busse & Or 1986; Or & Busse 1987). Hide (1971) has suggested that the momentum and heat transport in the core might be influenced significantly by so-called bumps, which are heterogeneities on the mantle-core boundary. To model the effect of surface topography on the transport of momentum and energy in the liquid outer core of the Earth, Bell & Soward (1996), Herrmann & Busse (1998) and Westerburg & Busse (2001) have studied the nonlinear dynamics of thermal Rossby waves in a cylindrical annulus with azimuthally modulated height.

Abarzhi, S. I.↗

Finite element thermo-viscoplastic analysis of aerospace structures

The time-dependent thermo-viscoplastic response of aerospace structures subjected to intense aerothermal loads is predicted using the finite-element method. The finite-element analysis uses the Bodner-Partom unified viscoplastic constitutive relations to determine rate-dependent nonlinear material behavior. The methodology is verified by comparison with experimental data and other numerical results for a uniaxially-loaded bar. The method is then used (1) to predict the structural response of a rectangular plate subjected to line heating along a centerline, and (2) to predict the thermal-structural response of a convectively-cooled engine cowl leading edge subjected to aerodynamic shock-shock interference heating. Compared to linear elastic analysis, the viscoplastic analysis results in lower peak stresses and regions of plastic deformations.

Pandey, Ajay K.↗

Finite-element thermo-viscoplastic analysis of aerospace structures

The time-dependent thermo-viscoplastic response of aerospace structures subjected to intense aerothermal loads is predicted using the finite-element method. The finite-element analysis uses the Bodner-Partom unified viscoplastic constitutive relations to determine rate-dependent nonlinear material behavior. The methodology is verified by comparison with experimental data and other numerical results for a uniaxially-loaded bar. The method is then used (1) to predict the structural response of a rectangular plate subjected to line heating along a centerline, and (2) to predict the thermal-structural response of a convectively-cooled engine cowl leading edge subjected to aerodynamic shock-shock interference heating. Compared to linear elastic analysis, the viscoplastic analysis results in lower peak stresses and regions of plastic deformations.

Pandey, Ajay↗

The spectrum of steady state turbulent convection.

Based on Heisenberg's statistical theory of turbulence, a model for steady state turbulent convection is herein proposed, and on the basis of this model, equations for the energy spectrum for steady state turbulent convection are derived. The spectrum is obtained from the solution of a nonlinear integral equation. After the integral equation is brought into a universally valid nondimensional form, it is transformed into a nonlinear first order differential equation to be solved numerically, with the Rayleigh number appearing as the only parameter. The energy spectrum has a substantial deviation from the Kolmogoroff law, as a result of the buoyancy force acting on the rising and falling eddies. The presented theory may be applicable to convection in planetary and stellar atmospheres wherein the radiative heat transport is small.

Winterberg, F.↗

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is paramterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations.

Hanna, A. F.↗

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is parameterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations. Previously announced in STAR as N83-14822

Hanna, A. F.↗

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is parameterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations. Previously announced in STAR as N83-14822

Hanna, A. F.↗

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

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↗

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

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

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

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↗

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗

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

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

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

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↗