Direct Computational Simulations for Internal Condensing Flows and Results on Attainability/Stability of Steady Solutions, Their Intrinsic Waviness, and Their Noise Sensitivity

[+] Author and Article Information
A. Narain, Q. Liang, G. Yu, X. Wang

Department of Mechanical Engineering–Engineering Mechanics, Michigan Technological University, 1400 Townsend Drive, Houghton, MI 49931

J. Appl. Mech 71(1), 69-88 (Mar 17, 2004) (20 pages) doi:10.1115/1.1641063 History: Received December 12, 2002; Revised June 09, 2003; Online March 17, 2004
Copyright © 2004 by ASME
Your Session has timed out. Please sign back in to continue.


Liang, Q., Wang, X., and Narain, A., 2004, “Effects of Gravity, Shear, and Surface Tension in Internal Condensing Flows—Results from Direct Computational Simulations,” (accepted for publication in the ASME Journal of Heat Transfer).
Liang, Q., 2003, “Unsteady Computational Simulations and Code Development for a Study of Internal Film Condensation Flows’ Stability, Noise-Sensitivity, and Waviness,” Ph.D. thesis, Michigan Technological University.
Krotiuk, W. J., 1990, Thermal-Hydraulics for Space Power, Propulsion, and Thermal Management System Design, American Institute of Aeronautics and Aeronautics, Washington, DC.
Faghri, A., 1995, Heat Pipe Science and Technology, Taylor and Francis, Washington, DC.
Nusselt,  W., 1916, “Die Oberflächenkondesation des Wasserdampfes,” Z. Ver. Dt. Ing., 60(27), pp. 541–546.
Rohsenow,  W. M., 1956, Heat Transfer and Temperature Distribution in Laminar Film Condensation, Trans. ASME, 78, pp. 1645–1648.
Sparrow,  E. M., and Gregg,  J. L., 1959, “A Boundary Layer Treatment of Laminar Film Condensation,” ASME J. Heat Transfer, 81, pp. 13–18.
Koh,  J. C. Y., Sparrow,  E. M., and Hartnett,  J. P., 1961, “The Two-Phase ry Layer in Laminar Film Condensation,” Int. J. Heat Mass Transfer, 2, pp. 69–82.
Dhir,  V. K., and Lienhard,  J. H., 1971, “Laminar Film Condensation on Plane and Axisymmetric Bodies in Nonuniform Gravity,” ASME J. Heat Transfer, 93, pp. 97–100.
Rose,  J. W., 1988, “Fundamentals of Condensation Heat Transfer: Laminar Film Condensation,” JSME Int. J., 31, pp. 357–375.
Tanasawa,  I., 1991, “Advances in Condensation Heat Transfer,” Adv. Heat Transfer, 21, pp. 55–131.
Cess,  R. D., 1960, “Laminar Film Condensation on a Flat Plate in the Absence of a Body Force,” ZAMP, 11, pp. 426–433.
Koh,  J. C. Y., 1962, “Film Condensation in a Forced Convection Boundary Layer Flow,” Int. J. Heat Mass Transfer, 5, pp. 941–954.
Kutateladze, S. S., 1963, Fundamentals of Heat Transfer, Academic Press, San Diego, CA.
Labuntsov,  D. A., 1957, “Heat Transfer in Film Condensation of Pure Steam on Vertical Surfaces and Horizontal Tubes,” Teploenergetica, 4, p. 72. Also see (in Russian) “Heat Transfer During Condensation of Steam on a Vertical Surface in Conditions of Turbulent Flow of a Condensate Film,” 1960, Inghenerno-Fizicheski Zhurnal, 3, pp. 3–12.
Incropera, F. P., and DeWitt, D. P., 1996, Fundamentals of Heat and Mass Transfer, 4th Ed., John Wiley and Sons, New York.
Yu, G., 1999, “Development of a CFD Code for Computational Simulations and Flow Physics of Annular/Stratified Film Condensation Flows,” Ph.D. thesis, ME-EM Department, Michigan Technological University.
Chow, L. C., and Parish, R. C., 1986, “Condensation Heat Transfer in Microgravity Environment,” Proceedings of the 24th Aerospace Science Meeting, AIAA, New York.
Narain,  A., Yu,  G., and Liu,  Q., 1997, “Interfacial Shear Models and Their Required Asymptotic Form for Annular/Stratified Film Condensation Flows in Inclined Channels and Vertical Pipes,” Int. J. Heat Mass Transfer, 40(15), pp. 3559–3575.
Henstock,  W. H., and Hanratty,  T. J., 1976, “The Interfacial Drag and the Height of the Wall Layer in Annular Flows,” AIChE J., 22, pp. 990–1000.
Lu, Q., 1992, “An Experimental Study of Condensation Heat Transfer With Film Condensation in a Horizontal Rectangular Duct,” Ph.D. thesis, Michigan Technological University.
Sussman,  M., Smereka,  P., and Osher,  S., 1994, “A Level Set Approach for Computing Solutions to Incompressible Two-Phase Flow,” J. Comput. Phys., 114, pp. 146–159.
Son,  G., and Dhir,  V. K., 1998, “Numerical Simulation of Film Boiling Near Critical Pressures with a Level Set Method,” ASME J. Heat Transfer, 120, pp. 183–192.
Hirt,  C. W., and Nichols,  B. D., 1981, “Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries,” J. Comput. Phys., 39, pp. 201–255.
Tezduyar, T. E., 2001, “Finite Element Interface-Tracking and Interface-Capturing Techniques for Flows with Moving Boundaries and Interfaces,” Proceedings of the ASME Symposium on Fluid Physics and Heat Transfer for Macro and Micro-Scale Gas-Liquid and Phase-Change Flows, HTD-Vol. 369-3, ASME, New York.
Abbott, M. B., and Basco, D. R., 1989, Computational Fluid Dynamics—An Introduction For Engineers, Longman, pp. 86–88, Ch. 3, pp. 151–155, Ch. 5.
Li,  J., and Renardy,  Y., 2000, “Numerical Study of Flows of Two Immiscible Liquids at Low Reynolds Number,” SIAM Rev., 42(3), pp. 417–439.
Lu,  Q., and Suryanarayana,  N. V., 1995, “Condensation of a Vapor Flowing Inside a Horizontal Rectangular Duct,” ASME J. Heat Transfer, 117, pp. 418–424.
Bhatt,  B. L., Wedekind,  G. L., and Jung,  K., 1989, “Effects of Two-Phase Pressure Drop on the Self-Sustained Oscillatory Instability in Condensing Flows,” ASME J. Heat Transfer, 111, pp. 538–545.
Liu,  J., and Gollub,  J. P., 1994, “Solitary Wave Dynamics of Film Flows,” Phys. Fluids, 6(5), pp. 1702–1712.
Alekseenko, S. V., Nakoryakov, V. E., and Pokusaev, B. G., 1994, Wave Flow of Liquid Films, Begell House, New York.
Miyara,  A., 2001, “Flow Dynamics and Heat Transfer of Wavy Condensate Film,” ASME J. Heat Transfer, 123, pp. 492–500.
Lighthill, J., 1979, Waves in Fluids, Cambridge University Press, Cambridge, UK.
Traviss,  D. P., Rohsenow,  W. M., and Baron,  A. B., 1973, “Forced Convection Condensation Inside Tubes: A Heat Transfer Equation for Condenser Design,” ASHRAE Trans., 79, Part 1, pp. 157–165.
Shah,  M. M., 1979, “A General Correlation for Heat Transfer during Film Condensation inside Pipes,” Int. J. Heat Mass Transfer, 22, pp. 547–556.
Hewitt, G. F., Shires, G. L., and Polezhaev, Y. V., eds., 1997, International Encyclopedia of Heat and Mass Transfer, CRC Press, Boca Raton, FL.
Carey, V. P., 1992, Liquid-Vapor Phase-Change Phenomena, Series in Chemical and Mechanical Engineering, Hemisphere, Washington, DC.
Palen, J. W., Kistler, R. S., and Frank Y. Z., 1993, “What We Still Don’t Know About Condensation in Tubes,” Condensation and Condenser Design J. Taborek, J. Rose, and I. Tanasawa, eds., United Engineering Trustees, Inc. for Engineering Foundation and ASME, New York, pp. 19–53.
Delhaye,  J. M., 1974, “Jump Conditions and Entropy Sources in Two Phase Systems: Local Instant Formulation,” Int. J. Multiphase Flow, 1, pp. 395–409.
Patankar, S. V., 1980, Numerical Heat Transfer and Fluid Flow, Hemisphere, Washington DC.
Ferziger, J. H., and Peric, M., 1997, Computational Methods for Fluid Dynamics, Springer, Berlin.
Greenberg, M. D., 1978, Foundations of Applied Mathematics, Prentice-Hall, Englewood Cliffs, Jersey.
Joseph, D. D., 1976, Stability of Fluid Motions, Vol. 1, Springer-Verlag, New York, pp. 7–8.
Plesset,  M. S., and Prosperetti,  A., 1976, “Flow of Vapor in a Liquid Enclosure,” J. Fluid Mech., 78(3), pp. 433–444.
ASHRAE Handbook, 1985, Fundamentals SI Edition, American Society of Heating, Refrigeration and Air-Conditioning Engineers, Atlanta, GA.


