No description available.

The world population is expected to increase to approximately 11 billion by 2100. The ageing population (aged 60 and over) is projected to exceed the number of children in 2047. This will be a situation without precedent. The number of citizens with disorders of old age like Dementia will rise to 115 million worldwide by 2050. The estimated cost of Dementia will also increase, from $604 billion in 2010, to $1,117 billion by 2030. At the same time, medical expertise, evidence-driven policymaking and commissioning of services are increasingly evolving the definitive architecture of comprehensive long-term care to account for these changes. Technological advances, such as those provided by computational science and biomedical engineering, will allow for an expansion in our ability to model and simulate an almost limitless variety of complex problems that have long defied traditional methods of medical practice. Numerical methods and simulation offer the prospect of improved clinically relevant predictive information, and of course optimisation, enabling more efficient use of resources for designing treatment protocols, risk assessment and urgently needed management of a long term care system for a wide spectrum of brain disorders. Within this paradigm, the importance of the relationship of senescence of cerebrospinal fluid transport to dementia in the elderly make the cerebral environment notably worthy of investigation through numerical and computational modelling. Hydrocephalus can be succinctly described as the abnormal accumulation (imbalance between production and circulation) of cerebrospinal fluid (CSF) within the brain. Using hydrocephalus as a test bed, one is able to account for the necessary mechanisms involved in the interaction between cerebral fluid production, transport and drainage. The current state of knowledge about hydrocephalus, and more broadly integrative cerebral dynamics and its associated constitutive requirements, advocates that poroelastic theory provides a suitable framework to better understand the disease. In this work, Multiple-network poroelastic Theory (MPET) is used to develop a novel spatio-temporal model of fluid regulation and tissue displacement in various scales within the cerebral environment. The model is discretised in a variety of formats, through the established finite difference method, finite difference – finite volume coupling and also the finite element method. Both chronic and acute hydrocephalus was investigated in a variety of settings, and accompanied by emerging surgical techniques where appropriate. In the coupled finite difference – finite volume model, a key novelty was the amalgamation of anatomically accurate choroid plexuses with their feeding arteries and a simple relationship relaxing the constraint of a unique permeability for the CSF compartment. This was done in order to account for Aquaporin-4 sensitisation. This model is used to demonstrate the impact of aqueductal stenosis and fourth ventricle outlet obstruction. The implications of treating such a clinical condition with the aid of endoscopic third (ETV) and endoscopic fourth ventriculostomy (EFV) are considered. It was observed that CSF velocity in the aqueduct, along with ventricular displacement, CSF pressure, wall shear stress and pressure difference between lateral and fourth ventricles increased with applied stenosis. The application of ETV reduced the aqueductal velocity, ventricular displacement, CSF pressure, wall shear stress and pressure difference within nominal levels. The greatest reversal of the effects of atresia come by opting for ETV rather than the more complicated procedure of EFV. For the finite difference model incorporating nonlinear permeability, qualitatively similar results were obtained in comparison to the pertinent literature, however, there was an overall amplification of ventriculomegaly and transparenchymal pressure difference using this model. A quantitative and qualitative assessment is made of hydrocephalus cases involving aqueductal stenosis, along with the effects to CSF reabsorption in the parenchyma and subarachnoid space. The finite element discretisation template produced for the n^th- dimensional transient MPET system allowed for novel insight into hydrocephalus. In the 1D formulation, imposing the breakdown of the blood-CSF barrier responsible for clearance resulted in an increase in ventricular displacement, transparenchymal venous pressure gradient and transparenchymal CSF pressure gradient, whilst altering the compliance proved to markedly alter the rate of change of displacement and CSF pressure gradient. The influence of Poisson's ratio was investigated through the use of the dual-grid solver in order to distinguish between possible over or under prediction of the ventricular displacement. In the 2D model based on linear triangles, the importance of the MPET boundary conditions is acknowledged, along with the quality of the underlying mesh. Interesting results include that the fluid content is highest in the periventricular region and the skull, whilst after longer time scales, the peak CSF content becomes limited to the periventricular region. Venous fluid content is heavily influenced by the Biot-Willis constant, whilst both the venous and CSF/ISF compartments show to be strongly influenced by breakdown in the blood-CSF barrier. Increasing the venous compliance effects the arterial, capillary and venous compartments. Decreasing the venous compliance shows an accumulation of fluid, possibly helping to explain why the ventricles can be induced to compress rather than expand under decreased compliance. Finally, a successful application of the 3D-MPET template is shown for simple geometries. It is envisaged that future observations into the biology of cerebral fluid flow (such as perivascular CSF-ISF fluid exchange) and its interaction with the surrounding parenchyma, will demand the evolution of the MPET model to reach a level of complexity that could allow for an experimentally guided exploration of areas that would otherwise prove too intricate and intertwined under conventional settings.

