SEARCH · Search NASA
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.
Penetrative thermal convection stability for statically unstable layer surmounted by statically stable layer from numerical integration of linearized equations of motion and heat conduction
Experiments on thermal convection in a rotating, differentially-heated spherical shell with a radial buoyancy force were conducted in an orbiting microgravity laboratory. A variety of convective structures, or planforms, were observed depending on the magnitude of the rotation and the nature of the imposed heating distribution. The results are in agreement with numerical simulations that can be conducted at modest parameter values, and suggest possible regimes of motion in rotating planets and stars.
Microgravity conditions offer an environment in which convection in a nonisothermal gas could be driven primarily by thermal stress. A direct examination of thermal stress flows would be invaluable in assessing the accuracy of the Burnett terms in the fluid stress tensor. We present a preliminary numerical investigation of the competing effects of thermal stress, thermal creep at the side walls, and buoyancy on gas convection in nonuniformly heated containers under normal and reduced gravity levels. Conditions in which thermal stress convection becomes dominant are identified, and issues regarding the experimental measurement of the flows are discussed.
This study includes the development and design of a space experiment for the investigation of the thermally-induced convection heat transfer process in an enclosure containing air. Various features of the design including the optical method used in the fluid temperature measurement are presented and discussed. A series of ground experiments were carried out in support of the in-space experiment design. Their results are discussed briefly and compared with numerical predictions.
The project has sought to ascertain the veracity of the Burnett relations, as applied to slow moving, highly nonisothermal gases, by comparison of convection and stress predictions with those generated by the DSMC method. The Burnett equations were found to provide reasonable descriptions of the pressure distribution and normal stress in stationary gases with a 1-D temperature gradient. Continuum/Burnett predictions of thermal stress convection in 2-D heated enclosures, however, are not quantitatively supported by DSMC results. For such situations, it appears that thermal creep flows, generated at the boundaries of the enclosure, will be significantly larger than the flows resulting from thermal stress in the gas.
Models describing the steady-state thermal diffusion in a pure system, the thermal-solutal diffusion in a binary system, and heat and momentum transverse in a pure system are presented. The geometry of the model is described by a 2D Cartesian coordinate system that is applicable for crystal sheets. The melting, solidifying, and melt/gas interfacial shapes as well as the thermal, flow, and solutal profiles are analytically evaluated as functions of the heat and ambient temperature profiles and material properties. The solution procedure involves a coupled asymptotic/numerical approach which reduces the coupled set of partial differential equations to ordinary type. The results should be applicable in situations where melt flows are not intense enough to change the thermal field in pure systems, or where the physical properties of the melt are such that the convective field is decoupled from the thermal field, the latter being established primarily by diffusion.
One objective of this research proposal is to develop a 3-D thermal history model for Venus. The basis of our study is a finite-element computer model to simulate thermal convection of fluids with highly temperature- and pressure-dependent viscosities in a three-dimensional spherical shell. A three-dimensional model for thermal history studies is necessary for the following reasons. To study planetary thermal evolution, one needs to consider global heat budgets of a planet throughout its evolution history. Hence, three-dimensional models are necessary. This is in contrasts to studies of some local phenomena or local structures where models of lower dimensions may be sufficient. There are different approaches to treat three-dimensional thermal convection problems. Each approach has its own advantages and disadvantages. Therefore, the choice of the various approaches is subjective and dependent on the problem addressed. In our case, we are interested in the effects of viscosities that are highly temperature dependent and that their magnitudes within the computing domain can vary over many orders of magnitude. In order to resolve the rapid change of viscosities, small grid spacings are often necessary. To optimize the amount of computing, variable grids become desirable. Thus, the finite-element numerical approach is chosen for its ability to place grid elements of different sizes over the complete computational domain. For this research proposal, we did not start from scratch and develop the finite element codes from the beginning. Instead, we adopted a finite-element model developed by Baumgardner, a collaborator of this research proposal, for three-dimensional thermal convection with constant viscosity. Over the duration supported by this research proposal, a significant amount of advancements have been accomplished.
Not provided.
It is well known that in the liquid state eutectic alloys are theoretically homogeneous under 1 g conditions. However, the homogeneous solidified structure of this alloy is not obtained because thermal convection and non-equilibrium solidification occur. The present investigators have clarified the solidification mechanisms of the eutectic system alloys under 1 g conditions by using the in situ observation method; in particular, the primary crystals of the eutectic system alloys never nucleated in the liquid, but instead did so on the mold wall, and the crystals separated from the mold wall by fluid motion caused by thermal convection. They also found that the equiaxed eutectic grains (eutectic cells) are formed on the primary crystals. In this case, the leading phase of the eutectic must agree with the phase of the primary crystals. In space, no thermal convection occurs so that primary crystals should not move from the mold wall and should not appear inside the solidified structure. Therefore no equiaxed eutectic grains will be formed under microgravity conditions. Past space experiments concerning eutectic alloys were classified into two types of experiments: one with respect to the solidification mechanisms of the eutectic alloys and the other to the unidirectional solidification of this alloy. The former type of experiment has the problem that the solidified structures between microgravity and 1 g conditions show little difference. This is why the flight samples were prepared by the ordinary cast techniques on Earth. Therefore it is impossible to ascertain whether or not the nucleation and growth of primary crystals in the melt occur and if primary crystals influence the formation of the equiaxed eutectic grains. In this experiment, hypo- and hyper-eutectic aluminum copper alloys which are near eutectic point are used. The chemical compositions of the samples are Al-32.4mass%Cu (Hypo-eutectic) and Al-33.5mass%Cu (hyper-eutectic). Long rods for the samples are cast by the Ohno Continuous Casting Process and they show the unidirectionally solidified structure. Each flight and ground sample was made of these same rods. The dimensions of all samples are 4.5 mm in diameter and 23.5 mm in length. Each sample is put in a graphite capsule and then vacuum sealed in a double silica ampoule. Then the ampoule is put in the tantalum cartridge and sealed by electron beam welding. For onbard experiments, a Continuous Heating Furnance (CHF) will be used for melting and solidifying samples under microgravity conditions. Six flight samples will be used. Four samples are hypo-eutectic and two are hyper-eutectic alloys. The surface of the two hypo-eutectic alloy samples are covered with aluminum oxide film to prevent Marangoni convection expected under microgravity conditions. Each sample will be heated to 700 C and held at that temperature for 5 min. After that the samples will be allowed to cool to 500 C in the furnace and they will be taken out of the furnace for He gas cooling. The heating and cooling diagrams for the flight experiments are shown. After collecting the flight samples, the solidified structures of the samples will be examined and the mechanisms of eutectic solidification under microgravity conditions will be determined. It is likely that successful flight experiment results will lead to production of high quality eutectic alloys and eutectic composite materials in space.
The onset of thermal convection due to heating from below in a system consisting of a fluid layer overlying a porous layer with anisotropic permeability and thermal diffusivity is studied. The linear perturbation equations are solved numerically. It is found that the effects of anisotropy on the onset of thermal convection are most profound for small values of the depth ratio xi (ratio of fluid layer thickness to porous layer thickness), since in that case, the onset of convection corresponds to significant motion in both layers. For fixed vlues of the vertical permeability in the porous medium, decreasing the value of zeta (ratio of horizontal to vertical permeability) leads to stabilization of the superposed layer configuration because of increased resistance to motion in the porous medium. For larger values of xi, the onset of motion is increasingly confined to the fluid layer, with the transport of heat through the porous layer occurring primarily by conduction. Accordingly, the influence of zeta on the stability characteristics for larger xi is less significant than the effects of an anisotropic thermal conductivity.
Convection and macrosegregation in directionally solidified hypoeutectic Pb-38 wt pct Sn and hypereutectic Pb-64.5 wt pct Sn have been examined during upward and downward growth. Temperature fluctuations are observed along the length of the melt column during downward growth. With increasing Rayleigh number, these fluctuations change from none, to cyclic, to time periodic having multiple harmonics, and finally to random. At the higher convective driving force of 350 K temperature inversion, the transverse magnetic field decreased convective levels, strong random temperature fluctuations (flows) becoming smaller and periodic. The maximum field of 0.45 T was unable to completely eliminate convection. For the lower convective driving force of 150 K temperature inversion, the 0.05 T magnetic field decreased flows, and at 0.15 T, the field caused a dramatic decrease in the characteristic frequency of the temperature fluctuations, indicating a change in the nature of the flow, the waveform of the temperature fluctuations changing from sinusoidal to a pulsed wave. Temperature fluctuations and time delays between thermocouples were used to estimate flow velocities. Irrespective of the convection in the bulk melt (ahead of the mushy zone), longitudinal macrosegregation occurs only if the interdendritic melt mixes with the bulk melt.
Convection signatures are inferred in Jupiter's upper troposphere using three tools that are applied to global temperature maps: density maps, three dimensional geostrophy and the mean residual circulation.
The effects of E x B convection on the distribution of plasma parameters in the inner magnetosphere have been examined. Analytical solutions describing density distributions along convective trajectories in the equatorial plane have been found. These solutions suggest the following dependence of plasma concentration n on the magnetic field B along convective trajectories: n proportional to B(exp alpha), where the parameter alpha varies between alpha is less than or = between 4/3 and 2. The alpha = 2 case corresponds to disregarding transport parallel to B. The lower bound alpha = 4/3 describes the situation when parallel transport dominates over convective motions perpendicular to B. A solution has also been obtained describing the effect of convection on diffusive equilibrium profiles in the dayside plasmasphere. In addition, ion temperature variations due to adiabatic effects associated with plasma convection have been analyzed in detail. Convective drifts lead to ion temperature anisotropies with the value and sign of the anisotropy contingent on density and temperature variations, local time, and the location of the convective trajectory with respect to the location of the plasmapause. We have also found that convective motions lead to a substantial exchange of energy between the ionosphere and the plasmasphere through electron heat fluxes.
Thermal convection driven by uniform volumetric energy sources was studied in a horizontal fluid layer bounded from above by a rigid, isothermal surface and from below by a rigid, zero heat-flux surface. The side walls of the fluid domain were assumed to be rigid and perfectly insulating. The computations were formally restricted to two-dimensional laminar convection but were carried out for a range of Rayleigh numbers which spans the regimes of laminar and turbulent flow. The results of the computations consists of streamline and isotherm patterns, horizontally averaged temperature distributions, and horizontally averaged Nusselt numbers at the upper surface. Flow and temperature fields do not exhibit a steady state, but horizontally averaged Nusselt numbers reach limiting, quasi-steady values for all Rayleigh numbers considered. Correlations of the Nusselt number in terms of the Rayleigh and Prandtl numbers were determined.
Consideration is given to the numerical simulation of the solution growth of triglycine sulfate in space, with heat transfer, mass transfer, buoyancy-driven convection, the dependence of solubility on temperature, and finite interface kinetics taken into account. In the bulk solution, weak thermal convection was established quickly. The thermal convection then becomes dominated by solutionally driven motion. The convection due to the steady background g is predicted to influence both the growth rate and the crystal morphology. The time required for buoyancy-driven convection to influence the growth rate depends on the magnitude of g, i.e., more time for smaller g. Under the same cooling condition, the average growth rate increases as the magnitude of g increases. A steady average growth rate of 1.0 mm/d can be obtained when a properly designed cooling rate is applied to the sting. The most uniform growth is obtained when the g vector is aligned such that the solution flows normal towards the center of the crystal surface.
Semiconductor crystals such as Hg(1-x)CD(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. The value of effective diffusion coefficient calculated from fitting of the data to 1D model varies with Hg(1-x)Cd(x)Te growth conditions. This indicates that the growth condition of the Hg(1-x)Cd(x)Te is not purely diffusion controlled. Because of the higher thermal conductivity in the melt than that in the crystal in the growth system, there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart, and Thorpe, Hutt and Soulsby have shown that under such conditon a fluid will become convectively unstable as a result of different diffusitivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositonal segregation in both axial and radial directions in the unidirectionally solidified crystals under various gravitational directions. To reach this goal, we start with a simplified problem to study the effects of thermal-solutal convection on the temperature and solutal fields under various gravitional orientations. We begin by reviewing model governing equations.
The influences of depth-dependent thermal expansivity on the structure of Newtonian thermal convection and lateral thermal anomalies are studied. In an aspect-ratio ten box the effects of the depth-dependent expansivity are to produce large aspect-ratio cells with a few relatively stationary upwellings. The lateral thermal anomalies in the lower mantle are greatly diminished by depth-dependent expansivity. This together with a decrease in the temperature derivative of seismic velocities in the lower mantle would support the low level of lateral heterogeneities found in recent tomographic studies.
Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.