Cette thèse est dédiée à la conception optimale de la machine synchrone à rotor bobiné modulaire sans balais pour les systèmes embarqués. Cette machine est basée sur une structure POKIPOKITM développée par Mitsubishi Electric Coopération avec les convertisseurs de puissance intégrée pour augmenter la capacité de tolérance aux défauts. L'analyse électromagnétique est utilisée pour étudier les différentes machines synchrones à rotor bobiné et donc, pour sélectionner la structure qui offre la meilleure tolérance aux défauts et les performances les plus élevées. D’abord, le choix des nombres de phases, d’encoches et de pôles est un point critique. Ensuite, quelques machines sont analysées et comparées selon les critères tels que la densité de couple, le rendement, l'ondulation de couple. La machine avec 7 phases, 7 encoches et 6 pôles est alors choisie. Cette machine est ensuite comparée à la machine synchrone à aimant permanent monté en surface. Le résultat démontre que la machine synchrone à rotor bobiné modulaire sans balais possède le potentiel de remplacer la machine synchrone à aimant permanent dans notre application parce qu’elle présente des performances similaires avec une capacité de tolérance aux défauts élevée. Dans un second temps, une fois la structure 7phases/7encoches/6pôles choisie, cette machine est optimisée en utilisant NOMAD (qui est un logiciel d'optimisation de boîte noire) afin de minimiser le volume externe sous les contraintes électromagnétiques, thermiques et mécaniques. Comme ce problème d'optimisation est extrêmement difficile à résoudre, quelques relaxations ont été effectuées pour tester les différents algorithmes d'optimisation : fmincon (de Matlab) et NOMAD. Nous remarquons que NOMAD est plus efficace que fmincon pour trouver des solutions à ce problème de conception où certaines contraintes sont calculées par des simulations numériques (ANSYS Maxwell ; code éléments finis). En utilisant la méthode NOMAD basée sur l’algorithme Mesh Adaptive Direct Search, nous obtenons des résultats optimaux qui satisfont toutes les contraintes proposées. Il est nécessaire de valider ce design optimisé en vérifiant toutes les contraintes par des simulations électromagnétiques et thermiques en 3D. Les résultats montrent que le couple moyen obtenu par la simulation en 3D est inférieur à la valeur souhaitée. Par conséquent, en augmentant la longueur de la machine, une nouvelle machine corrigée est ainsi obtenue. Nous observons que les pertes de fer obtenues en 3D sont plus élevées qu'en 2D en raison du flux de fuite dans la tête de bobinage. En prenant les valeurs des pertes analysées par la simulation en 3D, la température de surface de la nouvelle machine analysée par la méthode Computational Fluid Dynamics est plus élevée que celle calculée dans l’optimisation. Enfin, un prototype de machine est construit et quelques tests expérimentaux est réalisés. Le résultat montre que la force électromotrice à vide a une forme d'onde similaire par rapport à la prédiction numérique en 3D et la différence de couple statique maximum entre les tests expérimentaux et les simulations par éléments finis en 3D est faible. / This thesis is dedicated to the design and the optimization of modular brushless wound rotor synchronous machine for embedded systems. This machine is constructed based on POKIPOKITM structure with integrated drive electronics. Finite element analysis based optimization becomes more popular in the field of electrical machine design because analytical equations are not easily formalized for the machines which have complicate structures. Using electromagnetic analysis to comparatively study different modular brushless wound rotor synchronous machines and therefore, to select the structure which offers the best fault tolerant capability and the highest output performances. Firstly, the fundamental winding factor calculated by using the method based on voltage phasors is considered as a significant criterion in order to select the numbers of phases, stator slots and poles. After that, 2D finite element numerical simulations are carried out for a set of 15 machines to analyze their performances. The simulation results are then compared to find an appropriate machine according to torque density, torque ripple and machine efficiency. The 7phase/7-slot/6-pole machine is chosen and compared with a reference design surfacemounted permanent magnet synchronous machine in order to evaluate the interesting performance features of the wound rotor synchronous machine. In the second design stage, this machine is optimized by using derivative-free optimization. The objective is to minimize external volume under electromagnetic, thermal and mechanical constraints. Given that an accurate finite element analysis for machine performance takes a long time. Moreover, considering that the average torque can be obtained by simulating the model with only four rotor positions instead of one electric period, optimization strategy is proposed to reduce computational time and therefore, obtain a fast convergence ability by defining relaxed problems which enable minimizing the external volume of the machine under only several constraints such as average torque, torque ripple and copper losses. By testing relaxed problems, two different optimization methods (NOMAD and fmincon) are compared in order to select an appropriate method for our optimization problem. Using NOMAD method based on Mesh Adaptive Direct Search, we achieve optimal results which satisfy all of the constraints proposed. In the third design stage, all constraints are validated by 3D electromagnetic and thermal simulations using finite element and computational fluid dynamics methods. The 3D results show that the average torque obtained is lower than the desired value. By increasing the length of the machine, a new corrected machine is thus obtained. It can be observed that the iron losses obtained in 3D are higher than that in 2D due to the leakage flux in the end-winding. Then, the machine temperature is analyzed by using ANSYS Fluent. Note that the surface temperature is higher than that calculated in the optimization and the coil temperature is 8.48°C higher than the desired value (105°C). However, some dissipation by the shaft and the bearings of the machine are expected to reduce the machine temperature. Finally, a machine prototype is built and some experimental tests are carried out. The results show that the electromotive force has a similar waveform compared to 3D prediction and the difference of the measured and predicted maximum static torques is small.

The numerics of hyperbolic conservation laws, e.g., the Euler equations of gas dynamics, shallow water equations and MHD equations, is non-trivial due to the convective terms being highly non-linear and equations being coupled. Many numerical methods have been developed to solve these equations, out of which central schemes and upwind schemes (such as Flux Vector Splitting methods, Riemann solvers, Kinetic Theory based Schemes, Relaxation Schemes etc.) are well known. The majority of the above mentioned schemes give rise to very dissipative solutions. In this thesis, we propose novel low dissipative numerical algorithms for some hyperbolic conservation laws representing fluid flows. Four different and independent numerical methods which give low diffusive solutions are developed and demonstrated. The first idea is to regulate the numerical diffusion in the existing dissipative schemes so that the smearing of solution is reduced. A diffusion regulator model is developed and used along with the existing methods, resulting in crisper shock solutions at almost no added computational cost. The diffusion regulator is a function of jump in Mach number across the interface of the finite volume and the average Mach number across the surface. The introduction of the diffusion regulator makes the diffusive parent schemes to be very accurate and the steady contact discontinuities are captured exactly. The model is demonstrated in improving the diffusive Local Lax-Friedrichs (LLF) (or Rusanov) method and a Kinetic Scheme. Even when employed together with accurate methods of Roe and Osher, improvement in solutions is demonstrated for multidimensional problems. The second method, a Central Upwind-Biased Scheme (CUBS), attempts to reorganize a central scheme such that information from irrelevant directions is largely reduced and the upwind biased information is retained. The diffusion co-efficient follows a new format unlike the use of maximum characteristic speed in the Local Lax-Friedrichs method and the scheme results in improved solutions of the flow features. The grid-aligned steady contacts are captured exactly with the reorganized format of diffusion co-efficient. The stability and positivity of the scheme are discussed and the procedure is demonstrated for its ability to capture all the features of solution for different flow problems. Another method proposed in this thesis, a Central Rankine-Hugoniot Solver, attempts to integrate more physics into the discretization procedure by enforcing a simplified Rankine-Hugoniot condition which describes the jumps and hence resolves steady discontinuities very accurately. Three different variants of the scheme, termed as the Method of Optimal Viscosity for Enhanced Resolution of Shocks (MOVERS), based on a single wave (MOVERS-1), multiple waves (MOVERS-n) and limiter based diffusion (MOVERS-L) are presented. The scheme is demonstrated for scalar Burgers equation and systems of conservation laws like Euler equations, ideal Magneto-hydrodynamics equations and shallow water equations. The new scheme uniformly improves the solutions of the Local Lax-Friedrichs scheme on which it is based and captures steady discontinuities either exactly or very accurately. A Grid-Free Central Solver, which does not require a grid structure but operates on any random distribution of points, is presented. The grid-free scheme is generic in discretization of spatial derivatives with the location of the mid-point between a point and its neighbor being used to define a relevant coefficient of numerical dissipation. A new central scheme based on convective-pressure splitting to solve for mid-point flux is proposed and many test problems are solved effectively. The Rankine-Hugoniot Solver, which is developed in this thesis, is also implemented in the grid-free framework and its utility is demonstrated. The numerical methods presented are solved in a finite volume framework, except for the Grid-Free Central Solver which is a generalized finite difference method. The algorithms developed are tested on problems represented by different systems of equations and for a wide variety of flow features. The methods presented in this thesis do not need any eigen-structure and complicated flux splittings, but can still capture discontinuities very accurately (sometimes exactly, when aligned with the grid lines), yielding low dissipative solutions. The thesis ends with a highlight on the importance of developing genuinely multidimensional schemes to obtain accurate solutions for multidimensional flows. The requirement of simpler discretization framework for such schemes is emphasized in order to match the efficacy of the popular dimensional splitting schemes.

This thesis describes design and implementation of 3D space physical real-time simulation system as add-on to graphical engine. System contains processing of rigid solid simulation, collision detection and response, which are provided by physical laws. Document contains theoretical introduction of this course of study, design and implementation as multiplatform library and description of user interfaces. Output of this thesis, except library, is also example applications, which demonstrate functionality of designed simulation model.

[pt] Métodos computacionais em engenharia são usados na solução de problemas físicos que não possuem solução analítica ou sua perfeita representação matemática é inviável. Técnicas de métodos numéricos, incluindo o amplamente usado método dos elementos finitos, podem exigir a solução de sistemas lineares com centenas de milhares de equações, demandando altos recursos computacionais (memória e tempo). Nesta tese, nós apresentamos um sistema baseado em plugins para análise numérica em larga escala. O sistema é usado como uma ferramenta original na solução de problemas de otimização topológica usando o método dos elementos finitos com milhões de elementos. Nossa estratégia utiliza uma técnica elemento-por-elemento para implementar um código altamente paralelo para um solver iterativo com baixo consumo de memória. Além disso, a abordagem de plugin proporciona um ambiente completamente flexível e fácil de estender, onde diferentes aplicações, exigindo diferentes tipos de elementos finitos, materiais, solvers lineares e formulações podem ser desenvolvidos e melhorados. O kernel do sistema é mínimo, com apenas um módulo gerenciador de plugin, responsável por carregar os plugins desejados em tempo real usando um arquivo de configuração de entrada. Todas as funcionalidades necessárias para uma determinada aplicação são definidas dentro dos plugins, sem a necessidade de mudar o kernel. Plugins podem disponibilizar ou exigir interfaces adicionais especializadas, onde outros plugins podem ser conectados para compor um sistema mais complexo e completo. Nós apresentamos resultados para uma análise estrutural estática linear elástica e para uma análise estrutural de otimização topológica. As simulações utilizam elementos Q4, hexagonal (Brick8) e prisma hexagonal (Honeycomb), com solvers diretos e iterativos usando computação sequencial, paralela e distribuída. Nós investigamos o desempenho com relação ao uso de memória e escalabilidade da solução para problemas com diferentes tamanhos, de exemplos pequenos a muito grandes em apenas uma máquina e em um cluster. Foi simulado um exemplo de análise estática linear elástica com 500 milhões de elementos em 300 máquinas. / [en] Computational methods in engineering are used to solve physical problems that do not have analytical solution or their perfect mathematical representation is unfeasible. Numerical techniques, including the largely used finite element method, require the solution of linear systems with hundreds of thousands equations, demanding high computational resources (memory and time). In this thesis, we present a plugin-based framework for large-scale numerical analysis. The framework is used as an original tool to solve topology optimization problems using the finite element method with millions of elements. Our strategy uses an element-by-element technique to implement a highly parallel code for an iterative solver with low memory consumption. Besides, the plugin approach provides a fully flexible and easy to extend environment, where different types of applications, requiring different types of finite elements, materials, linear solvers, and formulations, can be developed and improved. The kernel of the framework is minimum with only a plugin manager module, responsible to load the desired plugins during runtime using an input configuration file. All the features required for a specific application are defined inside plugins, with no need to change the kernel. Plugins may provide or require additional specialized interfaces, where other plugins may be connected to compose a more complex and complete system. We present results for a structural linear elastic static analysis and for a structural topology optimization analysis. The simulations use elements Q4, hexahedron (Brick8), and hexagonal prism (Honeycomb), with direct and iterative solvers using sequential, parallel and distributed computing. We investigate the performance regarding the use of memory and the scalability of the solution for problems with different sizes, from small to very large examples on a single machine and on a cluster. We simulated a linear elastic static example with 500 million elements on 300 machines.

Магистерская диссертация посвящена разработке web-приложения решения задачи оптимизации состава многокомпонентной плавильной шихты. В ходе работы рассмотрены основные этапы разработки программного обеспечения: анализ предметной области; создание архитектуры программного обеспечения; разработка алгоритмического обеспечения и справочной документации; подготовка дистрибутива. В процессе выполнения работы реализованы все цели и задачи, стоящие перед проектом. Основными пользователями системы являются специалисты инженерно-технологического персонала доменных цехов, студенты. Научная новизна полученных в работе результатов заключается в применении нового метода эффективной организации и ведения специализированного алгоритмического и программного обеспечения решения задачи оптимизации расчета многокомпонентной плавильной шихты, ориентированного на повышение эффективности управления процессами получения качественных сплавов с использованием современных методов обработки информации: использование гибкой методологии разработки (Agile) и таск-трекера Atlassian JIRA для ведения проекта, взаимодействия с заказчиком во время разработки, отслеживания ошибок, визуального отображения задач и мониторинга процесса их выполнения; функциональное моделирование процессов для реализации web-приложения решения задачи оптимизации затрат на перевозку продукции на основе методологии IDEF0 и средства реализации Ramus Educational; использование методики коллективного владения программным кодом на основе сервиса (удаленного репозитория) Atlassian Bitbucket. Практическая значимость результатов заключается в том, что разработанное программное обеспечение позволит: производить расчёт оптимального состава многокомпонентной плавильной шихты; инженерно-технологическому персоналу литейных цехов металлургических предприятий сократить время на выполнение расчетов состава многокомпонентной плавильной шихты за счет реализации эргономичного web-интерфейса; специалистам отдела сопровождения информационных систем предоставляет условия для снижения трудозатрат на сопровождение, совершенствование и развитие системы с учетом пожеланий пользователей. Результаты работы могут быть использованы также в учебном процессе для обучения бакалавров и магистрантов по направлению «Информационные системы и технологии». / The master's thesis is devoted to the development of a web application for solving the problem of optimizing the composition of a multicomponent melting mixture. In the course of the work, the main stages of software development were considered: analysis of the subject area; creation of software architecture; development of algorithmic support and reference documentation; distribution kit preparation. In the process of performing the work, all the goals and objectives of the project have been realized. The main users of the system are specialists of the engineering and technological personnel of blast-furnace shops, students. The scientific novelty of the results obtained in the work lies in the application of a new method of effective organization and maintenance of specialized algorithmic and software for solving the problem of optimizing the calculation of a multicomponent melting charge, focused on increasing the efficiency of control of the processes of obtaining high-quality alloys using modern information processing methods: use of flexible development methodology (Agile) and the Atlassian JIRA task tracker for project management, interaction with the customer during development, tracking errors, visual display of tasks and monitoring the process of their implementation; functional modeling of processes for the implementation of a web-application for solving the problem of optimizing the costs of transportation of products based on the IDEF0 methodology and Ramus Educational tools; using the method of collective ownership of the program code based on the service (remote repository) Atlassian Bitbucket. The practical significance of the results lies in the fact that the developed software will allow: to calculate the optimal composition of the multicomponent melting mixture; to the engineering and technological personnel of foundries of metallurgical enterprises to reduce the time for performing calculations of the composition of a multicomponent melting mixture by implementing an ergonomic web interface; for specialists of the information systems support department, it provides conditions for reducing labor costs for maintaining, improving and developing the system, taking into account the wishes of users. The results of the work can also be used in the educational process for training bachelors and undergraduates in the direction of "Information systems and technologies".

Объектом исследования является процесс транспортных перевозок. Предметом исследования выступают пункты потребления и пункты производства, автоматизация системы расчета оптимальной стоимости перевозки. Поставленные задачи: 1. Возможность ввода, корректировки и сохранения вариантов расчёта по оптимизации. 2. Отображение результатов расчета в графическом виде на пользовательской форме. Целью данной работы является создание информационного Web-приложения, который позволит рассчитывать оптимальную стоимость перевозки продукции, предоставлять пользователю результаты расчета в графическом виде. Научная новизна полученных в работе результатов заключается в применении нового метода эффективной организации и ведения специализированного алгоритмического и программного обеспечения решения задачи оптимизации затрат на перевозку продукции, ориентированного на повышение эффективности управления процессами грузоперевозок с использованием современных методов обработки информации: использование гибкой методологии разработки (Agile) и таск-трекера Atlassian JIRA для ведения проекта, взаимодействия с заказчиком во время разработки, отслеживания ошибок, визуального отображения задач и мониторинга процесса их выполнения; функциональное моделирование процессов для реализации web-приложения решения задачи оптимизации затрат на перевозку продукции на основе методологии IDEF0 и средства реализации Ramus Educational; использование методики коллективного владения программным кодом на основе сервиса (удаленного репозитория) Atlassian Bitbucket. Практическая значимость результатов заключается в том, что разработанное программное обеспечение позволит: производить расчёт оптимальной себестоимости транспортных перевозок для любого количества пунктов производства; специалистам транспортно-логистического операционного отдела сократить время на формирование отчетных документов, сократить время поиска необходимой фактической отчетной информации за счет реализации эргономичного web-интерфейса; специалистам отдела сопровождения информационных систем предоставляет условия для снижения трудозатрат на сопровождение, совершенствование и развитие системы с учетом пожеланий пользователей. Результаты работы могут быть использованы также в учебном процессе для обучения бакалавров и магистрантов по направлению «Информационные системы и технологии». / The object of the research is the process of transportation. The subject of the research is points of consumption and points of production, automation of the system for calculating the optimal cost of transportation. Assigned tasks: 1. Possibility of entering, adjusting and saving options for the calculation of optimization. 2. Displaying the calculation results in a graphical form on the user form. The purpose of this work is to create an information Web-application that will allow you to calculate the optimal cost of transportation of products, provide the user with the results of the calculation in a graphical form. The scientific novelty of the results obtained in the work lies in the application of a new method of effective organization and maintenance of specialized algorithmic and software solutions for the optimization of the cost of transportation of products, focused on improving the efficiency of management of cargo transportation processes using modern information processing methods: the use of flexible development methodology (Agile) and the Atlassian JIRA task tracker for project management, interaction with the customer during development, tracking errors, visual display of tasks and monitoring the process of their implementation; functional modeling of processes for the implementation of a web-application for solving the problem of optimizing the costs of transportation of products based on the IDEF0 methodology and Ramus Educational tools; using the method of collective ownership of the program code based on the service (remote repository) Atlassian Bitbucket. The practical significance of the results lies in the fact that the developed software will allow: to calculate the optimal cost of transportation for any number of points of production; for specialists of the transport and logistics operations department, to reduce the time for the formation of reporting documents, to reduce the time to search for the necessary actual reporting information due to the implementation of an ergonomic web interface; for specialists of the information systems support department, it provides conditions for reducing labor costs for maintaining, improving and developing the system, taking into account the wishes of users. The results of the work can also be used in the educational process for training bachelors and undergraduates in the direction "Information systems and technologies".

<p>Combustion instability, a complex phenomenon observed in combustion chambers is due to the coupling between heat release and other unsteady flow processes. Combustion instability has long been a topic of interest to rocket scientists and has been extensively investigated experimentally and computationally. However, to date, there is no computational tool that can accurately predict the combustion instabilities in full-size combustors because of the amount of computational power required to perform a high-fidelity simulation of a multi-element chamber. Hence, the focus is shifted to reduced fidelity computational tools which may accurately predict the instability by using the information available from the high-fidelity simulations or experiments of single or few-element combustors. One way of developing reduced fidelity computational tools involves using a reduced fidelity solver together with the flame transfer functions that carry important information about the flame behavior from a high-fidelity simulation or experiment to a reduced fidelity simulation.</p> <p> </p> <p>To date, research has been focused mainly on premixed flames and using acoustic solvers together with the global flame transfer functions that were obtained by integrating over a region. However, in the case of rockets, the flame is non-premixed and distributed in space and time. Further, the mixing of propellants is impacted by the level of flow fluctuations and can lead to non-uniform mean properties and hence, there is a need for reduced fidelity solver that can capture the gas dynamics, nonlinearities and steep-fronted waves accurately. Nonlinear Euler equations have all the required capabilities and are at the bottom of the list in terms of the computational cost among the solvers that can solve for mean flow and allow multi-dimensional modeling of combustion instabilities. Hence, in the current work, nonlinear Euler solver together with the spatially distributed local flame transfer functions that capture the coupling between flame, acoustics, and hydrodynamics is explored.</p> <p> </p> <p>In this thesis, the approach to extract flame transfer functions from high-fidelity simulations and their integration with nonlinear Euler solver is presented. The dynamic mode decomposition (DMD) was used to extract spatially distributed flame transfer function (FTF) from high fidelity simulation of a single element non-premixed flame. Once extracted, the FTF was integrated with nonlinear Euler equations as a fluctuating source term of the energy equation. The time-averaged species destruction rates from the high-fidelity simulation were used as the mean source terms of the species equations. Following a variable gain approach, the local species destruction rates were modified to account for local cell constituents and maintain correct mean conditions at every time step of the nonlinear Euler simulation. The proposed reduced fidelity model was verified using a Rijke tube test case and to further assess the capabilities of the proposed model it was applied to a single element model rocket combustor, the Continuously Variable Resonance Combustor (CVRC), that exhibited self-excited combustion instabilities that are on the order of 10% of the mean pressure. The results showed that the proposed model could reproduce the unsteady behavior of the CVRC predicted by the high-fidelity simulation reasonably well. The effects of control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size are also studied. The reduced fidelity model could reproduce the limit cycle amplitude within a few percent of the mean pressure. The successful constraints on the model include good spatial resolution and FTF with all modes up to at least one dominant frequency higher than the frequencies of interest. Furthermore, the reduced fidelity model reproduced consistent mode shapes and linear growth rates that reasonably matched the experimental observations, although the apparent ability to match growth rates needs to be better understood. However, the presence of significant heat release near a pressure node of a higher harmonic mode was found to be an issue. This issue was rectified by expanding the pressure node of the higher frequency mode. Analysis of two-dimensional effects and coupling between the local pressure and heat release fluctuations showed that it may be necessary to use two dimensional spatially distributed local FTFs for accurate prediction of combustion instabilities in high energy devices such as rocket combustors. Hybrid RANS/LES-FTF simulation of the CVRC revealed that it might be necessary to use Flame Describing Function (FDF) to capture the growth of pressure fluctuations to limit cycle when Navier-Stokes solver is used.</p> <p> </p> <p>The main objectives of this thesis are:</p> <p>1. Extraction of spatially distributed local flame transfer function from the high fidelity simulation using dynamic mode decomposition and its integration with nonlinear Euler solver</p> <p>2. Verification of the proposed approach and its application to the Continuously Variable Resonance Combustor (CVRC).</p> <p>3. Sensitivity analysis of the reduced fidelity model to control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size.</p> <p> </p> <p>The goal of this thesis is to contribute towards a reduced fidelity computational tool which can accurately predict the combustion instabilities in practical systems using flame transfer functions, by providing a path way for reduced fidelity multi-element simulation, and by defining the limitations associated with using flame transfer functions and nonlinear Euler equations for non-premixed flames.</p> <p> </p><br>

Cette étude s’inscrit dans le domaine de la simulation électromagnétique des problèmes de grande taille tels que la diffraction d’ondes planes par de larges plateformes et le rayonnement d’antennes aéroportées. Elle consiste à développer une méthode combinant décomposition en sous-domaines et compression hiérarchique des équations intégrales de frontière. Pour cela, nous rappelons dans un premier temps les points importants de la méthode des équations intégrales de frontière et de leur compression hiérarchique par l’algorithme ACA (Adaptive Cross Approximation). Ensuite, nous présentons la formulation IE-DDM (Integral Equations – Domain Decomposition Method) obtenue à partir d’une représentation intégrale des sous-domaines. Les matrices résultant de la discrétisation de cette formulation sont stockées au format H-matrice (matricehiérarchique). Un solveur spécialement adapté à la résolution de la formulation IE-DDM et à sa représentation hiérarchique a été conçu. Cette étude met en évidence l’efficacité de la décomposition en sous-domaines en tant que préconditionneur des équations intégrales. De plus, la méthode développée est rapide pour la résolution des problèmes à incidences multiples ainsi que la résolution des problèmes basses fréquences / This thesis is about the electromagnetic simulation of large scale problems as the wave scattering from aircrafts and the airborne antennas radiation. It consists in the development of a method combining domain decomposition and hierarchical compression of the surface integral equations. First, we remind the principles of the boundary element method and the hierarchical representation of the surface integral equations with the Adaptive Cross Approximation algorithm. Then, we present the IE-DDM formulation obtained from a sub-domain integral representation. The matrices resulting of the discretization of the formulation are stored in the H-matrix format. A solver especially fitted with the hierarchical representation of the IE-DDM formulation has been developed. This study highlights the efficiency of the sub-domain decomposition as a preconditioner of the integral equations. Moreover, the method is fast for the resolution of multiple incidences and the resolution of low frequencies problems

La thèse vise à intégrer des critères de transition dans le solveur Navier-Stokes non structuré Aether utilisé chez Dassault Aviation. Une méthodologie de prévision de la transition laminaire/turbulent a été élaborée et implémentée dans le solveur RANS Aether. Deux stratégies de calcul de transition ont été testées. Soit Aether est couplé avec le code de couche limite de l’ONERA 3C3D. Soit la position de transition est calculée en utilisant directement les profils de vitesse RANS. Les deux méthodes ont été testées pour des écoulements subsoniques et transsoniques. L’influence des solveurs numériques, des critères de transition et du processus de couplage sont étudiés. L’utilisation de schémas numériques d’ordre élevé est également considérée. / This thesis present the state-of-the-art of the transition prediction numerical chain which has been developed at Dassault Aviation in the RANS solver AETHER. Two strategies for transition location estimations exist. First, AETHER is coupled with the ONERA boundary layer code 3C3D. Second, the transition location is computed by using directly the RANS velocity profiles. Both methods were preliminarily tested in subsonic and transonic. The issues of the influence of the numerical solvers, transition onset criteria and coupling process are studied. The influence of higher order numerical method, are also considered.

Cette thèse se concentre sur la résolution de systèmes linéaires creux dans le contexte d’applications massivement parallèles. Ce type de problèmes s’exprime sous la forme AX=B, où A est une matrice creuse d’ordre n x n, i.e. qui possède un nombre d’entrées nulles suffisamment élevé pour pouvoir être exploité, et B et X sont respectivement la matrice de seconds membres et la matrice de solution de taille n x nrhs. Cette résolution par des méthodes dites directes est effectuée grâce à une étape de factorisation qui réduit A en deux matrices triangulaires inférieure et supérieure L et U, suivie de deux résolutions triangulaires pour calculer la solution.Nous nous intéressons à ces résolutions avec une attention particulière apportée à la première, LY=B. Dans beaucoup d’applications, B possède un grand nombre de colonnes (nrhs >> 1) transformant la phase de résolution en un goulot d’étranglement. Elle possède souvent aussi une structure creuse, donnant l’opportunité de réduire la complexité de cette étape.Cette étude aborde sous des angles complémentaires la résolution triangulaire de systèmes linéaires avec seconds membres multiples et creux. Nous étudions dans un premier temps la complexité asymptotique de cette étape dans différents contextes (2D, 3D, facteurs compressés ou non). Nous considérons ensuite l’exploitation de cette structure et présentons de nouvelles approches s’appuyant sur une modélisation du problème par des graphes qui permettent d’atteindre efficacement le nombre minimal d’opérations. Enfin, nous donnons une interprétation concrète de son exploitation sur une application d’électromagnétisme pour la géophysique. Nous adaptons aussi des algorithmes parallèles aux spécificités de la phase de résolution.Nous concluons en combinant l'ensemble des résultats précédents et en discutant des perspectives de ce travail. / We consider direct methods to solve sparse linear systems AX = B, where A is a sparse matrix of size n x n with a symmetric structure and X and B are respectively the solution and right-hand side matrices of size n x nrhs. A is usually factorized and decomposed in the form LU, where L and U are respectively a lower and an upper triangular matrix. Then, the solve phase is applied through two triangular resolutions, named respectively the forward and backward substitutions.For some applications, the very large number of right-hand sides (RHS) in B, nrhs >> 1, makes the solve phase the computational bottleneck. However, B is often sparse and its structure exhibits specific characteristics that may be efficiently exploited to reduce this cost. We propose in this thesis to study the impact of the exploitation of this structural sparsity during the solve phase going through its theoretical aspects down to its actual implications on real-life applications.First, we investigate the asymptotic complexity, in the big-O sense, of the forward substitution when exploiting the RHS sparsity in order to assess its efficiency when increasing the problem size. In particular, we study on 2D and 3D regular problems the asymptotic complexity both for traditional full-rank unstructured solvers and for the case when low-rank approximation is exploited. Next, we extend state-of-the-art algorithms on the exploitation of RHS sparsity, and also propose an original approach converging toward the optimal number of operations while preserving performance. Finally, we show the impact of the exploitation of sparsity in a real-life electromagnetism application in geophysics that requires the solution of sparse systems of linear equations with a large number of sparse right-hand sides. We also adapt the parallel algorithms that were designed for the factorization to solve-oriented algorithms.We validate and combine the previous improvements using the parallel solver MUMPS, conclude on the contributions of this thesis and give some perspectives.