Les nouveaux enjeux environnementaux motivent la recherche par les acteurs de l'industrie aéronautique de méthodes de calcul de trajectoires optimales. Les contributions de cette thèse se déclinent selon trois axes. Dans un premier temps, plusieurs techniques d'optimisation de trajectoire avion sont comparées sur un cas simple traité dans la bibliographie universitaire. Puis, un modèle réduit pour prendre en compte le niveau des nuisances sonores dans un algorithme d'optimisation de trajectoire est proposé.Enfin, un problème d'optimisation de trajectoire de montée d'un avion civil est résolu par une approche directe. Les spécificités du problème consistent en la présence de plusieurs phases au sein de la trajectoire, la formulation de contraintes égalités à vérifier par des composantes du vecteur d'état sur des intervalles de temps et enfin la difficulté d'intégration numérique du modèle de bruit. / Forthcoming environmental challenges stimulate the development of trajectory optimization methods by aeronautical actors. This contribution consists in three parts. First, several trajectory optimization techniques are compared. The comparison is based on a simple academic problem. After that, a model is proposed for considering noise nuisance level in the framework of trajectory optimization. Finally, the optimization problem of an ascent phase of a civil aircraft is solved using a direct approach. The specific issues of the problem are tackled with a general formulation. They consist in the presence of several phases along the trajectory, running state equality constraints and tough numerical integration of the noise model.

Cet ouvrage aborde la génération de grilles de Voronoï sous contrainte pour réduire les erreurs liées à la géométrie des cellules lors de la simulation réservoir. Les points de Voronoï sont optimisés en minimisant des fonctions objectif correspondant à différentes contraintes géométriques. L'originalité de cette approche est de pouvoir combiner les contraintes simultanément : - la qualité des cellules, en plaçant les points de Voronoï aux barycentres des cellules ; - le raffinement local, en fonction d'un champ de densité [rho], correspondant à la perméabilité, la vitesse ou la vorticité ; - l'anisotropie des cellules, en fonction d'un champ de matrice M contenant les trois vecteurs principaux de l'anisotropie, dont l'un est défini par le vecteur vitesse ou par le gradient stratigraphique ; - l'orientation des faces des cellules, en fonction d'un champ de matrice M contenant les trois vecteurs orthogonaux aux faces, dont l'un est défini par le vecteur vitesse ; - la conformité aux surfaces du modèle structural, failles et horizons ; - l'alignement des points de Voronoï le long des puits. La qualité des grilles générées est appréciée à partir de critères géométriques et de résultats de simulation comparés à des grilles fines de référence. Les résultats indiquent une amélioration de la géométrie, qui n'est pas systématiquement suivie d'une amélioration des résultats de simulation / Voronoi grids are generated under constraints to reduce the errors due to cells geometry during flow simulation in reservoirs. The Voronoi points are optimized by minimizing objective functions relevant to various geometrical constraints. An original feature of this approach is to combine simultaneously the constraints: - Cell quality, by placing the Voronoi points at the cell barycenters. - Local refinement according to a density field rho, relevant to permeability, velocity or vorticity. - Cell anisotropy according to a matrix field M built with the three principal vectors of the anisotropy, which one is defined by the velocity vector or by the stratigraphic gradient. - Faces orientation according to a matrix field M built with the three vectors orthogonal to the faces, which one is defined by the velocity vector. - Conformity to structural features, faults and horizons. - Voronoï points alignment along well path. The quality of the generated grids is assessed from geometrical criteria and from comparisons of flow simulation results with reference fine grids. Results show geometrical improvements, that are not necessarily followed by flow simulation results improvements

Approved for public release, distribution is unlimited / The extreme strain, strain rate and temperature gradients during Friction Stir Processing (FSP) render measurement of key parameters in the stir zone infeasible with common methods. The objective of this research was to separate the effects that temperature and deformation in an experimental study of the microstructure and mechanical properties of Ni-AL bronze (NAB). This was accomplished by subjecting as-cast NAB material to several isothermal annealing and quenching treatments as well as isothermal hot rolling processes. Sufficient material was generated to provide results and data for subsequent optical microscopy, tensile, & hardness tests. All results were then compared to similar data collected from previous works completed here at Naval Postgraduate School and with other DARPA FSP program participants. During the course of this work correlations were drawn between FSP material and the material subjected to isothermal hotworking, which may enhance our understanding of the roles that various FSP process parameters have on the microstructural transformation sequence within this material. The hot-rolling study conducted here suggests that FSP process parameters leading to severe deformation at temperatures between 950-1000 C in the NAB material provides high ductility (elongation approximately 28%) with moderate strengths. / Lieutenant, United States Coast Guard