Grahic Jump Location
Flow geometry for simulations
Grahic Jump Location
The solenoid valve V1 (actively controlled by flow meter F1) and constant heat input Q̇in≈Ṁin⋅hfg(pB) to the boiler fixes inlet pressure p0 and inlet flow rate to the test section. The mass flow rate through pump P is adjusted to a value that matches the corresponding value at F1. A high flow rate of the coolant (water) flow around the test section fixes condensing surface temperature at a nearly uniform value of Tw while it still allows for different heat removed rates for different exit qualities Ze.
Grahic Jump Location
Computational grids for flow simulation. For chosen xui lines, yvj lines in grid A are first generated by points Pi on δ(x,t). Above the “highest” yvj line thus obtained, the remaining yvj lines are independently generated with suitable unequal spacings. Grid B lines at x=xudi are different from xui lines and are used for tracking the interface δ(x,t).
Grahic Jump Location
(a) The liquid domain calculations underneath δshift(x,t) with prescribed values of (u1si,v1si1si) on δshift(x,t) satisfy the shear and pressure condition on δ(x,t). (b) The vapor domain calculations above δ(x,t) with prescribed values of (u2i,v2i2i) on δ(x,t) satisfy ṁVK=ṁEnergy and the requirement of continuity of tangential velocities.
Grahic Jump Location
The above predictions are for vertical channel flows of saturated R-113 vapor. The flow cases are specified in Table 1 with α=90 deg, xe=50 and two different exit conditions, viz. Ze1=0.5 and Ze2=0.38.
Grahic Jump Location
For the flow situations specified in Table 1 with α=90 deg, xe=50, the figure shows the equivalence of specifying exit vapor quality Ze or exit pressure π̄e≡1/(1−δ)∫δ1π2dy to specify exit conditions. It is computationally more convenient to specify exit condition Ze.
Grahic Jump Location
With all remaining flow parameters specified as in Table 1 with α=90 deg, the above figure shows that exit condition specified by the number Ze at a given xe must lie within two well-defined values, viz. Ze|min(xe)≤Ze≤Ze|max(xe). This restriction, presumably, arises from the fact (see Carey 37) that the assumed annular/stratified flows only occur within certain parameter ranges.
Grahic Jump Location
(a) For flow situation specified in Table 1 with α=90 deg and xe=30, the figure depicts two sets of δ(x,t) predictions for t>0. One curve C1 starts at Ze=0.51 at t=0, and tends, as t→∞, to the solution for Ze|Na=0.47. The other curve C2 starts at Ze=0.44 at t=0 and tends, as t→∞, to the same steady Ze|Na solution. (b) For flow situations considered in Fig. 7(a), the above predictions for t>0 starts at t=0 from the same curves C1 and C2 in Fig. 7(a). However, at t>0, there is a condensing surface noise given by v1(x,0,t)=ε⋅ sin(2πx/λ)⋅ sin(2πt /T), with ε=0.3E-6,λ=10, and T=24. As t→∞, the mean part of wavy quasi-steady solutions coincides with the smooth solution, shown in Fig. 8(a) for Ze=Ze|Na=0.47.
Grahic Jump Location
Qualitative nature of the stable, steady/quasi-steady solutions
Grahic Jump Location
For flow situations specified in Table 1 with α=90 deg and xe=20, the above δ(x,t) predictions (Δt=2.5) are for initial data δ(x,0)=δsteady(x)+δ(x,0), where a nonzero disturbance δ(x,0) has been superposed at t=0 on the steady solution δsteady (shown as curve C1 above for t<0). The steady solution corresponds to Ze=Ze|Na=0.5. Here δ(x,0)=0 except in the interval x*<x<x*+10, where x*=3.5 and δ(x,0)=0.5⋅δsteady(x)⋅ sin(2πx/5). It is clear that even this large a disturbance damps at later times.
Grahic Jump Location
(a) For flow situation specified in Table 1 with α=90 deg and xe=40, the characteristics curve C1 denote curves along which infinitesimal initial disturbances naturally propagate on the stable steady solution. Curve C1 denotes characteristics along which finite disturbance arising from forced bottom wall noise actually propagate. On characteristics originating at x=0,δ(0,t)≅0 implies δ≅0. (b) For flow situations defined in Table 1 with α=90 deg and xe=48, the above ū(x,t) predictions for t≥0 are for (i) steady flow with Ze|Na=0.524, (ii) resonant case in Fig. 13, (iii) nonresonant case in Fig. 13, (iv) large initial disturbance of Fig. 10 and (v) small initial disturbance δ(x,0) which is one-fifth of δ(x,0) in Fig. 10. (c) For flow situations defined in Table 1 with α=90 deg and xe=50, the above v̿(t) values are along actual characteristics curves like C1 in Fig. 11(a). The predictions are for (i) the steady and stable flow with Ze|Na=0.578,(ii) the resonant case of Fig. 13, (iii) the nonresonant case of Fig. 13, and (iv) the small initial disturbance case (scaled up and shown in the lower figure) in Fig. 11(c).
Grahic Jump Location
For the above flow situation specified in Table 1 with α=90 deg and xe=50, the steady solutions are obtained for Ze=Ze|min=0.26 and Ze|Na=0.36
Grahic Jump Location
(a) For the flow situations specified in Table 1 with α=90 deg, xe=50 and Ze=0.578; the above δ(x,t) predictions compare the nonresonant noise with a resonant noise of the same amplitude (ε=0.9E-4). The noise is given by: v1(x,0,t)=ε⋅ sin(2πx/λ)⋅ sin(2πt /T), where (i) λ=10 and T=24 for the nonresonant case, and (ii) λ=10 and T=T(x)=λ/ūsteady(x). (b) For the flow situations considered in Fig. 13(a), the above depicts the wall heat flux qw(x,t), in W, at t=25 for the resonant case, and its time-averaged values q̄w(x), in W, for all other cases.
Grahic Jump Location
(a) The above δ(x,t) predictions for t>0 are for the steady solution curve C1=C2 at t=0 and initial noise specified in Fig. 10. The t=0 solutions are obtained on two grids I and II with (ni×nj|L×nj|V)I=(30×30×20) and (ni×nj|L×nj|V)II=(50×50×30). The t>0 solution are shown as curves C1 and C2 and are, respectively, obtained on grids that have: (ni×nj|L×nj|V)I×Δt=(30×30×20)×2.5 and (ni×nj|L×nj|V)II×Δt=(50×50×30)×5. At t>0, the number of grid lines (ni×nj|L×nj|V) changes somewhat from their value at t=0.
Grahic Jump Location
For the vertical plate situation specified in Table 1 with α=90 deg, xe=48 and Lc=0.004 m, Curve 1 is a plot of the analytical solution of δ(x) as given in Nusselt 5. Curve 2 is the computational solution under the Nusselt assumption for stagnant vapor and zero liquid inertia. Curve 3 is the computational solution under the assumptions of stagnant vapor while allowing for liquid inertia. Curve 4 is the computational solution that allows vapor motion and liquid inertia (the vapor/liquid velocity profiles are shown only for this case). Though not shown above, vapor velocity tends to zero as y→∞.
Grahic Jump Location
The above is a plot of natural values of Ze|Na for different xe values for a representative flow situation specified in Table 1 with α=90 deg and xe=48. The “Increased Rein” case just changes Rein to a new value of 1300. The “Increased Ja” case just changes Ja to a new value of 0.0443 (i.e., ΔT=65°C).
Grahic Jump Location
For the flow situations described in Fig. 16 and xe=25.0, the above figure reports the representative wall heat flux values q̄w(x), in W, as a function of x with 0≤x≤xe
Grahic Jump Location
For the flow situations described in Figs. 1617 and xe=30.0, the above figure reports the values of δsteady(x)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In