Dans cette thèse, nous proposons un modèle pour le calcul des écoulements le long des failles. Ce modèle, baptisé modèle double interface permet de traiter deux difficultés majeures rencontrées lors de la modélisation des failles. Tout d'abord, l'utilisation d'un modèle interface, dans lequel les failles sont représentées par des éléments de dimension inférieure permet de s'affranchir du problème d'échelle spatiale. Ensuite, l'utilisation de deux interfaces pour représenter chaque faille permet de traiter naturellement les maillages non-conformes apparaissant dans ce type de problème. Les questions de failles non-planes et de réseaux de failles sont aussi abordées. Ce modèle est validé numériquement sur différents cas tests académiques et un cas synthétique inspiré du stockage du CO2 a aussi été réalisé. Finalement, une étude théorique a été menée afin de confirmer mathématiquement l'approche retenue. / In this thesis, we are interested in the modelisation of fluid flow along conductive faults. This model, so-called double interface model tackles two majors difficulties encountered when modelising faults. First of all, the use of an interface model, in which the faults are represented by lower dimension elements allows to treat the problem of space scale. Then, the use of two interfaces to modelise each fault allows to handle quite naturally the non-matching grid problem arising from this kind of problem. The question of non-planar fault and fault networks is also addressed. This model is then validated on several academic test cases and a synthetic case inspire by CO2 storage is also performed. Finally, a theoric study is also conducted in order to validate our approach.

