Dissertations / Theses: 'Stokes Law (Fluid mechanics)'

7

Mushtaq, Jawaid. "Aspects of implementation and performance evaluation of a parallel multiblock Navier-Stokes solver on a multiprocessor network." Thesis, University of Hertfordshire, 1994. http://ethos.bl.uk/OrderDetails.do?uin=uk.bl.ethos.387191.

Full text

APA, Harvard, Vancouver, ISO, and other styles

8

Botte, Vincenzo. "A robust and accurate Navier-Stokes algorithm for three-dimensional applications adopting arbitrary modelling of the Reynolds stresses." Thesis, Cranfield University, 1999. http://hdl.handle.net/1826/3609.

Full text

Abstract:

In this thesis a new Navier-Stokes solver for complex three-dimensional geometries, adopting arbitrary modelling of the Reynolds stresses, is presented. Ak-s model, adopting a modelling of the turbulent transport not based on the eddy viscosity, has been written in generalised coordinates and solved with a finite volume approach, using both a GMRES solver and a direct solver for the solution of the linear systems of equations. The results presented show that the modification adopted for the modelling of'the turbulent transport also provides a more accurate value of the physical diffusion and, as a consequence, improves the increase in accuracy when using higher-order convection schemes. A simple non-linear modelling of the Reynolds stresses has been designed introducing an additional term, quadratic in the main strain rate, to the basic Boussinesq's form; the corresponding constant has been evaluated through comparison with experimental data. The computational procedure is implemented for the flow analysis in a 90 ° square section bend and the obtained results show that with the non- linear modelling a much better agreement with the measured data is obtained, both for the velocity and the pressure. The importance of the convection scheme is also discussed, showing how the effect of the non-linear correction added to the Reynolds stresses is effectively hidden by the additional numerical diffusion introduced by a low- order convection scheme as the first-order Upwind, thus making necessary the use of higher-order schemes. Some results for centrifugal turbo machinery are also presented, giving some initial indications on the effects of the proposed modification in the modelling of the turbulent diffusion on the prediction of the flow in rotating passages.

APA, Harvard, Vancouver, ISO, and other styles

9

Langeard, Olivier. "Numerical study of a Navier-Stokes flow through a fibrous porous medium." Thesis, Georgia Institute of Technology, 2002. http://hdl.handle.net/1853/15944.

Full text

APA, Harvard, Vancouver, ISO, and other styles

10

Lu, Dongqiang, and 盧東強. "Unsteady free-surface waves generated by bodies in a viscous fluid." Thesis, The University of Hong Kong (Pokfulam, Hong Kong), 2002. http://hub.hku.hk/bib/B29750787.

Full text

APA, Harvard, Vancouver, ISO, and other styles

11

Strodtbeck, Joshua. "A FILTER-FORCING TURBULENCE MODEL FOR LARGE EDDY SIMULATION INCORPORATING THE COMPRESSIBLE "POOR MAN'S" NAVIER--STOKES EQUATIONS." UKnowledge, 2012. http://uknowledge.uky.edu/me_etds/13.

Full text

Abstract:

A new approach to large-eddy simulation (LES) based on the use of explicit spatial filtering combined with backscatter forcing is presented. The forcing uses a discrete dynamical system (DDS) called the compressible ``poor man's'' Navier--Stokes (CPMNS) equations. This DDS is derived from the governing equations and is shown to exhibit good spectral and dynamical properties for use in a turbulence model. An overview and critique of existing turbulence theory and turbulence models is given. A comprehensive theoretical case is presented arguing that traditional LES equations contain unresolved scales in terms generally thought to be resolved, and that this can only be solved with explicit filtering. The CPMNS equations are then incorporated into a simple forcing in the OVERFLOW compressible flow code, and tests are done on homogeneous, isotropic, decaying turbulence, a Mach 3 compression ramp, and a Mach 0.8 open cavity. The numerical results validate the general filter-forcing approach, although they also reveal inadequacies in OVERFLOW and that the current approach is likely too simple to be universally applicable. Two new proposals for constructing better forcing models are presented at the end of the work.

APA, Harvard, Vancouver, ISO, and other styles

12

Majer, Clemens Philipp 1963. "Stability investigations of a laminar wall jet using the complete Navier-Stokes equations." Thesis, The University of Arizona, 1991. http://hdl.handle.net/10150/278052.

Full text

Abstract:

The hydrodynamic stability of a plane, two-dimensional, incompressible wall jet subjected to small disturbances is investigated by direct numerical integration of the complete Navier-Stokes equations. The numerical model allows growing or decaying of disturbances in the downstream direction as in physical experiments. In the past, various numerical investigations were published using the linear stability theory for the case of temporally growing disturbances. In this work, the investigations are made for the case of spatially growing disturbances. The neutral curves of the linear stability theory are displayed, and in addition, the downstream development of spatial growing disturbances is provided by using the complete Navier-Stokes equations. It is shown that the behavior of the disturbances is as predicted by the linear stability theory for a certain frequency using small disturbances. The changes in the downstream development of the flow subjected to large disturbances compared to the results using small disturbances is discussed. For large disturbance amplitudes, it was found that for the frequency of the disturbance waves used in the investigations the boundary layer mode clearly dominates the hydrodynamic stability.

APA, Harvard, Vancouver, ISO, and other styles

13

Cox, Christopher. "Development of a High-Order Navier-Stokes Solver Using Flux Reconstruction to Simulate Three-Dimensional Vortex Structures in a Curved Artery Model." Thesis, The George Washington University, 2018. http://pqdtopen.proquest.com/#viewpdf?dispub=10742679.

Full text

Abstract:

Low-order numerical methods are widespread in academic solvers and ubiquitous in industrial solvers due to their robustness and usability. High-order methods are less robust and more complicated to implement; however, they exhibit low numerical dissipation and have the potential to improve the accuracy of flow simulations at a lower computational cost when compared to low-order methods. This motivates our development of a high-order compact method using Huynh's flux reconstruction scheme for solving unsteady incompressible flow on unstructured grids. We use Chorin's classic artificial compressibility formulation with dual time stepping to solve unsteady flow problems. In 2D, an implicit non-linear lower-upper symmetric Gauss-Seidel scheme with backward Euler discretization is used to efficiently march the solution in pseudo time, while a second-order backward Euler discretization is used to march in physical time. We verify and validate implementation of the high-order method coupled with our implicit time stepping scheme using both steady and unsteady incompressible flow problems. The current implicit time stepping scheme is proven effective in satisfying the divergence-free constraint on the velocity field in the artificial compressibility formulation. The high-order solver is extended to 3D and parallelized using MPI. Due to its simplicity, time marching for 3D problems is done explicitly. The feasibility of using the current implicit time stepping scheme for large scale three-dimensional problems with high-order polynomial basis still remains to be seen. We directly use the aforementioned numerical solver to simulate pulsatile flow of a Newtonian blood-analog fluid through a rigid 180-degree curved artery model. One of the most physiologically relevant forces within the cardiovascular system is the wall shear stress. This force is important because atherosclerotic regions are strongly correlated with curvature and branching in the human vasculature, where the shear stress is both oscillatory and multidirectional. Also, the combined effect of curvature and pulsatility in cardiovascular flows produces unsteady vortices. The aim of this research as it relates to cardiovascular fluid dynamics is to predict the spatial and temporal evolution of vortical structures generated by secondary flows, as well as to assess the correlation between multiple vortex pairs and wall shear stress. We use a physiologically (pulsatile) relevant flow rate and generate results using both fully developed and uniform entrance conditions, the latter being motivated by the fact that flow upstream of a curved artery may not have sufficient straight entrance length to become fully developed. Under the two pulsatile inflow conditions, we characterize the morphology and evolution of various vortex pairs and their subsequent effect on relevant haemodynamic wall shear stress metrics.

APA, Harvard, Vancouver, ISO, and other styles

14

Sandro, Manservisi. "Optimal Boundary and Distributed Controls for the Velocity Tracking Problem for Navier-Stokes Flows." Diss., Virginia Tech, 1997. http://hdl.handle.net/10919/30608.

Full text

Abstract:

The velocity tracking problem is motivated by the desire to match a desired target flow with a flow which can be controlled through time dependent distributed forces or time dependent boundary conditions. The flow model is the Navier-Stokes equations for a viscous incompressible fluid and different kinds of controls are studied. Optimal distributed and boundary controls minimizing a quadratic functional and an optimal bounded distributed control are investigated. The distributed optimal and the bounded control are compared with a linear feedback control. Here, a unified mathematical formulation, covering several specific classes of meaningful control problems in bounded domains, is presented with a complete and detailed analysis of all these time dependent optimal control velocity tracking problems. We concentrate not only on questions such as existence and necessary first order conditions but also on discretization and computational aspects. The first order necessary conditions are derived in the continuous, in the semidiscrete time approximation and in the fully finite element discrete case. This derivation is needed to obtain an accurate meaningful numerical algorithm with a satisfactory convergence rate. The gradient algorithm is used and several numerical computations are performed to compare and understand the limits imposed by the theory. Some computational aspects are discussed without which problems of any realistic size would remain intractable. Ph. D.

APA, Harvard, Vancouver, ISO, and other styles

15

Reid, Francis John Edward School of Mathematics UNSW. "The weakly nonlinear stability of an oscillatory fluid flow." Awarded by:University of New South Wales. School of Mathematics, 2006. http://handle.unsw.edu.au/1959.4/33364.

Full text

Abstract:

