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.

140 records · Page 8

A numerical study of transient heat and mass transfer in crystal growth

A numerical analysis of transient heat and solute transport across a rectangular cavity is performed. Five nonlinear partial differential equations which govern the conservation of mass, momentum, energy and solute concentration related to crystal growth in solution, are simultaneously integrated by a numerical method based on the SIMPLE algorithm. Numerical results showed that the flow, temperature and solute fields are dependent on thermal and solutal Grashoff number, Prandtl number, Schmidt number and aspect ratio. The average Nusselt and Sherwood numbers evaluated at the center of the cavity decrease markedly when the solutal buoyancy force acts in the opposite direction to the thermal buoyancy force. When the solutal and thermal buoyancy forces act in the same direction, however, Sherwood number increases significantly and yet Nusselt number decreases. Overall effects of convection on the crystal growth are seen to be an enhancement of growth rate as expected but with highly nonuniform spatial growth variations.

Han, Samuel Bang-Moo↗

Fluid Flow in Crystal Growth: Analysis of the Verticle Bridgman and Floating Zone Process

Research is aimed at the fundamental understanding of the interactions of heat and mass transport and fluid mechanics in crystal growth from the melt, especially in how it pertains to experiments supported by the Microgravity Science and Applications program. Emphasis has been on studies of the vertical Bridgman system, on small-scale floating zone systems proposed for space flight, and on quantitative prediction of nonlinear transitions in melt/solid interface morphologies leading to the formation of dendrites in directional solidification. The analysis of vertical Bridgman growth has lead to the complete numerical analysis of transport phenomena in directional solidification and to a detailed comparison of calculations to the GaGe growth experiments of Wang and Witt. Calculations show the central role of the ampoule material in determining the radial temperature gradients in the melt and thus, the intensity of convection. Results indicate that diffusion-controlled growth can not be achieved on Earth with conventional growth systems, but are feasible in microgravity

Brown, R. A.↗

Long-Wavelength Rupturing Instability in Surface-Tension-Driven Benard Convection

A liquid layer with a free upper surface and heated from below is subject to thermocapillary-induced convective instabilities. We use very thin liquid layers (0.01 cm) to significantly reduce buoyancy effects and simulate Marangoni convection in microgravity. We observe thermocapillary-driven convection in two qualitatively different modes, short-wavelength Benard hexagonal convection cells and a long-wavelength interfacial rupturing mode. We focus on the long-wavelength mode and present experimental observations and theoretical analyses of the long-wavelength instability. Depending on the depths and thermal conductivities of the liquid and the gas above it, the interface can rupture downwards and form a dry spot or rupture upwards and form a high spot. Linear stability theory gives good agreement to the experimental measurements of onset as long as sidewall effects are taken into account. Nonlinear theory correctly predicts the subcritical nature of the bifurcation and the selection between the dry spot and high spots.

Swift, J. B.↗

LOFTID Heat Flux Gauge Calibration: What is Truth?

The Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) is a demonstration of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology, which may enable the delivery of heavy payloads to Mars, Venus, and Titan, as well as return to Earth. Unlike rigid aeroshells that are constrained by the size of the rocket’s shroud, inflatable aeroshells can be deployed to a much larger scale, thus allowing a spacecraft to begin its deceleration earlier and experience less heating. On LOFTID, there will be 4 total heat flux gauges (HFG) with a range of 70 W/cm2 and 1 radiometer with a range of 3 W/cm2, arranged as shown in Fig. 1. Both the radiometer and total HFGs are Schmidt-Boelter gauges purchased from an external vendor. Radiative calibrations were performed in-house at NASA Ames’ Sensors and TPS Advanced Research Laboratories (STAR Labs) before and after environmental testing to investigate how the testing affected the sensors' response. Additional rounds of radiative calibration at STAR Labs were also performed in order to investigate the large uncertainties associated with these tests. For example, a survey of multiple calibration facilities concluded that the uncertainty within a given facility was +/-3% [1]. An additional NIST study that calibrated heat flux gauges at 7 different facilities also found the variation in calibration coefficients to be up to ~3% within a given facility, but up to 15% between facilities, suggesting systematic differences between test setups [2]. Finally, the response of heat flux gauges to radiative versus convective heat flux has shown to differ by up to 20% [3],[4]. Because the heat flux gauges on LOFTID will predominantly experience convective heat flux during flight, a convective calibration study was performed at Boeing's Large-Core Arc Tunnel (LCAT) facility. Radiative Calibration Procedure The calibrations performed at STAR Labs utilize a quartz lamp bank (QLB) that provides a maximum heat flux of 50 W/cm2, which bounds the expected LOFTID flight environment. The calibration involves exposing a water-cooled Gardon gauge (reference) and then the unit-under-test (UUT) to 5 different heat fluxes multiple times for 10 seconds each, and then calculating a linear fit. The test setup is shown in Fig. 2. The total HFGs were calibrated at STAR Labs 3 times, denoted as STAR 1 (before environmental testing), STAR 2 (after protoflight vibration and thermal-vacuum testing), and STAR 3 (no change from previous test). All 8 flight-lot total HFGs showed a decrease in full-scale output from STAR 1 to STAR 2 by between 0.5% and 10. The first portion of this investigation was to determine whether the change could be due to differences in temperature between the two calibration runs. A typical linear fit to the calibration data was performed using Eq. 1 where q’ is the heat flux in W/cm2, c is the calibration coefficient, and mV is the sensor output. To account for temperature, the data were fit to a nonlinear function that included both the sensor output (mV) and the temperature from the thermocouple embedded inside the HFG near the surface (T): q'=mV/(c1* T + c0}. The residuals between the fits and the actual data points were calculated for every point, and proven to be much smaller for the temperature-compensated fits than for the linear fits for all sensors. An example is shown in Fig. 3. When the temperature-compensated fits from STAR 1 were applied to the STAR 2 data, the residuals did not improve, suggesting that the change in sensitivity between these two calibration runs was not due to temperature. A third round of calibration (STAR 3) was conducted to further address the temperature dependence of the total HFGs, and the resulting sensitivities matched closely to STAR 2 (within 2%). Temperature-compensated calibration curves were once again fit to the data. In this case, when the temperature-compensated fits from STAR 3 were applied to STAR 2 data, the residuals between the fits and STAR 2 data were much lower than the residuals due to the linear fits. This suggests that the changes seen between STAR 1 and STAR 2 were likely due to actual changes in the sensors caused by the environmental testing between the two calibrations. A modification of the original calibration process, in which the UUT was exposed to each heat flux for just 3 seconds (instead of 10) to reduce the temperature increase during the test, was additionally performed on several of the HFGs. In general, the sensitivities were 1-1.5% lower than from the 10-second tests, but the temperatures were also significantly lower. When the temperature-compensated fits from the 10-second tests were applied to the 3-second test data, the residuals were greatly improved than when just using the linear fits, further suggesting that the temperature-compensated fits may lead to better accuracy than the linear fits in flight. Convective Calibration The second portion of this study was to create a mapping between the radiative and convective calibration coefficients. The majority of the heating during flight will be convective, so it is important to understand how the HFG response differs under these conditions. However, there are no standardized methods for convective calibration [5]. Because the TPS aerothermal response models were validated at LCAT, the same facility was chosen for convective calibration of two of the total HFGs (Fig. 4). Preliminary results showed that the full-scale output was 3% and 8% higher in convective heat flux as compared to radiative heat flux. However, tunnel variation may have contributed to noise and uncertainty in the measurements, and more testing and analysis remains to be done. Scope of Presentation The presentation will include an overview of the changes seen in HFG calibration before and after environmental testing, differences between radiative and convective calibrations, the modeling work done to aid in understanding the sensor response to varying environments, and recommended future work.

H S Alpert↗

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

Pressure- and buoyancy-driven thermal convection in a rectangular enclosure

Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.

Spradley, L. W.↗

Effects of Convection During the Photodeposition of Polydiacetylene Thin Films

