Engineering topics

Duh, J. C.

Publications and source records attributed to Duh, J. C..

8 records

NASA NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗

NASA NTRS1995

Marangoni instability in a liquid sheet

We study the onset of a pure Marangoni convection in a liquid layer with two deformable interfaces in the no-gravity environment. Both oscillatory and stationary instabilities are considered for a wide range of parameters. It is shown that only stationary instability is possible when surface tension at the colder interface is lower than that at the hotter one. Oscillatory instability tends to disappear and to be replaced by the stationary instability with increase of the Prandtl number and decrease of surface tension at the colder interface.

Oron, A.↗

NASA NTRS1994

Marangoni-Benard convection in a low-aspect-ratio liquid layer

This paper presents the results of a study on the variation of the critical Marangoni number (Ma(sub c)) for the onset of Benard convection in a finite liquid layer bounded horizontally as well as from below. A direct-numerical-simulation procedure is devised to determine the Ma(sub c) for aspect ratios (Ar) ranging from 0.8 to 10. The results predict a strong increase of Ma(sub c) as Ar decreases to below 2. A dip of Ma(sub c) occurs between Ar = 1.45 and 1.3, which is accompanied by a pattern transition from a two-cell convection to a unicellular flow. For Ar above 4, the calculated Ma(sub c) shows little change and asymptotically approach a value of 116.15, with Biot number (Bi) equal to 1.

Duh, J. C.↗

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 NTRS1991

Numerical investigation on Benard instability in a finite liquid layer

A numerical procedure for directly simulating the Benard-Marangoni instabilities (B-M-I) in a bounded liquid layer is presented in this paper. The procedure consists of applying a finite amplitude disturbance to the basic static state, and then integrating the Navier-Stokes equations to determine whether the disturbance will die down or will reach a state of finite-strength steady convection. The critical Marangoni number (Mac) for the onset of B-M-I can thus be determined and can be correlated as a function of the aspect ratio (Ar), Prandtl number (Pr), and Rayleigh number (Ra). The Biot number (B) between the liquid and the air is analyzed to approximate the heat transfer condition along the free surface. A 2D calculation is performed to investigate the effect of various initial disturbances, and the Mac is determined for Ar = 2, Ra = 0, and Pr = 0.7. Current results show that disturbances of different nature and amplitude have little effect on Mac. The Mac determined in this study also clearly demonstrates the dominant effect of the sidewalls.

Duh, J. C.↗

NASA NTRS1989

Numerical modeling of enclosure convection

A numerical study on the steady and unsteady natural convection in two-dimensional rectangular enclosures has been performed by a time-accurate ADI finite difference scheme. The study covered a range of Rayleigh numbers between 1000 and 10 to the 7th, aspect ratios between 0.2 and 10.0, and tilt angles between -90 (heating from bottom) and +90 deg (heating from top). Various Prandtl numbers have been studied, but only the results of water (Pr = 7.0) are reported here due to space limitations. The physics revealed, however, includes the convection phenomena and the Rayleigh-Benard stability, as well as the combined mechanism of these two. The onset of secondary cells is determined by using a velocity map, which is simpler and cleaner, instead of a streamline plot. The critical Ra number for the occurrence of these secondary cells is shown to be lower than can be detected by experimental studies. On the Rayleigh-Benard stability part, a second transition from stable single-cell convection to periodic multicellular convection is disclosed.

Duh, J. C.↗

NASA NTRS1989

Numerical analysis of natural convection in liquid droplets by phase change

A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Duh, J. C.↗

NASA NTRS1989

Comparison between pressure gradient method and MAC method on high Re calculation

A cavity flow driven by shear and buoyancy forces is used as a test problem in the application of a nonstaggered pressure gradient (PG) method in solving the two-dimensional incompressible Navier-Stokes equations. Twelve finite differencing schemes are used to solve the cavity flow problem. The schemes consist of various combinations of grid arrangements, upwinding treatments, and conservativeness of convection terms. An artificial source term is introduced, and the solutions are compared with those obtained by the conventional marker-and-cell (MAC) method. The comparisons favor the PG method. Numerical results obtained by the twelve schemes are compared with exact solutions in order to assess the stability and accuracy of each scheme.

Tan, C.-H.↗