The marker and cell method is an efficient co-volume technique suitable for the solution of incompressible flows using a Cartesian mesh. For flows around complex geometries the use of an unstructured mesh is desirable. For geometric flexibility an unstructured mesh implementation is desirable. A co-volume technique requires a dual orthogonal mesh, in the triangular case the Delaunay-Voronoi dual provides the means for determining this dual orthogonal mesh in an unstructured mesh framework. Certain mesh criteria must be placed on the Delaunay-Voronoi to ensure it meets the dual orthogonal requirements. The two dimensional extension of the marker and cell method to an unstructured framework is presented. The requirements of the mesh are defined and methods in their production are discussed. Initially an explicit time stepping scheme is implemented which allows efficient simulation of incompressible fluid flow problems. Limitations of the explicit time stepping scheme that were discovered, mean that high Reynolds number flows that require the use of stretched meshes cannot produce solutions in a reasonable time period. A semi-implicit time stepping routine removes this limitation allowing these types of flows to be successfully modelled. To validate the solvers accuracy and demonstrate its performance, a number of test cases are presented. These include the lid driven cavity, flow over a backward facing step, inviscid flow around a circular cylinder, unsteady flow around a circular cylinder, flow around an SD7003 aerofoil, flow around a NACA0012 aerofoil and flow around a multi element aerofoil. The investigation although revealing a high dependence on the quality of the mesh still demonstrates that accurate results can be obtained efficiently. The efficiency is demonstrated by comparison to the in-house 2D incompressible finite volume solver for flow around a circular cylinder. For this case the unstructured MAC method produced a solution four times faster than the finite volume code.

No description available.

No description available.

Accurate finite volume solution algorithms for the two dimensional Navier Stokes equations and the three dimensional Euler equations for both structured and unstructured grid topologies are presented. Results for two dimensional quadrilateral and triangular elements and three dimensional tetrahedral elements will be provided. Fundamental to the solution algorithm is a technique for generating multidimensional polynomials which model the spatial variation of the flow variables. Cell averaged data is used to reconstruct pointwise distributions of the dependent variables. The reconstruction errors are evaluated on triangular meshes. The implementation of the algorithm is unique in that three reconstructions are performed for each cell face in the domain. Two of the reconstructions are used to evaluate the inviscid fluxes and correspond to the right and left interface states needed for the solution of a Riemann problem. The third reconstruction is used to evaluate the viscous fluxes. The gradient terms that appear in the viscous fluxes are formed by simply differentiating the polynomial. By selecting the appropriate cell control volumes, centered, upwind and upwind-biased stencils are possible. Numerical calculations in two dimensions include solutions to elliptic boundary value problems, Ringleb’s flow, an inviscid shock reflection, a flat plate boundary layer, and a shock induced separation over a flat plate. Three dimensional results include the ONERA M6 wing. All of the unstructured grids were generated using an advancing front mesh generation procedure. Modifications to the three dimensional grid generator were necessary to discretize the surface grids for bodies with high curvature. In addition, mesh refinement algorithms were implemented to improve the surface grid integrity. Examples studied include a Glasair fuselage, High Speed Civil Transport, and the ONERA M6 wing. The role of reconstruction as applied to adaptive remeshing is discussed and a new first order error estimator is presented. Numerical examples of the remeshing procedure include both smooth and discontinuous flows. / Ph. D.

A Petrov-Galerkin method for the solution of the compressible Euler and Navier-Stokes equations is presented. The method is based on the introduction of an anisotropic balancing diffusion in the local direction of the propogation of the scalar variables. The direction in which the diffusion is added and its magnitude are automatically calculated locally using a criterion that is optimal for one-dimensional transport equations. Algorithms are developed using bilinear quadrilateral and linear triangular elements. The triangular elements are used in conjunction with an adaptive scheme using unstructured meshes. Several applications are presented that show the exceptional stability and accuracy of the method, including the ram accelerator concept for the acceleration of projectiles to ultrahigh velocities. Both two-dimensional and axisymmetric models are employed to evaluate multiple projectile configurations and flow conditions.

