Search NASA⌕ Search

Engineering topics

Yao, Minwu

Publications and source records attributed to Yao, Minwu.

9 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Shear Thinning in Xenon

We measured shear thinning, a viscosity decrease ordinarily associated with complex liquids such as molten plastics or ketchup, near the critical point of xenon. The data span a wide range of dimensionless shear rate: the product of the shear rate and the relaxation time of critical fluctuations was greater than 0.001 and was less than 700. As predicted by theory, shear thinning occurred when this product was greater than 1. The measurements were conducted aboard the Space Shuttle Columbia to avoid the density stratification caused by Earth's gravity.

Bergm Robert F.↗
NASA NTRS2008

Shear Thinning Near the Critical Point of Xenon

We measured shear thinning, a viscosity decrease ordinarily associated with complex liquids, near the critical point of xenon. The data span a wide range of reduced shear rate: 10(exp -3) < gamma-dot tau < 700, where gamma-dot tau is the shear rate scaled by the relaxation time tau of critical fluctuations. The measurements had a temperature resolution of 0.01 mK and were conducted in microgravity aboard the Space Shuttle Columbia to avoid the density stratification caused by Earth's gravity. The viscometer measured the drag on a delicate nickel screen as it oscillated in the xenon at amplitudes 3 mu,m < chi (sub 0) >430 mu, and frequencies 1 Hz < omega/2 pi < 5 Hz. To separate shear thinning from other nonlinearities, we computed the ratio of the viscous force on the screen at gamma-dot tau to the force at gamma-dot tau approximates 0: C(sub gamma) is identical with F(chi(sub 0), omega tau, gamma-dot tau )/F)(chi(sub 0, omega tau, 0). At low frequencies, (omega tau)(exp 2) < gamma-dot tau, C(sub gamma) depends only on gamma-dot tau, as predicted by dynamic critical scaling. At high frequencies, (omega tau)(exp 2) > gamma-dot tau, C(sub gamma) depends also on both x(sub 0) and omega. The data were compared with numerical calculations based on the Carreau-Yasuda relation for complex fluids: eta(gamma-dot)/eta(0)=[1+A(sub gamma)|gamma-dot tau|](exp - chi(sub eta)/3+chi(sub eta)), where chi(sub eta) =0.069 is the critical exponent for viscosity and mode-coupling theory predicts A(sub gamma) =0.121. For xenon we find A(sub gamma) =0.137 +/- 0.029, in agreement with the mode coupling value. Remarkably, the xenon data close to the critical temperature T(sub c) were independent of the cooling rate (both above and below T(sub c) and these data were symmetric about T(sub c) to within a temperature scale factor. The scale factors for the magnitude of the oscillator s response differed from those for the oscillator's phase; this suggests that the surface tension of the two-phase domains affected the drag on the screen below T(sub c).

Zimmerli, Gregory A.↗
NASA NTRS2007

Hydrodynamic Force on a Cylinder Oscillating at Low Frequency

The hydrodynamic force on a cylinder oscillating transversely to its axis is a nonlinear function of the displacement amplitude x0. We report measurements and numerical calculations of the force at frequencies low enough that delta > R, where delta is the viscous penetration length and R is the cylinder radius. For small amplitudes, the numerically calculated Fourier transform of the force per unit length, F(sub small), agrees with Stokes' analytical calculation. For larger amplitudes, the force per unit length found by both calculation and measurement is F = F(sub small)C (x(sub 0)/delta,R/delta). The complex function C depends only weakly on R/delta, indicating that x0/delta is more appropriate as a scaling variable than the Keulegan-Carpenter number KC = pi*x(sub 0)/R. The measurements used a torsion oscillator driven at frequencies from 1 to 12 Hz while immersed in dense xenon. The oscillator comprised cylinders with an effective radius of R = 13.4 micron and oscillation amplitudes as large as x(sub 0)/delta = 4 (corresponding to KC as large as 71). The calculations used similar conditions except that the amplitudes were as large as x0/delta = 28.

Berg, Robert F.↗
NASA NTRS1999

Melt Stabilization of PbSnTe in a Magnetic Field

Both the experimental observation and numerical simulation indicate that the Bridgman growth of PbSnTe under the microgravity environment in space is still greatly influenced by buoyancy-induced convection. The application of a magnetic field during the semiconductor growth can dampen the convective flow in the metal-like melt. However, for Bridgman growth of PbSnTe on earth (with either vertical or horizontal configuration), both experimental observation and numerical modeling suggest that even with a strong magnetic furnace (5-Tesla constant axial magnetic field), the convective flow in the melt still cannot be sufficiently suppressed to reach the diffusion-controlled level. In order to completely dampen the buoyancy-induced convection on earth, estimates based on scaling analysis indicate that for common experimental conditions, an extremely high magnetic field is required, far beyond the capacity of the experimental apparatus currently available. Therefore, it is proposed that only the combination of microgravity environment and magnetic damping will produce the desired diffusion-controlled growth state for this particular material. The primary objectives of this study are to provide a quantitative understanding of the complex transport phenomena during solidification of non-dilute binarys, to furnish a numerical tool for furnace design and growth condition optimization, to provide estimates of the required magnetic field strength for low gravity growth, and to assess the role of magnetic damping for space and earth control of the double-diffusive convection. As an integral part of a NASA research program, our numerical simulation supports both the flight and ground-based experiments in an effort to bring together a complete picture of the complex physical phenomena involved in the crystal growth process. For Bridgman growth of PbSnTe under microgravity (with both vertical and horizontal configurations), the simulations suggest that a moderate axial magnetic field of only a few kilo-Gauss in strength could effectively eliminate buoyancy-induced convection in the melt and control solute segregation. Therefore, this work confirms the idea that the combination of microgravity environment and the magnetic damping will indeed be sufficient to produce the desired diffusion-controlled growth state for PbSnTe.

