Not available / text

Thesis (Ph. D.)--University of Texas at Austin, 2006. / Vita. Includes bibliographical references.

Thesis (Ph. D.)--University of Washington, 2005. / Vita. Includes bibliographical references (p. 214-234).

Currently standard first order single-point upstream weighting methods are employed in reservoir simulation for integrating the essentially hyperbolic system components. These methods introduce both coordinate-line numerical diffusion (even in 1-D) and cross-wind diffusion into the solution that is grid and geometry dependent. These effects are particularly important when steep fronts and shocks are present and for cases where flow is across grid coordinate lines. In this thesis, families of novel edge-based and cell-based truly multidimensional upwind formulations that upwind in the direction of the wave paths in order to minimise crosswind diffusion are presented for hyperbolic conservation laws on structured and unstructured triangular and quadrilateral grids in two dimensions. Higher resolution as well as higher order multidimensional formulations are also developed for general structured and unstructured grids. The schemes are coupled with existing consistent and efficient continuous CVD (MPFA) Darcy flux approximations. They are formulated using an IMPES (Implicit in Pressure Explicit in Saturation) strategy for solving the coupled elliptic (pressure) and hyperbolic (saturation) system of equations governing the multi-phase multi-component flow in porous media. The new methods are compared with single point upstream weighting for two-phase and three-component two-phase flow problems. The tests arc conducted on both structured and unstructured grids and involve full-tensor coefficient velocity fields in homogeneous and heterogeneous domains. The comparisons demonstrate the benefits of multidimensional and higher order multidimensional schemes in terms of improved front resolution together with significant reduction in cross-wind diffusion.

A new implicit and compact optimization-based method is presented for high order derivative calculation for finite-volume numerical method on unstructured meshes. Highorder approaches to gradient calculation are often based on variants of the Least-Squares (L-S) method, an explicit method that requires a stencil large enough to accommodate the necessary variable information to calculate the derivatives. The new scheme proposed here is applicable for an arbitrary order of accuracy (demonstrated here up to 3rd order), and uses just the first level of face neighbors to compute all derivatives, thus reducing stencil size and avoiding stiffness in the calculation matrix. Preliminary results for a static variable field example and solution of a simple scalar transport (advection) equation show that the proposed method is able to deliver numerical accuracy equivalent to (or better than) the nominal order of accuracy for both 2nd and 3rd order schemes in the presence of a smoothly distributed variable field (i.e., in the absence of discontinuities). This new Optimization-based Gradient REconstruction (herein denoted OGRE) scheme produces, for the simple scalar transport test case, lower error and demands less computational time (for a given level of required precision) for a 3rd order scheme when compared to an equivalent L-S approach on a two-dimensional framework. For three-dimensional simulations, where the L-S scheme fails to obtain convergence without the help of limiters, the new scheme obtains stable convergence and also produces lower error solution when compared to a third order MUSCL scheme. Furthermore, spectral analysis of results from the advection equation shows that the new scheme is better able to accurately resolve high wave number modes, which demonstrates its potential to better solve problems presenting a wide spectrum of wavelengths, for example unsteady turbulent flow simulations.

Computational thermo-fluid models of a free surface flow under the dominant radiative transfer have been developed for the design and control of a modern optic fiber drawing process. Although modeling of the fiber drawing process has been of interest for the past three decades, most of the previous studies were limited to low draw speeds and small preforms. Large preforms drawn at high speeds have been used in modern fiber drawing systems to improve production efficiency and reduce cost. Several assumptions commonly made in previous studies have to be relaxed to address the new challenges. In this study, instead of using the Rosseland approximation, the radiative transfer equation (RTE) was solved directly for the radiation fluxes using the finite volume method (FVM). The complete two-dimensional free surface flow was simulated along with the coupling of the radiative transfer. Unlike most of the previous studies that only considered the furnace domain and that assumed the glass velocity at the exit, we included the post-chamber in the computation domain and predicted the fiber solidification location. Furthermore, the mixed convection of the air in the post-chamber was also considered, and was shown to have significant effects on the fiber solidification. On the basis of the computational model, a reduced order model (ROM) was developed for a mixed HŁ /LQG controller designed to regulate the fiber diameter under the effects of disturbances. The ROM was derived on the basis of the computational model. Optimal numerical eigenfunctions were obtained through the Karhunen-Loeve expansion using the computational model. The GalerkinŁ s method was then applied to obtain the state space ROM. The numerical model was shown to be efficient and was verified experimentally. The ROM characterizes the dynamics of the system accurately as compared with the computational model. The simulations using the full computational model showed that the closed-loop system is robust and superior to the open-loop system in the regulation of fiber diameter. The modeling and control methods can be applied to the design optimization and parameter regulation of the high-speed large-preform draw processes as well as other manufacturing processes that involve similar thermal-fluid transports.

