Steady “Attractors” for Forced Flow of Vapor Condensing Over a Horizontal Plate (Problem in Cess and Koh) and Their Unsteady Responses to Initial Disturbances and Noise

Author(s):  
S. Kulkarni ◽  
A. Narain ◽  
S. Mitra

Accurate steady and unsteady numerical solutions of the full 2-D governing equations – that model the film condensation of saturated vapor flowing over a horizontal plate (the problem of Cess [1] and Koh [2]) – are obtained and new results on the solutions’ unsteady response to disturbances are presented. The computations reveal important features of this classical condensing flow problem. The results highlight the scope and limitations of the well-known similarity solution given by Koh [2]. For the steady problem formulation, the paper discusses the similarities and differences between the solution obtained by solving the full 2-D governing equations and the one obtained semi-analytically by the similarity solution approach of Koh [2]. It is shown that the pressure variations in the vapor domain near the leading edge, though small, are important in deciding condensation dynamics (steady and unsteady) and cannot, in general, be neglected, as is the case with the similarity solution. For this shear driven flow, by considering the unsteady solutions, the paper finds that any initial guess leads to an unsteady solution which is attracted to a long-term steady solution (which is same as the solution as the steady problem). However, the attraction rates gradually diminish with increasing distance from the leading edge and decreasing inlet speed. The steady solutions for this external flow problem are generally found to be stable to initial disturbances at the interface or in the interior of the flow domain. However, since these flows can only be physically realized on suitable finite length portions of the plate, the issue of their stability and sensitivity to exit pressure disturbances and ever-present noise (through exit pressure or bottom plate) is also considered. For example, for the finite domain realization of this problem, it is found that the flows are stable to small initial disturbances to the nearly uniform value of exit pressure. These finite domain realizations of the flow are unique in the sense that they allow different non-uniform steady pressure prescriptions leading to different steady solutions – particularly near the exit zone. As a result, near the exit of a long plate, large unsteadiness is expected due to sensitivity to small exit pressure noise/fluctuations. The exit pressure noise for finite domain realization of these flows is expected because of practical difficulties in maintaining constant uniform pressures at downstream locations of the top and exit boundaries. It is shown that the transverse component of gravity does not affect the solution or its dynamic response except for the expected changes in the nature of hydrostatic pressure variations.

2010 ◽  
Vol 132 (10) ◽  
Author(s):  
S. Kulkarni ◽  
A. Narain ◽  
S. Mitra

Accurate steady and unsteady numerical solutions of the full 2D governing equations—which model the forced film condensation flow of saturated vapor over a semi-infinite horizontal plate (the problem of Cess and Koh)—are obtained over a range of flow parameters. The results presented here are used to better understand the limitations of the well-known similarity solutions given by Koh. It is found that steady/quasisteady filmwise solution exists only if the inlet speed is above a certain threshold value. Above this threshold speed, steady/quasisteady film condensation solutions exist and their film thickness variations are approximately the same as the similarity solution given by Koh. However, these steady solutions differ from the Koh solution regarding pressure variations and associated effects in the leading part of the plate. Besides results based on the solutions of the full steady governing equations, this paper also presents unsteady solutions that characterize the steady solutions’ attainability, stability (response to initial disturbances), and their response to ever-present minuscule noise on the condensing-surface. For this shear-driven flow, the paper finds that if the uniform vapor speed is above a threshold value, an unsteady solution that begins with any reasonable initial-guess is attracted in time to a steady solution. This long time limiting solution is the same—within computational errors—as the solution of the steady problem. The reported unsteady solutions that yield the steady solution in the long time limit also yield “attraction rates” for nonlinear stability analysis of the steady solutions. The attraction rates are found to diminish gradually with increasing distance from the leading edge and with decreasing inlet vapor speed. These steady solutions are generally found to be stable to initial disturbances on the interface as well as in any flow variable in the interior of the flow domain. The results for low vapor speeds below the threshold value indicate that the unsteady solutions exhibit nonexistence of any steady limit of filmwise flow in the aft portion of the solution. Even when a steady solution exists, the flow attainability is also shown to be difficult (because of waviness and other sensitivities) at large downstream distances.


1985 ◽  
Vol 107 (3) ◽  
pp. 687-695 ◽  
Author(s):  
M. di Marzo ◽  
M. J. Casarella