Desde o importante trabalho de Barth e Frederickson (1990), um certo número de pesquisadores têm estudado o método de Volumes Finitos de alta-ordem k-exato, por exemplo o grupo do Prof. Ollivier-Gooch: Ollivier-Gooch e van Altena (2002), Nejat (2007), Michalak (2009), etc. Outras discretizações espaciais de alta-ordem bastante populares são o método Galerkin Descontínuo e o método de Diferença Espectral; processos iterativos que involucram estes esquemas tem sido acelerados, nos últimos anos, por métodos p-multigrid. Porém, esta aceleração não tem sido aplicada no contexto do método de Volumes Finitos de alta-ordem, pelo menos para conhecimento do autor desta tese. Por isso, o objetivo desta pesquisa é adaptar o p-multigrid desenvolvido por Liang et al. (2009b) no contexto da Diferença Espectral, para o ambiente dos Volumes Finitos estudado pelo Prof. Ollivier-Gooch. A pesquisa começa implementando o solver VF-RK, de Volumes Finitos com avanço Runge-Kutta, para resolver as equações de advecção-difusão e de Euler aplicados a problemas estacionários, por exemplo, o escoamento transônico ao redor do NACA 0012. Depois, estuda-se o método p-multigrid no contexto da Diferença Espectral; o p-multigrid acelera o processo iterativo comutando níveis polinomiais de alta e de baixa-ordem. Após esse estudo, a adaptação ao âmbito dos Volumes Finitos é realizada resultando num p-multigrid relativamente mais simples porque, em contraposição com o p-multigrid para Diferença Espectral, não precisa de operadores de restrição e prolongação para a comunicação entre diferentes níveis polinomiais. A pesquisa conclui com uma comparação com o método de Volumes Finitos de 4a ordem sem p-multigrid (solver VF-RK). Nesse sentido, implementa-se o solver pMG, baseado no p-multigrid proposto, para resolver os problemas estacionários considerados na primeira parte do trabalho; o smoother do p-multigrid é o esquema Runge-Kutta do código VF-RK, e cada problema estacionário é resolvido utilizando diferentes Vciclos procurando sempre soluções de 4a ordem. Os resultados indicam que o método p-multigrid proposto é mais eficiente que o método de Volumes Finitos de 4a ordem sem p-multigrid, isto é, os dois métodos oferecem a mesma precisão mas o primeiro pode levar menos de 50% do tempo de CPU do segundo. / Since Barth and Frederickson\'s important work (Barth e Frederickson, 1990), a number of researchers have studied high-order k-exact Finite Volume method, for example Prof. Ollivier-Gooch\'s group: Ollivier-Gooch e van Altena (2002), Nejat (2007), Michalak (2009), etc. Other quite popular high-order spatial discretizations are the Discontinuous Galerkin methods and the Spectral Difference methods; the iterative processes involving these schemes have been accelerated in recent years by p-multigrid methods. However, this acceleration has not been applied in the context of the high-order Finite Volume method, at least for the knowledge of the author of this thesis. Therefore, the objective of this research is to adapt the p-multigrid developed by Liang et al. (2009b) in the context of Spectral Difference methods, to the environment of Finite Volume studied by Prof. Ollivier-Gooch. This research begins by implementing the solver VF-RK, Finite Volume solver with Runge-Kutta advance, to compute the advection-diffusion equation and Euler equations applied to steady state problems, for example, the transonic flow around NACA 0012. Then, it is studied the p-multigrid method in the context of Spectral Difference schemes; p-multigrid accelerates the iterative process by switching polynomial levels of high- and low-order. After this study, the adaptation to the context of the Finite Volume scheme is performed resulting in a relatively simple p-multigrid because, in contrast to the p-multigrid for Spectral Difference schemes, it doesn\'t need restriction and prolongation operators for communication between different polynomial levels. The research concludes with a comparison with 4th order Finite Volume method without p-multigrid (solver VF-RK). Accordingly, the solver pMG, based on the proposed p-multigrid, is implemented to resolve the steady state problems considered in the first part of the work; the p-multigrid smoother is the Runge-Kutta scheme from VF-RK code, and each steady state problem is solved using different Vcycles, looking for 4th order solutions ever. The results indicate that the proposed p-multigrid method is more efficient than the 4th order Finite Volume method without p-multigrid: the two methods give the same accuracy but the first one can take less than 50% of second one\'s CPU time.