In this work, we describe a preliminary investigation of buoyancy-driven heat transfer during the growth of thin films from solution following exposure to ultraviolet (UV) light. Irradiation of the growth cell occurs at various directions relative to gravitational acceleration. Through numerical computations, the steady-state flow and temperature profiles are simulated during the course of light exposure. Light-induced polymerization accompanies a heat transfer process through a fairly complicated recirculating flow pattern. A scaling analysis shows that buoyancy-driven velocities only reduce by a factor of 10 for gravity levels as low as 10(exp -2)g(sub 0). Paley et al. observe what appears to be gravitationally sensitive particle development and inclusion in thin films using a photodeposition process. From this study it is clear that production of homogeneous thin films would have to occur in the environment of a complicated flow pattern of recirculation with a nonuniform temperature distribution. Indeed, even when irradiation occurs from the top of the cell, the most stable stratified cell orientation, defects remain in our films due to the persistence of buoyancy-driven convection. To achieve homogeneity, minimal scattering centers, and possible molecular order, photodeposition of polymer films by UV light exposure must proceed in a reduced-convection environment. Fluid mechanics simulations are useful for establishing gravitational sensitivity to this recently discovered process (patent # 5,451,433) for preparing thin films having quite promising nonlinear optical characteristics.

Frazier, D. O.↗

Effects of Convection during the Photodeposition of Polydiacetylene Thin Films

In this work, we describe a preliminary investigation of buoyancy-driven heat transfer during the growth of thin films from solution following exposure to ultraviolet (UV) light. Irradiation of the growth cell occurs at various directions relative to gravitational acceleration. Through numerical computations, the steady-state flow and temperature profiles are simulated during the course of light exposure. Light-induced polymerization accompanies a heat transfer process through a fairly complicated recirculating flow pattern. A scaling analysis shows that buoyancy-driven velocities only reduce by a factor of 10 for gravity levels as low as 10(exp -2) g(sub 0). Paley et al. observe what appears to be gravitationally sensitive particle development and inclusion in thin films using a photodeposition process. From this study, it is clear that production of homogeneous thin films would have to occur in the environment of a complicated flow pattern of recirculation with a nonuniform temperature distribution. Indeed, even when irradiation occurs from the top of the cell, the most stable stratified cell orientation, defects remain in our films due to the persistence of buoyancy-driven convection. To achieve homogeneity, minimal scattering centers, and possible molecular order, photodeposition of polymer films by UV light exposure must proceed in a reduced-convection environment. Fluid mechanics simulations are useful for establishing gravitational sensitivity to this recently discovered process (patent # 5,451,433) for preparing thin films having quite promising nonlinear optical characteristics.

Frazier, D. O.↗

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗

Buoyancy-driven flows of a radiatively participating fluid in a vertical cylinder heated from below

The effect of radiative energy transport on the onset and evolution of natural convective flows is studied in a Rayleigh-Benard system. Steady, axisymmetric flows of a radiatively participating fluid contained in a rigid-walled, vertical cylinder which is heated on the base, cooled on top, and insulated on the side wall are calculated by using the Galerkin FEM. Bifurcation analysis techniques are used to investigate the changes in the flow structure due to internal radiation. The results of this two-parameter study - where the Rayleigh number, Ra and optical thickness, tau, are varied - apply to fluids ranging from opaque to nearly transparent with respect to IR radiation. For any nonopaque fluid, internal radiation eliminates the static state that without radiation, exists for all values of the Rayleigh number. This heat transfer mechanism also destroys a symmetry of the system that relates clockwise and counterclockwise flows. The connectivity between characteristic flow families and the range of Ra where families are stable are found to depend greatly on tau. Results demonstrate the inadequacy of characterizing the behavior of this system using simple notions of radiative transfer in optically thick or thin media; the nonlinear interaction of radiation and flow are far more complicated than these asymptotic limits would imply.

Salinger, A. G.↗

Gravitational Effects on the Morphology and Kinetics of Photodeposition of Polydiacetylene Thin Films From Monomer Solutions

The goal of this proposed work is to study gravitational effects on the photodeposition of polydiacetylene thin films from monomer solutions onto transparent substrates. Polydiacetylenes have been an extensively studied class of organic polymers because they exhibit many unusual and interesting properties, including electrical conductivity and optical nonlinearity. Their long polymeric chains render polydiacetylenes readily conducive to thin film formation, which is necessary for many applications. These applications require thin polydiacetylene films possessing uniform thicknesses, high purity, minimal inhomogeneities and defects (such as scattering centers), etc. Also, understanding and controlling the microstructure and morphology of the films is important for optimizing their electronic and optical properties. The lack of techniques for processing polydiacetylenes into such films has been the primary limitation to their commercial use. We have recently discovered a novel method for the formation of polydiacetylene thin films using photo-deposition from monomer solutions onto transparent substrates with UV light. This technique is very simple to carry out, and can yield films with superior quality to those produced by conventional methods. Furthermore, these films exhibit good third-order properties and are capable of waveguiding. We have been actively studying the chemistry of diacetylene polymerization in solution and the photo-deposition of polydiacetylene thin films from solution. It is well-known that gravitational factors such as buoyancy-driven convection and sedimentation can affect chemical and mass transport processes in solution. One important aspect of polydiacetylene thin film photodeposition in solution, relevant to microgravity science, is that heat generated by absorption of UV radiation induces thermal density gradients that under the influence of gravity, can cause fluid flows (buoyancy-driven convection). Additionally, changes in the chemical composition of the solution during polymerization may cause solutal convection. These fluid flows affect transport of material to and from the film surface and thereby affect the kinetics of the growth process. This manifests itself in the morphology of the resulting films; films grown under the influence of convection tend to have less uniform thicknesses, and can possess greater inhomogeneities and defects. Specifically, polydiacetylene films photodeposited from solution, when viewed under a microscope, exhibit very small particles of solid polymer which get transported by convection from the bulk solution to the surface of the growing film and become embedded. Even when carried out under conditions designed to minimize unstable density gradients (i.e., irradiating the solution from the top), some fluid flow still takes place (particles remain present in the films). It is also possible that defect nucleation may be occurring within the films or on the surface of the substrate; this, too, can be affected by convection (as is the case with crystal growth). Hence films grown in 1-g will, at best, still possess some defects. The objective of this proposal is to investigate, both in 1-g and in low-g, the effects of gravitational factors (primarily convection) on the dynamics of these processes, and on the quality, morphology, and properties of the films obtained.

Paley, Mark S.↗

Material processing of convection-driven flow field and temperature distribution under oblique gravity

A set of mathematical formulation is adopted to study vapor deposition from source materials driven by heat transfer process under normal and oblique directions of gravitational acceleration with extremely low pressure environment of 10(exp -2) mm Hg. A series of time animation of the initiation and development of flow and temperature profiles during the course of vapor deposition has been obtained through the numerical computation. Computations show that the process of vapor deposition has been accomplished by the transfer of vapor through a fairly complicated flow pattern of recirculation under normal direction gravitational acceleration. It is obvious that there is no way to produce a homogeneous thin crystalline films with fine grains under such a complicated flow pattern of recirculation with a non-uniform temperature distribution under normal direction gravitational acceleration. There is no vapor deposition due to a stably stratified medium without convection for reverse normal direction gravitational acceleration. Vapor deposition under oblique direction gravitational acceleration introduces a reduced gravitational acceleration in vertical direction which is favorable to produce a homogeneous thin crystalline films. However, oblique direction gravitational acceleration also induces an unfavorable gravitational acceleration along horizontal direction which is responsible to initiate a complicated flow pattern of recirculation. In other words, it is necessary to carry out vapor deposition under a reduced gravity in the future space shuttle experiments with extremely low pressure environment to process vapor deposition with a homogeneous crystalline films with fine grains. Fluid mechanics simulation can be used as a tool to suggest most optimistic way of experiment with best setup to achieve the goal of processing best nonlinear optical materials.

Hung, R. J.↗

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗

Investigating the Response of Land-Atmosphere Interactions and Feedbacks to Spatial Representation of Irrigation in a Coupled Modeling Framework

The transport of water, heat, and momentum from the surface to the atmosphere is dependent in part on the 10 characteristics of the land surface. Together with the model physics, parameterization schemes, and parameters employed, land datasets determine the spatial variability in land surface states (i.e., soil moisture and temperature) and fluxes. Despite the importance of these datasets, they are often chosen out of convenience or regional limitations without due assessment of their impacts on model results. Irrigation is an anthropogenic form of land heterogeneity that has been shown to alter the land surface energy balance, ambient weather, and local circulations. As such, irrigation schemes are becoming more 15 prevalent in weather and climate models with rapid developments in dataset availability and parameterization scheme complexity. Thus, to address pragmatic issues related to modeling irrigation, this study uses a high-resolution, regional coupled modeling system to investigate the impacts of irrigation dataset selection on land-atmosphere (L-A) coupling using a case study from the Great Plains Irrigation Experiment (GRAINEX) field campaign. The simulations are assessed in the context of irrigated versus non-irrigated regions, subregions across the irrigation gradient, and sub-grid scale process 20 representation in coarser scale models. The results show that L-A coupling is sensitive to the choice of irrigation dataset and resolution and that the irrigation impact on surface fluxes and near surface meteorology can be dominant, conditioned on the details of the irrigation map (i.e., boundaries, heterogeneity, etc), or minimal. A consistent finding across several analyses was that even a low percentage of irrigation fraction (i.e., 4-16%) can have significant local and downstream atmospheric impacts (e.g., lower PBL height), suggesting that representation of boundaries and heterogeneous areas within irrigated 25 regions is particularly important for the modeling of irrigation impacts on the atmosphere in this model. When viewing the simulations presented here as a proxy for ‘ideal’ tiling in a Earth System Model scale gridbox, the results show that some ‘tiles’ will reach critical nonlinear moisture and planetary boundary layer (PBL) thresholds that could be important for clouds and convection, implying that heterogeneity resulting from irrigation should be taken into consideration in new sub-grid land-atmosphere exchange parameterizations.

Patricia Lawston-Parker↗