Mean-Field Free-Energy Lattice Boltzmann Method for Liquid-Vapor Interfacial Flows

dc.contributor.authorLi, Shi-Mingen
dc.contributor.committeechairTafti, Danesh K.en
dc.contributor.committeememberPuri, Ishwar K.en
dc.contributor.committeemembervon Spakovsky, Michael R.en
dc.contributor.committeememberPaul, Mark R.en
dc.contributor.committeememberEllis, Michael W.en
dc.contributor.departmentMechanical Engineeringen
dc.date.accessioned2014-03-14T20:18:38Zen
dc.date.adate2007-12-10en
dc.date.available2014-03-14T20:18:38Zen
dc.date.issued2007-11-02en
dc.date.rdate2008-12-10en
dc.date.sdate2007-11-16en
dc.description.abstractThis dissertation includes a theoretical and numerical development to simulate liquid-vapor flows and the applications to microchannels. First, we obtain a consistent non-local pressure equation for simulating liquid-vapor interfacial flows using mean-field free-energy theory. This new pressure equation is shown to be the general form of the classical van der Waals" square-gradient theory. The new equation is implemented in two-dimensional (2D) D2Q7, D2Q9, and three-dimensional (3D) D3Q19 lattice Boltzmann method (LBM). The three LBM models are validated successfully in a number of analytical solutions of liquid-vapor interfacial flows. Second, we have shown that the common bounceback condition in the literature leads to an unphysical velocity at the wall in the presence of surface forces. A few new consistent mass and energy conserving velocity-boundary conditions are developed for D2Q7, D2Q9, and D3Q19 LBM models, respectively. The three LBM models are shown to have the capabilities to successfully simulate different wall wettabilities, the three typical theories or laws for moving contact lines, and liquid-vapor channel flows. Third, proper scaling laws are derived to represent the physical system in the framework of the LBM. For the first time, to the best of the author's knowledge, we obtain a flow regime map for liquid-vapor channel flows with a numerical method. Our flow map is the first flow regime map so far for submicrochannel flows, and also the first iso-thermal flow regime map for CO₂ mini- and micro-channel flows. Our results show that three major flow regimes occur, including dispersed, bubble/plug, and liquid strip flow. The vapor and liquid dispersed flows happen at the two extremities of vapor quality. When vapor quality increases beyond a threshold, bubble/plug patterns appear. The bubble/plug regimes include symmetric and distorted, submerged and non-wetting, single and train bubbles/plugs, and some combination of them. When the Weber number<10, the bubble/plug flow regime turns to a liquid strip pattern at the increased vapor quality of 0.5~0.6. When the Weber number>10, the regime transition occurs around a vapor quality of 0.10~0.20. In fact, when an inertia is large enough to destroy the initial flow pattern, the transition boundary between the bubble and strip regimes depends only on vapor quality and exists between x=0.10 and 0.20. The liquid strip flow regimes include stratified strip, wavy-stratified strip, intermittent strip, liquid lump, and wispy-strip flow. We also find that the liquid-vapor interfaces become distorted at the Weber number of 500~1000, independent of vapor quality. The comparisons of our flow maps with two typical experiments show that the simulations capture the basic and important flow mechanisms for the flow regime transition from the bubble/plug regimes to the strip regimes and from the non-distorted interfaces to the distorted interfaces. Last, our available results show that the flow regimes of both 2D and 3D fall in the same three broad categories with similar subdivisions of the flow regimes, even though the 3D duct produces some specific 3D corner flow patterns. The comparison between 2D and 3D flows shows that the flow map obtained from 2D flows can be generally applied to a 3D situation, with caution, when 3D information is not available. In addition, our 3D study shows that different wettabilities generate different flow regimes. With the complete wetting wall, the flow pattern is the most stable.en
dc.description.degreePh. D.en
dc.identifier.otheretd-11162007-181247en
dc.identifier.sourceurlhttp://scholar.lib.vt.edu/theses/available/etd-11162007-181247/en
dc.identifier.urihttp://hdl.handle.net/10919/29621en
dc.publisherVirginia Techen
dc.relation.haspartPhDDissertationShiMingLi2007.pdfen
dc.rightsIn Copyrighten
dc.rights.urihttp://rightsstatements.org/vocab/InC/1.0/en
dc.subjectinterfacial flowen
dc.subjectliquid-vapor interfaceen
dc.subjectD3Q19en
dc.subjecttwo-phase flowen
dc.subjectD2Q9en
dc.subjectminichannelen
dc.subjectmicrochannelen
dc.subjectD2Q7en
dc.subjectlattice Boltzmann methoden
dc.subjectcapillarityen
dc.subjectcontact angleen
dc.subjectwettabilityen
dc.subjectCO2en
dc.subjectnear critical pointen
dc.subjectphase changeen
dc.subjectflow regimeen
dc.subjectflow regime mapen
dc.titleMean-Field Free-Energy Lattice Boltzmann Method for Liquid-Vapor Interfacial Flowsen
dc.typeDissertationen
thesis.degree.disciplineMechanical Engineeringen
thesis.degree.grantorVirginia Polytechnic Institute and State Universityen
thesis.degree.leveldoctoralen
thesis.degree.namePh. D.en
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
PhDDissertationShiMingLi2007.pdf
Size:
5.39 MB
Format:
Adobe Portable Document Format