A combination of theoretical, numerical, and experimental methods were utilized to determine that shock waves can form in brain tissue from smooth boundary conditions. The conditions that lead to the formation of shock waves were determined. The implication of this finding was that the high gradients of stress and strain that could occur at the shock wave front could contribute to mechanism of brain injury in blast loading conditions. The approach consisted of three major steps. In the first step, a viscoelastic constitutive model of bovine brain tissue under finite step-and-hold uniaxial compression with 10 1/s ramp rate and 20 s hold time has been developed. The assumption of quasi-linear viscoelasticity (QLV) was validated for strain levels of up to 35%. A generalized Rivlin model was used for the isochoric part of the deformation and it was shown that at least three terms (C_10, C_01 and C_11) are needed to accurately capture the material behavior. Furthermore, for the volumetric deformation, a linear bulk modulus model was used and the extent of material incompressibility was studied. The hyperelastic material parameters were determined through extracting and fitting to two isochronous curves (0.06 s and 14 s) approximating the instantaneous and steady-state elastic responses. Viscoelastic relaxation was characterized at five decay rates (100, 10, 1, 0.1, 0 1/s) and the results in compression and their extrapolation to tension were compared against previous models. In the next step, a framework for understanding the propagation of stress waves in brain tissue under blast loading was developed. It was shown that tissue nonlinearity and rate dependence are key parameters in predicting the mechanical behavior under such loadings, as they determine whether traveling waves could become steeper and eventually evolve into shock discontinuities. To investigate this phenomenon, the QLV material model developed based on finite compression results mentioned above was extended to blast loading rates, by utilizing the stress data published on finite torsion of brain tissue at high rates (up to 700 1/s). It was shown that development of shock waves is possible inside the head in response to compressive pressure waves from blast explosions. Furthermore, it was argued that injury to the nervous tissue at the microstructural level could be attributed to the high stress and strain gradients with high temporal rates generated at the shock front and this was proposed as a mechanism of injury in brain tissue. In the final step, the phenomenon of shock wave formation and propagation in brain tissue was further studied by developing a one-dimensional model of brain tissue using the Discontinuous Galerkin finite element method. This model is capable of capturing high-gradient waves with higher accuracy than commercial finite element software. The deformation of brain tissue was investigated under displacement input and pressure input boundary conditions relevant to blast over-pressure reported in the literature. It was shown that a continuous wave can become a shock wave as it propagates in the tissue when the initial changes in acceleration are beyond a certain limit. The high spatial gradients of stress and strain at the shock front cause large relative motions at the cellular scale at high temporal rates even when the maximum strains and stresses are relatively low. This gradient-induced local deformation occurs away from the boundary and can therefore contribute to the diffuse nature of blast-induced injuries.   / Mechanical Engineering

This work is devoted to the numerical solution of equations describing the dynamics of flocks of birds. Specifically, we pay attention to the Euler equations for compressible flow with a right-hand side correction. This model is based on the work Fornasier et al. (2010). Due to the complexity of the model, we focus only on the one-dimensional case. For the numerical solution we use a semi-implicit discontinuous Galerkin method. Discretization of the right-hand side is chosen so that we preserve the structure of the semi-implicit scheme for the Euler equations presented in the work Feistauer, Kučera (2007). The proposed numerical scheme was implemented and numerical experiments showing the robustness of the scheme were carried out. Powered by TCPDF (www.tcpdf.org)

The presented work deals with the discontinuous Galerkin method with the anisotropic mesh adaptation for stationary convection-diffusion problems. Basic definitions are included in an introduction where we also present the used method. The following parts describe various methods for evaluating a Riemann metric, which is necessary for anisotropic mesh adaptation. The most important part of work follows - numerical experiments carried out with ADGFEM and ANGENER software packages. In these experiments, we compare different approaches for the definition of Riemann metrics and compare their efficiency. The main output of this thesis are subroutines for evaluation of the Riemann metric including its source code.