A weakly nonlinear stability analysis was conducted for the flow induced in an incompressible, Newtonian, viscous fluid lying between two infinite parallel plates which form a channel. The plates are oscillating synchronously in simple harmonic motion. The disturbed velocity of the flow was written in the form of a series in powers of a parameter which is a measure of the distance away from the linear theory neutral conditions. The individual terms of this series were decomposed using Floquet theory and Fourier series in time. The equations at second order and third order in were derived, and solutions for the Fourier coefficients were found using pseudospectral methods for the spatial variables. Various alternative methods of computation were applied to check the validity of the results obtained. The Landau equation for the amplitude of the disturbance was obtained, and the existence of equilibrium amplitude solutions inferred. The values of the coefficients in the Landau equation were calculated for the nondimensional channel half-widths h for the cases h = 5, 8, 10, 12, 14 and 16. It was found that equilibrium amplitude solutions exist for points in wavenumber Reynolds number space above the smooth portion of the previously determined linear stability neutral curve in all the cases examined. Similarly, Landau coefficients were calculated on a special feature of the neutral curve (called a ???finger???) for the case h = 12. Equilibrium amplitude solutions were found to exist at points inside the finger, and in a particular region outside near the top of the finger. Traces of the x-components of the disturbance velocities have been presented for a range of positions across the channel, together with the size of the equilibrium amplitude at these positions. As well, traces of the x-component of the velocity of the disturbed flow and traces of the velocity of the basic flow have been given for comparison at a particular position in the channel.

APA, Harvard, Vancouver, ISO, and other styles

16

Erdmann, Robert Gerald. "Image-Based Numerical Simulation of Stokes Flow in Porous Media." Diss., The University of Arizona, 2006. http://hdl.handle.net/10150/195724.

Full text

Abstract:

Numerical models for the simulation of longitudinal and transverse Stokes flow in cylindrical periodic porous media are presented. The models, which are based on a finite-volume formulation in primitive variables, utilize digital image representations of the geometries to simulate, making them particularly well-suited for the rapid automated analysis of creeping flow in porous media with complex morphologies. Complete details of the model formulations are given, including extensive treatment of the pressure boundary conditions at the solid-liquid interface needed to guarantee convergence with all possible geometries. The convergence behavior of both models is tested, and the models are shown to be second-order accurate.The models are used to simulate flow over the whole range of volume fractions of liquid in several regular geometries. The longitudinal model is used to simulate flow in square arrays of circular and square ducts, and both models are used to simulate flow in square and hexagonal arrays of circular cylinders and square arrays of square cylinders rotated by varying amounts. For each of the geometries, accurate empirical expressions for the Darcy permeability as a function of volume fraction solid are presented. Where applicable, model predictions of permeability are compared to existing analytical results.Subsequently, the models are used to simulate Stokes flow in random domains over a wide range of fractions liquid. The sequential random packing algorithm is used to generate 1,000 random packings of circular cylinders at each of 14 fractions of liquid, and longitudinal and transverse flow simulations are performed for each geometry. Histograms and summary statistics are computed for the permeability for each fraction liquid, and empirical expressions for mean permeability as a function of fraction liquid are given. The autocorrelation structure of the geometry and of the fluid velocity is analyzed, and an analysis of the scaling of longitudinal permeability variance is presented. In transverse flow at high packing densities, it is found that lightning-like patterns emerge in the fluid velocity. It is also found that the details of flows in such geometries are strongly sensitive to the placement of individual solid obstacles.

APA, Harvard, Vancouver, ISO, and other styles

17

Durlofsky, Louis J. "Topics in fluid mechanics : I. flow between finite rotating disks II. simulation of hydrodynamically interacting particles in stokes flow." Thesis, Massachusetts Institute of Technology, 1986. http://hdl.handle.net/1721.1/15049.

Full text

Abstract:

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 1986. MICROFICHE COPY AVAILABLE IN ARCHIVES AND SCIENCE. Includes bibliographies. by Louis J. Durlofsky. I. Flow between finite rotating disks II. Simulation of hydrodynamically interacting particles in stokes flow. Ph.D.

APA, Harvard, Vancouver, ISO, and other styles

18

Tryggeson, Henrik. "Analytical vortex solutions to Navier-Stokes equation." Doctoral thesis, Växjö universitet, Matematiska och systemtekniska institutionen, 2007. http://urn.kb.se/resolve?urn=urn:nbn:se:vxu:diva-1282.

Full text

Abstract:

Fluid dynamics considers the physics of liquids and gases. This is a branch of classical physics and is totally based on Newton's laws of motion. Nevertheless, the equation of fluid motion, Navier-Stokes equation, becomes very complicated to solve even for very simple configurations. This thesis treats mainly analytical vortex solutions to Navier-Stokes equations. Vorticity is usually concentrated to smaller regions of the flow, sometimes isolated objects, called vortices. If one are able to describe vortex structures exactly, important information about the flow properties are obtained. Initially, the modeling of a conical vortex geometry is considered. The results are compared with wind-tunnel measurements, which have been analyzed in detail. The conical vortex is a very interesting phenomenaon for building engineers because it is responsible for very low pressures on buildings with flat roofs. Secondly, a suggested analytical solution to Navier-Stokes equation for internal flows is presented. This is based on physical argumentation concerning the vorticity production at solid boundaries. Also, to obtain the desired result, Navier-Stokes equation is reformulated and integrated. In addition, a model for required information of vorticity production at boundaries is proposed. The last part of the thesis concerns the examples of vortex models in 2-D and 3-D. In both cases, analysis of the Navier-Stokes equation, leads to the opportunity to construct linear solutions. The 2-D studies are, by the use of diffusive elementary vortices, describing experimentally observed vortex statistics and turbulent energy spectrums in stratified systems and in soapfilms. Finally, in the 3-D analysis, three examples of recent experimentally observed vortex objects are reproduced theoretically. First, coherent structures in a pipe flow is modeled. These vortex structures in the pipe are of interest since they appear for Re in the range where transition to turbulence is expected. The second example considers the motion in a viscous vortex ring. The model, with diffusive properties, describes the experimentally measured velocity field as well as the turbulent energy spectrum. Finally, a streched spiral vortex is analysed. A rather general vortex model that has many degrees of freedom is proposed, which also may be applied in other configurations.

APA, Harvard, Vancouver, ISO, and other styles

19

Li, Siran. "Analysis of several non-linear PDEs in fluid mechanics and differential geometry." Thesis, University of Oxford, 2017. https://ora.ox.ac.uk/objects/uuid:20866cbb-e5ab-4a6b-b9dc-88a247d15572.

Full text

Abstract:

In the thesis we investigate two problems on Partial Differential Equations (PDEs) in differential geometry and fluid mechanics. First, we prove the weak L p continuity of the Gauss-Codazzi-Ricci (GCR) equations, which serve as a compatibility condition for the isometric immersions of Riemannian and semi-Riemannian manifolds. Our arguments, based on the generalised compensated compactness theorems established via functional and micro-local analytic methods, are intrinsic and global. Second, we prove the vanishing viscosity limit of an incompressible fluid in three-dimensional smooth, curved domains, with the kinematic and Navier boundary conditions. It is shown that the strong solution of the Navier-Stokes equation in H r+1 (r &GT; 5/2) converges to the strong solution of the Euler equation with the kinematic boundary condition in H r, as the viscosity tends to zero. For the proof, we derive energy estimates using the special geometric structure of the Navier boundary conditions; in particular, the second fundamental form of the fluid boundary and the vorticity thereon play a crucial role. In these projects we emphasise the linkages between the techniques in differential geometry and mathematical hydrodynamics.

APA, Harvard, Vancouver, ISO, and other styles

20

Carter, E. J. "Flow of power law fluids with application to oil drilling." Thesis, University of Plymouth, 1986. http://hdl.handle.net/10026.1/458.

Full text

Abstract:

The thesis is concerned with a theoretical study of the flow behaviour of inelastic power law fluids in two different types of flow situation. These are: 1. The creeping motion of a sphere moving through an expanse of liquid. 2. The combined steady and oscillatory flow of a liquid through a straight tube of circular cross section. The first part of the work is devoted to the prediction of the drag correction factor for a sphere falling slowly through a bounded inelastic power law fluid. The analysis is carried out for the case when the outer spherical boundary has a finite or infinite radius. A perturbation technique is used to produce the resulting equations for a slightly power law fluid which are solved using the finite element method. An asymptotic expansion is used to provide an analytical far field solution for the infinite outer sphere case. The second part considers the combined steady and oscillatory flow of an inelastic power law liquid in a tube. The analysis is carried. out for the case when both the steady flow rate and the oscillatory flow rate are known. An expression for the pressure gradient reduction in the tube is then derived. The resulting partial differential equation is solved by finite difference techniques. An analytical solution for the pressure gradient is also obtained using a perturbation analysis for the case when the fluid inertial effects are small.

APA, Harvard, Vancouver, ISO, and other styles

21

Le, Anh Ha. "A posteriori error estimation for simulation of diffusion and fluid mechanics problems by finite volume techniques." Paris 13, 2011. http://www.theses.fr/2011PA132055.

Full text

Abstract:

Ce travail concerne les estimations a posteriori pour la simulation numérique, par des méthodes de volumes finis, de plusieurs modèles liés à la mécanique des fluides : le problème de Stokes stationnaire, l'équation de Darcy non-linéaire et l'équation de transport linéaire. De plus, un résultat concerant des inégalités de Poincaré discrètes sur des maillages bidimensionels quelconques est également présenté.

APA, Harvard, Vancouver, ISO, and other styles

22

Padhye, V. A. "A penalty finite element model for axisymmetric flows of power-law and white-metzner fluids." Thesis, Virginia Polytechnic Institute and State University, 1986. http://hdl.handle.net/10919/101358.

Full text

Abstract:

