The aim of this study is to create CFD tools and models capable of simulating pulsatile blood flow in abdominal aortic aneurysm (AAA) and stent graft. It helps to increase the current physiological understanding of rupture risk of AAA and stent graft fixation or migration. Firstly, in order to build a general solver for the AAA modeling with reasonable accuracy, a third/fourth order modified OCI scheme is originally developed for general numerical simulation. The modified OCI scheme has a wider cell Reynolds number limitation. This high order scheme performs well with general rectangular mesh for incompressible fluid. Second, a velocity based finite volume method is originally developed to calculate the stress field for solid in order to capture the transient changes of the blood vessel since the artery is a rubber like material. All one, two and three dimensional classical cases for solid are tested and good results are obtained. The velocity based finite volume method show good potential to calculate the stress field for solid and easy to blend with the finite volume fluid solver. It has been recognized that fluid structure interaction (FSI) is very crucial in biomechanics. In this regard, the velocity based finite volume method is then further developed for FSI application. A well known one dimensional piston problem is studied to understand the feasibility of the fluid structure coupling. The numerical prediction matches the analytical solution very well. The velocity based method introduces less numerical damping compared with a stagger method and a monolithic method. Finally, the work focuses on practical pulsatile boundary conditions, non-Newtonian blood viscous properties and bifurcating geometry, and provides an overview of the hemodynamic within the AAA model. A modified Womersley inlet and imbalance pressure outlet boundary conditions are originally used in this study. The Womersley inlet boundary represents better approximation for pulsatile flow compared with the parabolic inlet condition. Numerical results are presented providing comparison between different boundary conditions using different viscous models in both 2D and 3D aneurysms. Good agreement between the numerical predictions and the experimental data is achieved for 2D case. 3D stent models with different bifurcation angles are also tested. The Womersley inlet boundary condition improves the existing inlet conditions significantly and it can reduce the Aneurysm neck computation domain. The influence of the non-Newtonian model to the wall shear stress (WSS) and strain-rate is also studied. The non-Newtonian model tends to produce higher WSS at both proximal and distal end of the aneurysm as compared with the Newtonian model (both 2D and 3D cases). The computed strain-rate distribution at the centre of the aneurysm is different between these two models. The influence of imbalance outlet pressure at the iliac arteries to the blood flow is originally investigated. The imbalance outlet pressure boundary conditions affect the computed wall shear stress significantly near the bifurcation point. All the pulsatile Womersley inlet, non-Newtonian viscosity properties and the imbalance pressure outlet need to be considered in blood flow simulation of AAA.

It is well known that when a flexible or flexibly-mounted structure is placed perpendicular to the flow of a Newtonian fluid, it can oscillate due to the shedding of separated vortices at high Reynolds numbers. If the same flexible object is placed in non-Newtonian flows, however, the structure's response is still unknown. The main objective of this thesis is to introduce a new field of viscoelastic fluid-structure interactions by showing that the elastic instabilities that occur in the flow of viscoelastic fluids can drive the motion of a flexible structure placed in its path. Unlike Newtonian fluids, the flow of viscoelastic fluids can become unstable at infinitesimal Reynolds numbers due to the onset of a purely elastic flow instability. This instability occurs in the absence of nonlinear effects of fluid inertia and the Reynolds number of the flows studied here are in the order of 10-4. When such an elastic flow instability occurs in the vicinity of a flexible structure, the fluctuating fluid forces exerted on the structure grow large enough to cause a structural instability which in turn feeds back into the fluid resulting in a flow instability. Nonlinear periodic oscillations of the flexible structure are observed which have been found to be coupled to the time-dependent growth and decay of viscoelastic stresses in the wake of the structure. Presented in this thesis are the results of an investigation of the interaction occurring in the flow of a viscoelastic wormlike micelle solution past a flexible rectangular sheet. The structural geometries studied include: flexible sheet inclinations at 20°, 45° and 90° and flexible sheet widths of 5mm and 2.5mm. By varying the flow velocity, the response of the flexible sheet has been characterized in terms of amplitude and frequency of oscillations. Steady and dynamic shear rheology and filament stretching extensional rheology measurements are conducted in order to characterize the viscoelastic wormlike micelle solution. Bright field images show the deformation of the flexible sheet during an unstable oscillation while flow-induced birefringence images highlight the viscoleastic fluid stresses produced in the wake of the flexible sheet.

The main focus of this work is on the non-Newtonian effects on the inertial instabilities in shear flows. Both inelastic (Carreau) and elastic models (Oldroyd-B and FENE-P) have been employed to examine the main features of the non-Newtonian fluids; shear-thinning, shear-thickening and elasticity. Several classical configurations have been considered; flow past a circular cylinder, in a lid-driven cavity and in a channel. We have used a wide range of tools for linear stability analysis, modal, non-modal, energy and sensitivity analysis, to determine the instability mechanisms of the non-Newtonian flows and compare them with those of the Newtonian flows. Direct numerical simulations have been also used to prove the results obtained by the linear stability analysis. Significant modifications/alterations in the instability of the different flows have been observed under the action of the non-Newtonian effects. In general, shear-thinning/shear-thickening effects destabilize/stabilize the flow around the cylinder and in a lid driven cavity. Viscoelastic effects both stabilize and destabilize the channel flow depending on the ratio between the viscoelastic and flow time scales. The instability mechanism is just slightly modified in the cylinder flow whereas new instability mechanisms arise in the lid-driven cavity flow. We observe that the non-Newtonian effect can alter the inertial flow at both baseflow and perturbation level (e.g. Carreau fluid past a cylinder or in a lid driven cavity) or it may just affect the perturbations (e.g. Oldroyd-B fluid in channel). In all the flow cases studied, the modifications in the instability dynamics are shown to be strongly connected to the contribution of the different terms in the perturbation kinetic energy budget. / <p>QC 20140113</p>

