Although the flow of granular material in silos and the pressure acting on the silo walls have been studied for over a century, many challenges still remain in silo design. In particular, during the discharge process some dynamic phenomena in silos can often be observed to display large, self-induced and dynamic pulsations which may endanger the stability of the silo structure. The aim of this thesis is to study the flow and pressure in silos using numerical modelling and analytical methods, and to further understand the mechanical behaviour of granular material and mechanism of dynamic phenomena during silo discharge. The Finite Element (FE) method can be used to analyse the behaviour of the granular material in silos by considering the material as a continuum. In this thesis, FEM modelling of silo flow was developed using the Arbitrary Lagrangian-Eulerian (ALE) formulation in the Abaqus/Explicit program and the key parameters that affect the predictions of the flow and pressure during discharge were identified. Using the ALE technique, almost the entire silo discharge process can be simulated without mesh distortion problems. The mass flow rate and temporally averaged discharge pressure predicted by the FE model were first investigated in a conical hopper and were found to be in good agreement with those from the most commonly quoted theoretical solutions. The transient dynamic pressure fluctuations during incipient silo discharge were predicted and the causes for these dynamic events have been investigated which led to the conclusion that the stress wave propagation and the moving shear zone phenomena within the bulk solid were responsible for the dominant higher and lower frequencies effects respectively. A one-dimensional dynamic model of granular columns subject to Coulomb wall friction was developed to investigate the propagation of stress waves, focusing on the effect of geometry by examining converging and diverging tapered columns. The analytical solutions of this model are compared to the FE model based on the ALE formulation. This FE model was first validated using the known behaviour for cylindrical columns. In all cases, the stress impulse set off by incipient discharge at the silo outlet grew with the distance travelled up the column, however the rate was shown to depend on the halfangle of the taper. Over a range of small angles, the proposed analytical model was found to accurately predict this behaviour. After the successful application of the ALE technique for a conical hopper, the FE model was extended to simulate the granular flow in a flat-bottomed model silo. The FE predictions were compared with the silo pressure measurements in a model silo (Rotter et al, 2004). Pressure cells mounted along a vertical line on the silo walls were used to measure the pressure distribution in the silo tests using dry sand. The FE model was further extended to simulate the granular flow in a model silo consisting of a cylindrical section with a conical hopper. The prediction was compared with the experimental observations from a model silo (Munch-Andersen et al, 1992), together with the well-known theoretical solutions. Two numerical issues were addressed in some detail: one is the numerical treatment of the abrupt transition between the cylinder section and the conical hopper, the other is the interaction between the granular solid and the silo walls that was modelled using a dynamic friction model. In addition, the dynamic pressure events during discharge were examined and plausible explanations were given. Finally, this thesis deployed a non-coaxial elastoplastic constitutive model to explore the effect of non-coaxiality on silo phenomena. The non-coaxial FE modelling was performed on three problems: a simple shear test under various initial conditions, a steep hopper and a flat-bottomed silo. The results show that non-coaxiality did not influence the prediction of wall pressure during filling and storing, on the other hand, the discharge pressure was predicted to be larger when non-coaxiality is considered.

No description available.

No description available.