Cette thèse vise à étudier le couplage entre méthodes de Galerkine discontinue (DG) et méthodes de différences finies (DF) en maillages hybrides non structuré / cartésien, en vue d'applications en aéroacoustique numérique. L'idée d'une telle approche consiste à pouvoir tirer profit localement des avantages respectifs de ces méthodes, soit, en d'autres termes, à pouvoir prendre en compte la présence de géométries complexes par une méthode DG en maillage non structuré, et les zones qui en sont suffisamment éloignées par une méthode DF en maillage cartésien, moins coûteuse. Plus précisément, il s'agit de concevoir un algorithme d'hybridation de ces deux types de schémas pour l'approximation des équations d'Euler linéarisées, puis d'évaluer avec attention le comportement numérique des solutions qui en sont issues. De par le fait qu'aucun résultat théorique ne semble actuellement atteignable dans un cas général, cette étude est principalement fondée sur une démarche d'expérimentation numérique. Par ailleurs, l'intérêt d'une telle hybridation est illustré par son application à un calcul de propagation acoustique dans un cas réaliste / This thesis aims at studying coupling techniques between Discontinuous Galerkin (DG) and finite difference (FD) schemes in a non-structured / Cartesian hybrid-mesh context,in the framework of Aeroacoustics computations. The idea behind such an approach is the possibility to locally take advantage of the qualities of each method. In other words, the goal is to be able to deal with complex geometries using a DG scheme on a non-structured mesh in their neighborhood, while solving the rest of the domain using a FD scheme on a cartesian grid, in order to alleviate the needs in computational resources. More precisely, this work aims at designing an hybridization algorithm between these two types of numerical schemes, in the framework of the approximation of the solutions of the Linearized Euler Equations. Then, the numerical behaviour of hybrid solutions is cautiously evaluated. Due to the fact that no theoretical result seems achievable at the present time, this study is mainly based on numerical experiments. What's more, the interest of such an hybridization is illustrated by its application to an acoustic propagation computation in a realistic case

Estuda-se o problema de equilíbrio sem força de corpo de uma esfera anisotrópica sob compressão radial uniformemente distribuída sobre o seu contorno no contexto da teoria da elasticidade linear clássica. A solução deste problema prediz o fenômeno inaceitável da auto-intersecção em uma região próxima ao centro da esfera para uma dada faixa de parâmetros materiais. Sob o contexto de uma teoria de minimização do funcional de energia potencial total da elasticidade linear clássica com a restrição de que o determinante do gradiente da função mudança de configuração seja injetivo, este fenômeno é eliminado. Aplicam-se duas formulações do Método dos Elementos Finitos de Galerkin Descontínuo (MEFGD) para obter soluções aproximadas para o problema de equilíbrio da esfera sem restrição. A primeira formulação do MEFGD aproxima diretamente os campos de deslocamento e deformação infinitesimal. A consideração do campo adicional de deformação na formulação do MEFGD aumenta o número de graus de liberdade associados aos nós da malha de elementos finitos e, consequentemente, o custo computacional. Com o objetivo de reduzir o número de graus de liberdade, introduz-se neste trabalho uma formulação alternativa do MEFGD. Nesta formulação, o campo de deformação infinitesimal não é obtido diretamente da inversão do sistema de equações resultante, mas sim por pós-processamento, a partir do campo de deslocamento aproximado. As soluções aproximadas obtidas com ambas as formulações do MEFGD são comparadas com a solução exata do problema sem restrição e com soluções aproximadas obtidas com o Método dos Elementos Finitos de Galerkin Clássico (MEFGC). Ambas as formulações do MEFGD fornecem melhores aproximações para a solução exata do que as aproximações obtidas com o MEFGC. Os erros entre a solução exata e as soluções aproximadas obtidas com a formulação alternativa do MEFGD são um pouco maiores do que os erros correspondentes obtidos com a formulação original do MEFGD. Este aumento nos erros é compensado pelo menor esforço computacional exigido pela formulação alternativa. Este trabalho serve de base para o estudo de problemas com restrição de injetividade utilizando o método de Galerkin descontínuo. / The equilibrium problem without body force of an anisotropic sphere under radial compression that is uniformly distributed on the sphere\'s boundary is investigated in the context of the classical linear elasticity theory. The solution of this problem predicts the unacceptable phenomenon of self-intersection in a vicinity of the center of the sphere for a given range of material parameters. This phenomenon can be eliminated in the context of a theory that minimizes the total potential energy of classical linear elasticity subjected to the restriction that the deformation field be injective. Two formulations of the Finite Element Method using Discontinuous Galerkin (MEFGD) are used to obtain approximate solutions for the unconstrained problem. The first formulation of the MEFGD approximates both the displacement and the strain fields. The consideration of the strain as an additional field in the formulation of the MEFGD increases the number of degrees of freedom associated to the finite elements and, therefore, the computational cost. With the objective of reducing the number of degrees of freedom, an alternative formulation of the MEFGD is introduced in this work. In this formulation, the strain field is not obtained directly from the inversion of the resulting linear system of equations, but from a post-processing calculation using the approximate displacement field. The approximate solutions obtained with both formulations of the MEFGD are compared with the exact solution of the problem without restriction and with approximate solutions obtained with the Finite Element Method using Classical Galerkin (MEFGC). Both formulations of the MEFGD yield better approximations for the exact solution than the approximations obtained with the MEFGC. The errors between the exact solution and the approximate solutions obtained with the alternative formulation of the MEFGD are slightly higher than the corresponding errors obtained with the original formulation of the MEFGD. These errors are compensated by the fact that the alternative formulation requires less computational effort than the computational effort required by the original formulation. This work serves as a basis for the study of problems with the injectivity restriction using the discontinuous Galerkin method.