The problem of laminar film condensation of a saturated vapor flowing over a cold horizontal cylinder is investigated. A rigorous formulation of the governing equations for the vapor boundary layer and the condensed liquid film, including both the gravity-driven body forces and the imposed pressure gradient caused by the vapor flow, is presented. A generalized transformation of the governing equations allows a wide range of Froude numbers to be investigated. A unique value of the Froude number is defined which allows a distinction between the gravity-dominated flow (Fr→0) and the forced flow (Fr→∞) and basically defines the overlap region for the two solution domains. Numerical solutions are obtained in the merging flow regions controlled by both driving forces. The effects of density/viscosity ratio at the liquid-vapor interface, Prandtl number, Jakob number, and Froude number on the heat transfer characteristics are presented.


2016 ◽  
Vol 794 ◽  
pp. 595-638 ◽  
Author(s):  
M. J. Woodhouse ◽  
J. C. Phillips ◽  
A. J. Hogg

We model the unsteady evolution of turbulent buoyant plumes following temporal changes to the source conditions. The integral model is derived from radial integration of the governing equations expressing the evolution of mass, axial momentum and buoyancy in the plume. The non-uniform radial profiles of the axial velocity and density deficit in the plume are explicitly captured by shape factors in the integral equations; the commonly assumed top-hat profiles lead to shape factors equal to unity. The resultant model for unsteady plumes is hyperbolic when the momentum shape factor, determined from the radial profile of the mean axial velocity in the plume, differs from unity. The solutions of the model when source conditions are maintained at constant values are shown to retain the form of the well-established steady plume solutions. We demonstrate through a linear stability analysis of these steady solutions that the inclusion of a momentum shape factor in the governing equations that differs from unity leads to a well-posed integral model. Therefore, our model does not exhibit the mathematical pathologies that appear in previously proposed unsteady integral models of turbulent plumes. A stability threshold for the value of the shape factor is also identified, resulting in a range of its values where the amplitudes of small perturbations to the steady solutions decay with distance from the source. The hyperbolic character of the system of equations allows the formation of discontinuities in the fields describing the plume properties during the unsteady evolution, and we compute numerical solutions to illustrate the transient development of a plume following an abrupt change in the source conditions. The adjustment of the plume to the new source conditions occurs through the propagation of a pulse of fluid through the plume. The dynamics of this pulse is described by a similarity solution and, through the construction of this new similarity solution, we identify three regimes in which the evolution of the transient pulse following adjustment of the source qualitatively differs.


Author(s):  
Foad Hassaninejadfarahani ◽  
Scott Ormiston

Laminar film condensation is an important phenomenon which occurs in numerous industrial applications such as refrigeration, chemical processing, and thermal power generation industries. It is well known that film condensation heat transfer is greatly reduced in the presence of a non-condensing gas. The present work performs a numerical analysis of the steady-state, laminar film condensation from a vapour-gas mixture in vertical parallel plate channels to demonstrate a computer model that could assist engineering analysts designing systems involving these phenomena. The present model has three new aspects relative to other current work. First, the complete elliptic two-dimensional governing equations are solved in both phases. Thus, the entire channel domain is solved rather than using an approach that marches along the channel from inlet to a prescribed length. Second, a dynamically determined sharp interface is used between the phases. This sharp interface is determined during the solution on a non-orthogonal structured mesh. Third, the governing equations are solved in a fully-coupled approach. The equations for two velocities, pressure, temperature, and gas mass fraction are solved in a coupled method simultaneously for both phases. Discretisation has been done based on a finite volume method and a co-located variable storage scheme. An in-house computer code was developed to implement the numerical solution scheme. Detailed results are presented for laminar film condensation from steam-air mixtures flowing in vertical parallel-plate channels. The results include velocity and pressure profiles, as well as axial variations of film thickness, Nusselt number and interface gas mass fraction. Detailed comparisons are made with results from a parabolic solution approach.


1971 ◽  
Vol 93 (3) ◽  
pp. 297-304 ◽  
Author(s):  
V. E. Denny ◽  
A. F. Mills ◽  
V. J. Jusionis

An analytical study of the effects of noncondensable gas on laminar film condensation of vapor under going forced flow along a vertical surface is presented. Due to the markedly nonsimilar character of the coupled two-phase-flow problem, the set of parabolic equations governing conservation of momentum, species, and energy in the vapor phase was solved by means of finite-difference methods using a forward marching technique. Interfacial boundary conditions for the numerical solution were extracted from a locally valid Nusselt-type analysis of the liquid-film behavior. Locally variable properties in the liquid were treated by means of the reference-temperature concept, while those in the vapor were treated exactly. Closure of the numerical solution at each step was effected by satisfying overall mass and energy balances on the liquid film. A general computer program for solving the problem has been developed and is applied here to condensation from water-vapor–air mixtures. Heat-transfer results, in the form q/qNu versus x, are reported for vapor velocities in the range 0.1 to 10.0 fps with the mass fraction of air ranging from 0.001 to 0.1. The temperature in the free stream is in the range 100–212 deg F, with overall temperature differences ranging from 5 to 40 deg F. The influence of noncondensable gas is most marked for low vapor velocities and large gas concentrations. The nonsimilar character of the problem is especially evident near x = 0, where the connective behavior of the vapor boundary layer is highly position-dependent.