Le premier objectif de cette thèse est de générer et d'analyser une base de données pour l'écoulement intracardiaque dans des géométries réalistes. Dans ce but, une stratégie couplant simulation numérique et imagerie médicale est appliquée à un cœur gauche pathologique et à un cœur gauche sain. Le second objectif est d'illustrer comment cette base de données peut être analysée afin de mieux comprendre l'écoulement intracardiaque, en portant une attention particulière aux caractéristiques instationnaires de l'écoulement et à sa nature turbulente. Une chaîne numérique pour simuler l'écoulement dans des géométries spécifiques au patient est tout d'abord présentée. La cavité cardiaque et ses mouvements sont extraits à partir d'images médicales à l'aide d'un algorithme de recalage d'image afin d'obtenir le domaine de calcul. Les équations qui régissent l'écoulement sont écrites dans le cadre d'un maillage se déformant au cours du temps (approche arbitrairement Lagrangienne ou Eulérienne). Les valves cardiaques sont modélisées à l'aide de frontières immergées. L'application de cette chaîne numérique à deux cœurs gauches, l'un pathologique, l'autre sain est ensuite détaillée. L'écoulement sanguin est caractérisé par sa nature transitoire, donnant un écoulement complexe et cyclique. Il est montré que l'écoulement n'est ni laminaire, ni pleinement turbulent, justifiant a posteriori l'utilisation de simulation aux grandes échelles. Le développement instationnaire de la turbulence est analysé à l'aide de l'écoulement moyenné sur un nombre suffisant de cycles cardiaques. Les statistiques de l'écoulement, l'énergie turbulente, la production de turbulence et une analyse spectrale sont notamment présentées. Une étude Lagrangienne est aussi effectuée en utilisant des statistiques calculées à l'aide de particules ensemencées dans l'écoulement. En plus des caractéristiques habituellement rapportées, ce travail met en évidence le caractère perturbé et transitoire de l'écoulement, tout en identifiant les mécanismes de production de la turbulence. / The first objective of this thesis is to generate and analyse CFD-based databases for the intracardiac flow in realistic geometries. To this aim, an image-based CFD strategy is applied to both a pathological and a healthy human left hearts. The second objective is to illustrate how the numerical database can be analysed in order to gain insight about the intracardiac flow, mainly focusing on the unsteady and turbulent features. A numerical framework allowing insight in fluid dynamics inside patient-specific human hearts is first presented. The heart cavities and their wall dynamics are extracted from medical images, with the help of an image registration algorithm, in order to obtain a patient-specific moving numerical domain. Flow equations are written on a conformal moving computational domain, using an Arbitrary Lagrangian-Eulerian framework. Valves are modelled using immersed boundaries.Application of this framework to compute flow and turbulence statistics in both a realistic pathological and a realistic healthy human left hearts is presented. The blood flow is characterized by its transitional nature, resulting in a complex cyclic flow. Flow dynamics is analysed in order to reveal the main fluid phenomena and to obtain insights into the physiological patterns commonly detected. It is demonstrated that the flow is neither laminar nor fully turbulent, thus justifying a posteriori the use of Large Eddy Simulation.The unsteady development of turbulence is analysed from the phase averaged flow, flow statistics, the turbulent stresses, the turbulent kinetic energy, its production and through spectral analysis. A Lagrangian analysis is also presented using Lagrangian particles to gather statistical flow data. In addition to a number of classically reported features on the left heart flow, this work reveals how disturbed and transitional the flow is and describes the mechanisms of turbulence production.

Dans certaines simulations numériques exigeantes de mécanique des fluides, ilest nécessaire de simuler des écoulements multiphasiques impliquant de nombreuses contraintes simultanées : nombre de fluides important, évolutions compressibles à la fois isentropes et fortement choquées, équations d’états variables et contrastées, déformations importantes et transport surdes longues distances. Afin de remplir ces objectifs de manière robuste, il est nécessaire que la cohérence thermodynamique du schéma numérique soit vérifiée.Dans le premier chapitre, un schéma de type Lagrange plus projection est proposé pour la simulation d’écoulements diphasiques avec un modèle squelette à six équations et sans termes de dissipation. L’importance de la propriété de préservation des écoulements isentropiques est mise en évidence à l’aide d’une comparaison avec des résultats issus de la littérature pour le test deRansom. Ce chapitre souligne aussi certaines limitations de l’approche Lagrange plus projection pour simuler des modèles multiphasiques.Afin de pallier à ces limitations, une nouvelle procédure de dérivation est proposée afin de construire un schéma mimétique pour la simulation d’écoulements instationnaires compressibles dans un formalisme ALE direct (Arbitrary Lagrangian–Eulerian). La possibilité de choisir a prioriles degrés de liberté permet de s’inscrire dans une continuité avec les schémas historiques décalés, tout en imposant les conservations au niveau discret. L’équation de quantité de mouvement discrèteest obtenue par application d’un principe variationnel, assurant par construction la cohérence thermodynamique des efforts de pression. Cette approche est appliquée au cas d’écoulements monofluides comme preuve de concept au Chapitre 3, puis elle est étendue au cas d’écoulements à Nphasescompressibles au Chapitre 4. Des tests mono et multiphasiques montrent un comportement satisfaisant en terme de conservativité, versatilité aux mouvements de grilles et robustesse. / In some highly demanding fluid dynamics simulations, it appears necessary tosimulate multiphase flows involving numerous constraints at the same time : large numbers of fluids, both isentropic and strongly shocked compressible evolution, highly variable and contrasted equations of state, large deformations, and transport over large distances. Fulfilling such a challengein a robust and tractable way demands that thermodynamic consistency of the numerical scheme be carefully ensured.In the first chapter, a Lagrange plus remap scheme is proposed for the simulation of two-phase flows with a dissipation-free six-equation bakcbone model. The importance of the property of isentropic flow preservation is highlighted with a comparison with Ransom test results fromthe literature. This chapter also also point out certain limitations of the Lagrange plus remap approach for multiphase simulations.In order to overcome these limitations, a novel derivation procedure is proposed to construct a mimetic scheme for the simulation of unsteady and compressible flows in a direct ALE (ArbitraryLagrangian-Eulerian) formalism. The possibility to choose a priori the degrees of freedom allows to obtain a continuity with historical staggered scheme, while imposing conservativity at discretelevel. The discrete momentum evolution equation is obtained by application of a variational principle, thus natively ensuring the thermodynamic consistency of pressure efforts. This approach is applied to single-fluid flows as a proof of concept in Chapter 3, then it is extended to N-phasecompressible flows in Chapter 4. Single- and multi-phase tests show satisfactory behavior in terms on conservation, versatility to grid motions, and robustness.