This thesis is focused on direct numerical simulation (DNS) of compressible and incompressible fully developed and developing turbulent flows between isothermal walls using a discontinuous Galerkin method (DGM). Three cases (Ma = 0.2, 0.7 and 1.5) of DNS of turbulent channel flows between isothermal walls with Re ~ 2800, based on bulk velocity and half channel width, have been carried out. It is found that a power law seems to scale mean streamwise velocity with Ma slightly better than the more usual log-law. Inner and outer scaling of second-order and higher-order statistics have been analyzed. The linkage between the pressure gradient and vorticity flux on the wall has been theoretically derived and confirmed and they are highly correlated very close to the wall. The correlation coefficients are influenced by Ma, and viscosity when Ma is high. The near-wall spanwise streak spacing increases with Ma. Isosurfaces of the second invariant of the velocity gradient tensor are more sparsely distributed and elongated as Ma increases. DNS of turbulent isothermal-wall bounded flow subjected to favourable and adverse pressure gradient (FPG, APG) at Ma ~ 0.2 and Reref ~ 428000, based on the inlet bulk velocity and the streamwise length of the bottom wall, is also investigated. The FPG/APG is obtained by imposing a concave/convex curvature on the top wall of a plane channel. The flows on the bottom and top walls are tripped turbulent and laminar boundary layers, respectively. It is observed that the first and second order statistics are strongly influenced by the pressure gradients. The cross-correlation coefficients of the pressure gradients and vorticity flux remain constant across the FPG/APG regions of the flat wall. High correlations between the streamwise/wallnormal pressure gradient and the spanwise vorticity are found near the separation region close to the curved top wall. The angle of inclined hairpin structure to streamwise direction of the bottom wall is smaller (flatter) in the FPG region than the APG region. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2009-12-21 13:59:53.084