A finite element model based on the penalty function formulation of the equations governing unsteady axisymmetric flows of viscous incompressible fluids obeying power-law and White-Metzner constitutive relations is developed. The formulation accounts for inertial (or convective) terms. For power-law fluids, two different finite element models are developed: one based on the velocity formulation (i.e., only velocities as nodal variables) and the other based on a mixed formulation involving velocities and stress components. For the White-Metzner model, only the mixed model involving the velocities and extra stress components can be developed. The pressure variable does not enter the finite element model because of the application of the penalty method to introduce the incompressibility constraint. However, the pressure can be post-computed once the velocities are obtained. The finite element models are used to analyze several plane and axisymmetric flows of power-law and White-Metzner fluids. The effect of boundary conditions, power-law model, inertia terms on the velocity profiles is investigated. The numerical solutions agree, qualitatively, with the known experimental and numerical results. The finite element models developed here can be easily modified to include thermal effects and other constitutive models. M.S.

APA, Harvard, Vancouver, ISO, and other styles

23

Doradoux, Adrien. "Simulation numérique d’écoulements diphasiques autour d’un solide mobile." Thesis, Bordeaux, 2017. http://www.theses.fr/2017BORD0582/document.

Full text

Abstract:

Les méthodes de domaines fictifs permettent de simuler numériquement des écoulements autour de structures complexes et/ou mobiles à l’aide de maillages simples. L’objet solide est alors « immergé » dans un domaine de calcul englobant le fluide et le solide. Dans un premier temps, on étudie une méthode de pénalisation, qui consiste à ajouter un terme dans l’équation de conservation de la quantité de mouvement du fluide afin d’imposer la vitesse du solide. Grâce à des développements asymptotiques, on obtient une estimation de l’erreur induite par cette approche lorsque le solide est en mouvement. Ce procédé est ensuite couplé avec un schéma de projection vectorielle permettant d’imposer la contrainte d’incompressibilité. La convergence du schéma ainsi obtenu, vers les équations de Navier-Stokes, est établie. Dans un second temps,une approche originale capable de traiter des écoulements multiphasiques est développée : la méthode de porosité variable. L’idée principale est de considérer le solide comme un milieu sans masse. La discrétisation des bilans massiques de chaque phase est alors modifiée, de sorte que le volume total occupé par l’ensemble des phases fluides soit égal au volume laissé libre parle solide. Cette méthode est validée numériquement sur un ensemble varié de cas test comprenantd es écoulements monophasiques incompressibles et compressibles ainsi que des écoulements diphasiques Fictitious domain methods allow to simulate flows around complex and/or moving bodies with simple meshes. The object is "immersed" in a domain that contains fluid and solid volumes. The penalization method, which consists in adding a term in the momentum balance equation, in order to impose the solid velocity, is studied in a first part. Thanks to asymptotic expansions, the order of the error induced by this method is computed for moving bodies. This approach is then coupled with a Vector Penalty Projection scheme that permits to impose the incompressibility constraint. The convergence of the penalized scheme towards the Navier-Stokes equations is established. In a second part, an original approach, able to treat multiphase flowsis presented: the Time and Space Dependent Porosity method. The key idea is to consider the solid as a medium without mass. The discretization of the mass balance equation is modified,so that the total volume occupied by all fluid phases and the solid is equal to the total volume.This method is numerically validated on a set of various test cases including incompressible or compressible single phase flows and two-phase flows

APA, Harvard, Vancouver, ISO, and other styles

24

Kong, Dali. "Analytical and numerical studies of several fluid mechanical problems." Thesis, University of Exeter, 2012. http://hdl.handle.net/10036/3651.

Full text

Abstract:

In this thesis, three parts, each with several chapters, are respectively devoted to hydrostatic, viscous and inertial fluids theories and applications. In the hydrostatics part, the classical Maclaurin spheroids theory is generalized, for the first time, to a more realistic multi-layer model, which enables the studies of some gravity problems and direct numerical simulations of flows in fast rotating spheroidal cavities. As an application of the figure theory, the zonal flow in the deep atmosphere of Jupiter is investigated for a better understanding of the Jovian gravity field. High viscosity flows, for example Stokes flows, occur in a lot of processes involving low-speed motions in fluids. Microorganism swimming is such typical a case. A fully three dimensional analytic solution of incompressible Stokes equation is derived in the exterior domain of an arbitrarily translating and rotating prolate spheroid, which models a large family of microorganisms such as cocci bacteria. The solution is then applied to the magnetotactic bacteria swimming problem and good consistency has been found between theoretical predictions and laboratory observations of the moving patterns of such bacteria under magnetic fields. In the analysis of dynamics of planetary fluid systems, which are featured by fast rotation and very small viscosity effects, three dimensional fully nonlinear numerical simulations of Navier-Stokes equations play important roles. A precession driven flow in a rotating channel is studied by the combination of asymptotic analyses and fully numerical simulations. Various results of laminar and turbulent flows are thereby presented. Computational fluid dynamics requires massive computing capability. To make full use of the power of modern high performance computing facilities, a C++ finite-element analysis code is under development based on PETSc platform. The code and data structures will be elaborated, along with the presentations of some preliminary results.

APA, Harvard, Vancouver, ISO, and other styles

25

Li, Zhiyong. "Data-Driven Adaptive Reynolds-Averaged Navier-Stokes k - ω Models for Turbulent Flow-Field Simulations". UKnowledge, 2017. http://uknowledge.uky.edu/me_etds/93.

Full text

Abstract:

The data-driven adaptive algorithms are explored as a means of increasing the accuracy of Reynolds-averaged turbulence models. This dissertation presents two new data-driven adaptive computational models for simulating turbulent flow, where partial-but-incomplete measurement data is available. These models automatically adjust (i.e., adapts) the closure coefficients of the Reynolds-averaged Navier-Stokes (RANS) k-ω turbulence equations to improve agreement between the simulated flow and a set of prescribed measurement data. The first approach is the data-driven adaptive RANS k-ω (D-DARK) model. It is validated with three canonical flow geometries: pipe flow, the backward-facing step, and flow around an airfoil. For all 3 test cases, the D-DARK model improves agreement with experimental data in comparison to the results from a non-adaptive RANS k-ω model that uses standard values of the closure coefficients. The second approach is the Retrospective Cost Adaptation (RCA) k-ω model. The key enabling technology is that of retrospective cost adaptation, which was developed for real-time adaptive control technology, but is used in this work for data-driven model adaptation. The algorithm conducts an optimization, which seeks to minimize the surrogate performance, and by extension the real flow-field error. The advantage of the RCA approach over the D-DARK approach is that it is capable of adapting to unsteady measurements. The RCA-RANS k-ω model is verified with a statistically steady test case (pipe flow) as well as two unsteady test cases: vortex shedding from a surface-mounted cube and flow around a square cylinder. The RCA-RANS k-ω model effectively adapts to both averaged steady and unsteady measurement data.

APA, Harvard, Vancouver, ISO, and other styles

26

Baggio, Gustavo Aparecido Pita [UNESP]. "Solução das equações de Navier-Stokes e da equação de transporte por um método de elementos finitos estabilizado pela técnica de separação baseada na característica." Universidade Estadual Paulista (UNESP), 2015. http://hdl.handle.net/11449/132106.

Full text

Abstract:

Made available in DSpace on 2015-12-10T14:23:55Z (GMT). No. of bitstreams: 0 Previous issue date: 2015-07-08. Added 1 bitstream(s) on 2015-12-10T14:30:00Z : No. of bitstreams: 1 000853535.pdf: 3830279 bytes, checksum: 808e22162a974efb82537b4b5a583580 (MD5) O método clássico de elementos finitos denominado de método de Bubnov-Galerkin ou Método de Galerkin (GFEM) é ótimo para solução de problemas difusivos, entretanto, para o caso de problemas com convecção dominante como ocorre em muitos problemas de escoamentos de fluidos, o GFEM produz soluções oscilantes para altos números de Reynolds. Diversas técnicas têm sido associadas à discretização dos termos convectivos nas equações de Navier-Stokes para remediar a deficiência do GFEM. Uma destas técnicas é denominada na literatura de separação baseada na característica (em inglês CBS - Characteristic Based Split), e tem sido aplicada com sucesso para estabilização das soluções de problemas convectivos usando o GFEM. Neste trabalho, o GFEM associado ao esquema CBS será aplicado para simulação de casos de escoamentos incompressíveis bidimensionais e tridimensionais The classical finite element method called the method of Bubnov-Galerkin or Galerkin Method (GFEM) is optimum for solving diffusion problems, however, in case of problems with dominant convection as in many problems of fluid flow, the GFEM produces oscillating solutions for high Reynolds numbers. Several techniques have been associated with the discretization of the convective terms in the Navier-Stokes equations to remedy the deficiency of GFEM. One of such techniques is denominated in literature as Characteristic Based Split (CBS) and has been applied successfully to stabilize the solutions of convective problems using GFEM. In this work, the GFEM associated with the Characteristic Based Split (CBS) scheme will be applied to simulate two and three-dimensional incompressible fluid flows

APA, Harvard, Vancouver, ISO, and other styles

27

Baggio, Gustavo Aparecido Pita. "Solução das equações de Navier-Stokes e da equação de transporte por um método de elementos finitos estabilizado pela técnica de separação baseada na característica /." Ilha Solteira, 2015. http://hdl.handle.net/11449/132106.

Full text

Abstract:

Orientador: João Batista Campos Silva Banca: João Batista Aparecido Banca: Estaner Claro Romão Resumo: O método clássico de elementos finitos denominado de método de Bubnov-Galerkin ou Método de Galerkin (GFEM) é ótimo para solução de problemas difusivos, entretanto, para o caso de problemas com convecção dominante como ocorre em muitos problemas de escoamentos de fluidos, o GFEM produz soluções oscilantes para altos números de Reynolds. Diversas técnicas têm sido associadas à discretização dos termos convectivos nas equações de Navier-Stokes para remediar a deficiência do GFEM. Uma destas técnicas é denominada na literatura de separação baseada na característica (em inglês CBS - Characteristic Based Split), e tem sido aplicada com sucesso para estabilização das soluções de problemas convectivos usando o GFEM. Neste trabalho, o GFEM associado ao esquema CBS será aplicado para simulação de casos de escoamentos incompressíveis bidimensionais e tridimensionais Abstract: The classical finite element method called the method of Bubnov-Galerkin or Galerkin Method (GFEM) is optimum for solving diffusion problems, however, in case of problems with dominant convection as in many problems of fluid flow, the GFEM produces oscillating solutions for high Reynolds numbers. Several techniques have been associated with the discretization of the convective terms in the Navier-Stokes equations to remedy the deficiency of GFEM. One of such techniques is denominated in literature as Characteristic Based Split (CBS) and has been applied successfully to stabilize the solutions of convective problems using GFEM. In this work, the GFEM associated with the Characteristic Based Split (CBS) scheme will be applied to simulate two and three-dimensional incompressible fluid flows Mestre

