The first ocean general circulation models developed in the late sixties were based on finite differences schemes on structured grids. Many improvements in the fields of engineering have been achieved since three decades with the developments of new numerical methods based on unstructured meshes. Some components of the first models may now seem out of date and new second generation models are therefore under study, with the aim of taking advantage of the potential of modern numerical techniques such as finite elements. In particular, unstructured meshes are believed to be more efficient to resolve the large range of time and space scales present in the ocean. Besides the classical continuous finite element or finite volume methods, another popular new trend in engineering applications is the Discontinuous Galerkin (DG) method, i.e. discontinuous finite elements presenting many interesting numerical properties in terms of dispersion and dissipation, errors convergence rates, advection schemes, mesh adaptation, etc. The method is especially efficient at high polynomial orders. The motivation for this PhD research is therefore to investigate the use of the high-order DG method for geophysical flow modeling. A first part of the thesis is devoted to the mesh adaptation using the DG method. The inter-element jumps of the fields are used as error estimators. New mesh size fields or polynomial orders are then derived and local h- or p-adaptation is performed. The technique is applied to standard benchmarks and computations in more realistic domains as the Gulf of Mexico. A second part deals with the use of the high order DG method with high-order representation of geometrical features. On one hand, a method is proposed to deal with complex representations of the coastlines. Computations are performed using high-order mappings around the Rattray island, located in the Great Barier Reef. Numerical results are then compared to in-situ measurements. On the other hand, a new method is proposed to deal with curved manifolds in order to represents oceanic or atmospheric flows on the sphere. The approach is based on the use of a local high-order non-orthogonal basis, and is equivalent to the use of vectorial shape and test functions to represent the vectorial conservation laws on the manifold's surface. A method is finally proposed to analyze the dispersion and dissipation properties of any numerical scheme on any kind of grid, possibly unstructured. The DG method is then compared to other techniques as the mixed non-conforming linear elements, and the impact of unstructured meshes is studied.

L’adaptation des surfaces pour des fonctions prédéterminées par le choix des matériaux métalliques ou des couches minces ayant des propriétés mécaniques avancées peut potentiellement permettre de réaliser des nouvelles applications à petites échelles. Concevoir de telles applications utilisant des nouveaux matériaux nécessite en premier lieu la connaissance des propriétés mécaniques des matériaux ciblés à l’échelle microscopique et nanoscopique. Une méthode souvent appliquée pour caractériser les matériaux à petites échelles est la nanoindentation, qui peut être vue comme une mesure de dureté à l’échelle nanoscopique. Ce travail présente une contribution relative à l'interprétation des résultats de la nanoindentation, qui fait intervenir un grand nombre de phénomènes physiques couplés à l'aide de simulations numériques. A cette fin une approche interdisciplinaire, adaptée aux phénomènes apparaissant à petites échelles, et située à l’intersection entre la physique, la mécanique et la science des matériaux a été utilisée. Des modèles numériques de la nanoindentation ont été conçus à l'échelle atomique (modèle discret) et à l'échelle des milieux continus (méthode des éléments finis), pour étudier le comportement du nickel pur. Ce matériau a été choisi pour ses propriétés mécaniques avancées, sa résistance à l'usure et sa bio-compatibilité, qui peuvent permettre des applications futures intéressantes à l'échelle nanoscopique, particulièrement dans le domaine biomédical. Des méthodes avancées de mécanique du solide ont été utilisées pour prendre en compte les grandes déformations locales du matériau (par la formulation corotationelle), et pour décrire les conditions de contact qui évoluent au cours de l'analyse dans le modèle à l'échelle des milieux continus (traitement des conditions de contact unilatérales et tangentielles par une forme de Lagrangien augmenté). L’application des modèles numériques a permis de contribuer à l’identification des phénomènes qui gouvernent la nanoindentation du nickel pur. Le comportement viscoplastique du nickel pur pendant nanoindentation a été identifié dans une étude expérimentale-numérique couplée, et l'effet cumulatif de la rugosité et du frottement sur la dispersion des résultats de la nanoindentation a été montré par une étude numérique (dont les résultats sont en accord avec des tendances expérimentales). Par ailleurs, l’utilisation de l’outil numérique pour une autre application à petites échelles, la manipulation des objets par contact, a contribué à la compréhension de la variation de l’adhésion électrostatique pendant micromanipulation. La déformation plastique des aspérités de surface sur le bras de manipulateur (en nickel pur) a été identifiée comme une source potentielle d’augmentation importante de l'adhésion pendant la micromanipulation, qui peut potentiellement causer des problèmes de relâche et de précision de positionnement, observés expérimentalement. Les résultats présentés dans cette thèse montrent que des simulations numériques basées sur la physique du problème traité peuvent expliquer des tendances expérimentales et contribuer à la compréhension et l'interprétation d'essais couramment utilisé pour la caractérisation aux petites échelles. Le travail réalisé dans cette thèse s’inscrit dans un projet de recherche appelé "mini-micro-nano" (mµn), financé par la Communauté Française de Belgique dans le cadre de "l'Action de Recherche Concertée", convention 04/09-310.