Die Dynamik der Atmosphäre der Erde umfasst einen Bereich von mikrophysikalischer Turbulenz über konvektive Prozesse und Wolkenbildung bis zu planetaren Wellenmustern. Für Wettervorhersage und zur Betrachtung des Klimas über Jahrzehnte und Jahrhunderte ist diese Gegenstand der Modellierung mit numerischen Verfahren. Mit voranschreitender Entwicklung der Rechentechnik sind Neuentwicklungen der dynamischen Kerne von Klimamodellen, die mit der feiner werdenden Auflösung auch entsprechende Prozesse auflösen können, notwendig. Der dynamische Kern eines Modells besteht in der Umsetzung (Diskretisierung) der grundlegenden dynamischen Gleichungen für die Entwicklung von Masse, Energie und Impuls, so dass sie mit Computern numerisch gelöst werden können. Die vorliegende Arbeit untersucht die Eignung eines unstetigen Galerkin-Verfahrens niedriger Ordnung für atmosphärische Anwendungen. Diese Eignung für Gleichungen mit Wirkungen von externen Kräften wie Erdanziehungskraft und Corioliskraft ist aus der Theorie nicht selbstverständlich. Es werden nötige Anpassungen beschrieben, die das Verfahren stabilisieren, ohne sogenannte „slope limiter” einzusetzen. Für das unmodifizierte Verfahren wird belegt, dass es nicht geeignet ist, atmosphärische Gleichgewichte stabil darzustellen. Das entwickelte stabilisierte Modell reproduziert eine Reihe von Standard-Testfällen der atmosphärischen Dynamik mit Euler- und Flachwassergleichungen in einem weiten Bereich von räumlichen und zeitlichen Skalen. Die Lösung der thermischen Windgleichung entlang der mit den Isobaren identischen charakteristischen Kurven liefert atmosphärische Gleichgewichtszustände mit durch vorgegebenem Grundstrom einstellbarer Neigung zu(barotropen und baroklinen)Instabilitäten, die für die Entwicklung von Zyklonen wesentlich sind. Im Gegensatz zu früheren Arbeiten sind diese Zustände direkt im z-System(Höhe in Metern)definiert und müssen nicht aus Druckkoordinaten übertragen werden.Mit diesen Zuständen, sowohl als Referenzzustand, von dem lediglich die Abweichungen numerisch betrachtet werden, und insbesondere auch als Startzustand, der einer kleinen Störung unterliegt, werden verschiedene Studien der Simulation von barotroper und barokliner Instabilität durchgeführt. Hervorzuheben ist dabei die durch die Formulierung von Grundströmen mit einstellbarer Baroklinität ermöglichte simulationsgestützte Studie des Grades der baroklinen Instabilität verschiedener Wellenlängen in Abhängigkeit von statischer Stabilität und vertikalem Windgradient als Entsprechung zu Stabilitätskarten aus theoretischen Betrachtungen in der Literatur. / The dynamics of the Earth’s atmosphere encompass a range from microphysical turbulence over convective processes and cloud formation up to planetary wave patterns. For weather forecasting and the investigation of climate over decades and centuries, these are subject to modelling with numerical methods. With progressing development of computer technology, re-development of the dynamical cores of climate models is in order to properly handle processes covered by the increasing resolution. The dynamical core of a model consists of the adaptation(discretization)of the basic equations for the dynamics of mass, energy and momentum for solving them numerically employing computers. The presented work investigates the applicability of a low-order Discontinuous Galerkin (DG) method for atmospheric applications. With equations that include external forces like gravitation and the Coriolis force, that is not given by theory. Necessary changes for stabilizing the method without resorting to slope limiters are presented. For the unmodified method, the basic inability to properly keep atmospheric balances is demonstrated. The developed stabilized model reproduces a set of standard test cases in a wide range of spatial and temporal scales. The solution of the termal wind equation along its characteristics curves, those being identical to the isobars, produces balanced atmospheric states with tunable (barotropic and baroclinic) instability via a prescribed zonal wind field. The constructed instability directly relates to the generation of cyclones. In contrast to earlier works, these balanced states are directly given in the z system (height in meters), without need for elaborate conversion from pressure coordinates. With these constructed states, both as reference state, the deviations from which being considered numerically, and as especially as initial condition subject to a small perturbation, several studies of barotropic and baroclinic instability are conducted via simulations. Particularily, the construction of steady states with configurable zonal flows of certain baroclinity facilitates a simulation-based study of baroclinic instability of differing wavelengths, depending on static stability and vertical wind gradient, in correspondence with stability maps from theoretical considerations in the literature.