Ce sujet de thèse s’inscrit dans le cadre des études actuellement menées au CEA/DAM concernant des schémas numériques ALE (Arbitrary-Lagrangian-Eulerian)de type « Lagrange + Projection », dans le contexte des simulations hydrodynamiques mutli-matériaux en grandes déformations. Ces schémas doivent respecter les équations de conservation de la masse, de la quantité de mouvement et de l’énergie totale.Les schémas décalés en temps et en espace sont très utilisés dans les codes industriels. Ils sont robustes et permettent une bonne approximation des comportements complexes, mais sont connus pour ne pas conserver exactement l’énergie totale. Cela pose un problème dans le traitement des chocs, sur maillages raffinés ou dans la simulation des milieux réactifs.En 2016, des travaux originaux on été proposés par A. Llor et. al. pour rendre conservatif ce type de schéma dans un contexte lagrangien (sans projection), notamment en proposant une correction pour retrouver la conservation de l’énergie totale.Le travail de cette thèse a été d’étendre ce schéma lagrangien dans un contexte ALE multi-matériaux (avec interface), en garantissant la conservation de toutes les quantités, le respect du second principe de la thermodynamique et la robustesse. De nombreux cas tests ont été menés (en 2D plan et en 2D axisymétrique) et comparés aux méthodes existantes afin de montrer la pertinence de cette approche. / This PhD subject comes within actual studies managed by CEA/DAM about ALE (Arbitrary-Lagrangian-Eulerian) schemes (with a splitting of Lagrangian and Remapping steps) in the context of hydrodynamic simulations. These numerical schemes have to respect mass, momentum and total energy conservation, which are the fundamental equations of the studied systems.Space- and Time-Staggered are widely used in industrial codes for their simplicity androbustness despite their known lack of exact energy conservation. This is a major drawbackin presence of strong shocks. Among all existing schemes, none of them meet the expectations of robustness, conservation,thermodynamic consistency (both shocks and relaxations capture), accuracy andadaptibility to complex behaviors. Recently, some novel works have been proposed by A.Llor et. al. in order to make conservative this type of scheme in a Lagrangian context (without remapping step). Current remap methods, necessary in large deformations, donot guarantee simultaneously total energy conservation and thermodynamic consistency.This work aims at extending this conservative Lagrangian space- and time-staggeredscheme to a multi-material ALE methodology, keeping its good properties (conservation,accuracy, thermodynamic consistency, robustness) intact. Classical, but demanding, test cases have been performed (both in plane and axisymmetric 2D geometries) and have been compared to existing numerical methods in order to assess the relevance of our approach.