Finite amplitude perturbations in compressible media are ubiquitous in scientific and engineering applications such as gas-turbine engines, rocket propulsion systems, combustion instabilities, inhomogeneous solids, and traffic flow prediction models, to name a few. Small amplitude waves in compressible fluids propagate as sound and are very well described by linear theory. On the other hand, the theory of nonlinear acoustics, concerning high-amplitude wave propagation (Mach<2) is relatively underdeveloped. Most of the theoretical development in nonlinear acoustics has focused on wave steepening and has been centered around the Burgers' equation, which can be extended to nonlinear acoustics only for purely one-way traveling waves. In this dissertation, theoretical and computational developments are discussed with the objective of advancing the multi-fidelity modeling of nonlinear acoustics, ranging from quasi one-dimensional high-amplitude waves to combustion-induced detonation waves. <br> <br> We begin with the theoretical study of spectral energy cascade due to the propagation of high amplitude sound in the absence of thermal sources. To this end, a first-principles-based system of governing equations, correct up to second order in perturbation variables is derived. The exact energy corollary of such second-order system of equations is then formulated and used to elucidate the spectral energy dynamics of nonlinear acoustic waves. We then extend this analysis to thermoacoustically unstable waves -- i.e. amplified as a result of thermoacoustic instability. We drive such instability up until the generation of shock waves. We further study the nonlinear wave propagation in geometrically complex case of waves induced by the spark plasma between the electrodes. This case adds the geometrical complexity of a curved, three-dimensional shock, yielding vorticity production due to baroclinic torque. Finally, detonation waves are simulated by using a low-order approach, in a periodic setup subjected to high pressure inlet and exhaust of combustible gaseous mixture. An order adaptive fully compressible and unstructured Navier Stokes solver is currently under development to enable higher fidelity studies of both the spark plasma and detonation wave problem in the future. <br>

The thesis presents a comprehensive study of a posteriori error estimation in the adaptive solution to some classical elliptic partial differential equations. Several new error estimators are proposed for diffusion problems with discontinuous coefficients and for convection-reaction-diffusion problems with dominated convection/reaction. The robustness of the new estimators is justified theoretically. Extensive numerical results demonstrate the robustness of the new estimators for challenging problems and indicate that, compared to the well-known residual-type estimators, the new estimators are much more accurate.

Cette thèse s'intéresse à la conception et à l'analyse de méthodes de discrétisation hybrides pour les inégalités variationnelles non linéaires apparaissant en mécanique des fluides et des solides. Les principaux avantages de ces méthodes sont la conservation locale au niveau des mailles, la robustesse par rapport à différents régimes de paramètres et la possibilité d’utiliser des maillages polygonaux / polyédriques avec des nœuds non coïncidants, ce qui est très intéressant dans le contexte de l’adaptation de maillage. Les méthodes de discrétisation hybrides sont basées sur des inconnues discrètes attachées aux faces du maillage. Des inconnues discrètes attachées aux mailles sont également utilisées, mais elles peuvent être éliminées localement par condensation statique. Deux applications principales des discrétisations hybrides sont abordées dans cette thèse. La première est le traitement par la méthode de Nitsche du problème de contact de Signorini (dans le cas scalaire) avec une non-linéarité dans les conditions aux limites. Nous prouvons des estimations d'erreur optimales conduisant à des taux de convergence d'erreur d'énergie d'ordre (k + 1), si des polynômes de face de degré k >= 0 sont utilisés. La deuxième application principale concerne les fluides à seuil viscoplastiques. Nous concevons une méthode de Lagrangien augmenté discrète appliquée à la discrétisation hybride. Nous exploitons la capacité des méthodes hybrides d’utiliser des maillages polygonaux avec des nœuds non coïncidants afin d'effectuer l’adaptation de maillage local et mieux capturer la surface limite. La précision et la performance des schémas sont évaluées sur des cas tests bidimensionnels, y compris par des comparaisons avec la littérature / This thesis is concerned with the devising and the analysis of hybrid discretization methods for nonlinear variational inequalities arising in computational mechanics. Salient advantages of such methods are local conservation at the cell level, robustness in different regimes and the possibility to use polygonal/polyhedral meshes with hanging nodes, which is very attractive in the context of mesh adaptation. Hybrid discretizations methods are based on discrete unknowns attached to the mesh faces. Discrete unknowns attached to the mesh cells are also used, but they can be eliminated locally by static condensation. Two main applications of hybrid discretizations methods are addressed in this thesis. The first one is the treatment using Nitsche's method of Signorini's contact problem (in the scalar-valued case) with a nonlinearity in the boundary conditions. We prove optimal error estimates leading to energy-error convergence rates of order (k+1) if face polynomials of degree k >= 0 are used. The second main application is on viscoplastic yield flows. We devise a discrete augmented Lagrangian method applied to the present hybrid discretization. We exploit the capability of hybrid methods to use polygonal meshes with hanging nodes to perform local mesh adaptation and better capture the yield surface. The accuracy and performance of the present schemes is assessed on bi-dimensional test cases including comparisons with the literature