Este trabalho apresenta um estudo sobre a solução numérica da equação geral de advecção-difusão usando uma metodologia numérica conservativa. Para a discretização espacial, é usado o Método de Volumes Finitos devido à natureza conservativa da equação em questão. O método é configurado de modo a ter suas variáveis centradas em centro de célula e, para as variáveis, como a velocidade, centradas nas faces um método de interpolação de segunda ordem é utilizado para um ajuste numérico ao centro. Embora a implementação computacional tenha sido feita de forma paramétrica de maneira a acomodar outros esquemas numéricos, a discretização temporal dá ênfase ao Método de Crank-Nicolson. Tal método numérico, sendo ele implícito, dá origem a um sistema linear de equações que, aqui, é resolvido empregando-se o Método Multigrid-Multinível. A corretude do código implementado é verificada a partir de testes por soluções manufaturadas, de modo a checar se a ordem de convergência prevista em teoria é alcançada pelos métodos numéricos. Um jato laminar é simulado, com o acoplamento entre a equação de Navier-Stokes e a equação geral de advecção-difusão, em um domínio computacional tridimensional. O jato é uma forma de vericar se o algoritmo de geração de malhas adaptativas funciona corretamente. O módulo produzido neste trabalho é baseado no código computacional AMR3D-P desenvolvido pelos grupos de pesquisa do IME-USP e o MFLab/FEMEC-UFU (Laboratório de Dinâmica de Fluidos da Universidade Federal de Uberlândia). A linguagem FORTRAN é utilizada para o desenvolvimento da metodologia numérica e as simulações foram executadas nos computadores do LabMAP(Laboratório da Matemática Aplicada do IME-USP) e do MFLab/FEMEC-UFU. / This work presents a study about the numerical solution of variable coecients advectiondi usion equation, or simply, general advection-diusion equation using a conservative numerical methodology. The Finite Volume Method is choosen as discretisation of the spatial domain because the conservative nature of the focused equation. This method is set up to have the scalar variable in a cell centered scheme and the vector quantities, such velocity, are face centered and they need a second order interpolation to get adjusted to the cell center. The computational code is parametric, in which, any implicit temporal discretisation can be choosen, but the emphasis relies on Crank-Nicolson method, a well-known second order method. The implicit nature of aforementioned method gives a linear system of equations which is solved here by the Multilevel-Multigrid method. The correctness of the computational code is checked by manufactured solution method used to inspect if the theoretical order of convergence is attained by the numerical methods. A laminar jet is simulated, coupling the Navier-Stokes equation and the general advection-diusion equation in a 3D computational domain. The jet is a good way to check the corectness of adaptative mesh renement algorithm. The module designed here is based in a previous implemented code AMR3D-P designed by IME-USP and MFLab/FEMEC-UFU (Fluid Dynamics Laboratory, Federal University of Uberlândia). The programming language used is FORTRAN and the simulations were run in LabMAP(Applied Mathematics Laboratoy at IME-USP) and MFLab/FEMEC-UFU computers.

Dans cette thèse, on présente une étude numérique de la convection forcée turbulente dans un dissipateur thermique muni d'une chicane transversale dans le by-pass. Le premier modèle est composé d’ailettes planes et le second consiste à ajouter des ailettes broches entre les ailettes planes. Les équations gouvernantes basées sur le modèle de turbulences k- SSt sont discrétisées et résolues par la méthode des volumes finis et l'algorithme SIMPLE. Les résultats dynamiques sont présentés en en termes de champs de vitesse, des profils de vitesse axiales dans des sections choisies ainsi que la perte de charge. L'étude thermique est présentée en terme de champs de température et de distribution du nombre de Nusselt. Un rapport entre les performances thermique et dynamique est présenté pour évaluer les différents dissipateurs thermiques. / In this thesis, we present an numerical study of turbulent forced convection in a heat sink provided with a transverse baffle in the bypass. The first model is composed of plates fins and the second consists of adding pin fins between the plates fins. The governing equations, based on the k- SSt turbulence model, are disscredized and solved by the finite volume method and the SIMPLE algorithm. Dynamic results are presented in terms of velocity fields, profiles of the axial velocities in selected sections and pressure drop. The thermal study is presented in terms of temperature fields and the distribution of Nusselt number. A ratio between the thermal and dynamic performances is presented to evaluate the different heat sinks.

The main objective of this dissertation is to present an adaptation of some finite volume methods used in the resolution of problems arising in sedimentation processes of flocculated suspensions (or sedimentation with compression).<br />This adaptation is based on the utilization of multiresolution techniques, originally designed to reduce the computational cost incurred in solving using high resolution schemes in the numerical solution of hyperbolic systems of conservation laws.

