Browsing by Subject "Fluid dynamics--Mathematical models"
Now showing 1 - 9 of 9
Results Per Page
Sort Options
Item Analysis of a Darcy-Stokes system modeling flow through vuggy porous media(2004) Lehr, Heather Lyn; Arbogast, Todd James, 1957-Our goal is to accurately model flow through subsurface systems composed of vuggy porous media. A vug is a small cavity in a porous medium which is large relative to the intergranular pore size. A vuggy porous medium is a porous medium with vugs scattered throughout it. While the vugs are often small, they can have a tremendous effect on the flow of fluid through the medium. We first introduce our microscale mathematical model for flow of an incompressible, viscous fluid in vuggy porous media. Our next step is to obtain a homogenized macroscale model. In order to do so, we assume periodicity of the medium. We obtain necessary existence and uniqueness results, error estimates, and slight generalizations of two-scale convergence results for bi-modal media. First using formal homogenization and then the rigorous two-scale convergence method, we show that our microscale model homogenizes to give a much simpler modified Darcy’s law macroscale model. In this homogenized model, the permeability tensor is modified to capture the effects of the vugs on the flow through the medium. In order to compute the homogenized permeability tensor, we essentially compute our microscale system on a (much smaller) representative cell. Toward this end, we introduce two numerical methods for the microscale model. We combine a discontinuous Galerkin method with a low order RaviartThomas element and obtain suboptimal convergence rates for the first method. The second method differs only slightly from the first, but yields optimal convergence rates. Unfortunately, it is less efficient in practical implementations.Item A computational fluid dynamics simulation model for flare analysis and control(2006) Castiñeira Areas, David; Edgar, Thomas F.Industrial flares are units designed to safely dispose of waste hydrocarbon gases from chemical and petrochemical plants by burning gases to carbon dioxide and steam, which are then released to the atmosphere. There is still great uncertainty about flare efficiency and the resultant gas emissions under different operating conditions. For this reason, environmental agencies have encouraged the development of predictive models for flare gas combustion systems, so effective control and mitigation strategies can be implemented. The principal focus of this dissertation is to develop mathematical models of industrial flares that predict the efficiency of these industrial combustion systems. For this purpose, a computational fluid dynamics (CFD) simulation model is implemented to analyze the effects of variables such as ambient wind velocity, gas heating value, and steam injection on flare combustion efficiency. Some advanced chemistry and turbulence submodels are also implemented to describe the complex flare flow phenomena. Simulation results show that flares may represent an important source of gas emissions due to inefficient operation at high crosswinds and large steam/fuel ratios. The predictive models presented in this work will allow for better estimation of the resulting gas emissions from industrial plants. Use of these simulation models will also yield economic savings for environmental studies compared to setting up expensive flare experiments. In addition, these predictive models allow for a detailed analysis of species concentration profiles and turbulent flow patterns within the flames, data which is not available experimentally. Furthermore, several instrumentation and control strategies for industrial flares are analyzed in this dissertation. A new approach for flare monitoring based on multivariate image analysis is proposed so that flare combustion efficiency can be measured in real-time.Item A discontinuous Galerkin method for two- and three-dimensional shallow-water equations(2004) Aizinger, Vadym; Dawson, Clinton N.Item Numerical modeling of Stokesian emulsions(2002) Overfelt, James Robert; Rodin, G. J.; Van de Geijn, Robert A.The principle objective of this dissertation was to develop a fast solution method for the accurate large scale modeling of fluid dynamics problems involving emulsions. The objective has been achieved by combining boundary element methods with fast iterative solvers based on the fast multipole method. Our solution method required extending the fast multipole method to generalized electrostatics problems to periodic domains. In addition, we developed an efficient hierarchical time stepping scheme and adopted accurate integration schemes, particularly useful in simulation of dense emulsions. Our solution method was implemented on parallel computers and used to model a creaming experiment of an emulsion.Item A priori prediction of macroscopic properties of sedimentary rocks containing two immiscible fluids(2005) Gladkikh, Mikhail Nikolaevich; Bryant, Steven L.The processes in porous media governed by capillary forces, such as drainage and imbibition cycles, infiltration from surface water, flow, transport, adsorption and dissolution of contaminants, are of great importance in soil science, subsurface geochemistry, petroleum engineering, and hydrology. The methodology proposed in this work is devoted to the pore-scale modeling of such processes. The goal is to be able to make a priori predictions of the macroscopic properties of model sedimentary rocks, such as capillary pressure curves, interfacial area, relative permeabilities and electrical properties. The idea is to conduct a theoretical investigation in a simple but physically representative model porous medium. The model is a random packing of equal spheres for which the coordinates of the centers have been measured. Knowledge of the coordinates determines the grain space and the void space in the packing, thereby overcoming a long-standing difficulty for theoretical approaches to pore-level modeling. Geological processes, such as isopachous and pore-filling cementation, are simulated in the sphere pack, thus creating simple models of sedimentary rocks with predetermined pore space geometry. Despite the simplicity of the model porous medium, it is a powerful tool for investigation of the flow and transport in soils and sedimentary rocks. In particular it allows a priori predictions of macroscopic behavior. This capability is the most important aspect of the approach. Because there are no arbitrarily prescribed parameters, the predictions can be compared directly to experiments, providing a much stronger test than is possible with many other modeling approaches. For example, knowledge of the pore space geometry and wettability conditions allows computing the exact configuration of liquid phases in porous media. This capability enables the simulation of imbibition of a wetting phase into the model porous medium using a physically consistent dynamic criterion for the imbibition of individual pores. This approach allows a priori, quantitative prediction of the configuration of fluid phases during imbibition. This in turn allows a quantitative understanding of how different macroscopic processes and parameters (e.g. relative permeabilities or resistivity index) depend on the geologic (e.g. type and amount of cement) and physical (e.g. wettability) features of porous media.Item Simulating fluid flow in vuggy porous media(2005) Brunson, Dana Sue; Arbogast, Todd James, 1957-Item A three dimensional finite element method and multigrid solver for a Darcy-Stokes system and applications to vuggy porous media(2007-05) San Martin Gomez, Mario, 1968-; Arbogast, Todd James, 1957-Item A three dimensional finite element method and multigrid solver for a Darcy-Stokes system and applications to vuggy porous media(2007) San Martin Gomez, Mario; Arbogast, Todd James, 1957-Item The transport of suspensions in geological, industrial and biomedical applications(2008-12) Oguntade, Babatunde Olufemi; Bonnecaze, R. T. (Roger T.)Suspension flows in varied settings and at different concentrations of particles are studied theoretically using various modeling techniques. Particulate suspension flows are dispersion of particles in a continuous medium and their properties are a consequence of the interplay among hydrodynamic, buoyancy, interparticle and Brownian forces. The applicability of continuum modeling techniques to suspension flows at different particle concentration was assessed by studying systems at different time and length scales. The first two studies involve the use of modeling techniques that are valid in systems where the forces between particles are negligible, which is the case in dilute suspension flows. In the first study, the growth and progradation of deltaic geologic bodies from the sedimentation of particles from dilute turbidity currents is modeled using the shallow water equations or vertically averaged equations of motions coupled with a particle conservation equation. The shallow water model provides a basis for extracting grain size and depositional history information from seismic data. Next, the Navier-Stokes equations of motion and the convection-diffusion equation are used to model suspension flow in a biomedical application involving the flow and reaction of drug laden nanovectors in arteries. Results from this study are then used prescribe the best design parameters for optimal nanovector uptake at the desired sites within an artery. The third study involves the use of macroscopic two phase models to describe concentrated suspension flows where interparticle hydrodynamic forces cannot be neglected. The isotropic form of both the diffusion-flux and the suspension balance models are solved for a buoyant bidisperse pressure-driven flow system. The model predictions are found to compare fairly well with experimental results obtained previously in our laboratory. Finally, the power of discrete type models in connecting macroscopic observations to structural details is demonstrated by studying a system of aggregating colloidal particles via Brownian dynamics. The results from the simulations match experimental shear rheology and also provide a structural explanation for the observed macroscopic behavior of aging.