Sem resumo / Sem abstract

Neste trabalho desenvolveu-se uma formulação adaptativa do método de elementos de contorno (MEC) para a análise de problemas de fratura elástica linear. Foi utilizado o método da colocação para a formulação das equações integrais de deslocamento e de tensão. Para a discretização das equações integrais foram utilizados elementos lineares que possibilitaram a obtenção das expressões exatas das integrais (integração analítica) sobre elementos de contorno e fratura. Para a montagem do sistema de equações algébricas foram utilizadas apenas equações de deslocamento, apenas equações de forças de superfície, ou as duas escritas para nós opostos da fratura levando, portanto ao método dos elementos de contorno dual usualmente empregado na análise de fratura. Para o processo de crescimento da trinca foi desenvolvido um procedimento especial objetivando a correta determinação da direção de crescimento da trinca. Os fatores de intensidade de tensão são calculados por meio da conhecida técnica de correlação de deslocamentos a qual relaciona os deslocamentos atuantes nas faces da fissura. Após a determinação dos fatores de intensidade de tensão é utilizada a teoria da máxima tensão circunferencial para a determinação do ângulo de propagação. O modelo adaptativo empregado é do tipo h onde apenas a sub-divisão dos elementos é feita com base em erros estimados. O erro a ser considerado foi estimado a partir de normas onde se consideraram: a variação aproximada dos deslocamentos, a variação das forças de superfície e a variação da energia de deformação do sistema, calculada com a sua integração sobre o contorno. São apresentados exemplos numéricos para demonstrar a eficiência dos procedimentos propostos. / In this work, an adaptative formulation of the boundary element method is developed to analyze linear elastic fracture problems. The collocation point method was used to formulate the integral equations for the displacements and stresses (or tractions). To discretize the integral equations, linear elements were used to obtain the exact expressions of the integrals over boundary elements and fracture. To construct the linear system of equations were used only displacement equations, traction equations or both of them written for opposite nodes of the fracture, leading to the dual boundary element formulation usually employed in the fracture analyses. For the process of growth of the crack a special procedure was developed aiming at the correct determination of the direction of growth of the crack. The stress intensity factors, to calculate he crack growth angle, are calculated through of correlation displacements technique which relates the displacements actuants in the faces of the crack. The employed adaptative model is the h-type where only the sub-division of the elements is done based on error estimate. The error estimates considered in this work are based on the following norms: displacement, traction and strain energy variations, this last considered from the integration over the boundary. Numerical examples are presented to demonstrate the efficiency of the proposed procedures.