Underwater explosions are often modeled with Arbitrary Lagrangian-Eulerian (ALE) Finite Element Methods. The objective of this thesis is to validate the simulation method, with respect to the propagating shock wave. A two-dimensional axisymmetric model of a spherical TNT charge submerged in water is simulated using LS-DYNA. The explosive is modeled with the Burn Fraction technique and the Jones-Wilkins-Lee equation of state. Water is modeled as a non-viscous fluid, with the Grüneisen equation of state. The convergence for different mesh resolutions, the effect of different advection methods, and varied constants in the artificial viscosity are examined. Generally, the simulations agree well with empirical results, but the maximum pressure diminishes more rapidly with distance compared to experiments. The excessive dampening is most notable in the early stages of the propagation. Also, unexpected oscillations are observed near the discontinuity. The choice of advection scheme and constants in the artificial viscosity do not resolve the issues suggesting that other numerical techniques for treating the discontinuity should be considered. / Undervattensexplosioner simuleras ofta med ALE-baserade finita elementmetoder. Detta examensarbete avser att validera simuleringsmetoden med hänsyn till stötvågens utbredning i vattnet. En tvådimensionell axisymmetrisk modell av en sfärisk TNT-laddning nedsänkt i vatten simuleras i LS-DYNA. Laddningen modelleras med hjälp av brinnfraktioner och Jones-Wilkins-Lee tillståndsekvation. Vattnet modelleras som en inviskös fluid tillsammans med Grüneisens tillståndsekvation. Nätkonvergens, val av advektionsmetod och ändring av konstanter i den artificiella viskositeten studeras. Övergripande resultat stämmer väl överens med empirisk data, men stötvågens topptryck avtar fortare än väntat. Denna dämpning är tydligast i utredningens tidiga skeden. Dessutom observeras oväntade oscillationer kring stötvågens diskontinuerliga tryckprofil. Val av advektionsmetod och konstanter tillhörande artificiella viskositeten verkar ha liten betydelse för resultaten. En alternativ numerisk metod för behandling av stötvågens diskontinuitet bör implementeras.

A general framework is described for the solution of two-phase fluid-object interaction problems on the basis of coupling a distributed-Lagrange-multiplier fictitious domain method and a level-set method, intended for application to the problem of seabed scour by ice ridges. The resulting equations are discretized in space using stabilized finite-element methods and integrated in time using the generalized-α method. This approach is simple to implement and applicable to both structured and unstructured meshes in two and three dimensions. By means of examples, it is shown that despite the simplicity of the approach, good results are obtained in comparison with other more computationally demanding methods. A robust approach is utilized for constructing signed-distance functions on arbitrary meshes by introducing artificial numerical diffusivity to improve the robustness of classical signed-distance construction approaches without resorting to common pseudo-time relaxation. Under this approach, signed-distance functions can be rapidly constructed while preserving the numerical convergence properties and, generally, having minimal interfacial perturbation. The method is then applied with a modified deformation procedure for fast and efficient mesh adaptivity, including a discussion how it may be used in computational fluid dynamics. The two-phase fluid-object interaction approach is then customized for modeling of the seabed scour and soil-pipe interaction. In this approach, complex history-dependent soil constitutive models are replaced with a simple strain-rate dependent model. Utilization of this constitutive model along with the framework developed earlier leads to the treatment of seabed scour as a two-phase fluid-object interaction, and the soil-pipe interaction as a fluid-structure interaction problem without the need for remeshing. Good agreement with past experimental and numerical studies are obtained using our approach. The dissertation is concluded by conducting a parametric study of seabed scour in two- and three-dimensions. / text

N/A / Thesis / Master of Applied Science (MASc)

<p>Finite element analysis (FEA) employs a science-based approach in which the complete machining process can be simulated and optimized before resorting to costly and time-consuming experimental trials. In this work, cutting coefficient of AISI 1045 steel will be estimated using finite element modelling using Arbitrary Lagrangian Eulerian formulation (ALE). The estimated values are then experimentally validated. A parametric study is carried out after in order to investigate how some cutting parameters can affect the cutting coefficients. The process parameters to be varied include feed rate, cutting speed, and cutting edge radius.</p> / Master of Applied Science (MASc)