A suspension flow model based on the "suspension balance" approach has been developed. This work modifies the model to allow the solution of suspension flows under general flow conditions. This requires the development of a frame-invariant constitutive model for the particle stress which can take into account the spatially-varying local kinematic conditions. The mass and momentum balances for the bulk suspension and particle phase are solved numerically using a finite volume method. The particle stress is based upon the computed rate of strain and the local kinematic conditions. A nonlocal stress contribution corrects the continuum approximation of the particle phase for finite particle size effects. Local kinematic conditions are accounted through the local ratio of rotation to extension in the flow field. The coordinates for the stress definition are the local principal axes of the rate of strain field. The developed model is applied to a range of problems. (i) Axially-developing conduit flows are computed using both the full two-dimensional solution and the more computationally efficient "marching" method. The model predictions are compared to experimental results for cross-stream particle concentration profiles and axial development lengths. (ii) Model predictions are compared to experiments for wide-gap circular Couette flow of a concentrated suspension in a shear-thinning liquid. With minor modification, the suspension flow model predicts the major trends and results observed in this flow. (iii) Comparisons are made to experiments for an axisymmetric contraction-expansion. Model predictions for a two-dimensional planar contraction flow test the influence of model formulation. The variation of the magnitude of an isotropic particle normal stress with local kinematic conditions and anisotropy in the in-plane normal stresses are both explored. The formulation of the particle phase stress is found to have significant effects on the solid fraction and velocity. (iv) Finally, for a rectangular piston-driven flow and an obstructed channel flow, a "computational suspension dynamics" study explores the effect of particle migration on the bulk flow field, system pressure drop and particle phase composition.

Dielectric Barrier Discharge (DBD) excimer lamp is a very attractive source for Ultraviolet (UV) radiation. It has a number of advantages compared to the mercury lamp which is the main lamp used in the industry for UV production. Some of these advantages are instant UV radiation (no warm-up period), narrow UV spectrum, longer life times and simple construction. The DBD UV lamp can be used in number of applications like water disinfection, Plasma Display Panels (PDP) and surface treatment in the semiconductor industry. Yet, the full industrial application of this lamp still faces some problems mainly related to finding the optimum electrical driver to maximize the efficiency of such a lamp. This includes the type of the electrical waveform to generate and the power electronic driver to produce it. In this thesis, firstly a physically based circuit model for the DBD lamp using the Finite Volume Method (FVM) is developed. This model provides the electrical and optical characteristics of the lamp. Using this model the sensitivity of the lamp efficiency to the proposed electrical waveform has been determined. Secondly, the order of this FVM model has been reduced to obtain a model which is used in the design procedure of the proposed driver. Since the DBD lamp has a capacitive nature, a current controlled driver is proposed in this thesis as opposed to most of the published drivers which are voltage controlled drivers. The design of this driver is intended to enhance the electrical to optical efficiency of the lamp and therefore enhancing the overall efficiency of the system. The driver topology permits direct control of the peak lamp current and the operating frequency of the supplied current to the DBD lamp. The width of the current pulses is determined by the transformer magnetizing inductance and the lamp capacitance. Experimental results of the proposed driver connected to a XeCl DBD lamp are presented to validate the performance of the driver and to prove the concept of such a current controlled driver. The proposed driver performance is compared to a voltage source driver which was also implemented. The proposed driver produced higher overall system efficiency but at the expense of a reduction in the driver efficiency as compared to the voltage source driver. The complete system, which consists of the developed FVM based model and the equivalent circuit of the proposed driver, was simulated and the results were compared to the experimental results to validate the accuracy of the developed model for the DBD lamp.

