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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2020

Convectively Forced Diurnal Gravity Waves in the Maritime Continent

Long-lived, zonally propagating diurnal rainfall disturbances are a highly pronounced and common feature in the Maritime Continent (MC). A recent study argues that these disturbances can be explained as diurnally phase-locked gravity waves. Here we explore the origins of these waves through regional cloud-permitting numerical model experiments. The gravity waves are reproduced and isolated in the model framework through the combined use of realistic geography and diurnally cyclic lateral boundary conditions representative of both characteristic easterly and westerly background zonal flow regimes. These flow regimes are characteristic of the Madden–Julian oscillation (MJO) suppressed and active phase in the MC, respectively. Tests are conducted wherein Borneo, Sumatra, or both islands and/or their orography are removed. These tests imply that the diurnal gravity waves are excited and maintained directly by latent heating from the vigorous mesoscale convective systems (MCSs) that form nocturnally in both Borneo and Sumatra. Removing orography has only a secondary impact on both the MCSs and the gravity waves, implying that it is not critical to these waves. We therefore hypothesize that diurnal gravity waves are fundamentally driven by mesoscale organized deep convection, and are only sensitive to orography to the measure that the convection is affected by the orography and its mesoscale flows. Factor separation further reveals that the nonlinear interaction of synchronized diurnal cycles in Sumatra and Borneo slightly amplifies this gravity wave mode compared to if either island existed in isolation. This nonlinear feedback appears most prominently at longitudes directly between the two islands.

Ruppert, Jr., James H.↗
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.↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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.↗
NASA NTRS1982

Toward a unified theory of atmospheric convective instability

A nonlinear three-dimensional truncated spectral model of shallow and moist Boussinesq convection indicates that parallel instability and thermal forcing are linked, in view of the fact that only one convective mode exists in which either or both mechanisms are operating to generate convection in the planetary boundary layer. It is also established that the wind field causes two-dimensional roll convection formation, an alignment of the convection with the wind in a preferred manner, and a propagation speed that is related to the wind component perpendicular to the roll axis. Latent heating is responsible for the decrease of the critical value of the environmental lapse rate in accordance with the slice method stability criterion. When only the upper part of the upward branch is moist and all of the downward branch is dry, latent heating also causes a finite-amplitude convective solution for Rayleigh number values lower than the critical value of linear analysis.

Shirer, H. N.↗
NASA NTRS1990

Theoretical and experimental studies of baroclinic processes

The geophysical fluid flow simulation code, GEOSIM, is being used to study the phenomenon of vacillation in the baroclinic annulus. Having verified that the code predicts vacillation for the same points as the experiments, the work is aiming toward explaining the mechanics of vacillation and pointing out some of the sensitivities of the results to the numerical method. Researchers are finding that there is a structural change associated with amplitude vacillation, where the structural changes are in the vertical. The results disagree with the premise of Lindzen et al, that the vacillation is due to constructive and destructive interference of neutral modes with different phase speeds. The researchers are continuing to study the Spacelab 3 Geophysical Fluid Flow Cell (GFFC) results with horizontal temperature gradients and heating from below. GEOSIM has been used to compute a wide range of cases, and these are being compared with the observations. The computations and observations compare well, and the model is being used to extend the results beyond cases studied in the experiments and to study the mechanics and predictability of the flows. The study of fully nonlinear baroclinic instability using the GFFC apparatus is proceeding with the numerical code. While the first instability that occurs is of planetary scale, secondary instabilities consisting of small-scale, penetrative convection occurs where cold fluid flows over a warm surface. The simultaneous modeling of the planetary scale and the convective scale is possible because of the nonhydrostatic formulation of the model. Some of these results have been animated on the Stardent computer, which shows the explosive nature of the small-scale convection.

Miller, Tim↗
NASA NTRS2001

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗
DOE OSTI2025

Nonlinearity of the post-spinel transition and its expression in slabs and plumes worldwide

Phase transitions in the mantle control its internal dynamics and structure. The post-spinel transition marks the upper–lower mantle boundary, where ringwoodite dissociates into bridgmanite plus ferropericlase, and its Clapeyron slope regulates mantle flow across it. This interaction has previously been assumed to have no lateral spatial variations, based on the assumption of a linear post-spinel boundary in pressure and temperature. Here we present laser-heated diamond anvil cell experiments with synchrotron X-ray diffraction to better constrain this boundary, especially at higher temperatures. Combining our data with results from the literature, and using a global analysis based on machine learning, we find a pronounced nonlinearity in the post-spinel boundary, with its slope ranging from –4 MPa/K at 2100 K, to –2 MPa/K at 1950 K, and to 0 MPa/K at 1600 K. Changes in temperature over time and space can therefore cause the post-spinel transition to have variable effects on mantle convection and the movement of subducting slabs and upwelling plumes.

58 GEOSCIENCES↗
NASA NTRS1987

Nonlinear heat transfer and structural analyses of SSME turbine blades

Three-dimensional nonlinear finite-element heat transfer and structural analyses were performed for the first stage high-pressure fuel turbopump blade of the space shuttle main engine (SSME). Directionally solidified (DS) MAR-M 246 material properties were considered for the analyses. Analytical conditions were based on a typical test stand engine cycle. Blade temperature and stress-strain histories were calculated using MARC finite-element computer code. The study was undertaken to assess the structural response of an SSME turbine blade and to gain greater understanding of blade damage mechanisms, convective cooling effects, and the thermal-mechanical effects.

Abdul-Aziz, A.↗
NASA NTRS1988

Thermal-structural analyses of Space Shuttle Main Engine (SSME) hot section components

Three dimensional nonlinear finite element heat transfer and structural analyses were performed for the first stage high pressure fuel turbopump (HPFTP) blade of the space shuttle main engine (SSME). Directionally solidified (DS) MAR-M 246 and single crystal (SC) PWA-1480 material properties were used for the analyses. Analytical conditions were based on a typical test stand engine cycle. Blade temperature and stress strain histories were calculated by using the MARC finite element computer code. The structural response of an SSME turbine blade was assessed and a greater understanding of blade damage mechanisms, convective cooling effects, and thermal mechanical effects was gained.

Abdul-Aziz, Ali↗
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