The traditional first-order analysis in Euclidean spaces relies on the Sobolev spaces W1,p(Ω), where Ω ⊂ Rn is open and p ∈ [1, ∞].The Sobolev norm is then defined as the sum of Lp norms of a function and its distributional gradient.We generalize the notion of Sobolev spaces in two different ways. First, the underlying function norm will be replaced by the “norm” of a quasi-Banach function lattice. Second, we will investigate functions defined on an abstract metric measure space and that is why the distributional gradients need to be substituted. The thesis consists of two papers. The first one builds up the elementary theory of Newtonian spaces based on quasi-Banach function lattices. These lattices are complete linear spaces of measurable functions with a topology given by a quasinorm satisfying the lattice property. Newtonian spaces are first-order Sobolev-type spaces on abstract metric measure spaces, where the role of weak derivatives is passed on to upper gradients. Tools such asmoduli of curve families and the Sobolev capacity are developed, which allows us to study basic properties of the Newtonian functions.We will see that Newtonian spaces can be equivalently defined using the notion of weak upper gradients, which increases the number of techniques available to study these spaces. The absolute continuity of Newtonian functions along curves and the completeness of Newtonian spaces in this general setting are also established. The second paper in the thesis then continues with investigation of properties of Newtonian spaces based on quasi-Banach function lattices. The set of all weak upper gradients of a Newtonian function is of particular interest.We will prove that minimalweak upper gradients exist in this general setting.Assuming that Lebesgue’s differentiation theoremholds for the underlyingmetricmeasure space,wewill find a family of representation formulae. Furthermore, the connection between pointwise convergence of a sequence of Newtonian functions and its convergence in norm is studied.

No description available.

No description available.

[pt] Neste trabalho foi analisado o escoamento de fluidos não Newtonianos através de canais axisimétricos convergentes divergentes. A solução da conservação de massa e de conservação de momento foi obtida numericamente via volumes finitos utilizando o programa de computador Fluent. A equação constitutiva de fluidos Newtonianos generalizados foi utilizada para modelar o comportamento não Newtoniano, utilizando a equação constitucional de Shunk-Scriven para cálculo da viscosidade, que assume como sendo a média geométrica ponderada pelo classificador de escoamento R entre a viscosidade de cisalhamento e a viscosidade de extensão. Os resultados de perda de pressão e vazão são comparados com os resultados calculados pela relação simplificada proposta por Souza Mendes e Naccache, 2002 entre a perda de carga e vazão de fluidos viscoelásticos fluindo através do meio poroso, para analisar a sua performance. / [en] In this work, the flow of non-Newtonian fluids through axisimmetric convergingdiverging channels is analyzed. The solution of mass and momentum conservation equations is obtained numerically via finite volume technique using the Fluent software. The Generalized Newtonian Fluid constitutive equation was used to model the non- Newtonian fluid behavior, using the Shunk-Scriven model for the viscosity, where a weighted geometric mean by the flow classifier R between shear and extensional viscosities is assumed. The results of pressure drop and flow rate are compared to the ones predicted by a previously proposed simplified relation (Souza Mendes and Naccache, 2002) between pressure drop and flow rate, for viscoelastic fluids flow through porous media, in order to analyze its performance.

A novel method of calculating viscosity from molecular dynamics simulations is developed, benchmarked, and tested. The technique is a transient method which has the potential to reduce CPU requirements for many conditions. An initial sinusoidal velocity profile is overlaid upon the peculiar velocities of the individual molecules in an equilibrated simulation. The transient relaxation of this initial velocity profile is then compared to the corresponding analytical solution of the momentum equation by adjusting the viscosity-related parameters in the constitutive equation that relate the shear rate to the stress tensor. The newly developed Transient Molecular Dynamics (TMD) method was tested for a Lennard-Jones (LJ) fluid over a wide range of densities and temperatures. The simulated values were compared to an analytical solution of the boundary value problem for a Newtonian fluid. The resultant viscosities agreed well with those published for Equilibrium Molecular Dynamics (EMD) simulations up to a dimensionless density of 0.7. Application of a linear viscoelastic Maxwell constitutive equation was required to achieve good agreement at dimensionless densities greater than 0.7. When the Newtonian model is used for densities in the range of 0.1 to 0.3 and the Maxwell model is used for densities higher than 0.3, the TMD method was able to predict viscosities with an uncertainty of 10% or better. Application of the TMD method to multi-site molecules required the Jeffreys constitutive equation to adequately fit the simulation responses. TMD simulations were performed on model fluids representing n-butane, isobutane, n-hexane, water, methanol, and hexanol. Molecules with strong hydrogen bonding and Coulombic interactions agreed well with NEMD simulated values and experimental values. Simulated viscosities for nonpolar and larger molecules agreed with NEMD simulations at low to moderate densities, but deviated from these values at higher densities. These deviations are explainable in terms of potential model inaccuracies and the shear-rate dependence of both NEMD and TMD viscosity values. Results show that accurate viscosity predictions can be made for multi-site molecules as long as the shear-rate dependence of the viscosity is not too large or is adequately addressed.

No description available.

No description available.