APA, Harvard, Vancouver, ISO, and other styles

28

Finnie, John I. "An Application of the Finite Element Method and Two Equation (K and E) Turbulence Model to Two and Three Dimensional Fluid Flow Problems Governed by the Navier-Stokes Equations." DigitalCommons@USU, 1987. https://digitalcommons.usu.edu/etd/7350.

Full text

Abstract:

Finite Element computer codes in two and three dimensions were written that solve both laminar and turbulent flow. These codes use the two equation (k and E) turbulence model to evaluate turbulent viscosity. They were tested with 29 different flow problems. The largest two dimensional turbulent problem solved is flow under a sluice gate. A three dimensional vortex flow problem was attempted but was not feasible due to the size of the available computer. The Harwell sparse matrix subroutines of the United Kingdom Atomic Energy Authority were used to solve the set of simultaneous equations. The performance of these subroutines is evaluated. The importance of defining adequate finite element grids and setting proper boundary and initial conditions is discussed.

APA, Harvard, Vancouver, ISO, and other styles

29

Oliver, Todd A. "A High-Order, Adaptive, Discontinuous Galerkin Finite Element Method for the Reynolds-Averaged Navier-Stokes Equations." Ft. Belvoir : Defense Technical Information Center, 2008. http://handle.dtic.mil/100.2/ADA488912.

Full text

APA, Harvard, Vancouver, ISO, and other styles

30

Meslameni, Mohamed. "Équations de Stokes et d'Oseen en domaine extérieur avec diverses conditions aux limites." Thesis, Pau, 2013. http://www.theses.fr/2013PAUU3002/document.

Full text

Abstract:

On s’intéresse aux équations stationnaires de Navier-Stokes linéarisées, il s'agit ici des équations d'Oseen et des équations de Stokes posées dans des domaines infinis, comme les domaines extérieurs, en dimension trois et l'espace tout entier. Le but est d'étudier l'existence de solutions généralisés et de solutions fortes dans un cadre général non nécessairement hilbertien. On s'intéresse aussi au cas des solutions très faibles. Dans ce travail, on considère aussi bien des conditions aux limites classiques de type Dirichlet que des conditions aux limites non standard portant sur certaines composantes du champ de vitesses, du tourbillon, voir du champ de pression. Les espaces de Sobolev classiques ne sont pas adaptés à l'étude de ces problèmes pour une telle géométrie. Pour une bonne analyse mathématique, nous avons choisi de travailler dans le cadre des espaces de Sobolev avec poids, ce qui permet en particulier de mieux contrôler le comportement à l'infini de la solution In this work, we study the linearized Navier-Stokes equations in an exterior domain or in the whole space at the steady state, that is, the Stokes equations and the Oseen equations. We give existence, uniqueness and regularity of solutions. The case of very weak solutions is also treated. We consider not only the Dirichlet boundary conditions but also the Non Standard boundary conditions, on some components of the velocity field, vorticity and also on the pressure. Since the domain is not bounded, the classical Sobolev spaces are not adequate. Therefore, a specific functional framework is necessary which also has to take into account the behaviour of the functions at infinity. Our approach rests on the use of weighted Sobolev spaces

APA, Harvard, Vancouver, ISO, and other styles

31

McCraney, Joshua Thomas. "Analysis of Capillary Flow in Interior Corners : Perturbed Power Law Similarity Solutions." PDXScholar, 2015. https://pdxscholar.library.pdx.edu/open_access_etds/2725.

Full text

Abstract:

The design of fluid management systems requires accurate models for fluid transport. In the low gravity environment of space, gravity no longer dominates fluid displacement; instead capillary forces often govern flow. This thesis considers the redistribution of fluid along an interior corner. Following a rapid reduction of gravity, fluid advances along the corner measured by the column length z = L(t), which is governed by a nonlinear partial differential equation with dynamical boundary conditions. Three flow types are examined: capillary rise, spreading drop, and tapered corner. The spreading drop regime is shown to exhibit column length growth L ~ t2/5, where a closed form analytic solution exists. No analytic solution is available for the capillary rise problem. However, a perturbed power law similarity solution is pursued to approximate an analytic solution in the near neighborhood of the exact solution for the spreading drop. It is recovered that L ~ t1/2 for the capillary rise problem. The tapered corner problem is not analytically understood and hence its corresponding L is undocumented. Based on the slender corner geometry, it is natural to hypothesize the tapered corner column length initially behaves like the capillary rise regime, but after sufficient time has elapsed, it transitions into the spreading drop regime. This leads to a conjecture that its column length growth L is restricted to t2/5 < L < t1/2. To verify this conjecture an explicit finite difference numerical solution is developed for all three regimes. As will be shown, the finite difference scheme converges towards the analytic solutions for the spreading drop and capillary rise regimes. From this we assume the finite difference scheme is accurate for corner flows of similar geometries, and thus apply this scheme the more onerous criteria of the tapered corner. Numerical results support the conjectured L behavior for the tapered corner. Understanding the dynamics of such flows and responses to various geometries offers design advantages for spacecraft waste-management systems, fuel control, hydration containment, cryogenic flows, and a myriad of other fluid applications.

APA, Harvard, Vancouver, ISO, and other styles

32

Rahimian, Abtin. "Parallel algorithms for direct blood flow simulations." Diss., Georgia Institute of Technology, 2012. http://hdl.handle.net/1853/43611.

Full text

Abstract:

Fluid mechanics of blood can be well approximated by a mixture model of a Newtonian fluid and deformable particles representing the red blood cells. Experimental and theoretical evidence suggests that the deformation and rheology of red blood cells is similar to that of phospholipid vesicles. Vesicles and red blood cells are both area preserving closed membranes that resist bending. Beyond red blood cells, vesicles can be used to investigate the behavior of cell membranes, intracellular organelles, and viral particles. Given the importance of vesicle flows, in this thesis we focus in efficient numerical methods for such problems: we present computationally scalable algorithms for the simulation of dilute suspension of deformable vesicles in two and three dimensions. Our method is based on the boundary integral formulation of Stokes flow. We present new schemes for simulating the three-dimensional hydrodynamic interactions of large number of vesicles with viscosity contrast. The algorithms incorporate a stable time-stepping scheme, high-order spatiotemporal discretizations, spectral preconditioners, and a reparametrization scheme capable of resolving extreme mesh distortions in dynamic simulations. The associated linear systems are solved in optimal time using spectral preconditioners. The highlights of our numerical scheme are that (i) the physics of vesicles is faithfully represented by using nonlinear solid mechanics to capture the deformations of each cell, (ii) the long-range, N-body, hydrodynamic interactions between vesicles are accurately resolved using the fast multipole method (FMM), and (iii) our time stepping scheme is unconditionally stable for the flow of single and multiple vesicles with viscosity contrast and its computational cost-per-simulation-unit-time is comparable to or less than that of an explicit scheme. We report scaling of our algorithms to simulations with millions of vesicles on thousands of computational cores.

APA, Harvard, Vancouver, ISO, and other styles

33

ROBALINHO, ERIC. "Desenvolvimento de um modelo numerico computacional aplicado a uma celula a combustivel unitaria de 144 CMsup(2) tipo PEM." reponame:Repositório Institucional do IPEN, 2009. http://repositorio.ipen.br:8080/xmlui/handle/123456789/9396.

Full text

Abstract:

Made available in DSpace on 2014-10-09T12:26:29Z (GMT). No. of bitstreams: 0 Made available in DSpace on 2014-10-09T14:04:32Z (GMT). No. of bitstreams: 0 Tese (Doutoramento) IPEN/T Instituto de Pesquisas Energeticas e Nucleares - IPEN-CNEN/SP

APA, Harvard, Vancouver, ISO, and other styles

34

Bagge, Joar. "Numerical simulation of an inertial spheroidal particle in Stokes flow." Thesis, KTH, Numerisk analys, NA, 2015. http://urn.kb.se/resolve?urn=urn:nbn:se:kth:diva-180290.

Full text

Abstract:

Particle suspensions occur in many situations in nature and industry. In this master’s thesis, the motion of a single rigid spheroidal particle immersed in Stokes flow is studied numerically using a boundary integral method and a new specialized quadrature method known as quadrature by expansion (QBX). This method allows the spheroid to be massless or inertial, and placed in any kind of underlying Stokesian flow. A parameter study of the QBX method is presented, together with validation cases for spheroids in linear shear flow and quadratic flow. The QBX method is able to compute the force and torque on the spheroid as well as the resulting rigid body motion with small errors in a short time, typically less than one second per time step on a regular desktop computer. Novel results are presented for the motion of an inertial spheroid in quadratic flow, where in contrast to linear shear flow the shear rate is not constant. It is found that particle inertia induces a translational drift towards regions in the fluid with higher shear rate. Partikelsuspensioner förekommer i många sammanhang i naturen och industrin. I denna masteruppsats studeras rörelsen hos en enstaka stel sfäroidisk partikel i Stokesflöde numeriskt med hjälp av en randintegralmetod och en ny specialiserad kvadraturmetod som kallas quadrature by expansion (QBX). Metoden fungerar för masslösa eller tröga sfäroider, som kan placeras i ett godtyckligt underliggande Stokesflöde. En parameterstudie av QBX-metoden presenteras, tillsammans med valideringsfall för sfäroider i linjärt skjuvflöde och kvadratiskt flöde. QBX-metoden kan beräkna kraften och momentet på sfäroiden samt den resulterande stelkroppsrörelsen med små fel på kort tid, typiskt mindre än en sekund per tidssteg på en vanlig persondator. Nya resultat presenteras för rörelsen hos en trög sfäroid i kvadratiskt flöde, där skjuvningen till skillnad från linjärt skjuvflöde inte är konstant. Det visar sig att partikeltröghet medför en drift i sidled mot områden i fluiden med högre skjuvning.