Intervertebral disc (IVD) degeneration is primarily involved in back pain, a morbidity that strongly affects the quality of life of individuals nowadays. Lumbar IVDs undergo stressful mechanical loads while being the largest avascular tissues in our body: Mechanical principles alone cannot unravel the intricate phenomena that occur at the cellular scale which are fundamental for the IVD regeneration. The present work aimed at coupling biomechanical and relevant molecular transport processes for disc cells to provide a mechanobiological finite element framework for a deeper understanding of degenerative processes and the planning of regenerative strategies. Given the importance of fluid flow within the IVD, the influence of poroelastic parameters such as permeabilities and solid-phase stiffness of the IVD subtissues was explored. A continuum porohyperelastic material model was then implemented. The angles of collagen fibers embedded in the annulus fibrosus (AF) were calibrated. The osmotic pressure of the central nucleus pulposus (NP) was also taken into account. In a parallel study of the human vertebral bone, microporomechanics was used together with experimental ultrasonic tests to characterize the stiffness of the solid matrix, and to provide estimates of poroelastic coefficients. Fluid dynamics analyses and microtomographic images were combined to understand the fluid exchanges at the bone-IVD interface. The porohyperelastic model of a lumbar IVD with poroelastic vertebral layers was coupled with a IVD transport model of three solutes - oxygen, lactate and glucose - interrelated to reproduce the glycolytic IVD metabolism. With such coupling it was possible to study the effect of deformations, fluid contents, solid-phase stiffness, permeabilities, pH, cell densities of IVD subtissues and NP osmotic pressure on the solute transport. Moreover, cell death governed by glucose deprivation and lactate accumulation was included to explore the mechanical effect on cell viability. Results showed that the stiffness of the AF had the most remarkable role on the poroelastic behavior of the IVD. The permeability of the thin cartilage endplate and the NP stiffness were also relevant. The porohyperelastic model was shown to reproduce the local AF mechanics, provided the fiber angles were calibrated regionally. Such back-calculation led to absolute values of fibers angles and to a global IVD poromechanical behavior in agreement with experiments in literature. The inclusion of osmotic pressure in the NP also led to stress values under confined compression comparable to those measured in healthy and degenerated NP specimens. For the solid bone matrix, axial and transverse stiffness coefficients found experimentally in the present work agreed with universal mass density-elasticity relationships, and combined with continuum microporomechanics provided poroelastic coefficients for undrained and drained cases. The effective permeability of the vertebral bony endplate calculated with fluid dynamics was highly correlated with the porosity measured in microtomographic images. The coupling of transport and porohyperelastic models revealed a mechanical effect acting under large volume changes and high compliance, favored by healthy rather than degenerated IVD properties. Such effect was attributed to strain-dependent diffusivities and diffusion distances and was shown to be beneficial for IVD cells due to the load-dependent increases of glucose levels. Cell density, NP osmotic pressure and porosity were the most important parameters affecting the coupled mechano-transport of metabolites. This novel study highlights the restoration of both cellular and mechanical factors and has a great potential impact for novel designs of treatments focused on tissue regeneration. It also provides methodological features that could be implemented in clinical image-based tools and improve the multiscale understanding of the human spine mechanobiology.

