Dans cette thèse, nous étudions les propriétés des écoulements de fluides qui interagissent avec un corps rigide ou avec une source et un puits. Dans le cas d'un fluide visqueux incompressible qui satisfait les équations de Navier Stokes dans un domaine borné 2D, les solutions faibles de Leray-Hopf sont bien comprises. L'existence et l'unicité sont prouvées. De plus, les solutions sont continues en temps `a valeurs dans L 2 (Omega) et satisfont l’égalité d'énergie classique. Plus récemment, le problème d'un corps rigide en mouvement dans un fluide visqueux incompressible modélisé par les équations de Navier-Stokes couplées aux lois de Newton qui décrivent le mouvement du solide a également été abordé dans le cas où des conditions aux limites sans glissement ont été prescrites. Des résultats analogues concernant les solutions de Leray-Hopf ont également été démontrés dans ce contexte. Dans ce manuscrit, nous étudions le cas de conditions aux limites de Navier-Slip. Dans ce cadre, le résultat d'existence pour le système couplé a été prouvée par G'erard-Varet et Hillairet en 2014. Ici, nous montrons que les solutions sont continues en temps, qu'elles satisfont l’égalité d'énergie et qu’elles sont uniques. De plus, nous montrons un résultat d'existence des solutions faibles dans le cas d'un fluide incompressible visqueux auquel s'ajoute un corps rigide dans le cas où la vitesse du fluide a une partie orthonormale d'énergie infinie.Pour un fluide incompressible non visqueux modélisé par les équations d'Euler dans un domaine borné 2D, le cas où le fluide est autorisé à entrer et à sortir de la frontière a été traité par Judovic qui a introduit certaines conditions limites consistant à prescrire la composante normale de la vitesse et de la vorticité entrante. Dans ce manuscrit, nous considérons un domaine borné qui possède deux trous. L'un d'eux est une source, ce qui signifie que le fluide est autorisé à entrer dans le domaine et l'autre est un puits où le fluide peut sortir. En particulier, nous établissons les équations limites vérifiées par le fluide lorsque la source et le puits se contractent en deux points différents. Le système limite est caractérisé par un point source/puits et un point vortex en chacun des deux points où les trous se sont contractés. / In this thesis, we investigate properties of incompressible flows that interact with a rigid body or a source and a sink. In the case of an incompressible viscous fluid that satisfies the Navier Stokes equations in a 2D bounded domain well-posedness of Leray-Hopf weak solutions is well-understood. Existence and uniqueness are proved. Moreover solutions are continuous in time with values in L 2 (Omega) and they satisfy the energy equality. Recently the problem of a rigid body moving in a viscous incompressible fluid modeled by the Navier-Stokes equations coupled with the Newton laws that prescribe the motion of the solid, was also tackled in the case where the no-slip boundary conditions were imposed. And the correspondent well-posedness result for Leray-Hopf type weak solutions was proved. In this manuscript we consider the case of the Navier-slip boundary conditions. In this setting, the existence result for the coupled system was proved by G'erard-Varet and Hillairet in 2014. Here, we prove that solutions are continuous in time, that they satisfy the energy equality and that they are unique. Moreover we show an existence result for weak solutions of a viscous incompressible fluid plus rigid body system in the case where the fluid velocity has an orthoradial part of infinite energy.For an inviscid incompressible fluid modelled by the Euler equations in a 2D bounded domain, the case where the fluid is allowed to enter and to exit from the boundary was tackled by Judovic who introduced some conditions which consist in prescribing the normal component of the velocity and the entering vorticity. In this manuscript we consider a bounded domain with two holes, one of them is a source which means that the fluid is allowed to enter in the domain and the other is a sink from where the fluid can exit. In particular we find the limiting equations satisfied by the fluid when the source and the sink shrink to two different points. The limiting system is characterized by a point source/sink and a point vortex in each of the two points where the holes shrunk.