For the latter half of the last century renormalization methods have played an important part in tackling problems in fundamental physics and in providing a deeper understanding of systems with many interacting scales or degrees of freedom with strong coupling. The study of turbulence is no exception, and this thesis presents an investigation of renormalization techniques available in the study of the statistical theory of homogeneous and isotropic turbulence. The thesis consists of two parts which assess the two main renormalization approaches available in modeling turbulence. In particular we will be focusing on the renormalization procedures developed by McComb and others. The first part of this thesis will discuss Renormalization Group (RG) approaches to turbulence, with a focus on applications to reduce the degrees of freedom in a large-eddy simulation. The RG methods as applied to classical dynamical systems will be reviewed in the context of the Navier-Stokes equations describing fluid flow. This will be followed by introducing a functional based formalism of a conditional average first introduced by McComb, Roberts and Watt [Phys. Rev A 45, 3507 (1992)] as a tool for averaging out degrees of freedom needed in an RG calculation. This conditional average is then used in a formal RG calculation applied to the Navier-Stokes equations, originally done by McComb and Watt [Phys. Rev. A 46, 4797 (1992)], and later revised by Mc- Comb and Johnston [Physica A 292, 346 (2001)]. A correction to the summing of the time-integral detailed in the latter work is shown to introduce an extra viscous life-time term to the denominator of the increment to the renormalized viscosity and is shown to have a negligible effect in the numerical calculations. We follow this study by outlining some problems with the previous approach. In particular it is shown that a cross-term representing the interaction between high and low wavenumber modes which was neglected in the previous studies on the grounds that it does not contribute to energy dissipation, does in fact contribute significantly. A heuristic method is then put forward to include the effects of this term in the RG calculation. This leads to results which agree qualitatively with numerical calculations of eddy-viscosities. We finish this part of the thesis with an application of the RG method to the modeling of a passive scalar advected by a turbulent velocity field. The second part of this thesis will begin by reviewing Eulerian renormalized perturbation theory attempts in closing the infinite moment hierarchy introduced by averaging the Navier-Stokes equations. This is followed by presenting a new formulation of the local energy transfer theory (LET) of McComb et. al. [J. Fluid Mech. 245, 279 (1992)] which resolves some problems of previous derivations. In particular we show by the introduction of time-ordering that some previous problems with the exponential representation of the correlator can be overcome. Furthermore, we show that the singularity in the LET propagator equation cancels by way of a counter-term. We end this study by introducing a single-time Markovian closure based on LET which, unlike other Markovian closures, does not rely on any arbitrary parameters being introduced in the theory.

A dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, in fulfilment of requirements for the degree of Master of Science. February 16, 2015. / The Adomian decomposition method is an effective procedure for the analytical solution of a wide class of dynamical systems without linearization or weak nonlinearity assumptions, closure approximations, perturbation theory, or restrictive assumptions on stochasticity. Our aim here is to apply the Adomian decomposition method to steady two-dimensional blood flow through a constricted artery in the presence of a uniform transverse magnetic field. Blood flow is the study of measuring blood pressure and determining flow through arteries. Blood flow is assumed to be Newtonian and is governed by the equation of continuity and the momentum balanced equation (which are known as the Navier-Stokes equations). This model is consistent with the principles of ferro-hydrodynamics and magnetohydrodynamics and takes into account both magnetization and electrical conductivity of blood. We apply the Adomian decomposition method to the equations governing blood flow through arteries in the presence of an external transverse magnetic field. The results show that the e ect of a uniform external transverse magnetic field applied to blood flow through arteries favors the physiological condition of blood. The motion of blood in stenosed arteries can be regulated by applying a magnetic field externally and increasing/decreasing the intensity of the applied field.

Chu Shun Yin. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2000. / Includes bibliographical references (leaves 63-65). / Abstracts in English and Chinese. / Acknowledgments --- p.i / Abstract --- p.ii / Introduction --- p.3 / Chapter 1 --- Preliminaries --- p.8 / Chapter 2 --- Singular Integrals --- p.15 / Chapter 2.1 --- Marcinkiewicz Integral --- p.15 / Chapter 2.2 --- Decomposition in cubes of open sets in Rn --- p.17 / Chapter 2.3 --- Interpolation Theorem for Lp --- p.18 / Chapter 2.4 --- Singular Integrals on homogeneous of degree 0 --- p.25 / Chapter 3 --- Solutions to the Euler Equations --- p.36 / Chapter 3.1 --- Existence and Uniqueness of smooth solutions for Euler Equations --- p.36 / Chapter 3.2 --- Rate of Convergence and Decay in Time --- p.43 / Chapter 3.2.1 --- Rate of Convergence --- p.43 / Chapter 3.2.2 --- Lp Decay for Solutions of the Navier-Stokes Equations --- p.46 / Chapter 3.3 --- Weak Solution to the Euler Equations --- p.48 / Chapter 3.3.1 --- Weak Solution to the Velocity Formulation --- p.49 / Chapter 3.3.2 --- Weak Solution to the Vorticity Formulation --- p.52 / Bibliography --- p.63