Cette thèse se concentre sur le développement d'une méthode précise et efficace pour la simulation des grandes échelles (LES) des écoulements turbulents. Une approche de la LES basée sur la méthode variationnelle multi-échelles (VMS) est considérée. La VMS applique aux équations de la dynamique des fluides une séparation d'échelles a priori sans recours à des hypothèses sur les conditions aux limites ou sur l'uniformité du maillage. Afin d'assurer effectivement une séparation d'échelles dans l'espace des nombres d'onde associé, nous choisissons d'utiliser les ondelettes de deuxième génération (SGW), une base polynomiale qui présente des propriétés de localisation spatiale-fréquence optimales. A partir de la séparation d'échelles ainsi réalisée, l'action du modèle sous-maille est limitée à un intervalle de nombres d'onde proche de la coupure spectrale. Cette approche VMS-LES basée sur les ondelettes est désignée par WAVVMS-LES. Elle est incorporée dans un solveur d'ordre élevé pour la simulation des écoulements incompressibles sur la base d'une méthode de Galerkin discontinue (DG-FEM) stabilisée pour la pression. La méthode est évaluée par réalisation de LES sur des maillages fortement sous-résolus pour le cas test du tourbillon de Taylor-Green 3D à deux nombres de Reynolds différents. / This thesis focuses on the development of an accurate and efficient method for performing Large-Eddy Simulation (LES) of turbulent flows. An LES approach based upon the Variational Multiscale (VMS) method is considered. VMS produces an a priori scale-separation of the governing equations, in a manner which makes no assumptions on the boundary conditions and mesh uniformity. In order to ensure that scale-separation in wavenumber is achieved, we have chosen to make use of the Second Generation Wavelets (SGW), a polynomial basis which exhibits optimal space-frequency localisation properties. Once scale-separation has been achieved, the action of the subgrid model is restricted to the wavenumber band closest to the cutoff. We call this approach wavelet-based VMS-LES (WAV-VMS-LES). This approach has been incorporated within the framework of a high-order incompressible flow solver based upon pressure-stabilised discontinuous Galerkin FEM (DG-FEM). The method has been assessed by performing highly under-resolved LES upon the 3D Taylor-Green Vortex test case at two different Reynolds numbers.

L'étude de la susceptibilité électromagnétique de circuits électroniques nécessite l'utilisation d'un outil de simulation rapide, précis et suffisamment flexible pour intégrer les dernières innovations technologiques. La méthode itérative basée sur le concept d'onde (notée WCIP pour Wave Concept Iterative Procedure) initialement proposée par H. Baudrand est particulièrement adaptée pour la modélisation numérique de circuits multicouches à plusieurs niveaux de métallisation. Pour ce type de circuits, elle se révèle être l'une des méthodes qui utilise le plus petit nombre d'inconnues pour atteindre une précision donnée. Néanmoins, la WCIP n'est pas adaptée à la prise en compte des diélectriques inhomogènes et des trous d'interconnexion. L'objectif de la thèse est de s'affranchir de ces limitations par un couplage avec des méthodes numériques volumiques. En premier lieu, l'hybridation a été mise en œuvre avec une méthode basée sur la théorie des lignes de transmission pour des raisons de correspondance de maillages. Par la suite, le couplage avec une technique d'éléments finis de type Galerkin Discontinu (notée GD) Hybridée permet d'atteindre des objectifs de précision et de rapidité car GD apporte une flexibilité dans la discrétisation. En effet, c'est une méthode d'éléments finis non conforme qui permet notamment de faire varier d'un élément à l'autre l'ordre polynomial d'approximation. On a ainsi développé une nouvelle méthode numérique hybride couplant la WCIP avec des méthodes volumiques qui offrent plus de souplesse pour la prise en compte des milieux complexes. Enfin, une stratégie de résolution par décomposition de domaines est également abordée à la fin du manuscrit. / Electromagnetic susceptibility study of electronic circuits requires the use of a simulation tool which is fast, accurate and flexible enough to incorporate the latest technological innovations. The Wave Concept Iterative Procedure (WCIP) initially proposed by H. Baudrand is particularly adapted for numerical modeling of multilayered circuits with multilevel metallization. For this kind of circuits, it turns out to be one of the methods that uses the smallest number of unknowns to reach a given accuracy. However, the WCIP is not appropriate for inhomogeneous dielectric substrates and metallized via holes. The aim of this PhD thesis is to overcome these limitations coupling the WCIP with volume numerical methods. First, hybridization is carried out with the Frequency Domain Transmission Line Matrix (denoted FDTLM) assuming matching meshes at the interface between computational domains of both methods. Subsequently, the coupling with a finite element technique like a Hybridized Discontinuous Galerkin (denoted DG) method is considered to achieve the objectives of accuracy and speed because DG brings flexibility in the discretization. Indeed, it is a nonconforming finite element method which allows in particular changing the polynomial approximation order from one element to another. Therefore, a new hybrid method is developed coupling the WCIP with volume numerical methods which offer more flexibility for dealing with complex environments. Finally, a domain decomposition solution strategy is also discussed at the end of the manuscript.