This dissertation is concerned with the mathematical modeling of the flow in a porous cylinder with a focus on applications to solid rocket motors. After discussing the historical development and major contributions to the understanding of wall injected flows, we present an inviscid rotational model for solid and hybrid rockets with arbitrary headwall injection. Then, we address the problem of pressure integration and find that for a given divergence free velocity field, unless the vorticity transport equation is identically satisfied, one cannot find an analytic expression for the pressure by direct integration of the Navier-Stokes equations. This is followed by the application of a variational procedure to seek novel solutions with varying levels of kinetic energies. These are found to cover a wide spectrum of admissible motions ranging from purely irrotational to highly rotational fields. Subsequently, a second law analysis as well as an extension of Kelvin's energy theorem to open boundaries are presented to verify and corroborate the variational model. Finally, the focus is shifted to address the problem of laminar viscous flow in a porous cylinder with regressing walls. This is tackled using two different analytical techniques, namely, perturbation and decomposition. Comparisons with numerical Runge--Kutta solutions are also provided for a variety of wall Reynolds numbers and wall regression speeds.

In this thesis study, a least-squares finite element method (LSFEM) based conjugate heat transfer solver was developed. In the mentioned solver, fluid flow and heat transfer computations were performed separately. This means that the calculated velocity values in the flow calculation part were exported to the heat transfer part to be used in the convective part of the energy equation. Incompressible Navier-Stokes equations were used in the flow simulations. In conjugate heat transfer computations, it is required to calculate the heat transfer in both flow field and solid region. In this study, conjugate behavior was accomplished in a fully coupled manner, i.e., energy equation for fluid and solid regions was solved simultaneously and no boundary conditions were defined on the fluid-solid interface. To assure that the developed solver works properly, lid driven cavity flow, backward facing step flow and thermally driven cavity flow problems were simulated in three dimensions and the findings compared well with the available data from the literature. Couette flow and thermally driven cavity flow with conjugate heat transfer in two dimensions were modeled to further validate the solver. Finally, a microchannel conjugate heat transfer problem was simulated. In the flow solution part of the microchannel problem, conservation of mass was not achieved. This problem was expected since the LSFEM has problems related to mass conservation especially in high aspect ratio channels. In order to overcome the mentioned problem, weight of continuity equation was increased by multiplying it with a constant. Weighting worked for the microchannel problem and the mass conservation issue was resolved. Obtained results for microchannel heat transfer problem were in good agreement in general with the previous experimental and numerical works. In the first computations with the solver / quadrilateral and triangular elements for two dimensional problems, hexagonal and tetrahedron elements for three dimensional problems were tried. However, since only the quadrilateral and hexagonal elements gave satisfactory results, they were used in all the above mentioned simulations.

This thesis is concerned with the development of numerical techniques to simulate compressible multi-phase flows, in particular a high-accuracy numerical approach with mesh adaptivity. The Discontinuous Galerkin (DG) method was chosen as the framework for this work for being characterized for its high-order of accuracy -thus low numerical diffusion- and being compatible with mesh adaptivity due to its locality. A DG solver named DiGGIT (Discontinuous Galerkin at the Georgia Institute of Technology) has been developed and several aspects of the method have been studied. The Local Discontinuous Galerkin (LDG) method -an extension of DG for equations with high-order derivatives- was extended to solve multiphase flows using Diffused Interface Methods (DIM). This multi-phase model includes the convection of the volume fraction, which is treated as a Hamilton-Jacobi equation. This is the first study, to the author's knowledge, in which the volume fraction of a DIM is solved using the DG and the LDG methods. The formulation is independent of the Equation of State (EOS) and it can differ for each phase. This allows for a more accurate representation of the different fluids by using cubic EOSs, like the Peng-Robinson and the van der Waals models. Surface tension is modeled with a new numerical technique appropriate for LDG. Spurious oscillations due to surface tension are common to all the capturing schemes, and this new approach presents oscillations comparable in magnitude to the most common schemes. The moment limiter (ML) was generalized for non-uniform grids with hanging nodes that result from adaptive mesh refinement (AMR). The effect of characteristic, primitive, or conservative decomposition in the limiting stage was studied. The characteristic option cannot be used with the ML in multi-dimensions. In general, primitive variable decomposition is a better option than with conservative variables, particularly for multiphase flows, since the former type of decomposition reduces the numerical oscillations at material discontinuities. An additional limiting technique was introduced for DIM to preserve positivity while minimizing the numerical diffusion, which is especially important at the interface. The accuracy-preserving total variation diminishing (AP-TVD) marker for ``troubled-cell' detection, which uses an averaged-derivative basis, was modified to use the Legendre polynomial basis. Given that the latest basis is generally used for DG, the new approach avoids transforming to the averaged-derivative basis, what results in a more efficient technique. Furthermore, a new error estimator was proposed to determine where to refine or coarsen the grid. This estimator was compared against other estimator used in the literature and it showed an improved performance. In order to provide equal order of accuracy in time as in space, the commonly used 3rd-order TVD Runge-Kutta (RK) scheme in the DG method was replaced in some cases by the Spectral Deferred Correction (SDC) technique. High orders in time were shown to only be required when the error in time is significant. For instance, convection-dominated compressible flows require for stability a time step much smaller than is required for accuracy, so in such cases 3rd-order TVD RK resulted to be more efficient than SDC with higher orders. All these new capabilities were included in DiGGIT and have provided a generalized approach capable of solving sub- and super-critical flows at sub- and super-sonic speeds, using a high-order scheme in space and time, and with AMR. Canonical test cases are presented to verify and validate the formulation in one, two, and three dimensions. Finally, the solver is applied to practical applications. Shock-bubble interaction is studied and the effect of the different thermodynamic closures is assessed. Interaction between single-drops and a wall is simulated. Sticking and the onset of splashing are observed. In addition, the solver is used to simulate turbulent flows, where the high-order of accuracy clearly shows its benefits. Finally, the methodology is challenged with the simulation of a liquid jet in cross flow.