APA, Harvard, Vancouver, ISO, and other styles

35

Bravin, Marco. "Dynamics of a viscous incompressible flow in presence of a rigid body and of an inviscid incompressible flow in presence of a source and a sink." Thesis, Bordeaux, 2019. http://www.theses.fr/2019BORD0192.

Full text

Abstract:

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

APA, Harvard, Vancouver, ISO, and other styles

36

van, den Bremer T. S. "The induced mean flow of surface, internal and interfacial gravity wave groups." Thesis, University of Oxford, 2014. http://ora.ox.ac.uk/objects/uuid:e735afe7-a77d-455d-a560-e869a9941f69.

Full text

Abstract:

Although the leading-order motion of waves is periodic - in other words backwards and forwards - many types of waves including those driven by gravity induce a mean flow as a higher-order effect. It is the induced mean flow of three types of gravity waves that this thesis examines: surface (part I), internal (part II) and interfacial gravity waves (part III). In particular, this thesis examines wave groups. Because they transport energy, momentum and other tracers, wave-induced mean flows have important consequences for climate, environment, air traffic, fisheries, offshore oil and other industries. In this thesis perturbation methods are used to develop a simplified understanding of the physics of the induced mean flow for each of these three types of gravity wave groups. Leading-order estimates of different transport quantities are developed. For surface gravity wave groups (part I), the induced mean flow consists of two compo- nents: the Stokes drift dominant near the surface and the Eulerian return flow acting in the opposite direction and dominant at depth. By considering subsequent orders in a separation of scales expansion and by comparing to the Fourier-space solutions of Longuet-Higgins and Stewart (1962), this thesis shows that the effects of frequency dis- persion can be ignored for deep-water waves with realistic bandwidths. An approximate depth scale is developed and validated above which the Stokes drift is dominant and below which the return flow wins: the transition depth. Results are extended to include the effects of finite depth and directional spreading. Internal gravity wave groups (part II) do not display Stokes drift, but a quantity analogous to Stokes transport for surface gravity waves can still be developed, termed the “divergent- flux induced flow” herein. The divergent-flux induced flow it itself a divergent flow and induces a response. In a three-dimensional geometry, the divergent-flux induced flow and the return flow form a balanced circulation in the horizontal plane with the former transporting fluid through the centre of the group and the latter acting in the opposite direction around the group. In a two-dimensional geometry, stratification inhibits a balanced circulation and a second type of waves are generated that travel far ahead and in the lee of the wave group. The results in the seminal work of Bretherton (1969b) are thus validated, explicit expressions for the response and return flow are developed and compared to numerical simulations in the two-dimensional case. Finally, for interfacial wave groups (part III) the induced mean flow is shown to behave analogously to the surface wave problem of part I. Exploring both pure interfacial waves in a channel with a closed lid and interacting surface and interfacial waves, expressions for the Stokes drift and return flow are found for different configurations with the mean set-up or set-down of the interface playing an important role.

APA, Harvard, Vancouver, ISO, and other styles

37

Simon, Helene Anne. "Numerical simulations of the micro flow field in the hinge region of bileaflet mechanical heart valves." Diss., Georgia Institute of Technology, 2009. http://hdl.handle.net/1853/34861.

Full text

Abstract:

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.

APA, Harvard, Vancouver, ISO, and other styles

38

Abbasi, Baharanchi Ahmadreza. "Development of a Two-Fluid Drag Law for Clustered Particles Using Direct Numerical Simulation and Validation through Experiments." FIU Digital Commons, 2015. http://digitalcommons.fiu.edu/etd/2489.

Full text

Abstract:

This dissertation focused on development and utilization of numerical and experimental approaches to improve the CFD modeling of fluidization flow of cohesive micron size particles. The specific objectives of this research were: (1) Developing a cluster prediction mechanism applicable to Two-Fluid Modeling (TFM) of gas-solid systems (2) Developing more accurate drag models for Two-Fluid Modeling (TFM) of gas-solid fluidization flow with the presence of cohesive interparticle forces (3) using the developed model to explore the improvement of accuracy of TFM in simulation of fluidization flow of cohesive powders (4) Understanding the causes and influential factor which led to improvements and quantification of improvements (5) Gathering data from a fast fluidization flow and use these data for benchmark validations. Simulation results with two developed cluster-aware drag models showed that cluster prediction could effectively influence the results in both the first and second cluster-aware models. It was proven that improvement of accuracy of TFM modeling using three versions of the first hybrid model was significant and the best improvements were obtained by using the smallest values of the switch parameter which led to capturing the smallest chances of cluster prediction. In the case of the second hybrid model, dependence of critical model parameter on only Reynolds number led to the fact that improvement of accuracy was significant only in dense section of the fluidized bed. This finding may suggest that a more sophisticated particle resolved DNS model, which can span wide range of solid volume fraction, can be used in the formulation of the cluster-aware drag model. The results of experiment suing high speed imaging indicated the presence of particle clusters in the fluidization flow of FCC inside the riser of FIU-CFB facility. In addition, pressure data was successfully captured along the fluidization column of the facility and used as benchmark validation data for the second hybrid model developed in the present dissertation. It was shown the second hybrid model could predict the pressure data in the dense section of the fluidization column with better accuracy.

APA, Harvard, Vancouver, ISO, and other styles

39

Polavieja, Gonzalo Garcia de. "Geometric phase and angle for noncyclic adiabatic change, revivals and measures of quantal instability." Thesis, University of Oxford, 1999. http://ethos.bl.uk/OrderDetails.do?uin=uk.bl.ethos.325986.

Full text

APA, Harvard, Vancouver, ISO, and other styles

40

Lustri, Christopher Jessu. "Exponential asymptotics in unsteady and three-dimensional flows." Thesis, University of Oxford, 2013. http://ora.ox.ac.uk/objects/uuid:9fe9517a-1733-496f-914b-b6739b2dce59.

Full text

Abstract:

The behaviour of free-surface gravity waves on small Froude number fluid flow past some obstacle cannot be determined using ordinary asymptotic power series methods, as the amplitude of the waves is exponentially small. An exponential asymptotic method is used by Chapman and Vanden-Broeck (2006) to consider the problem of two-dimensional, steady flow past a submerged obstacle in the small Froude number limit, finding that a steady downstream wavetrainis switched on rapidly across a curve known as a Stokes line. Here, equivalent wavetrains on three-dimensional and unsteady flow configurations are considered, and Stokes switching causedby the interaction between exponentially small free-surface components is shown to play an important role in both cases. The behaviour of free-surface gravity waves is introduced by considering the problem of steady free-surface flow due to a line source. A steady wavetrain is shown to exist in the far field, and the behaviour of these waves is compared to existing numerical results. The problem of unsteady flow over a step is subsequently investigated, with the flow behaviour formulated in terms of Lagrangian coordinates so that the position of the free surface is fixed. Initially, the problem is linearized in the step-height, and the steady wavetrain is shown to spread downstream over time. The position of the wavefront is determined by considering the full Stokes structure present in the problem. The equivalent fully-nonlinear problem is then considered, with the position of the Stokes lines, and hence the wavefront, being determined numerically. Finally, linearized three-dimensional free-surface flow past an obstacle is considered in both the steady and unsteady case. The surface is shown to contain downstream longitudinal and transverse waves. These waves are shown to propagate downstream in the unsteady case, with the position of the wavefront again determined by considering the full Stokes structure of the problem.

APA, Harvard, Vancouver, ISO, and other styles

41

Dhifaoui, Anis. "Équations de Stokes en domaine extérieur avec des conditions aux limites de type Navier." Thesis, Bourgogne Franche-Comté, 2020. http://www.theses.fr/2020UBFCD009.

Full text

Abstract:

On s'intéresse aux équations stationnaires de Stokes posées dans un domaine extérieur tridimensionnel décrivant l'écoulement d'un fluide visqueux et incompressible autour d'obstacle supposé borné. La particularité ici réside dans les conditions au bord de l'obstacle que nous avons imposées. En effet, nous supposons que l'obstacle a une certaine rugosité et par conséquent, le fluide n'adhère pas au bord de l'obstacle mais, au contraire, il existe une friction dont on suppose décrite par les conditions aux limites de type Navier. Ces dernières modélisent d'une part l'imperméabilité de l'obstacle et d'autre part le fait que la composante tangentielle de la vitesse du fluide sur l'obstacle est proportionnelle au tenseur des déformations. Ce problème a été bien étudié lorsqu'il est posé dans un domaine borné. Les espaces de Sobolev classiques fournissent, dans ce cas, un cadre fonctionnel adéquat pour une étude complète. Cependant lorsque le domaine n'est pas borné, ces espaces ne sont plus adaptés car il est nécessaire de décrire le comportement à l'infini des solutions. On choisit alors de poser le problème dans des espaces de Sobolev avec des poids polynomiaux qui précisent la croissance ou la décroissance des fonctions à l'infini. Dans ce travail, nous commençons par effectuer une analyse hilbertienne du problème. Le point-clé ici est d'établir des inégalités de type Korn avec poids afin d'obtenir la coercivité de la forme bilinéaire associée à la formulation variationnelle. Nous continuons par démontrer des résultats d'existence, d'unicité et de régularité de solutions fortes et très faibles. Enfin, nous étudions l'extension de certains résultats en théorie L^p In this manuscript, we study the three-dimensional stationary Stokes equations set in a exterior domain. The problem describes the flow of a viscous and incompressible fluid past a bounded obstacle. The distinctif feature here relies on the fact that the obstacle is assumed to a rough boundary. As a result, the fluid may slip on the boundary of the obstacle and, to take into account this property, we use the Navier boundary conditions. On the one hand, They model the impermeability of the obstacle, and on the other hand, the fact that the tangential component of the fluid velocity on the obstacle is proportional to the stress tensor. This problem has been well studied when set in a bounded domain. The standard Sobolev spaces provides, in this case, an adequate functional framework for a complete study. Since in our case, the domain is unbounded, these spaces are not adapted since it is necessary to describe the behaviour of the solutions to infinity. Therefore, we choose to set the problem in weighted Sobolev spaces where the weights describe the behaviour at infinity of the function (growth or decay).In this work, we first start by performing the mathematical analysis in the Hilbert setting. The key point here is to establish variant weighted Korn’s inequalities in order to get the coercivity of the bilinear form associated to the variational formulation. Next, we proved the existence, uniqueness of strong and very weak solutions. Finally, we study the extension of some of thses results to a weightedL^p-theory