Fripp, Archibald L.↗
NASA NTRS1995

Numerical modeling of Bridgman growth of PbSnTe in a magnetic field

In this work we study heat and mass transport, fluid motion, and solid/liquid phase change in the process of steady Bridgman growth of Pb(.8)Sn(.2)Te (LTT) in an axially-imposed uniform magnetic field under terrestrial and microgravity conditions. In particular, this research is concerned with the interrelationships among segregation, buoyancy-driven convection, and magnetic damping in the LTT melt. The main objectives are to provide a quantitative understanding of the complex transport phenomena during solidification of the nondilute binary of LTT, to provide estimates of the strength of magnetic field required to achieve the desired diffusion-dominated growth, and to assess the role of magnetic damping for space and earth based control of the buoyancy-induced convection. The problem was solved by using FIDAP and numerical results for both vertical and horizontal growth configurations with respect to the acceleration of gravity vector are presented.

Yao, Minwu↗
NASA NTRS1995

Numerical simulation of heat and mass transport during space crystal growth with MEPHISTO

The MEPHISTO space experiments are collaborative United States and French investigations aimed at understanding the fundamentals of crystal growth. Microgravity experiments were conducted aboard the USMP-1 and -2 missions on STS-52 and 62 in October 1992 and March 1994 respectively. MEPHISTO is a French designed and built Bridgman type furnace which uses the Seebeck technique to monitor the solid/liquid interface temperature and Peltier pulsing to mark the location and shape of the solid/liquid interface. In this paper the Bridgman growth of Sn-Bi and Bi-Sn under terrestrial and microgravity conditions is modeled using the finite element code, FIDAP*. The numerical model considers fully coupled heat and mass transport, fluid motion and solid/liquid phase changes in the crystal growth process. The primary goals of this work are: to provide a quantitative study of the thermal buoyancy-induced convection in the melt for the two flight experiments; to compare the vertical and horizontal growth configurations and systematically evaluate the effects of various gravity levels on the solute segregation. Numerical results of the vertical and horizontal Bridgman growth configurations are presented.

Yao, Minwu↗
NASA NTRS1994

Numerical and Experimental Study of Transport Phenomena in Directional Solidification of Succinonitrile

A numerical and experimental study of the growth of succinonitrile (SCN) using a horizontal Bridginan furnace and transparent glass ampoule was conducted. Two experiments were considered: one in which the temperature profile was fixed relative to the ampoule (no-growth case); and a second in which the thermal profile was translated at a constant rate (steady growth case). Measured temperature profiles on the outer surface of the ampoule were used as thermal boundary conditions for the modelling. The apparent heat capacity formulation combined with the variable viscositymeth~ was used to model the phase change in SeN. Both 2-D and 3-D models were studied and numerical solutions obtained using the commercial finite element code, FIDAP1. Comparison of the numerical results to experimental data showed excellent agreement. The complex 3-D shallow-cavity flow in the melt, differences between 2-D and 3-D models, effects of natural convection on the thermal gradient and shape of the solid/liquid interface, and the sensitivity of simulations to specific assumptions, are also discussed.

de Groh, Henry C., III↗
NASA NTRS1992

Asymptotic integration algorithms for first-order ODEs with application to viscoplasticity

When constructing an algorithm for the numerical integration of a differential equation, one must first convert the known ordinary differential equation (ODE), which is defined at a point, into an ordinary difference equation (O(delta)E), which is defined over an interval. Asymptotic, generalized, midpoint, and trapezoidal, O(delta)E algorithms are derived for a nonlinear first order ODE written in the form of a linear ODE. The asymptotic forward (typically underdamped) and backward (typically overdamped) integrators bound these midpoint and trapezoidal integrators, which tend to cancel out unwanted numerical damping by averaging, in some sense, the forward and backward integrations. Viscoplasticity presents itself as a system of nonlinear, coupled first-ordered ODE's that are mathematically stiff, and therefore, difficult to numerically integrate. They are an excellent application for the asymptotic integrators. Considering a general viscoplastic structure, it is demonstrated that one can either integrate the viscoplastic stresses or their associated eigenstrains.

Freed, Alan D.↗