Many important factors affect the spread of childhood infectious disease. To understand better the fundamental drivers of infectious disease spread, several researchers have estimated seasonal transmission coefficients using discrete-time models. This research addresses several shortcomings of the discrete-time approaches, including removing the need for the reporting interval to match the serial interval of the disease using infectious disease data from three major cities: New York City, London, and Bangkok. Using a simultaneous approach for optimization of differential equation systems with a Radau collocation discretization scheme and total variation regularization for the transmission parameter profile, this research demonstrates that seasonal transmission parameters can be effectively estimated using continuous-time models. This research further correlates school holiday schedules with the transmission parameter for New York City and London where previous work has already been done, and demonstrates similar results for a relatively unstudied city in childhood infectious disease research, Bangkok, Thailand.

This study aims to develop a computer program to perform finite elements based fracture mechanics analyses of three dimensional surface cracks in T-welded joints of construction machinery. The geometrical complexity of the finite elements models and the requirement of large computer resources for the analyses necessitate the use of shell elements for general stress distribution optimization. A sub-modeling technique, together with a shell to solid conversion method, enables the user to model a local region and analyze it by defining the weld and crack parameters. It is assumed that the weld material is the same with the sheet metal material and the surface cracks are considered to occur on two weld toes and weld root. The surface cracks are assumed to have a semi elliptical crack front profile. In order to simulate the square-root strain singularity around the crack front, collapsed 20-node three dimensional brick elements are utilized. The rest of the local model is modeled by using 20-node three dimensional brick elements. The main objective of this work is to calculate the mixed mode energy release rates around the crack front for a sub-model of a global shell model by using J-integral method.

In an effort to achieve lighter and more economical space structures, a new patented steel composite member has been suggested and used in the construction of some steel roof structures. This special element has a sandwich construction composed of some strips of steel plates placed longitudinally along a helical wire core. The function of the helical core is to transfer the shear between the flange plates and increase the sectional inertia of the resulting composite member by keeping the flange plates at a desired distance from each other. Because of the lack of research, design engineers usually treat such elements as a solid member as if it has a full shear transfer between the flanges. However, a detailed analysis shows that this is not a valid assumption and leads to very unsafe results. In this context, the purpose of this study is to investigate the behavior of such members under axial compression and determine their effective sectional flexural rigidity by taking into account the shear deformations. This study applies an analytical investigation to a specific form of such elements with four flange plates placed symmetrically around a helical wire core. Five independent parameters of such a member are selected for this purpose. These are the spiral core and core wire diameters, the pitch of the spiral core, and the flange plate dimensions. Elements with varying combinations of the selected parameters are first analyzed in detail by finite element method, and some design charts are generated for the determination of the effective sectional properties to be used in the structural analysis and the buckling loads. For this purpose, an alternative closed-form approximate analytical solution is also suggested.

Failure of metallic structures operating under shock loading is a common occurrence in engineering applications. It is difficult to estimate the response of complicated systems analytically, due to structure&rsquo / s dynamic characteristics and varying loadings. Therefore, experimental, numerical or a combination of both methods are used for evaluations. The experimental analysis of the shocks due to firing is done for 12.7mm Gatling gun and 25mm cannon. During the tests, the Gatling gun and the cannon are located on military Stabilized Machine Gun Platform and Stabilized Cannon Platform respectively. For the firing tests, ICP (integrated circuit piezoelectric) accelerometers are attached to obtain the loading history for corresponding points. Shock Response Spectrum (SRS) analysis (nCode Glypworks) is done to define the equivalent shock profiles created on test pieces and the mount of 25mm cannon by means of the gun and the cannon firing. Transient shock analysis of the test pieces and the mount are done by applying the obtained shock profiles on the parts in a finite element model (ANSYS). Furthermore, experimental stress analysis due to shock loading is performed for two different types of material and different thicknesses of the test pieces. The input data for the analysis is obtained through measurements from strain rosette precisely located at the critical location of the test pieces. As a result of the thesis, a proposal is tried to be introduced where strain energy density theory is applied to predict the shock failure at military structures.

A novel and efficient methodology is developed for computer simulations of realistic two-dimensional (2D) and three-dimensional (3D) microstructures. The simulations incorporate realistic 2D and 3D complex morphologies/shapes, spatial patterns, anisotropy, volume fractions, and size distributions of the microstructural features statistically similar to those in the corresponding real microstructures. The methodology permits simulations of sufficiently large 2D as well as 3D microstructural windows that incorporate short-range (on the order of particle/feature size) as well as long-range (hundred times the particle/feature size) microstructural heterogeneities and spatial patterns at high resolution. The utility of the technique has been successfully demonstrated through its application to the 2D microstructures of the constituent particles in wrought Al-alloys, the 3D microstructure of discontinuously reinforced Al-alloy (DRA) composites containing SiC particles that have complex 3D shapes/morphologies and spatial clustering, and 3D microstructure of boron modified Ti-6Al-4V composites containing fine TiB whiskers and coarse primary TiB particles. The simulation parameters are correlated with the materials processing parameters (such as composition, particle size ratio, extrusion ratio, extrusion temperature, etc.), which enables the simulations of rational virtual 3D microstructures for the parametric studies on microstructure-properties relationships. The simulated microstructures have been implemented in the 3D finite-elements (FE)-based framework for simulations of micro-mechanical response and stress-strain curves. Finally, a new unbiased and assumption free dual-scale virtual cycloids probe for estimating surface area of 3D objects constructed by 2D serial section images is also presented.

In several practical applications, hot-finished steel pipe that exhibits Lüders bands is bent to strains of 2-3%. Lüders banding is a material instability that leads to inhomogeneous plastic deformation in the range of 1-4%. This work investigates the influence of Lüders banding on the inelastic response and stability of tubes under rotation controlled pure bending. It starts with the results of an experimental study involving tubes of several diameter-to-thickness ratios in the range of 33.2 to 14.7 and Lüders strains of 1.8% to 2.7%. In all cases, the initial elastic regime terminates at a local moment maximum and the local nucleation of narrow angled Lüders bands of higher strain on the tension and compression sides of the tube. As the rotation continues, the bands multiply and spread axially causing the affected zone to bend to a higher curvature while the rest of the tube is still at the curvature corresponding to the initial moment maximum. With further rotation of the ends, the higher curvature zone(s) gradually spreads while the moment remains essentially unchanged. For relatively low D/t tubes and/or short Lüders strains, the whole tube eventually is deformed to the higher curvature entering the usual hardening regime. Subsequently it continues to deform uniformly until the usual limit moment instability is reached. For high D/t tubes and/or materials with longer Lüders strains, the propagation of the larger curvature is interrupted by collapse when a critical length is Lüders deformed leaving behind part of the structure essentially undeformed. The higher the D/t and/or the longer the Lüders strain is, the shorter the critical length. This class of problems is analyzed using 3D finite elements while the material is modeled as an elastic-plastic solid with an “up-down-up” response over the extent of the Lüders strain, followed by hardening. The analysis reproduces the main features of the mechanical behavior provided the unstable part of the response is suitably calibrated. The uniform curvature elastic regime terminates with the nucleation of localized banded deformation. The bands appear in pockets on the most deformed sites of the tube and propagate into the hitherto intact part of the structure while the moment remains essentially unchanged. The Lüders-deformed section has a higher curvature, ovalizes more than the rest of the tube, and develops wrinkles with a characteristic wavelength. For every tube D/t there exists a threshold of Lüders strain separating the two types of behavior. This bounding value of Lüders strain was studied parametrically. / text

Term of remote-lab is certain types of laboratories which practical experiments are directedfrom a separate area by remote controller devices. This study is part of developing andupgrading advanced vibration remote laboratory. In the new remote lab, users have theability to measure the dynamic characteristics of the test object similar to the current existingremote lab. But in addition to current existing remote lab, they are capable to modifydynamic properties of the test object remotely by attaching vibration test instruments; such asa block of mass, spring-mass or non-linear spring. Doing several accurate experimental testsremotely on the test object are the toughest issues we faced as designers. In creating anddeveloping of this remote-lab, number of different approaches was adopted for producingwell-defined tests. Also, instead of implementing routine devices and techniques for regularvibration laboratories, the new prototypes were designed by finite elements method (FEM)and LABVIEW. For instance, the desirable test object, the attachment mechanism, usefulapplications, and proper software for managing via internet were prepared.