Cette thèse se penche sur la propagation d’ondes au sein du système coupléTerre-océan-atmosphère. La compréhension de ces phénomènes a une importance majeure pour l’étude de perturbations sismiques et d’explosions atmosphériques notamment dans le cadre de missions spatiales planétaires. Les formes d’ondes issues du couplage fluide-solide permettent d’obtenir de précieuses informations sur la source du signal ou les propriétés des milieux de propagation. On développe donc deux outils numériques d’ordre élevé pour la propagation d’ondes acoustiques et de gravité. L'u en différences finies et se concentre sur le milieu atmosphérique et la propagation d’ondes linéaires dans un milieu stratifié visqueux et avec du vent. Cette méthode linéaire est validée par des solutions quasi-analytiques reposant sur les équations de dispersion dans une atmosphère stratifiée. Elle est aussi appliquée à deux cas d’études : la propagation d’ondes liée à l’impact d’une météorite à la surface de Mars (mission NASA INSIGHT), et la propagation d’ondes atmosphériques liées au tsunami de Sumatra en 2004. La seconde méthode résout la propagation non-linéaire d’ondes gravito-acoustiques dans une atmosphère couplée, avec topographie, à la propagation d’ondes élastiques dans un solide visco-élastique. Cette méthode repose sur sur le couplage d’une formulation en éléments finis discontinus, pour résoudre les équations de Navier-Stokes la partie fluide, par éléments finis continus pour résoudre les équations de l’élastodynamique dans la partie solide. Elle a été validée grâce à des solutions analytiques ainsi que par des comparaisons avec les résultats de la méthode par différences finies. / This thesis deals with the wave propagation problem within the Earth-ocean-atmosphere coupled system. A good understanding of the these phenomena has a major importance for seismic and atmospheric explosion studies, especially for planetary missions. Atmospheric wave-forms generated by explosions or surface oscillations can bring valuable information about the source mechanism or the properties of the various propagation media. We develop two new numerical full-wave high-order modeling tools to model the propagation of acoustic and gravity waves in realistic atmospheres. The first one relies on a high-order staggered finite difference method and focus only on the atmosphere. It enables the simultaneous propagation of linear acoustic and gravity waves in stratified viscous and windy atmosphere. This method is validated against quasi-analytical solutions based on the dispersion equations for a stratified atmosphere. It has also been employed to investigate two cases : the atmospheric propagation generated by a meteor impact on Mars for the INSIGHT NASA mission and for the study of tsunami-induced acoutic and gravity waves following the 2004 Sumatra tsunami. The second numerical method resolves the non-linear acoustic and gravity wave propagation in a realistic atmosphere coupled, with topography, to the elastic wave propagation in a visco-elastic solid. This numerical tool relies on a discontinuous Galerkin method to solve the full Navier-Stokes equations in the fluid domain and a continuous Galerkin method to solve the elastodynamics equations in the solid domain. It is validated against analytical solutions and numerical results provided by the finite-difference method.