1995 ◽  
Vol 2 (3/4) ◽  
pp. 186-193 ◽  
Author(s):  
A. Stegner ◽  
V. Zeitlin

Abstract. The problem of the large-scale quasi-geostrophic anticyclonic vortices is studied in the framework of the baratropic rotating shallow- water equations on the β-plane. A systematic approach based on the multiplescale asymptotic expansions is used leading to a hierarchy of governing equations for the large-scale vortices depending on their characteristic size, velocity and a free surface elevation. Among them are the Charney-Obukhov equation, the intermediate geostrophic model equation, the frontal dynamics equation and some new nonlinear quasi-geostrophic equation. We are looking for steady-drifting axisymmetric anticyclonic solutions and find them in a consistent way only in this last equation. These solutions are soliton-like in the sense that the effects of weak non-linearity and dispersion balance each other. The same regimes on the paraboloidal β-plane are studied, all giving a negative result in what concerns the axisymmetric steady solutions, except for a strong elevation case where any circular profile is found to be steadily propagating within the accuracy of the approximation.


2018 ◽  
Vol 02 (04) ◽  
pp. 1850008
Author(s):  
Ting Yi

The attempt of this paper is to suggest a new theoretical medium to support the efforts of building particle models as mechanical field structures. The medium model is based on two simple fundamental assumptions. All the properties of the medium, as counterparts in this medium of the well-known physics laws, including energy and momentum conservations, Lorentz transformation and special relativity, electromagnetic interaction, etc., are derived from these two assumptions. The governing equations are established based on the derived properties. Some steady solutions as well as the feasibility of interpreting these solutions as particles are also briefly discussed.


1988 ◽  
Vol 110 (2) ◽  
pp. 390-394 ◽  
Author(s):  
P. H. Oosthuizen

Two-dimensional flow over a horizontal plate in a saturated porous medium mounted near an impervious adiabatic horizontal surface and subjected to a horizontal forced flow has been numerically investigated. The plate is heated to a uniform temperature that is higher than the temperature of the flowing fluid. The conditions considered are such that the buoyancy forces have an effect on the flow and, therefore, on the heat transfer rate from the plate. The full governing equations, written in dimensionless form, have been solved for a range of values of the governing parameters using the finite element method. The heat transfer rate from the plate is influenced both by the dimensionless depth of the plate below the surface and the importance of the buoyancy forces, the latter having been characterized by a parameter which is equal to the ratio of the Darcy–Rayleigh number to Peclet number. The conditions under which these parameters have a negligible effect on the heat transfer rate are discussed.


1969 ◽  
Vol 39 (1) ◽  
pp. 173-192 ◽  
Author(s):  
Zeev Rotem ◽  
Lutz Claassen

The paper discusses free convective flows above a horizontal plate, both theoretically and on the basis of experiments which yield quantitative data. The theory is applicable to the semi-infinite plate and is extended to cover the complete range of Prandtl number values including Pr → 0 and Pr → ∞. Experiments were carried out to demonstrate the existence of a laminar boundary layer above a horizontal plate at intermediate Grashof (respectively Rayleigh) numbers, and its extent along the plate. This layer breaks down into large-eddy instability some distance from the leading edge. The value of the critical Rayleigh number for this to occur, obtained experimentally using semi-focusing colour-Schlieren photography is in reasonable qualitative agreement with previously known data (Tritton 1963a,b).


2014 ◽  
Vol 30 (4) ◽  
pp. 339-347 ◽  
Author(s):  
S. Munawar ◽  
A. Ali ◽  
N. Saleem ◽  
A. Naqeeb

AbstractIn this work a numerical investigation has been conducted to study the unsteady oscillatory flow of a viscous fluid induced by a swirling disk. The disk stretches radially with the time-based sinusoidal oscillations. The governing equations for the three-dimensional boundary layer-flow are normalized using a suitable set of similarity transformations. The normalized partial differential equations are then solved numerically using a finite difference scheme by altering the semi-infinite domain to a finite domain. The effects of various imperative parameters on the oscillatory flow are discussed with graphs and tables.


Sign in / Sign up

Export Citation Format

Share Document