This study solves the problem of unsteady free convection from an inclined heated tube both numerically and analytically. The tube is taken to have an elliptic cross-section having a constant heat flux applied to its surface. The surrounding fluid is viscous and incompressible and infinite in extent. The Boussinesq approximation is used to describe the buoyancy force driving the flow. The underlying assumptions made in this work are that the flow remains laminar and two-dimensional for all time. This enables the Navier-Stokes and energy equations to be formulated in terms of the streamfunction, and vorticity. We assume that initially an impulsive heat flux is applied to the surface and that both the tube and surrounding fluid have the same initial temperature. The problem is solved subject to the no-slip and constant heat flux conditions on the surface together with quiescent far-field and initial conditions. An approximate analytical-numerical solution was derived for small times, t and large Grashof numbers, Gr. This was done by expanding the flow variables in a double series in terms of two small parameters and reduces to solving a set of differential equations. The first few terms were solved exactly while the higher-order terms were determined numerically. Flow characteristics presented include average surface temperature plots as well as surface vorticity and surface temperature distributions. The results demonstrate that the approximate analytical-numerical solution is in good agreement with the fully numerical solution for small t and large Gr.

This study solves the problem of unsteady free convection from an inclined heated tube both numerically and analytically. The tube is taken to have an elliptic cross-section having a constant heat flux applied to its surface. The surrounding fluid is viscous and incompressible and infinite in extent. The Boussinesq approximation is used to describe the buoyancy force driving the flow. The underlying assumptions made in this work are that the flow remains laminar and two-dimensional for all time. This enables the Navier-Stokes and energy equations to be formulated in terms of the streamfunction, and vorticity. We assume that initially an impulsive heat flux is applied to the surface and that both the tube and surrounding fluid have the same initial temperature. The problem is solved subject to the no-slip and constant heat flux conditions on the surface together with quiescent far-field and initial conditions. An approximate analytical-numerical solution was derived for small times, t and large Grashof numbers, Gr. This was done by expanding the flow variables in a double series in terms of two small parameters and reduces to solving a set of differential equations. The first few terms were solved exactly while the higher-order terms were determined numerically. Flow characteristics presented include average surface temperature plots as well as surface vorticity and surface temperature distributions. The results demonstrate that the approximate analytical-numerical solution is in good agreement with the fully numerical solution for small t and large Gr.

An axisymmetric, Navier-Stokes finite volume flow solver, which uses Harten, Lax and van Leer (HLL) and Harten, Lax and van Leer&ndash / Contact (HLLC) upwind flux differencing scheme for spatial and uses Runge-Kutta explicit multi-stage time stepping scheme for temporal discretization on unstructured meshe is developed. Developed solver can solve the compressible axisymmetric flow. The spatial accuracy of the solver can be first or second order accurate. Second order accuracy is achieved by piecewise linear reconstruction. Gradients of flow variables required for piecewise linear reconstruction are calculated by Green-Gauss theorem. Baldwin-Lomax turbulent model is used to compute the turbulent viscosity. Approximate Riemann solver of HLL and HLLC implemented in solver are validated by solving a cylindrical explosion case. Also the solver&rsquo / s capability of solving unstructured, multi-zone domain is investigated by this problem. First and second order results of solver are compared by solving the flow over a circular bump. Axisymmetric flow in solid propellant rocket motor is solved in order to validate the axisymmetric feature of solver. Laminar flow over flat plate is solved for viscous terms validation. Turbulent model is studied in the flow over flat plate and flow with mass injection test cases.