Este trabalho tem como propósito investigar como a variação geométrica de determinados parâmetros envolvidos na construção de uma geometria bifurcada de seção circular, em forma de Y, afeta a resistência ao escoamento, tanto de fluidos newtonianos como não newtonianos. As geometrias estudadas foram construídas utilizando-se o princípio do Design Construtal. Os parâmetros variados foram a relação entre os comprimentos dos dutos pais e filhos, a relação entre os diâmetros dos mesmos dutos, e o ângulo central da estrutura em forma de Y. Para as relações geométricas lineares foram utilizados os valores de 0,5; 0,6; 0,7; 0,8; 0,9 e 1, enquanto para os ângulos, foram utilizados os valores de 155°, 135°, 115°, 95°, 75°, 45°, 25° e 10°. Os fluidos utilizados foram do tipo newtoniano e não newtoniano, dentre estes últimos, foram estudados fluidos pseudoplásticos e dilatantes. O trabalho foi realizado através de simulações numéricas, implementadas com a utilização do software comercial Ansys Fluent, o qual resolve as equações governantes através do método dos volumes finitos. As malhas utilizadas foram do tipo poliédrica. Os resultados indicam que há uma diferença em relação ao que se espera da literatura para as relações entre os diâmetros e os comprimentos. A Lei Hess-Murray indica que estas relações ótimas seriam de 2-1/3 para as relações entre os diâmetros e comprimentos. No presente trabalho, foram determinadas relações entre os diâmetros próximas de 0,6, e entre os comprimentos, iguais a 1. Os ângulos ótimos ficaram localizados no intervalo entre 100° e 135°. / This work aims to investigate how the geometric variation of certain parameters involved in the construction of a bifurcated Y-shaped circular cross-section geometry affects the flow resistance of both Newtonian and non-Newtonian fluids. The geometries studied were constructed using the Constructal Design principle. The parameters were the relationship between the lengths of the daughter and parent ducts, the relationship between the diameters of the same ducts, and the central angle of the Y-shaped structure. For the linear geometric relations, values of 0.5; 0.6; 0.7; 0.8; 0.9 and 1 where used, for the angles, the values of 155 °, 135 °, 115 °, 95°, 75 °, 45 °, 25 ° and 10 ° were used. The fluids used were of the Newtonian and non-Newtonian type, among the latter, pseudo plastic and dilatant fluids were studied. The work was carried out through numerical simulations, implemented with the commercial software Ansys Fluent, which solves the governing equations through the finite volume method. The meshes used were of the polyhedral type. The results indicate that there is a difference in relation to what is expected from the literature for the relationships between diameters and lengths. The Hess-Murray Law indicates that these optimal relations would be 2-1/3 for the relationships between diameters and lengths. In the present work, relationships between the diameters close to 0,6 were found and s equal to 1 between the lengths. The optimum angles were located in the range between 100 ° and 135 °.

This work investigates numerically the process of Teflon ablation using a finite-volume discretization, implicit time integration and a domain decomposition method in three-dimensions. The interest in Teflon stems from its use in Pulsed Plasma Thrusters and in thermal protection systems for reentry vehicles. The ablation of Teflon is a complex process that involves phase transition, a receding external boundary where the heat flux is applied, an interface between a crystalline and amorphous (gel) phase and a depolymerization reaction which happens on and beneath the ablating surface. The mathematical model used in this work is based on a two-phase model that accounts for the amorphous and crystalline phases as well as the depolymerization of Teflon in the form of an Arrhenius reaction equation. The model accounts also for temperature-dependent material properties, for unsteady heat inputs and boundary conditions in 3D. The model is implemented in 3D domains of arbitrary geometry with a finite volume discretization on unstructured grids. The numerical solution of the transient reaction-diffusion equation coupled with the Arrhenius-based ablation model advances in time using implicit Crank-Nicolson scheme. For each time step the implicit time advancing is decomposed into multiple sub-problems by a domain decomposition method. Each of the sub-problems is solved in parallel by Newton-Krylov non-linear solver. After each implicit time-advancing step, the rate of ablation and the fraction of depolymerized material are updated explicitly with the Arrhenius-based ablation model. After the computation, the surface of ablation front and the melting surface are recovered from the scalar field of fraction of depolymerized material and the fraction of melted material by post-processing. The code is verified against analytical solutions for the heat diffusion problem and the Stefan problem. The code is validated against experimental data of Teflon ablation. The verification and validation demonstrates the ability of the numerical method in simulating three dimensional ablation of Teflon.