APA, Harvard, Vancouver, ISO, and other styles

42

Conti, Thadeu das Neves. "Aplicação do método da expansão em funções hierárquicas na solução das equações de Navier-Stokes em duas dimensões para fluidos compressíveis em alta velocidade." Universidade de São Paulo, 2006. http://www.teses.usp.br/teses/disponiveis/3/3132/tde-05092006-154143/.

Full text

Abstract:

O trabalho desenvolvido nesta tese propõe a aplicação do método da expansão em funções hierárquicas elaborado por Zienkiewics e Morgan (1983), para a solução das equações de conservação da massa (continuidade), conservação da quantidade de movimento (Navier-Stokes) e conservação da energia, para fluidos compressíveis em duas dimensões e em alta velocidade. Esse método consiste no emprego do método de elementos finitos utilizando a formulação Petrov-Galerkin, conhecido como SUPG (Streamline Upwind Petrov-Galerkin), desenvolvido por Brooks e Hughes (1982), aplicado em conjunto com uma expansão das variáveis em funções hierárquicas. A fim de testar e validar o método numérico proposto, assim como o programa computacional elaborado, foram simulados alguns casos conhecidos da literatura. Os casos estudados foram os seguintes: teste de Continuidade; teste de convergência e estabilidade; problema do degrau de temperatura e problema do choque oblíquo, onde o objetivo desse último caso era, basicamente, verificar a captura da onda de choque pelo método numérico desenvolvido. Através dos casos estudados e em função dos resultados obtidos nas simulações realizadas, conclui-se que o objetivo desse trabalho foi alcançado de maneira satisfatória, pois os resultados obtidos com o método desenvolvido nesse trabalho foram qualitativamente e quantitativamente bons, quando comparados com os resultados teóricos. The Thesis develops a new application for the Hierarchical Function Expansion Method, proposed by Zienkiewics and Morgan (1983), for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as, SUPG (Streamline Upwind Petrov-Galerkin) developed by Brooks and Hughes (1982), and applied in conjunction with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed some cases whose theoretical solution are known simulated. These cases are the following: continuity test; stability and convergence test; temperature step problem; and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations of the cases performed with the proposed method were good both qualitatively and quantitatively when compared with the teorethical solutions. This allows us to conclude that the objective of this Thesis was satisfactorily reached.

APA, Harvard, Vancouver, ISO, and other styles

43

Turuban, Régis. "Études numérique et expérimentales du mélange en milieux poreux 2D et 3D." Thesis, Rennes 1, 2017. http://www.theses.fr/2017REN1S051/document.

Full text

Abstract:

Le mélange de solutés par les écoulements en milieux poreux contrôle les réactions chimiques dans un grand nombre d'applications souterraines, dont le transport et la remédiation des contaminants, le stockage et l'extraction souterrains d'énergie, et la séquestration du CO2. Nous étudions les mécanismes du mélange à l'échelle du pore et plus précisément comment la topologie de l'écoulement est reliée à la dynamique du mélange d'espèces conservatives; en particulier, l'émergence d'un mélange chaotique est-elle possible dans un milieu poreux tridimensionnel (3D) ? Nous calculons donc numériquement ou mesurons expérimentalement les vitesses d'écoulement et l'évolution temporelle des champs de concentration afin de caractériser la déformation et le mélange à l'échelle du pore. Une première étude, expérimentale, permet de caractériser le mélange dans un fluide s’écoulant à travers un milieu poreux bidimensionnel (2D). Nous mesurons les vitesses par suivi de microparticules solides (''PTV''). L’évolution temporelle de la distance séparant deux particules permet de caractériser la dynamique de la déformation lagrangienne. Des mesures de transport conservatif dans le même milieu fournissent l'évolution temporelle du gradient de concentration moyen (une mesure du mélange). À partir de ces résultats expérimentaux nous proposons la première validation expérimentale à l'échelle du pore de la théorie lamellaire du mélange, reliant les propriétés de la déformation du fluide à la dynamique du mélange. Dans une deuxième étude nous examinons les conditions d'apparition du mélange chaotique dans l’écoulement dans des milieux poreux 3D granulaires ordonnés. Nous effectuons des calculs numériques hautement résolus de d'écoulement de Stokes entre des sphères empilées selon une structure cristalline (cubique simple ou cubique centrée), périodique. La déformation lagrangienne, obtenue à partir des champs de vitesse à l'aide d'outils numériques développés spécifiquement, met en lumière une large variété de dynamiques de la déformation dans ces milieux 3D, selon l'orientation de l'écoulement. Quand la direction de l'écoulement n'est pas normale à l'un des plans de symétrie de réflection du cristal, l'évolution temporelle de la déformation est exponentielle, traduisant une advection chaotique. L’émergence (ou non) du chaos est contrôlée par un mécanisme similaire à la ''transformation du boulanger'': les particules fluides se déplaçant autour d'un grain solide se retrouvent séparées par une surface virtuelle (appelée “variété”) qui émerge de la surface du grain. De multiples variétés existent dans l’écoulement, et la façon dont elles s'intersectent contrôle la nature - chaotique ou non - du mélange, et l'intensité du chaos. En particulier, l'exposant de Lyapunov (une mesure du chaos), est contrôlé par la fréquence spatiale des intersections appropriées à la génération du chaos, nommées ''connections hétéroclinines'' entre variétés. L'image conventionnelle, 2D, des mécanismes du mélange, impose des contraintes topologiques qui ne permettent pas le développement de ces mécanismes 3D. Elle pourrait donc être inadaptée aux milieux poreux naturels. La troisième étude a deux objectifs: (i) fournir une preuve expérimentale de la nature chaotique de l'advection, par la visualisation des variétés et par l'obtention d'une mesure de l'exposant de Lyapunov; et (ii), évaluer si nos résultats numériques obtenus pour des milieux granulaires ordonnés peuvent être généralisés à des milieux désordonnés, plus proches des milieux naturels. L’expérience est fondée sur un empilement désordonné de sphères rendu transparent par l'ajustement optique du liquide avec les sphères. La fluorescence induite par laser (''LIF'') permet de détecter les variétés au sein de l'écoulement, et des techniques PTV de mesurer les vitesses d'écoulement et quantifier l'exposant de Lyapunov. Les premiers résultats expérimentaux sont prometteurs Solute mixing in porous media flows plays a central role in driving chemical reactions in a number of subsurface applications, including contaminant transport and remediation, subsurface energy storage and extraction, and CO2 sequestration. We study the mechanisms of solute mixing, in particular how the pore scale flow topology is related to the mixing dynamics of conservative solutes, with a particular emphasis on the possible emergence of chaotic mixing processes in three-dimensional (3D) porous media. To do so, we perform numerical computations or experimental measurements of the flow velocities and temporal evolution of the concentration fields, and characterize fluid deformation and mixing at the pore scale. This PhD work consists of three main studies. In the first study, we experimentally characterize mixing in a fluid flowing through a two-dimensional (2D) porous medium built by lithography. We measure the velocity distributions from Particle Tracking Velocimetry (PTV). The time evolution of the separation distance between two particles is analyzed to characterize the Lagrangian deformation dynamics. In parallel we perform conservative transport experiments with the same porous media, and quantify the temporal evolution of the mean concentration gradient, which is a measure of the mixing rate. From these experimental results we obtain the first experimental pore scale validation of the lamella mixing theory, which relates the fluid deformation properties to the mixing dynamics. In the second study, we investigate the conditions of emergence of chaotic mixing in the flow through 3D ordered granular porous media. In these periodic cubic crystalline packings (Simple Cubic - SC - and Body-Centered Cubic - BCC) of spheres, we are able to perform highly resolved computations of the 3D Stokes flow. Using custom-developed numerical tools to measure the Lagrangian deformation from the computed velocity fields, we uncover the existence of a rich array of Lagrangian deformation dynamics in these 3D media, depending on the flow orientation. When the flow direction is not normal to one of the reflection symmetry planes of the crystalline lattice, we find that the Lagrangian deformation dynamics follow an exponential law, which indicates chaotic advection. This chaotic behavior is controlled by a mechanism akin to the baker's transformation: fluid particles traveling around a solid grain along different paths end up either separated by, or on the same side of, a virtual surface projecting from the grain surface and called a manifold. Multiple such manifolds exist within the flow, and the way they intersect controls the nature of mixing (that is, either non-chaotic or chaotic), and the strength of chaos. We show in particular that the magnitude of the Lyapunov exponent (a measure of the vigor of chaos) is controlled by the spatial frequency of transverse connections between the manifolds (called heteroclinic intersections). We thus demonstrate that the conventional 2D picture of the mechanisms of mixing may not be adapted for natural porous media because that picture imposes topological constraints which cannot account for these important 3D mechanisms. The third study has two objectives: (i) provide experimental evidence of the chaotic nature of pore scale advection/mixing, both by visualizing the manifolds and by obtaining a quantitative estimate of the Lyapunov exponent; and (ii) assess if the results obtained numerically in ordered packings of spheres extend to random packings, which are closer to natural porous media. The experiment features a random packing of glass beads rendered transparent by optical index-matching between the fluid and solid grains. We use Laser Induced Fluorescence (LIF) to detect the manifolds, and PTV techniques to measure flow velocities and subsequently quantify Lyapunov exponent. The first experimental results are promising