<p>The MUSCL-Hancock finite volume method with different slope limiters has been analyzed in the context of a fast transient flow problem. A derivation and analysis of the axial forces inside a pipe system due to a flow transient is also performed. </p><p>The following slope limiters were implemented and compared: MC, van Leer, van Albada, Minmod and Superbee. The comparison was based on the method's ability to calculate the forces due to a flow transient inside a pipe system.</p><p>The tests and comparisons in this thesis show that the MC, van Leer, van Albada and Minmod limiters behave very much the same for the flow transient problem. If one would rank these four limiters with respect to the numerical error, the order would be the one presented above, the MC limiter being the most accurate. The error the four limiters produce is mainly of diffusive nature and it is just the magnitude of the diffusion that seems to differ between the methods. One should also note that the workload rank of the four limiters is the same as the order presented above. The MC limiter being the least efficient of the four and the Minmod limiter the most efficient.</p><p>In most of the tests performed the Superbee limiter display a rather negative unpredictable behavior. For some relatively simple cases this particular approach shows big difficulties maintaining the dynamical properties of the force. However, the upside of the Superbee limiter is its remarkable ability to maintain the maximum value of the forces present in the pipe system, preventing underestimation of the maximum magnitude of the force.</p>

In this thesis study, Detached Eddy Simulation turbulence model is studied in two dimension mainly for flow over single element airfoils in high Reynolds numbers to gain experience with model before applying it to a three dimensional simulations. For this aim, Spalart-Allmaras and standard DES ,DES97, turbulence models are implemented to parallel, viscous, hybrid grid flow solver. The flow solver ,Set2d, is written in FORTRAN language. The Navier-Stokes equations are discretized by first order accurately cell centered finite volume method and solved explicitly by using Runge-Kutta dual time integration technique. Inviscid fluxes are computed using Roe flux difference splitting method. The numerical simulations are performed in parallel environment using domain decomposition and PVM library routines for inter-process communications. To take into account the effect of unsteadyness after the convergence is ensured by local time stepping technique for four order magnitude drop in density residual, global time stepping is applied for 20000 iterations. The solution algorithm is validated aganist the numerical and experimental studies for single element airfoils in subsonic and transonic flows. It is seen that Spalart-Allmaras and DES97 turbulence models give the same results in the non-seperated flows. Grey area is investigated by changing $C_{DES}$ coefficient. Modeled Stress Depletion which cause reduction of eddy viscosity is observed.

Conservation laws with uncertain initial and boundary conditions are approximated using a generalized polynomial chaos expansion approach where the solution is represented as a generalized Fourier series of stochastic basis functions, e.g. orthogonal polynomials or wavelets. The stochastic Galerkin method is used to project the governing partial differential equation onto the stochastic basis functions to obtain an extended deterministic system. The stochastic Galerkin and collocation methods are used to solve an advection-diffusion equation with uncertain viscosity. We investigate well-posedness, monotonicity and stability for the stochastic Galerkin system. High-order summation-by-parts operators and weak imposition of boundary conditions are used to prove stability. We investigate the impact of the total spatial operator on the convergence to steady-state.  Next we apply the stochastic Galerkin method to Burgers' equation with uncertain boundary conditions. An analysis of the truncated polynomial chaos system presents a qualitative description of the development of the solution over time. An analytical solution is derived and the true polynomial chaos coefficients are shown to be smooth, while the corresponding coefficients of the truncated stochastic Galerkin formulation are shown to be discontinuous. We discuss the problematic implications of the lack of known boundary data and possible ways of imposing stable and accurate boundary conditions. We present a new fully intrusive method for the Euler equations subject to uncertainty based on a Roe variable transformation. The Roe formulation saves computational cost compared to the formulation based on expansion of conservative variables. Moreover, it is more robust and can handle cases of supersonic flow, for which the conservative variable formulation fails to produce a bounded solution. A multiwavelet basis that can handle  discontinuities in a robust way is used. Finally, we investigate a two-phase flow problem. Based on regularity analysis of the generalized polynomial chaos coefficients, we present a hybrid method where solution regions of varying smoothness are coupled weakly through interfaces. In this way, we couple smooth solutions solved with high-order finite difference methods with non-smooth solutions solved for with shock-capturing methods.