Recentemente introduzida no mercado, a liga à base de titânio &beta;-metaestável Ti-5Al-5V-5Mo3Cr-1Zr (Ti-55531) é um material estrutural para aplicações nos setores automotivo e aéreo, além de compor peças para centrífugas. A presente Dissertação de Mestrado tem como objetivo o estudo in situ das transformações de fase em estado sólido da liga Ti-55531 durante o envelhecimento a diferentes temperaturas e condições iniciais por meio de refinamento Rietveld. A ênfase do projeto será a análise qualitativa e quantitativa das principais fases presentes nesta liga, em especial aquelas &beta; de estrutura cúbica de corpo centrado, alfa (&alpha;) e ômega (&omega;) de estrutura hexagonal e martensita (&alpha;\'\') de estrutura ortorrômbica. Os padrões de difração desta liga para os diferentes tempos de envelhecimento foram obtidas com o auxílio da técnica de difração de raio-X de alta energia in situ utilizando luz síncrotron, enquanto que as respectivas análises das imagens dos anéis de difração a partir destes experimentos foram realizadas com o auxílio do programa de integração e calibração instrumental FIT2D e do programa de análise de materiais usando difração (do inglês MAUD). / Recently introduced, the new Ti-55531 titanium based alloy has been developed aiming structural applications such as the automotive and aerospace ones. The present Thesis is focused on the assessment of Ti-55531 alloy in situ solid state transformations during ageing at diferente temperatures as well as initial conditions. In this sense, the study emphasis will deal with both quantitative and qualitative analysis of the main constitutive phases found in the alloy, namely the body-centered cubic &beta;, both hexagonal &alpha; and &omega; and the martensitic &alpha;\'\', this one orthorhombic, by means of Rietveld refinement. Concerning the data, it was obtained by high energy X-ray diffraction (HEXRD) technique in a synchrotron light source facility. All the Debye-Scherrer rings resulting from these in situ experiments were analyzed by both FIT2D and MAUD software.

In dieser Arbeit geht es um die schrittweise Verfeinerung verteilter Algorithmen. Dabei wird ein einfacher Algorithmus, der einige gewünschte Eigenschaften hat, Schritt für Schritt zu einem komplexen Algorithmus verfeinert, der konkrete Implementationsanforderungen erfüllt, so daß in jedem Schritt die gewünschten Eigenschaften erhalten bleiben. Wir stellen einen neuen eigenschaftserhaltenden Verfeinerungsbegriff vor, der auf der kausalen Ordnung der Aktionen eines Algorithmus basiert. Diesen Begriff definieren wir als Transitionsverfeinerung für elementare Petrinetze und diskutieren Beweiskriterien. Danach definieren und diskutieren wir die simultane Verfeinerung mehrerer Transitionen. Zur Modellierung komplexer verteilter Algorithmen sind elementare Petrinetze oft nicht adäquat. Wir benutzen deshalb algebraische Petrinetze. Wir definieren Transitionsverfeinerung für algebraische Petrinetze und stellen einen Zusammenhang zur simultanen Verfeinerung von Transitionen in elementaren Petrinetzen her. Transitionsverfeinerung ist besonders für Verfeinerungsschritte geeignet, in denen synchrone Kommunikation zwischen Agenten durch asynchronen Nachrichtenaustausch ersetzt wird. Wir zeigen dies am Beispiel eines komplexen verteilten Algorithmus, zur Berechnung des minimalen spannenden Baumes in einem gewichteten Graphen. Wir zeigen die Korrektheit dieses Algorithmus in mehreren Schritten, von denen einige Schritte Transitionsverfeinerungen sind. In anderen Schritten sind klassische Verfeinerungsbegriffe ausreichend. Wir übertragen deshalb auch einen klassischen Verfeinerungsbegriff in unser formales Modell. / The topic of this PhD thesis is the stepwise refinement of distributed algorithms. Stepwise refinement starts with a simple algorithm with certain desired properties. This algorithm is refined step by step such that the desired properties are preserved in each refinement step. The result is a complex distributed algorithm which satisfies concrete implementation requirements and which still has the desired properties. We propose a new property preserving notion of refinement which is based on the causal ordering of actions of an algorithm. We call this notion transition refinement and we define it first for elementary Petri nets. Furthermore, we discuss proof criteria. Then, we define and discuss the simultaneous refinement of several transitions. For modelling complex distributed algorithms, we use algebraic Petri nets instead of elementary Petri nets. We define transition refinement for algebraic Petri nets, and we show its relationship to simultaneous transition refinement in elementary Petri nets. Transition refinement is particularly suitable for refinement steps in which synchronous communication between agents is replaced by asynchronous message passing. We show this by means of a complex distributed algorithm for determining the minimal spanning tree of a weighted graph. We prove the correctness of this algorithm in several steps. Some of these steps are transition refinements. For other steps, well-known notions of refinement are sufficient. Therefore, we also carry over a well-known notion of refinement into our formal model.