A fully automated Cartesian/Quad grid generator and laminar flow solver have been developed for external flows by using C++. After defining the input geometry by nodal points, adaptively refined Cartesian grids are generated automatically. Quadtree data structure is used in order to connect the Cartesian cells to each other. In order to simulate viscous flows, body-fitted quad cells can be generated optionally. Connectivity is provided by cut and split cells such that the intersection points of Cartesian cells are used as the corners of quads at the outmost row. Geometry based adaptation methods for cut, split cells and highly curved regions are applied to the uniform mesh generated around the geometry. After obtaining a sufficient resolution in the domain, the solution is achieved with cellcentered approach by using multistage time stepping scheme. Solution based grid adaptations are carried out during the execution of the program in order to refine the regions with high gradients and obtain sufficient resolution in these regions. Moreover, multigrid technique is implemented to accelerate the convergence time significantly. Some tests are performed in order to verify and validate the accuracy and efficiency of the code for inviscid and laminar flows.

Generating natural phenomena in a virtual world has a number of practical applications. Thanks to the rich and complicated details in the real world, the goal of realistically and efficiently reproducing natural phenomena is well known as an open problem for graphics researchers. In this dissertation, three different issues in modeling liquid animations have been addressed. First, a virtual surface method is proposed to account for surface tension effects and their interactions with solid surfaces in physically based fluid simulation. This allows us to generate various surface tension behaviors in small scale liquid. The second issue that is addressed is how to make small scale fluid simulation more efficient. The proposed solution is a general shallow wave equation model, extended from the original shallow wave equations. By simplifying 3D incompressible fluid dynamics into 2D, small scale liquid can be stably and efficiently simulated over arbitrarily curved surfaces using implicit numerical schemes. The third contribution is a novel hybrid framework that combines image based reconstruction techniques with physically based fluid simulation. While image based methods cannot correctly generate fluid animations alone frame by frame, physics is used as a refinement tool to enforce physical soundness by propagating shape information back and forth in space and time. In this way, water animations can be realistically and faithfully generated from images without error accumulation or stability issues.

Native heart valves with limited functionality are commonly replaced by a bileaflet mechanical heart valve (BMHV). However, despite their widespread use, BMHVs still cause major complications, including hemolysis, platelet activation, and thromboembolic events. These complications are believed to be due to the non-physiologic hemodynamic stresses imposed on blood elements by the hinge flows. Three-dimensional characterization of the hinge flows is therefore crucial to ultimately design BMHVs with lower complication rates. This study aims at simulating the pulsatile 3D hinge flows of various BMHVs placed and estimating the thromboembolic potential associated with each hinge. The hinge and leaflet geometries of clinical BMHVs are reconstructed from micro-computed tomography scans. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Physiologic flow boundary conditions and leaflet motion are extracted from the Fluid-Structure-Interaction simulations of the BMHV bulk flow. The accuracy of the solver is assessed by comparing the results with experimental data. The numerical results are analyzed using a particle tracking approach coupled with existing blood damage models to relate the flow structures to the risk for blood damage. Calculations reveal complex, unsteady, and highly 3D flow fields. Zones of flow stagnation and recirculation, favorable to thrombosis and regions of elevated shear stresses, which may induce platelet activation, are identified throughout the hinge and cardiac cycle. The hinge gap width and, more importantly, the shape of the hinge recess and leaflet are found to impact the flow distribution. Avoiding sharp corners or sudden shape transitions appear as key geometrical design parameters to minimize flow disturbances and thromboembolic potential. The computed flows underscore the need to perform full 3D pulsatile simulations throughout the cardiac cycle to fully capture the complexity and unsteadiness of the hinge flows. Though based only on three different designs, this study provides general guidelines to optimize the hinge design based on hemodynamic performance and thromboembolic potential. The developed framework enables rapid and cost-efficient pre-clinical evaluation of prototype BMHV designs prior to valve manufacturing. Application to a wide range of hinges with varying design parameters will eventually help in determining the optimal hinge design.