APA, Harvard, Vancouver, ISO, and other styles

44

Barker, Tobias. "Uniqueness results for viscous incompressible fluids." Thesis, University of Oxford, 2017. https://ora.ox.ac.uk/objects/uuid:db1b3bb9-a764-406d-a186-5482827d64e8.

Full text

Abstract:

First, we provide new classes of initial data, that grant short time uniqueness of the associated weak Leray-Hopf solutions of the three dimensional Navier-Stokes equations. The main novelty here is the establishment of certain continuity properties near the initial time, for weak Leray-Hopf solutions with initial data in supercritical Besov spaces. The techniques used here build upon related ideas of Calderón. Secondly, we prove local regularity up to the at part of the boundary, for certain classes of solutions to the Navier-Stokes equations, provided that the velocity field belongs to L∞(-1; 0; L3, β(B(1) &xcap; &Ropf;3 +)) with 3 ≤ β < ∞. What enables us to build upon the work of Escauriaza, Seregin and &Scaron;verák [27] and Seregin [100] is the establishment of new scale-invariant estimates, new estimates for the pressure near the boundary and a convenient new ϵ-regularity criterion. Third, we show that if a weak Leray-Hopf solution in &Ropf;3 +×]0,∞[ has a finite blow-up time T, then necessarily limt↑T||v(·, t)||L3,β(&Ropf;3 +) = ∞ with 3 < β < ∞. The proof hinges on a rescaling procedure from Seregin's work [106], a new stability result for singular points on the boundary, suitable a priori estimates and a Liouville type theorem for parabolic operators developed by Escauriaza, Seregin and &Scaron;verák [27]. Finally, we investigate a notion of global-in-time solutions to the Navier- Stokes equations in &Ropf;3, with solenoidal initial data in the critical Besov space ?-1/44,∞(&Ropf;3), which has certain continuity properties with respect to weak&ast; convergence of the initial data. Such properties are motivated by the strategy used by Seregin [106] to show that if a weak Leray-Hopf solution in &Ropf;3×]0,∞[ has a finite blow-up time T, then necessarily limt↑T ||v(·, t)||L3(&Ropf;3) = ∞. We prove new decomposition results for Besov spaces, which are key in the conception and existence theory of such solutions.

APA, Harvard, Vancouver, ISO, and other styles

45

McCormick, Matthew. "Ventricular function under LVAD support." Thesis, University of Oxford, 2012. http://ora.ox.ac.uk/objects/uuid:0d49ba30-b508-4c69-9ba6-b398d4338c01.

Full text

Abstract:

This thesis presents a finite element methodology for simulating fluid–solid interactions in the left ventricle (LV) under LVAD support. The developed model was utilised to study the passive and active characteristics of ventricular function in anatomically accurate LV geometries constructed from normal and patient image data. A non–conforming ALE Navier–Stokes/finite–elasticity fluid–solid coupling system formed the core of the numerical scheme, onto which several novel numerical additions were made. These included a fictitious domain (FD) Lagrange multiplier method to capture the interactions between immersed rigid bodies and encasing elastic solids (required for the LVAD cannula), as well as modifications to the Newton–Raphson/line search algorithm (which provided a 2 to 10 fold reduction in simulation time). Additional developments involved methods for extending the model to ventricular simulations. This required the creation of coupling methods, for both fluid and solid problems, to enable the integration of a lumped parameter representation of the systemic and pulmonary circulatory networks; the implementation and tuning of models of passive and active myocardial behaviour; as well as the testing of appropriate element types for coupling non–conforming fluid– solid finite element models under high interface tractions (finding that curvilinear spatial interpolations of the fluid geometry perform best). The behaviour of the resulting numerical scheme was investigated in a series of canonical test problems and found to be convergent and stable. The FD convergence studies also found that discontinuous pressure elements were better at capturing pressure gradients across FD boundaries. The ventricular simulations focused firstly on studying the passive diastolic behaviour of the LV both with and without LVAD support. Substantially different vortical flow features were observed when LVAD outflow was included. Additionally, a study of LVAD cannula lengths, using a particle tracking algorithm to determine recirculation rates of blood within the LV, found that shorter cannulas improved the recirculation of blood from the LV apex. Incorporating myocardial contraction, the model was extended to simulate the full cardiac cycle, converging on a repeating pressure–volume loop over 2 heart beats. Studies on the normal LV geometry found that LVAD implementation restricts the recirculation of early diastolic inflow, and that fluid–solid coupled models introduce greater heterogeneity of myocardial work than was observed in equivalent solid only models. A patient study was undertaken using a myocardial geometry constructed using image data from an LVAD implant recipient. A series of different LVAD flow regimes were tested. It was found that the opening of the aortic valve had a homogenising effect on the spatial variation of work, indicating that the synchronisation of LVAD outflow with the cardiac cycle is more important if the valve remains shut. Additionally, increasing LVAD outflow during systole and decreasing it during diastole led to improved mixing of blood in the ventricular cavity – compared with either the inverse, or holding outflow constant. Validation of these findings has the potential to impact the treatment protocols of LVAD patients.

APA, Harvard, Vancouver, ISO, and other styles

46

Djati, Nabil. "Study of interface capturing methods for two-phase flows." Thesis, Lyon, 2017. http://www.theses.fr/2017LYSEI052/document.

Full text

Abstract:

Cette thèse est consacrée au développement et à la comparaison des méthodes de suivi d'interface pour les écoulements diphasiques incompressibles. Elle s'intéresse à la sélection de méthodes robustes de suivi d'interface, puis à leur couplage avec le solveur des équations de Navier-Stokes. La méthode level-set est en premier lieu étudiée, en particulier l'influence du schéma d'advection et de l'étape de réinitialisation sur la qualité des résultats du suivi d'interface. Il a été montré que la méthode de réinitialisation avec contrainte de volume est robuste et précise en combinaison avec des schémas conservatifs WENO d'ordre 5 pour l'advection. Il a été constaté que les erreurs du suivi d'interface augmentent de manière abrupte lorsque la condition CFL est trop petite. Comme remède, la réinitialisation du champ level-set effectuée moins souvent réduit la diffusion numérique et le déplacement non-physique de l'interface. La conservation de la masse n'est pas assurée avec les méthodes level-set. Les méthodes VOF (volume-of-fluid) qui conservent naturellement la masse du fluide de référence sont alors étudiées. Une résolution géométrique avec un schéma consistent et conservatif est alors adoptée, ainsi qu'une autre technique alternative plus aisément extensible en 3D. Il a été trouvé que ces deux dernières méthodes donnent des résultats très proches. La méthode MOF (moment-of-fluid), qui reconstruit l'interface en utilisant le centre de masse du fluide de référence, est plus précise que les méthodes VOF. Différentes méthodes couplées entre level-set et VOF sont alors étudiées, notamment: CLSVOF, MCLS, VOSET et CLSMOF. Il a été observé que la méthode level-set tend à épaissir les filaments minces, tandis que VOF et les méthodes couplées les fragmentent en petites particules. Finalement, on a couplé les méthodes level-set et VOF avec le solveur incompressible des équations de Navier-Stokes. On a comparé différentes manières de prise en compte des conditions de saut à l'interface (lisse et raide). Il a été montré que les méthodes VOF sont plus robustes, et donnent d'excellents résultats pour quasiment toutes les simulations. Deux méthodes level-set donnant de très bons résultats, comparables à ceux de VOF, sont aussi identifiées This thesis is devoted to the development and comparison of interface methods for incompressible two-phase flows. It focuses on the selection of robust interface capturing methods, then on the manner of their coupling with the Navier-stokes solver. The level-set method is first investigated, in particular the influence of the advection scheme and the reinitialization step on the accuracy of the interface capturing. It is shown that the volume constraint method for reinitialization is robust and accurate in combination with the conservative fifth-order WENO schemes for the advection. It is found that interface errors increase drastically when the CFL number is very small. As a remedy, reinitializing the level-set field less often reduces the amount of numerical diffusion and non-physical interface displacement. Mass conservation is, however, not guaranteed with the level-set methods. The volume-of-fluid (VOF) method is then investigated, which naturally conserves the mass of the reference fluid. A geometrical consistent and conservative scheme is adopted, then an alternative technique more easily extended to 3D. It is found that both methods give very similar results. The moment-of-fluid (MOF) method, which reconstructs the interface using the reference fluid centroid, is found to be more accurate than the VOF methods. Different coupled level-set and VOF methods are then investigated, namely: CLSVOF, MCLS, VOSET and CLSMOF. It is observed that the level-set method tends to thicken thin filaments, whereas the VOF and coupled methods break up thin structures in small fluid particles. Finally, we coupled the level-set and volume-of-fluid methods with the incompressible Navier-Stokes solver. We compared different manners (sharp and smoothed) of treating the interface jump conditions. It is shown that the VOF methods are more robust, and provide excellent results for almost all the performed simulations. Two level-set methods are also identified that give very good results, comparable to those obtained with the VOF methods

APA, Harvard, Vancouver, ISO, and other styles

47

Lissoni, Giulia. "Méthode DDFV : applications en mécanique des fluides et décomposition des domaines." Thesis, Université Côte d'Azur (ComUE), 2019. http://www.theses.fr/2019AZUR4060.

Full text

Abstract:

L’objectif de cette thèse est d'étudier et développer des schémas numériques du type volume finis pour des problèmes provenant de la mécanique des fluides, notamment le problème de Stokes et Navier-Stokes. Les schémas choisis sont du type dualité discrète, dénotés DDFV ; cette méthode travaille sur des grilles décalées, où les inconnus de vitesse sont placés aux centres des volumes de contrôle et aux sommets du maillage, et les inconnus de pression aux arêtes du maillage. Ce type de construction a deux avantages principaux : elle permet de considérer des maillages généraux (qui ne vérifient pas nécessairement la condition d’orthogonalité classique des maillages volumes finis) et de reconstruire à niveau discret les propriétés de dualité des opérateurs différentiels continus. On commence par l'étude de la discrétisation du problème de Stokes avec des conditions aux bords mixtes de type Dirichlet/Neumann ; le caractère bien posé de ce problème est strictement lié à l'inégalité Inf-sup, qui doit être vérifiée. Dans le cadre DDFV, cette inégalité a été prouvée pour des maillages particuliers ; on peut éviter cette hypothèse, en ajoutant des termes de stabilisation dans l’équation de conservation de masse. Dans un premier temps, on étudie un schéma stabilisé pour le problème de Stokes en forme de Laplace, en montrant son caractère bien posé, des estimations d'erreur et des tests numériques. On étudie ensuite le même problème en forme divergence, où le tenseur des contraintes remplace le gradient ; ici, on suppose que l'inégalité Inf-sup est vérifiée, et on écrit un schéma bien posé suivi des tests numériques. On considère ensuite le problème de Navier-Stokes incompressible. Initialement, on étude ce problème couplé avec des conditions aux bords « ouvertes » en sortie ; ce type de conditions apparaissent lors qu'on veut introduire une frontière artificielle, qui peut arriver pour des raisons de coût de calcul ou physiques. On écrit un schéma bien posé et des estimations d’énergie, validés par des simulations numériques. Deuxièmement, on s'intéresse à la méthode de décomposition de domaines sans recouvrement pour le problème de Navier-Stokes incompressible, en écrivant un algorithme de Schwarz discret. On discrétise le problème avec un schéma de type Euler semi-implicite en temps, et à chaque itération on applique l’algorithme de Schwarz au système linéaire résultant. Nous montrons également la convergence de cet algorithme et nous terminons par des expériences numériques. Cette thèse se termine par un cinquième chapitre issu d’une collaboration lors du CEMRACS 2019, où le but est d'étendre DPIR (une technique récente pour la reconstruction d'interfaces entre deux matériaux) au cas d'interfaces courbes et de trois matériaux. Des simulations numériques montrent les résultats The goal of this thesis is to study and develop numerical schemes of finite volume type for problems arising in fluid mechanics, namely Stokes and Navier-Stokes problems. The schemes we choosed are of discrete duality type, denoted by DDFV; this method works on staggered grids, where the velocity unknowns are located at the centers of control volumes and at the vertices of the mesh, and the pressure unknowns are on the edges of the mesh. This kind of construction has two main advantages: it allows to consider general meshes (that do not necessarily verify the classical ortogonality condition required by finite volume meshes) and to reconstruct and mimic at the discrete level the dual properties of the continuos differential operators. We start by the study of the discretization of Stokes problem with mixed boundary conditions of Dirichlet/Neumann type; the well-posed character of this problem is strictly relied to Inf-sup inequality, that has to be verified. In the DDFV setting, this inequality has been proven for particular meshes; we can avoid this hypothesis, by adding some stabilization terms in the equation of conservation of mass. In the first place, we study a stabilized scheme for Stokes problem in Laplace form, by showing its well-posedness, some error estimates and numerical tests. We study the same problem in divergence form, where the strain rate tensor replaces the gradient; here, we suppose that the Inf-sup inequality is verified, and we design a well-posed scheme followed by some numerical tests. We consider then the incompressible Navier-Stokes problem. At first, we study this problem coupled with « open » boundary conditions on the outflow; this kind of conditions arises when an artificial boundary is introduced, to save computational ressources or for physical reasons. We write a well-posed scheme and some energy estimates, validated by numerical simulations. Secondly, we address the domain decomposition method without overlap for the incompressible Navier-Stokes problem, by writing a Schwarz algorithm. We discretize the problem with a semi-implicit Euler scheme in time, and at each time iteration we apply Schwarz algorithm to the resulting linear system. We show the convergence of this algorithm and we end by some numerical experiments. This thesis ends with a last chapter concerning the work done during CEMRACS 2019, where the goal is to extend DPIR (a recent technique for interface reconstruction between two materials) to the case of curved interfaces and of three materials. Some numerical simulations show the results

APA, Harvard, Vancouver, ISO, and other styles

48

Paillere, Henri J. "Multidimensional upwind residual distribution schemes for the Euler and Navier-Stokes equations on unstructured grids." Doctoral thesis, Universite Libre de Bruxelles, 1995. http://hdl.handle.net/2013/ULB-DIPOT:oai:dipot.ulb.ac.be:2013/212553.

Full text

Abstract:

Une approche multidimensionelle pour la résolution numérique des équations d'Euler et de Navier-Stokes sur maillages non-structurés est proposée. Dans une première partie, un exposé complet des schémas de distribution, dits de "fluctuation-splitting" ,est décrit, comprenant une étude comparative des schémas décentrés, positifs et de 2ème ordre, pour résoudre l'équation de convection à coefficients constants, ainsi qu'une étude théorique et numérique de la précision des schémas sur maillages réguliers et distordus. L'extension à des lois de conservation non-linéaires est aussi abordée, et une attention particulière est portée au problème de la linéarisation conservative. Dans une deuxième partie, diverses discrétisations des termes visqueux pour l'équation de convection-diffusion sont développées, avec pour but de déterminer l'approche qui offre le meilleur compromis entre précision et coût. L'extension de la méthode aux systèmes des lois de conservation, et en particulier à celui des équations d'Euler de la dynamique des gaz, représente le noyau principal de la thèse, et est abordée dans la troisième partie. Contrairement aux schémas de distribution classiques, qui reposent sur une extension formelle du cas scalaire, l'approche développée ici repose sur une décomposition du résidu par élément en équations scalaires, modélisant le transport de variables caracteristiques. La difficulté vient du fait que les équations d'Euler instationnaires ne se diagonalisent pas, et admettent une infinité de solutions élémentaires (ondes simples) se propageant dans toutes les directions d'espace. En régime stationnaire, en revanche, les équations se diagonalisent complètement dans le cas des écoulements supersoniques, et partiellement dans le cas des écoulements subsoniques. Ainsi, les équations sous forme conservative peuvent être remplacées par un système équivalent comprenant deux équations totalement découplées, exprimant l'invariance de l'entropie et de l'enthalpie totale le long des lignes de courant, et deux autres équations, modélisant les effets purement acoustiques. En régime supersonique, celles-ci se découplent aussi, et expriment la convection le long des lignes de Mach d'invariants de Riemann généralisés. La discrétisation de ces équations par des schémas scalaires décentrés permet de simuler des écoulements continus et discontinus avec une grande précision et sans oscillations. Finalement, dans une dernière partie, l'extension aux équations de Navier-Stokes est abordée, et la discrétisation des termes visqueux par une approche éléments finis est proposée. Les résultats numériques confirment la précision et la robustesse de la méthode. Doctorat en sciences appliquées info:eu-repo/semantics/nonPublished

APA, Harvard, Vancouver, ISO, and other styles

49

Ben-Aïm, Laurence. "Applications des methodes particulaires en mecanique des fluides et en physique des plasmas." Paris 6, 1988. http://www.theses.fr/1988PA066060.

Full text

Abstract:

Etude des possibilites numeriques offertes par les methodes particulaires. La methode repose sur la possibilite de recreer l'evolution d'une quantite (continue) a partir de mouvements ponctuels. Application aux problemes d'ecoulement en mecanique des fluides et a celui du mouvement des electrons dans un plasma soumis a une difference de potentiel

APA, Harvard, Vancouver, ISO, and other styles

50

Huang, Yu-Lin. "Ecoulement d'un film de biomolécules tensioactives : expériences et modélisation." Phd thesis, Université Joseph Fourier (Grenoble), 2009. http://tel.archives-ouvertes.fr/tel-00429430.

Full text

Abstract:

L'objectif de la thèse est de contrôler, à partir d'un écoulement recirculant, la densification en 2-D de biomolécules amphiphiles à une surface liquide. Pour cela, une expérience basée sur une cuve annulaire est développée et l'on étudie le comportement d'une monocouche d'acide pentadecanoïque (PDA) soumise à des cisaillements volumique et surfacique. Une organisation mésoscopique avec deux phases 2-D (liquide expansée et liquide condensée) est mise en évidence par microscopie Brewster. A partir de la fraction aréolaire de phase condensée, on montre que l'écoulement centrifuge le long du fond tournant engendre une compression radiale centripète le long de l'interface et compacte les domaines condensés vers les petits rayons. Pour un niveau suffisant de centrifugation, le film de Langmuir expérimente une transition morphologique due à une forte atomisation des domaines condensés et pilotée par l'équilibre entre cisaillement de surface et tension de ligne effective. A vitesse de rotation élevée (nombre de Reynolds de l'ordre de 1000 ou plus), lorsque la phase condensée est organisée sous forme de films stratifiés, le cisaillement volumique engendre une fragmentation, en bon accord avec la littérature. Finalement, une formulation intégrale de l'écoulement annulaire cisaillé est proposée afin d'estimer au mieux le nombre de Boussinesq et la stratification visqueuse.

APA, Harvard, Vancouver, ISO, and other styles

Dissertations / Theses on the topic 'Stokes Law (Fluid mechanics)'

Create a spot-on reference in APA, MLA, Chicago, Harvard, and other styles