Chan Man-fai. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2002. / Includes bibliographical references (leaves 40-41). / Abstracts in English and Chinese. / Acknowledgement --- p.i / Abstract --- p.ii / Chapter 1 --- Introduction --- p.1 / Chapter 2 --- Equation of motion and main results --- p.3 / Chapter 3 --- The space K(x) --- p.9 / Chapter 4 --- Proof of the main theorem --- p.17 / Chapter 4.1 --- The passage to the limit as ε →0 --- p.18 / Chapter 4.2 --- The passage to the limit as δ→ 0 --- p.26 / Chapter 4.3 --- Properties of the solution --- p.29 / Chapter 5 --- Conclusion and comments on future works --- p.36 / Appendix --- p.38 / Bibliography --- p.40

Neste trabalho estudamos a existência de soluções de equações do tipo Navier-Stokes com retardo na força externa e no termo n~ao linear. Usando a teoria de semigrupos estudamos a existência de soluções para um problema da forma \'d. SUP. dt u(t) - v\'delta\'u(t) + (F(t, \'u IND.t\'). abla)u(t) + abla p = g(t, \'u IND.t\'), em \'OMEGA\' x (0, T), div u(t) = 0 em \'OMEGA\' x (0, T), u(0, x) = \'u POT.0 (x) x PERTENCE a \' OMEGA\', u(t, x) = 0 t > 0, X \'PERTENCE A\' \' PARTIAL\' \'OMEGA\', u(t, x) =\\phi (t, x) t \'PERTENCE A\' (- \'INFINITO\', 0) x \'PERTENCE A\' \'OMEGA\', onde F9t, \'uIND.t) = INT.IND.t SUP. -\' INFINITO\' \' ALFA1(s-t)u(s)ds + u(t-r), g(t, \'u IND.t\') = INT. SUP. t IND. - INFINITO \'BETA\' (s-t)u(s)ds. Similarmente, usando a tecnica de aproximac~oes de Galerkin, estudamos o problema anterior com F(.) e g(.) dadas por f(t; \'u INDS.t\') = u(t-r(t)); e g(t; \'u IND.t\') = G(u(t-\'rô\' (t))), para alguma função G apropriada. Neste caso, também estudamos a estabilidade de soluções estacionarias / In this work we stuy the existence of solutions for a Navier-Stokes typt equations with delay in the external force and in the nonlinear term. Using the semi-group theory we study the existence of solution for a problem in the form \'d. SUP. dt u(t) - v\'delta\'u(t) + (F(t, \'u IND.t\'). abla)u(t) + abla p = g(t, \'u IND.t\'), ijn \'OMEGA\' x (0, T), div u(t) = 0 in \'OMEGA\' x (0, T), u(0, x) = \'u POT.0 (x) x \'IT BELONGS \' OMEGA\', u(t, x) = 0 t > 0, X \'IT BELONGS\' \'PARTIAL\' \'OMEGA\', u(t, x) =\\phi (t, x) t \'IT BELONGS\' (- \'INFINITY\', 0) x \'IT BELONGS\' \'OMEGA\', where F(t, \'u .t) = INT.IND.t SUP. -\' INFINITY\' \' ALFA(s-t)u(s)ds + u(t-r), g(t, \'u IND.t\') = INT. SUP. t IND. - INFINITY \'BETA\' (s-t)u(s)ds. On another hand using the Galerkin appreoximations method we study the same with F(.) e g(.) given by f(t; \'u INDS.t\') = u(t-r(t)); and g(t; \'u IND.t\') = G(u(t-\'rô\' (t))), for some G appropriated. In thiis case, we study also the stability of stanionary solutions