Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems

Convective Instability of a Gravity Modulated Fluid Layer with Surface Tension Variation

Gravity modulation of an unbounded fluid layer with surface tension variations along its free surface is investigated. In parameter space of (wavenumber, Marangoni number) modulation has a destabilizing effect on the unmodulated neutral stability curve for large Prandtl number, Pr, and small modulation frequency, Omega, while a stabilizing effect is observed for small Pr and large Omega. As Omega yields infinity, the modulated neutral stability curves approach the unmodulated neutral stability curve. At certain values of Pr and L2 multiple minima are observed and the neutral stability curves become highly distorted. Closed regions of subharmonic instability are also observed. Alternating regions of synchronous and subharmonic instability separated by very thin stable regions are observed in (1/Omega,g(sub 1)) space for the singly diffusive cases. Quasiperiodic behavior in addition to the synchronous and subharmonic responses, are observed for the case of a double diffusive fluid layer. Minimum acceleration amplitudes were observed to closely correspond with a subharmonic response, Lambda(sub im) = Omega/2

Thermal modeling with solid/liquid phase change of the thermal energy storage experiment

A thermal model which simulates combined conduction and phase change characteristics of thermal energy storage (TES) materials is presented. Both the model and results are presented for the purpose of benchmarking the conduction and phase change capabilities of recently developed and unvalidated microgravity TES computer programs. Specifically, operation of TES-1 is simulated. A two-dimensional SINDA85 model of the TES experiment in cylindrical coordinates was constructed. The phase change model accounts for latent heat stored in, or released from, a node undergoing melting and freezing

Scaling analysis applied to the NORVEX code development and thermal energy flight experiment

A scaling analysis is used to study the dominant flow processes that occur in molten phase change material (PCM) under 1 g and microgravity conditions. Results of the scaling analysis are applied to the development of the NORVEX (NASA Oak Ridge Void Experiment) computer program and the preparation of the Thermal Energy Storage (TES) flight experiment. The NORVEX computer program which is being developed to predict melting and freezing with void formation in a 1 g or microgravity environment of the PCM is described. NORVEX predictions are compared with the scaling and similarity results. The approach to be used to validate NORVEX with TES flight data is also discussed. Similarity and scaling show that the inertial terms must be included as part of the momentum equation in either the 1 g or microgravity environment (a creeping flow assumption is invalid). A 10(exp -4) environment was found to be a suitable microgravity environment for the proposed PCM

The Benard problem: A comparison of finite difference and spectral collocation eigen value solutions

The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes