Le contexte général de cette thèse est la Simulation Numérique Directe des écoulements diphasiques dilués anisothermes. Un accent particulier est mis sur la détermination précise de la dispersion des particules et du transfert de chaleur entre la phase porteuse et dispersée. Cette dernière est décrite à l’aide d’une approche Eulérienne aux moments : le Formalisme Eulérien Mésoscopique (FEM) [41, 123], récemment étendu aux écoulements anisothermes [78]. Le principal objectif de ce travail est de déterminer si ce formalisme est capable de prendre en compte de manière précise l’inertie dynamique et thermique des particules dans un écoulement turbulent, et particulièrement dans une configuration avec un gradient moyen. Le code de calcul utilisé est AVBP. La simulation numérique d’un spray dilué avec une approche Eulerienne soulève des questions supplémentaires sur les méthodes numériques et les modèles employés. Ainsi, les méthodes numériques spécifiques aux écoulements diphasiques implémentées dans AVBP [69, 103, 109] ont été testées et revisitées. L’objectif est de proposer une stratégie numérique précise et robuste qui résiste aux forts gradients de fraction volumique de particule provoqués par la concentration préférentielle [132], tout en limitant la diffusion numérique. Ces stratégies numériques sont comparées sur une série de cas tests de complexité croissante et des diagnostics pertinents sont proposés. Par exemple, les dissipations dues `a la physique et au numérique sont extraites des simulations et quantifiées. Le cas test du tourbillon en deux dimensions chargé en particules est suggéré comme une configuration simple pour mettre en évidence l’impact de l’inertie des particules sur leur champ de concentration et pour discriminer les stratégies numériques. Une solution analytique est aussi proposée pour ce cas dans la limite des faibles nombres de Stokes. Finalement, la stratégie numérique qui couple le schéma centré d’ordre élevé TTGC et une technique de stabilisation, aussi appelée viscosité artificielle, est celle qui fournit les meilleurs résultats en terme de précision et de robustesse. Les paramètres de viscosité artificielle (c'est-à-dire les senseurs) doivent néanmoins être bien choisis. Ensuite, la question des modèles nécessaires pour d´écrire correctement la dispersion des particules dans une configuration avec un gradient moyen est abordée. Pour ce faire, un des modèles RUM (appel´e AXISY-C), proposé par Masi [78] et implémenté dans AVBP par Sierra [120], est validé avec succès dans deux configurations: un jet plan diphasique anisotherme 2D et 3D. Contrairement aux anciens modèles RUM, les principales statistiques de la phase dispersée sont désormais bien prédites au centre et aux bords du jet. Finalement, l’impact de l’inertie thermique des particules sur leur température est étudié. Les résultats montrent un effet important de cette inertie sur les statistiques mettant en évidence la nécessité pour les approches numériques de prendre en compte ce phénomène. Ainsi, l’extension du FEM aux écoulements anisothermes, c’est-à-dire les flux de chaleur RUM (notés RUM HF), est implémentée dans AVBP. L’impact des RUM HF sur les statistiques de température des particules est ensuite évalué sur les configurations des jets 2D et 3D. Les champs Eulériens sont comparés à des solutions Lagrangiennes de référence calculées par B. Leveugle au CORIA et par E. Masi à l’IMFT pour les jets 2D et 3D, respectivement. Les résultats montrent que les RUM HF améliorent la prédiction des fluctuations de température mésoscopique, et dans une moindre mesure la température moyenne des particules en fonction de la configuration. Les statistiques Lagrangiennes sont retrouvées lorsque les RUM HF sont pris en compte alors que les résultats sont dégradés dans le cas contraire. / This work addresses the Direct Numerical Simulation of non-isothermal turbulent flows laden with solid particles in the dilute regime. The focus is set on the accurate prediction of heat transfer between phases and of particles dispersion. The dispersed phase is described by an Eulerian approach : the Mesoscopic Eulerian Formalism [41, 123], recently extended to non-isothermal flows [78]. The main objective of this work is to assess the ability of this formalism to accurately account for both dynamic and thermal inertia of particles in turbulent sheared flows. The CFD code used in this work is AVBP. The numerical simulation of dilute sprays with an Eulerian approach calls for specific modelling and raises additional numerical issues. First, the numerical methods implemented in AVBP for two-phase flows [69, 103, 109] were tested and revisited. The objective was to propose an accurate and robust numerical strategy that withstands the steep gradients of particle volume fraction due to preferential concentration [132] with a limited numerical diffusion. These numerical strategies have been tested on a series of test cases of increasing complexity and relevant diagnostics were proposed. In particular, the two-dimensional vortex laden with solid particles was suggested as a simple configuration to illustrate the effect of particle inertia on their concentration profile and to test numerical strategies. An analytical solution was also derived in the limit of small inertia. Moreover, dissipations due to numerics and to physical effects were explicitly extracted and quantified. Eventually, the numerical strategy coupling the highorder centered scheme TTGC with a stabilization technique –the so called artificial viscosity– proved to be the most accurate and robust alternative in AVBP if an adequate set-up is used (i.e. sensors). Then, the issue of the accurate prediction of particle dispersion in configurations with a mean shear was adressed. One of the RUM model (denoted AXISY-C), proposed by Masi [78] and implemented by Sierra [120], was successfully validated in a two-dimensional and a three-dimensional non-isothermal jet laden with solid particles. Contrary to the former RUM models [63, 103], the main statistics of the dispersed phase were recovered at both the center and the edges of the jet. Finally, the impact of the thermal inertia of particles on their temperature statistics has been investigated. The results showed a strong dependency of these statistics to thermal inertia, pinpointing the necessity of the numerical approaches to account for this phenomenon. Therefore, the extension of the MEF to non isothermal conditions, i.e. the RUM heat fluxes, has been implemented in AVBP. The impact of the RUM HF terms on the temperature statistics was evaluated in both configurations of 2D and 3D jets. Eulerian solutions were compared with Lagrangian reference computations carried out by B. Leveugle at CORIA and by E. Masi at IMFT for the 2D and 3D jets, respectively. Results showed a strong positive impact of the RUM HF on the fluctuations of mesoscopic temperature, and to a lesser extent on the mean mesoscopic temperature depending of the configuration. Neglecting the RUM HF leads to erroneous results whereas the Lagrangian statistics are recovered when they are accounted for.

Este estudo se dedica a uma parte pouco conhecida dos trabalhos de Leonhard Euler (1707 - 1783) relacionado à música enquanto uma ciência matemática. Tais trabalhos mostram, em certo sentido, um lado pitagórico do pensador e também algumas contribuições do mesmo à teoria musical. O interesse deste matemático pelo assunto permeia várias obras e épocas de sua vida, mas neste trabalho focamos em três delas: um conjunto de nove cartas que compõem a obra Lettres a une princesse dAllemagne sur divers sujets de physique e de philosophie (1768) e dois artigos, Conjecture sur la raison de quelques dissonances generalement reçues dans la musique (1766) e De harmoniae veris principiis per speculum musicum repraesentatis (1774). Para possibilitar uma melhor compreensão desses textos, faremos uma revisão histórica do tratamento aritmético dado à música desde os tempos clássicos até o Renascimento com enfoque nas principais contribuições que levam à construção da escala da entonação pura ou justa, trabalhada por Euler nos originais estudados. Após a apresentação e análise desses trabalhos eulerianos, especialmente seus diagramas para representar sons e acordes e do seu expoente de um acorde para medir consonâncias, terminaremos essa dissertação refletindo sobre as implicações pedagógicas e históricas, bem como as potencialidades e limitações do uso de fontes originais de determinados mestres do pensamento matemático ocidental, nomeadamente o próprio Euler, na formação de matemáticos, professores e licenciados. / This present research takes a journey into a little know part of Leonhard Eulers works about music as a mathematical science. Those works show, at certain level, his Pythagorean thoughts and also his contributions to music theory. Eulers interest for this subject permeates many of his works through his lifetime (1707 - 1783), but here we focus on three of them: a set of nine letters from the book Lettres a une princesse dAllemagne sur divers sujets de physique e de philosophie (1768), and two articles, Conjecture sur la raison de quelques dissonances generalement reçues dans la musique (1766) and De harmoniae veris principiis per speculum musicum repraesentatis (1774). To allow a better understanding of his ideas, first we put on a historical review of the arithmetic treatment of music since the ancient Greece to the Renaissance Era pointing out the main contributions to the development of the pure just intonation scale, the one used by Euler on the works just mentioned. After presenting and analyzing these Eulers contributions, specially his diagrams to represent sounds and chords to the eyes and the concept of exponent of a chord to measure consonances, we conclude this dissertation thinking about the pedagogical and historical impact, and also the potential and limitations concerning the use of historical sources of the masters of mathematics, to the training of future mathematics and teachers.

[pt] Resposta a uma preocupação importante na resolução de equações, a impossibilidade de determinar, no conjunto dos números reais, as raízes quadradas de números negativos, o conjunto dos números complexos é estudado em detalhes. O rigor nas demonstrações dos teoremas e propriedades é prioritário para um bom entendimento do assunto. A história dos números complexos permite ao leitor visualizar como e porque surgiram estes números. Seu aparecimento não foi por acaso, e sim, pela necessidade de resolução de equações cúbicas. / [en] Response to a major concern in solving equations, the impossibility of determining, in the set of real numbers, the square roots of negative numbers, the set of complex numbers is studied in details. The rigour in statements of theorems and properties is a priority for a good understanding of the subject. The history of complex numbers allow the reader to see how and why these figures emerged. His appearance was no accident, and yes, by the necessity of solving cubic equations.

Ces travaux de recherche s'intéressent à la caractérisation des transferts thermiques dans les milieux fluide-particules, et en particulier, les lits fluidisés au sein desquels un solide divisé est mis en suspension par un fluide. La grande diversité d'échelles spatiales et temporelles dans ces procédés nécessite d'étudier les interactions hydrodynamiques, thermiques et/ou chimiques entre les particules et le fluide à l'aide d'une approche multi-échelles. Une étude des transferts thermiques dans des lits fixes puis fluidisés, est réalisée à deux échelles : locale (Particle Resolved Simulation) et moyennée (Discrete Element Method-Computional Fluids Dynamics). L'étude PRS permet de caractériser les couplages locaux des transferts thermiques entre particules ainsi que la dynamique de ces transferts dans les configurations fluidisées. Une étude comparative entre les échelles met en évidence les limites du modèle DEM-CFD à capter les fluctuations des transferts thermiques observées dans les simulations PRS. Dans un dernier temps, les fermetures du modèle DEM-CFD sont améliorées de manière à réintroduire les fluctuations perdues par le changement d'échelles. / This work aims at characterizing heat transfer into fluid-solid flows, and more particularly fluidized beds, into which a solid phase is suspended by a flowing fluid. The wide range of spatial and temporal scales present in such processes encourage to study hydrodynamic, thermal and/or chemical interactions between the particles and the fluid through a multi-scale strategy. The analysis of thermal interactions was first carried out for fixed bed configurations and then, fluidized beds at two overlapping scales: local (PRS; Particle Resolved Simulation) and mesoscopic (DEMCFD; Discrete Element Method-Computional Fluids Dynamics). The PRS approach accounts for the local coupling of heat transfer between the particles and its dynamics into fluidized beds. A comparative study of the two scales indicated the limits of the DEM-CFD model to capture the heat transfer fluctuations observed into PRS. In a last step, the closure laws for DEM-CFD were improved to reintroduce the fluctuations lost at this scale.

Les robots manipulateurs sont généralement des machines rigides, conçues pour que leurflexibilité ne perturbe pas leurs mouvements. En effet, des flexibilités mécaniques importantesdans la structure d’un système introduisent des degrés de liberté supplémentaires dont le comportementest complexe et difficile à maîtriser. Cependant, la réduction de la masse d’un systèmeest bénéfique du point de vue des coûts, de la performance énergétique, de la sécurité et des performancesdynamiques. Afin de faciliter l’accès aux nombreux avantages d’une structure légèremalgré la présence de fortes flexibilités, cette thèse porte sur la modélisation, la conception et lacommande de robots manipulateurs flexibles. Elle est motivée par le projet YAKA, dont l’applicationest le lancement et la récupération de drones à voilure fixe depuis un navire faisant route.Cette application nécessite une importante dynamique sur un vaste espace de travail, bien au-delàdes spécifications des robots rigides classiques. Les outils de modélisation, de conception et decommande proposés prennent en compte la flexibilité des segments et des articulations, pour unnombre quelconque de degrés de liberté et de segments flexibles. Le modèle dynamique flexibleest obtenu par le formalisme de Lagrange, les poutres flexibles sont représentées par le modèled’Euler-Bernoulli. Le schéma de commande proposé se décompose en une inversion de modèledynamique rigide et un bloc de précommande par Input Shaping adapté aux robots manipulateursflexibles. Les outils de conception proposés permettent de baser le processus de conceptionsur des performances prédites du système complet muni de ses actionneurs et de son contrôleuravec une simulation réaliste. Les validations expérimentales effectuées sur le robot YAKA permettentde valider la pertinence de la démarche suivie. Les résultats du projet YAKA confirment lafaisabilité de la mise en oeuvre d’un robot flexible de grande envergure et à forte dynamique dansun contexte industriel, en particulier pour le lancement et la récupération d’un drone à voilurefixe depuis un navire faisant route. / Robot manipulators are generally stiff machines, designed in a way that flexibility does not affecttheir movements. Indeed, significant flexibility introduces additional degrees of freedom witha complex behavior. However, reducing the mass of a system allows for costs, performance, andsafety improvements. In order to allow those benefits despite important flexibility, this thesis focuseson modeling, design and control of flexible robot manipulators. It is motivated by the YAKAproject, which aims at developing a robot to launch and recover fixed wing UAVs from a movingship. It implies reaching very high dynamics on a large workspace, way beyond the specificationsof common rigid robots. The proposed tools for modeling, design and control allow for taking intoaccount both joint and link flexibility, for any number of degrees of freedom and flexible links.The elastodynamic model is obtained with Lagrange principle, each flexible link being representedwith one ormany Euler-Bernouilli beams. The proposed control scheme uses a nonlinear rigiddynamic inversion and extends classical Input Shaping techniques to flexible robot manipulators.The proposed design tools allow for performance prediction of the system including its actuatorsand controllers thanks to a realistic simulation. Experiments conducted with the YAKA robot validatedthe proposed approach. The results of the YAKA project confirmed the feasibility of usinga large scale, highly dynamic flexible robot in an industrial context, in particular for UAVs launchand recovery operations from amoving ship.

Este trabalho trata do estudo de robôs de arquitetura paralela, focando na modelagem e otimização dos mesmos. Não foi construído nenhum tipo de protótipo físico, contudo os modelos virtuais poderão, no futuro, habilitar tal façanha. Após uma busca por uma aplicação que se beneficie do uso de um robô de arquitetura paralela, fez-se uma pesquisa por arquiteturas viáveis já existentes ou relatadas na literatura. Escolheu-se a mais apta e prosseguiu-se com os estudos e modelagem cinemática e dinâmica, dando uma maior ênfase na cinemática e dinâmica inversa, esta última utilizando a formulação de Newton - Euler. Foi construído um simulador virtual em ambiente MATLAB 6.5, dotado de várias capacidades como interpolação linear e circular, avanço e uso de múltiplos eixos coordenados. Seu propósito principal é o de demonstrar a funcionalidade e eficácia dos métodos utilizados. Depois foi incorporado ao simulador um algoritmo de cálculo do volume de trabalho da máquina que utiliza alguns dados do usuário para calcular o volume, que pode ser aquele atrelado a uma postura em particular ou o volume de trabalho de orientação total. Algoritmos para medir o desempenho da máquina quanto à uniformidade e utilização da força dos atuadores foram construídos e também incorporados ao simulador, que consegue mostrar o elipsóide de forças ao longo de quaisquer movimentos executados pela plataforma móvel. Quanto à otimização, parte do ferramental previamente construído foi utilizado para que se pudesse chegar a um modelo de uma máquina que respeitasse restrições mínimas quanto ao tamanho e forma de seu volume de trabalho, mas ainda mantendo o melhor desempenho possível dentro deste volume. / This work is about the study of parallel architecture robots, focusing in modeling and optimization. No physical prototypes were built, although the virtual models can help those willing to do so. After searching for an application that could benefit from the use of a parallel robot, another search was made, this time for the right architecture type. After selecting the architecture, the next step was the kinematics and dynamics analysis. The dynamics model is developed using the Newton ? Euler method. A virtual simulator was also developed in MATLAB 6.5 environment. The simulator?s main purpose was to demonstrate that the methods applied were correct and efficient, so it has several features such as linear and circular interpolations, capacity to use multiple coordinate systems and others. After finishing the simulator, an algorithm to calculate the machine workspace was added. The algorithm receives as input some desired requirements regarding the manipulator pose and then calculates the workspace, taking into consideration imposed constraints. Lastly, algorithms capable to measure the manipulator?s performance regarding to its actuator and end-effector force relationship were also incorporated into the simulator that calculates the machine?s force ellipsoid during any movement, for each desired workspace point. For the optimization procedures, some previously developed tools were used, so that the resulting model was capable to respect some workspace constraints regarding size and shape, but also maintaining the best performance possible inside this volume.

Les Structures atomiques Quasi-périodiques (QP) possèdent des propriétés particulières, notamment dans le domaine vibrationnel. Il pourrait être intéressant de pouvoir transférer ces propriétés à des méta-matériaux macroscopiques. Des réseaux de poutres quasi-périodiques 2D sont étudiés dans cette thèse dans le cadre du modèle élément finis (EF) poutre Euler Bernoulli. Ces réseaux de poutres peuvent facilement être produits par fabrication additive ou par découpe laser. Il est possible de faire varier l'élancement des poutres (le ratio hauteur sur longueur) qui est un paramètre intéressant pour modifier la réponse mécanique du réseau. En utilisant la méthode EF, l'influence de l'élancement des poutres sur la réponse vibratoire des réseaux de poutres QP est étudiée. La méthode numérique Kernel Polynomial (KPM) est adaptée avec succès de la dynamique moléculaire aux réseaux de poutres pour étudier leurs modes de vibration sans avoir à diagonaliser complètement la matrice dynamique. Les réseaux de poutres QP présentent des propriétés similaires à leur compère atomique: en particulier la localisation de modes sur des sous-structures et une relation de dispersion hiérarchisée. Le comportement à la fracture est aussi étudié étant donné que les symétries présentes dans les QP pourraient permettre des réseaux de poutres ne présentant pas de plans faibles pour la propagation de fissures. Cela a été démontré d'après des calculs EF statiques avec un critère de fracture fragile sur l'énergie de déformation. Les simulations statiques ne suffisent pas car elles ne peuvent pas capturer les phénomènes dynamiques complexes qui apparaissent lors de la fissuration fragile. Les propriétés de vibration du QP pourraient aussi avoir un impact sur la propagation dynamique de fissure. Un modèle dynamique de fissuration est développé afin d'étudier l'impact de l'élancement sur la capacité des réseaux de poutres QP à dissiper de l'énergie par fissuration. Finalement une méthode Coarse Graining est développée pour identifier un milieu Cosserat continu équivalant au réseau de poutres QP pour différentes échelles. Cette méthode permet d'identifier la densité, les déformations, les contraintes et donc les modules d'élasticité du milieu Cosserat équivalent, permettant ainsi une meilleure compréhension du rôle des sous structures précédemment identifiées. / Quasi periodic (QP) structures have shown peculiar properties in the atomistic domain, especially the vibrational one. It could be interesting to be able to transpose these properties in macroscopic meta-materials. Quasi periodic 2D beam lattices are studied in this thesis due to the simplicity of the Euler Bernoulli finite element (FE) model. These beam lattices can easily be produced by additive manufacturing or by laser cutting. It is possible to vary the beam slenderness (i.e the ratio of height over length) that is a interesting parameter to modify the mechanical response of the lattice. Using finite element method, the influence of the beam slenderness over the vibration behavior of the QP beam lattices will be studied. The Kernel Polynomial numerical Method (KPM) is successfully adapted from molecular dynamics simulations in order to study vibrational modes of FE beam lattices without having to fully diagonalize the dynamical matrix. The QP lattices show similar properties as their atomic counterpart e.g mode localization over sub-stuctures and hierarchical dispersion relation. The fracture behavior is also studied, as the special symmetries allowed by the quasi periodicity could result in beam lattices without weak planes for crack propagation. It was proved to be true from static FE simulations with a brittle strain energy breaking criterion. Static simulations were not enough and do not grasp the complex dynamical phenomena taking place in brittle fracture. A dynamic crack propagation model was thus developed. The vibrational properties of quasi periodic structures could also have an impact on the dynamic crack propagation. Several simulations are run in order to study the impact of the slenderness on the energy dissipated by fracture of QP lattices. Finally, a coarse graining method (CG) was developed to identify a continuous Cosserat medium at different scales from the FE beam model. This CG method allows to identify, density, strain, stress and elastic moduli of an equivalent continuous Cosserat. This allows a better understanding of the role of previously identified characteristic sub structures.

In der vorliegenden Arbeit werden kontinuumsmechanische Probleme mit nichtmateriellen Randbedingungen untersucht. Randbedingungen gelten dabei als nichtmateriell, wenn sie im Zeitverlauf nicht ein und demselben materiellen Punkt zugeordnet werden können. Die Erweiterung der klassischen kontinuumsmechanischen Feldgleichungen um solche Randbedingungen erfolgt unter Anwendung einer Arbitrary-LAGRANGE-EULER-Kinematik. Hierbei wird eine Notation entwickelt, bei der Feldgrößen und Operatoren ihre jeweilige Platzierung eindeutig zugeordnet wird. Insbesondere in Hinblick auf eine konsistente Darstellung von Ableitungsoperatoren werden die Vorteile dieser Schreibweise dargelegt. Zur Ermittlung und Untersuchung (semi-)analytischer Lösungen dienen Beispiele eindimensionaler Kontinua, die sich zwei unterschiedlichen Problemklassen zuordnen lassen. In der ersten Problemklasse gelingen analytische Lösungen mit Hilfe eines Integrations- und eines Separationsansatzes für das Modell einer axial unbewegten, schwingenden Saite. Als nichtmaterielle Randbedingungen werden dabei die transversalen Verschiebungen an zwei zeitabhängigen Positionen zu null vorgeschrieben. In der zweiten Problemklasse sind eine Saite sowie ein Seil, die einer vorgegebenen axialen Führungsbewegung unterliegen, Gegenstand der Untersuchung. In diesem Fall sind die zwei vorgegebenen, räumlich festen Verschiebungsrandbedingungen nichtmateriell. Es finden (semi-)analytische Verfahren Anwendung. Die Relativgeschwindigkeit zwischen den Randbedingungen und dem jeweils betrachteten Kontinuum wird dabei als beliebig zeitabhängig angenommen. Eine experimentelle Studie zum Schwingungsverhalten eines Monochords mit nichtmateriellen Randbedingungen vervollständigt die Analyse eindimensionaler Kontinua. Aus den ermittelten (semi-)analytischen Lösungen werden Rückschlüsse auf das Transformationsverhalten der Bewegungsgleichungen dreidimensionaler Kontinua gezogen. Damit sind die entwickelten Methoden in vielen technischen Anwendungen einsetzbar. Als ein wirtschaftlich bedeutendes Beispiel ist die Schwingungsanalyse axial bewegter Papierbahnen in Papierproduktionsmaschinen zu nennen. / Within this work, problems of continuum mechanics with non-material boundary conditions are investigated. Boundary conditions are classified as non-material if they can not be assigned to one and only one material particle over time. The extension of the classical field-equations of continuum mechanics by such boundary conditions is realized by application of Arbitrary-LAGRANGE -E ULER -Kinematics. Therefore a notation, which assigns the particular placement to field quantities and operators, is developed. The advantages of this notation can be identified particularly with regard to a consistent representation of derivative operators. Examples of one-dimensional continua, which can be assigned to different problem categories, are used to determine and investigate (semi-)analytical solutions. In the first category, analytical solutions can be found using an integral and a separation formulation for the model of an axially non-moving, vibrating string. As non-material boundary conditions the transverse displacements at two time-dependent positions are prescribed to zero. A string and a wire, which are moved axially, are investigated within the second problem category. In this case, the prescribed, spatially fixed displacement conditions are non-material. The applied methods are (semi-)analytical. The relative velocity between the boundary conditions and the considered continuum is assumed to be arbitrary time-dependent. An experimental study on the vibration behaviour of a monochord with non-material boundary conditions completes the analysis of one-dimensional continua. Conclusions on the transformation of the equations of motion of three-dimensional continua are derived from the determined (semi-)analytical solutions. For this reason the developed methods are usable in many technical applications. The vibration analysis of axially moving paper sheets in papermaking machines can be stated as an economical important example.

A finite-volume procedure, comprising a gradient-reconstruction technique and a multidimensional limiter, has been proposed for upwind algorithms on unstructured grids. The high-resolution strategy, with its inherent dependence on a wide computational stencil, does not suffer from a catastrophic loss of accuracy on a grid with poor connectivity as reported recently is the case with many unstructured-grid limiting procedures. The continuously-differentiable limiter is shown to be effective for strong discontinuities, even on a grid which is composed of highly-distorted triangles, without adversely affecting convergence to steady state. Numerical experiments involving transient computations of two-dimensional scalar convection to steady-state solutions of Euler and Navier-Stokes equations demonstrate the capabilities of the new procedure.

Analysis of unsteady flows is a very challenging topic of research. A decade ago, potential flow equations were used to predict unsteady pressures on oscillating bodies. Recognising the fact that nonlinear aerodynamics is essential to analyse unsteady flows accurately, particularly in transonic and supersonic flows, different Euler formulations operating on moving grids have emerged recently as important CFD tools for unsteady aerodynamics. Numerical solution of Euler equations on moving grids based on upwind schemes such as the ones due to van Leer and Roe have been developed for the purpose of numerical simulation of unsteady transonic and supersonic flows. In the present work, Euler computations based on yet another recent robust upwind scheme (for steady flows) namely Kinetic Flux Vector Splitting (KFVS) scheme due to Deshpande and Mandal is chosen for further development of a time accurate Euler solver to operate on problems involving moving boundaries. The development of an Euler code based on this scheme is likely to be highly useful to analyse problems of unsteady aerodynamics and computational aeroelasiticity especially when it is noted that KFVS has been found to be an extremely robust scheme for computation of subsonic, transonic, supersonic and hypersonic flows. The KFVS scheme, basically exploits the connection between the linear scalar Boltzmann equation of kinetic theory of gases and the nonlinear vector conservation law, that is, Euler equations of fluid dynamics through moment method strategy. The KFVS scheme has inherent simplicity in splitting the flux even on moving grids due to underlying particle model. The inherent simplicity of KFVS for moving grid problems is due to its relationship with the Boltzmann equation. If a surface is moving with velocity w and a particle has velocity v, then it is quite reasonable to do the splitting based on (v-w)<0 or >0. Only particles having velocity v greater than w will cross the moving surface from left to right and similar arguments hold good for particles moving in opposite direction. It is therefore quite natural to extend KFVS by splitting the Maxwellian velocity distribution at Boltzmann level based on the sign of the normal component of the relative velocity. The relative velocity is the difference between the molecular velocity (v) and the velocity of the moving surface(w). This inherent simplicity of the Kinetic Flux Vector Splitting scheme on Moving Grids (KFMG) method has prompted us to extend the same ideas to 2-D and 3-D problems leading to the present KFMG method. If w is set to zero then KFMG formulation reduces to the one corresponding to KFVS. Thus KFMG formulations axe generalisation of the KFVS formulation. In 2-D and 3-D cases, in addition to the KFMG formulation, the method to move the grids, the appropriate boundary conditions for treating moving surfaces and techniques to improve accuracy in space and time are required to be developed. The 2-D and 3-D formulations based on Kinetic Flux Vector Splitting scheme on Moving Grids method have been developed for computing unsteady flows. Between two successive time steps, the body changes its orientation in case of an oscillation or it deforms when subjected to, aerodynamic loads. In either of these cases the grid corresponding to the first time step has to be moved or regenerated around the displaced or deformed body. There are several approaches available to generate grids around moving bodies. In the present work, the 'spring analogy method' is followed to obtain grid around deflected geometries within the frame work of structured grid. Using this method, the grids are moved from previous time to the current time. This method is capable of tackling any kind of aeroelastic deformation of the body. For oscillating bodies, a suitable boundary condition enforcing the flow tangency on the body needs to be developed. As a first attempt, the body surface has been treated as an 1-D piston undergoing compression and expansion. Then, a more general Kinetic Moving Boundary Condition(KMBC) has been developed. The KMBC uses specular reflection model of kinetic theory of gases. In order to treat fixed outer boundary, Kinetic Outer Boundary Condition(KOBC) has been applied. The KOBC is more general in the sense that, it can treat different type of boundaries (subsonic, supersonic, inflow or out flow boundary). A 2-D cell-centered finite volume KFMG Euler code to operate on structured grid has been developed. The time accuracy is achieved by incorporating a fourth order Runge-Kutta time marching method. The space accuracy has been enhanced by using high resolution scheme as well as second order scheme using the method of reconstruction of fluxes. First, the KFMG Euler code has been applied to standard test cases for computing steady flows around NACA 0012 and NACA 64AQ06 airfoils in transonic flow. For these two airfoils both computational and experimental results are available in literature. It is thus possible to verify (that is, prove the claim that code is indeed solving the partial differential equations + boundary conditions posed to the code) and validate(that is, comparison with experimental results) the 2-D KFMG Euler code. Having verified and validated the 2-D KFMG Euler code for the standard test cases, the code is then applied to predict unsteady flows around sinusoidally oscillating NACA 0012 and NACA 64A006 airfoils in transonic flow. The computational and experimental unsteady results are available in literature for these airfoils for verification and validation of the present results. The unsteady lift and normal force coefficients have been predicted fairly accurately by all the CFD codes. However there is some difficulty about accurate prediction of unsteady pitching moment coefficient. Even Navier-Stokes code could not predict pitching moment accurately. This issue needs further in depth study and probably intensive computation which have not been undertaken in the present study. Next, a two degrees of £reedom(2-DOF) structural dynamics model of an airfoil undergoing pitch and plunge motions has been coupled with the 2-D KFMG Euler code for numerical simulation of aeroelastic problems. This aeroelastic analysis code is applied to NACA 64A006 airfoil undergoing pitch and plunge motions in transonic flow to obtain aeroelastic response characteristics for a set of structural parameters. For this test case also computed results are available in literature for verification. The response characteristics obtained have showed three modes namely stable, neutrally stable and unstable modes of oscillations. It is interesting to compare the value of airfoil-to-air mass ratio (Formula) obtained by us for neutrally stable condition with similar values obtained by others and some differences between them are worth mentioning here. The values of \i for neutral stability are different for different authors. The differences in values of (Formula) predicted by various authors are primarily due to differences which can be due to grid as well as mathematical model used. For example, the Euler calculations, TSP calculations and full potential calculations always show differences in shock location for the same flow problem. Changes in shock location will cause change in pressure distribution on airfoil which in turn will cause changes in values of \L for conditions of neutral stability. The flutter speed parameter(U*) has also been plotted with free stream Mach number for two different values of airfoil - to - air mass ratio. These curves shown a dip when the free stream Mach number is close to 0.855. This is referred as "Transonic Dip Phenomenon". The shock waves play a dominant role in the mechanism of transonic dip phenomenon. Lastly, cell-centered finite volume KFMG 3-D Euler code has been developed to operate on structured grids. The time accuracy is achieved by incorporating a fourth order Runge-Kutta method. The space accuracy has been enhanced by using high resolution scheme. This code has 3-D grid movement module which is based on spring analogy method. The KMBC to treat oscillating 3-D configuration and KOBC for treating 3-D outer boundary have also been formulated and implemented in the code. The 3-D KFMG Euler code has been first verified and validated for 3-D steady flows around standard shapes such as, transonic flow past a hemisphere cylinder and ONERA M6 wing. This code has also been used for predicting hypersonic flow past blunt cone-eylinder-flare configuration for which experimental data are available. Also, for this case, the results are compared with a similar Euler code. Then the KFMG Euler code has been used for predicting steady flow around ogive-cylinder-ogive configuration with elliptical cross section. The aerodynamic coefficients obtained have been compared with those of another Euler code. Thus, the 3-D KFMG Euler code has been verified and validated extensively for steady flow problems. Finally, the 3-D KFMG based Euler code has been applied to an oscillating ogive-cylinder-ogive configuration in transonic flow. This test case has been chosen as it resembles the core body of a flight vehicle configuration of interest to DRDO,India. For this test case, the unsteady lift coefficients are available in literature for verifying the present results. Two grid sizes are used to perform the unsteady calculations using the present KFMG 3-D Euler code. The hysteresis loops of lift and moment coefficients confirmed the unsteady behaviour during the oscillation of the configuration. This has proved that, the 3-D formulations are capable of predicting the unsteady flows satisfactorily. The unsteady results obtained for a grid with size of 45x41x51 which is very close to the grid size chosen in the reference(Nixon et al.) are considered for comparison. It has been mentioned in the reference that, a phase lag of (Formula) was observed in lift coefficients with respect to motion of the configuration for a free stream Mach number of 0.3 with other conditions remaining the same. The unsteady lift coefficients obtained using KFMG code as well as those available in literature are plotted for the same flow conditions. Approximately the same phase lag of (Formula) is present (for (Formula)) between the lift coefficient curves of KFMG and due to Nixon et al. The phase lag corrected plot of lift coefficient obtained by Nixon et al. is compared with the lift coefficient versus time obtained by 3-D KFMG Euler code. The two results compare well except that the peaks are over predicted by KFMG code. It is nut clear at this stage whether our results should at all match with those due to Nixon et al. Further in depth study is obviously required to settle the issue. Thus the Kinetic Flux Vector Splitting on Moving Grids has been found to be a very good and a sound method for splitting fluxes and is a generalisation of earlier KFVS on fixed grids. It has been found to be very successful in numerical simulation of unsteady aerodynamics and computational aeroelasticity.