Dans cette thèse nous mettons en oeuvre une méthode vortex hybride pénalisée (HVP) afin desimuler des écoulements incompressibles autour de corps non profilés dans des milieux complexessolides-fluides-poreux. Avec cette approche particulaire hybride, le phénomène de convection estmodélisé à l’aide d’une méthode vortex afin de bénéficier du caractère peu diffusif et naturel desméthodes particulaires. Un remaillage des particules est alors réalisé systématiquement sur unegrille cartésienne sous-jacente afin d’éviter les phénomènes de distorsion. D’autre part, les effetsdiffusifs et d’étirement ainsi que le calcul de la vitesse sont traités sur la grille cartésienne, àl’aide de schémas eulériens. Le traitement des conditions de bords aux parois de l’obstacle esteffectué à l’aide d’une technique de pénalisation, particulièrement bien adaptée au traitementde milieux solides-fluides-poreux.Cette méthode HVP est appliquée au contrôle passif d’écoulement. Cette étude de contrôleest effectuée respectivement en 2D et en 3D autour d’un demi-cylindre et d’un hémisphère parl’ajout d’un revêtement poreux à la surface de l’obstacle. La présence de cette couche poreusemodifiant la nature des conditions aux interfaces, permet de régulariser l’écoulement global etde diminuer la traînée aérodynamique de l’obstacle contrôlé. A travers des études paramétriquessur la perméabilité, l’épaisseur et la position du revêtement poreux, ce travail vise à identifier desdispositifs de contrôles efficaces pour des écoulements autour d’obstacles comme des rétroviseursautomobiles. / In this work we use a hybrid vortex penalization method (HVP) to simulate incompressibleflows past bluff bodies in complex solid-fluid-porous media. In this hybrid particle approach,the advection phenomenon is modeled through a vortex method in order to benefit from thenatural description of the flow supplied by particle methods and their low numerical diffusionfeatures. A particle remeshing is performed systematically on an underlying Cartesian grid inorder to prevent distortion phenomena. On the other hand, the viscous and stretching effects aswell as the velocity calculation are discretized on the mesh through Eulerian schemes. Finally,the treatment of boundary conditions is handled with a penalization method that is well suitedfor the treatment of solid-fluid-porous media.The HVP method is applied to passive flow control. This flow control study is realized pasta 2D semi-circular cylinder and a 3D hemisphere by adding a porous layer on the surface of thebody. The presence of such porous layer modifies the characteristics of the conditions at theinterfaces and leads to a regularization of the wake and to a decrease of the aerodynamic dragof the controlled obstacle. Through parametric studies on the permeability, the thickness andthe position of the porous coating, this works aims to identify efficient control devices for flowsaround obstacles like the rear-view mirrors of a ground vehicle.

O presente trabalho tem como objetivo realizar um estudo comparativo entre os três métodos para obtenção das características aerodinâmicas de tabuleiros de pontes: analítico, experimental e numérico. O tema é motivado pelo crescimento das necessidades de conhecimento científico e tecnológico na área de desempenho aerodinâmico de seções típicas de pontes nas etapas mais iniciais de projeto. A precisão de uma ferramenta de previsão certamente acelera a convergência ao modelo final; porém, qualquer metodologia numérica deve ser extensivamente testada antes de ser utilizada nos escritórios de cálculo. Inicialmente mostra-se a lógica da análise de instabilidades provocadas pelo vento, resume-se a modelagem de tabuleiros de pontes para ensaios em túnel de vento e descreve-se o método numérico implementado no programa comercial de análise de pontes RM. Comparam-se dados disponíveis de testes de modelos reduzidos no túnel de vento com as estimativas de aproximação do Método da Vorticidade em Partículas (MVP). Confrontam-se resultados para casos clássicos, um estudo paramétrico, dados dos relatórios de três pontes Brasileiras já ensaiadas no Túnel de Vento Joaquim Blessmann do LAC – PPGEC/UFRGS: Paulicéia, Guamá e Roberto Marinho. Finalmente, considera-se o método experimental como preciso e julga-se o método numérico comparativamente, analisando a confiabilidade dos resultados e procurando definir a sua faixa de aplicabilidade. Examinando a variabilidade dos resultados, bastante baixa, definiu-se sua faixa de aplicação considerando os resultados obtidos como satisfatórios para projetos básicos e executivos de estruturas; necessitando, porém, para os casos especiais, de comprovação experimental em túnel de vento antes da sua execução. Exemplifica-se então o uso do método para alguns problemas de análises de seções. Os resultados e comparações corroboram verificações anteriores do método dos vórtices discretos para verdadeiras seções de pontes e mostram porque o seu uso vem se tornando cada vez mais aceitável. / The work has the objective to realize a comparative study between three methods to obtain bridge aerodynamic coefficients: analytic, experimental and numerical. The theme is motivated by the increasing need of knowledge about the aerodynamic behaviour of bridge decks in the early stages of the design process. The accuracy of a predictive tool will certainly accelerate the convergence to the final design. However, any numerical methodology must be extensively tested ascertain and hence validated before it becomes of current use by designers. Collaboration is made comparing wind tunnel data of reduced models and numerical estimations by a simplified approach. After a brief overview of wind loading and wind tunnel modeling, the bridge wind resistance design in the commercially available bridge design software package RM is shown. The Discrete Vortex Method and the implementation are shortly described. Initially, classical fluid problems are indicated together with the comparison of parametric cross sectional bridge shapes. After, three Brazilian large cablestayed bridge projects are presented as instance of analysis: Paulicéia, Guamá and Roberto Marinho. Bridges aerodynamic coefficients were experimentally obtained at the Boundary Layer Wind Tunnel Joaquim Blessmann at Universidade Federal do Rio Grande do Sul (UFRGS). Finally, experimental results are considered precise and those results evaluated with CFD are compared analyzing the liability and applicability. According to the results variability, quite low, we considered them appropriated for basic and final structural designs; needing, although, for special cases, of experimental probation in wind tunnel tests before the construction phase. Some CFD analysis problems of cross sections are then presented as application examples. The results and similitude corroborate earlier verifications of the method and this implementation, demonstrating why its use is becoming increasingly acceptable.

O presente trabalho tem como objetivo realizar um estudo comparativo entre os três métodos para obtenção das características aerodinâmicas de tabuleiros de pontes: analítico, experimental e numérico. O tema é motivado pelo crescimento das necessidades de conhecimento científico e tecnológico na área de desempenho aerodinâmico de seções típicas de pontes nas etapas mais iniciais de projeto. A precisão de uma ferramenta de previsão certamente acelera a convergência ao modelo final; porém, qualquer metodologia numérica deve ser extensivamente testada antes de ser utilizada nos escritórios de cálculo. Inicialmente mostra-se a lógica da análise de instabilidades provocadas pelo vento, resume-se a modelagem de tabuleiros de pontes para ensaios em túnel de vento e descreve-se o método numérico implementado no programa comercial de análise de pontes RM. Comparam-se dados disponíveis de testes de modelos reduzidos no túnel de vento com as estimativas de aproximação do Método da Vorticidade em Partículas (MVP). Confrontam-se resultados para casos clássicos, um estudo paramétrico, dados dos relatórios de três pontes Brasileiras já ensaiadas no Túnel de Vento Joaquim Blessmann do LAC – PPGEC/UFRGS: Paulicéia, Guamá e Roberto Marinho. Finalmente, considera-se o método experimental como preciso e julga-se o método numérico comparativamente, analisando a confiabilidade dos resultados e procurando definir a sua faixa de aplicabilidade. Examinando a variabilidade dos resultados, bastante baixa, definiu-se sua faixa de aplicação considerando os resultados obtidos como satisfatórios para projetos básicos e executivos de estruturas; necessitando, porém, para os casos especiais, de comprovação experimental em túnel de vento antes da sua execução. Exemplifica-se então o uso do método para alguns problemas de análises de seções. Os resultados e comparações corroboram verificações anteriores do método dos vórtices discretos para verdadeiras seções de pontes e mostram porque o seu uso vem se tornando cada vez mais aceitável. / The work has the objective to realize a comparative study between three methods to obtain bridge aerodynamic coefficients: analytic, experimental and numerical. The theme is motivated by the increasing need of knowledge about the aerodynamic behaviour of bridge decks in the early stages of the design process. The accuracy of a predictive tool will certainly accelerate the convergence to the final design. However, any numerical methodology must be extensively tested ascertain and hence validated before it becomes of current use by designers. Collaboration is made comparing wind tunnel data of reduced models and numerical estimations by a simplified approach. After a brief overview of wind loading and wind tunnel modeling, the bridge wind resistance design in the commercially available bridge design software package RM is shown. The Discrete Vortex Method and the implementation are shortly described. Initially, classical fluid problems are indicated together with the comparison of parametric cross sectional bridge shapes. After, three Brazilian large cablestayed bridge projects are presented as instance of analysis: Paulicéia, Guamá and Roberto Marinho. Bridges aerodynamic coefficients were experimentally obtained at the Boundary Layer Wind Tunnel Joaquim Blessmann at Universidade Federal do Rio Grande do Sul (UFRGS). Finally, experimental results are considered precise and those results evaluated with CFD are compared analyzing the liability and applicability. According to the results variability, quite low, we considered them appropriated for basic and final structural designs; needing, although, for special cases, of experimental probation in wind tunnel tests before the construction phase. Some CFD analysis problems of cross sections are then presented as application examples. The results and similitude corroborate earlier verifications of the method and this implementation, demonstrating why its use is becoming increasingly acceptable.

O presente trabalho tem como objetivo realizar um estudo comparativo entre os três métodos para obtenção das características aerodinâmicas de tabuleiros de pontes: analítico, experimental e numérico. O tema é motivado pelo crescimento das necessidades de conhecimento científico e tecnológico na área de desempenho aerodinâmico de seções típicas de pontes nas etapas mais iniciais de projeto. A precisão de uma ferramenta de previsão certamente acelera a convergência ao modelo final; porém, qualquer metodologia numérica deve ser extensivamente testada antes de ser utilizada nos escritórios de cálculo. Inicialmente mostra-se a lógica da análise de instabilidades provocadas pelo vento, resume-se a modelagem de tabuleiros de pontes para ensaios em túnel de vento e descreve-se o método numérico implementado no programa comercial de análise de pontes RM. Comparam-se dados disponíveis de testes de modelos reduzidos no túnel de vento com as estimativas de aproximação do Método da Vorticidade em Partículas (MVP). Confrontam-se resultados para casos clássicos, um estudo paramétrico, dados dos relatórios de três pontes Brasileiras já ensaiadas no Túnel de Vento Joaquim Blessmann do LAC – PPGEC/UFRGS: Paulicéia, Guamá e Roberto Marinho. Finalmente, considera-se o método experimental como preciso e julga-se o método numérico comparativamente, analisando a confiabilidade dos resultados e procurando definir a sua faixa de aplicabilidade. Examinando a variabilidade dos resultados, bastante baixa, definiu-se sua faixa de aplicação considerando os resultados obtidos como satisfatórios para projetos básicos e executivos de estruturas; necessitando, porém, para os casos especiais, de comprovação experimental em túnel de vento antes da sua execução. Exemplifica-se então o uso do método para alguns problemas de análises de seções. Os resultados e comparações corroboram verificações anteriores do método dos vórtices discretos para verdadeiras seções de pontes e mostram porque o seu uso vem se tornando cada vez mais aceitável. / The work has the objective to realize a comparative study between three methods to obtain bridge aerodynamic coefficients: analytic, experimental and numerical. The theme is motivated by the increasing need of knowledge about the aerodynamic behaviour of bridge decks in the early stages of the design process. The accuracy of a predictive tool will certainly accelerate the convergence to the final design. However, any numerical methodology must be extensively tested ascertain and hence validated before it becomes of current use by designers. Collaboration is made comparing wind tunnel data of reduced models and numerical estimations by a simplified approach. After a brief overview of wind loading and wind tunnel modeling, the bridge wind resistance design in the commercially available bridge design software package RM is shown. The Discrete Vortex Method and the implementation are shortly described. Initially, classical fluid problems are indicated together with the comparison of parametric cross sectional bridge shapes. After, three Brazilian large cablestayed bridge projects are presented as instance of analysis: Paulicéia, Guamá and Roberto Marinho. Bridges aerodynamic coefficients were experimentally obtained at the Boundary Layer Wind Tunnel Joaquim Blessmann at Universidade Federal do Rio Grande do Sul (UFRGS). Finally, experimental results are considered precise and those results evaluated with CFD are compared analyzing the liability and applicability. According to the results variability, quite low, we considered them appropriated for basic and final structural designs; needing, although, for special cases, of experimental probation in wind tunnel tests before the construction phase. Some CFD analysis problems of cross sections are then presented as application examples. The results and similitude corroborate earlier verifications of the method and this implementation, demonstrating why its use is becoming increasingly acceptable.

La simulation des grandes échelles (SGE ou LES pour large eddy simulation) commence à être très utilisée dans l’industrie. Par résolution directe des structures turbulents de grande tailles, le calcul LES est capable de calculer le bruit générée par la voilure ou de prédire avec précision le décollement de coin dans une configuration très simplifiée du compresseur. L’un des problèmes les plus importants pour effectuer un calcul LES est de fournir des conditions d’entrée avec des champs turbulents.Pour une approche hybride RANS/LES (RANS pour Reynolds Averaged Navier-Stokes), les conditions d’entrée turbulentes pour un calcul LES sont générées à l’aide des solutions fournies par le calcul RANS en amont. Il existe plusieurs méthodes pour générer les conditions d’entrée pour LES. Elles peuvent principalement être classées en deux catégories : 1) simulation avec pré-calcul ; 2) la méthode de turbulence synthétique. La simulation avec pré-calcul consiste à effectuer un calcul LES indépendant pour générer un champ turbulent comme conditions d’entrée pour alimenter le calcul principal. Cette méthode peut obtenir des turbulences de haute qualité, mais elle augmente considérablement le temps de calcul et le stockage des données. Le champ turbulent généré par la méthode de turbulence synthétique exige une « distance de adaptation », pendante laquelle le champ turbulent devient pleinement développé. L’objectif principal pour améliorer ce genre de méthodes est donc de diminuer cette distance nécessaire.Dans cette thèse, la méthode de vortex, qui est une approche de turbulence synthétique, est présentée et améliorée. A travers des expériences numériques, les paramètres de la méthode de vortex améliorée sont systématiquement optimisés. L’application à l’écoulement en canal plan et à couche limite en plaque plane, montrent que la méthode de vortex améliorée génère de manière efficace pour fournir des conditions d’entrée pour LES. Dans le cas de l’écoulement en canal plan, la distance d’adaptation nécessaire pour le rétablissement de la turbulence est d’environ 6 fois la demi-hauteur du canal. Pour le cas de l’écoulement en plaque plane, cette distance est environ 21 fois l’épaisseur de la couche limite. Enfin, dans le but de qualifier la turbulence obtenue par des calculs LES, nous utilisons les coefficients de dissymétrie des dérivées des fluctuations de vitesse, et, nous les introduisons comme un nouveau critère pour la qualité de LES. / Large eddy simulation is becoming an important numerical tool in industry recently. Resolving large scale turbulent motions directly, LES is capable to compute the aeroacoustic noise generated by the airfoil or to precisely capture the corner separation in a linear compressor cascade. The main challenge to perform a LES calculation is to prescribe a realistic unsteady inflow field. For hybrid RANS/LES approaches, inflow conditions for downstream LES region must be generated from the upstream RANS solutions. There exist several methods to generate inflow conditions for LES. They can mainly be divided into two categories: 1) Precursor simulation; 2) Synthetic turbulence methods. Precursor simulation requires to run a separate calculation to generate a turbulent ow or a database to feed the main computation. This kind of methods can generate high quality turbulence. However, it requires heavy extra computing load. Synthetic turbulence methods consist in generating a fluctuating velocity field, and within a short “adaptation distance”, the field get fully developed. So main goal of synthetic turbulence methods is to decrease the required adaptation distance. The vortex method which is a synthetic turbulence method is presented and improved here. Parameters of the improved vortex method are optimized systematically with a series of calculations in this thesis. Applications on channel and flat-plate flows show that the improved vortex method is effective in generating the LES inflow conditions. The adaptation distance required for turbulence recovery is about 6 times the half channel height for channel flow, and 21 times the boundary-layer thickness (at the inlet of vortex) for at-plate ow. The velocity-derivative skewness is used to qualify the generated turbulence, and is introduced as a new criterion of LES calculation.

Neste trabalho, a resposta dinâmica de um riser marítimo devido à geração e desprendimento alternado de vórtices é investigada numericamente. O riser é dividido em seções bidimensionais ao longo de seu comprimento. O Método dos Vórtices Discretos é empregado para a determinação das forças hidrodinâmicas que agem nestas seções bidimensionais. As seções hidrodinâmicas são resolvidas independentemente, e o acoplamento entre as mesmas é feito através da solução da estrutura no domínio do tempo pelo Método dos Elementos Finitos. Os resultados numéricos são comparados com resultados obtidos experimentalmente. Processamento paralelo é empregado para melhorar a performance do método. As simulações são realizadas através de uma metodologia mestre-escravo, utilizando MPI Message Passing Interface para explorar o paralelismo. A escalabilidade do algoritmo é mostrada e discutida. Este trabalho representa o desenvolvimento de um simulador que permite, efetivamente, a análise dinâmica de um riser com características e dimensões representativas das condições reais encontradas em campo, a um custo computacional factível para seu uso como uma ferramenta de engenharia. Isto é obtido por meio da técnica de processamento paralelo, aliada à solução do escoamento através de um método eficiente de CFD Método dos Vórtices Discretos e à solução da estrutura através do Método dos Elementos Finitos. / In this work the dy6namic response of a marine riser due to vortex shedding is numerically investigated. The riser is divided in two-dimensional sections along the riser length. The Discrete Vortex Method is employed for the assessment of the hydrodynamic forces acting on these two-dimensional sections. The hydrodynamic sections are solved independently, and the coupling among the sections is taken into account by the solution of the structure in the time domain by the Finite Element Method. The numerical results are compared with results obtained experimentally. Parallel processing is employed to improve the performance of the method. The simulations are carried out through a master-slave approach using MPI Message Passing Interface to exploit the parallelism. Scalability of the algorithm is shown and discussed. This work represents the development of a simulator that effectively allows the dynamic analysis of a riser with representative characteristics and dimensions of real field conditions, with a feasible computational cost for its use as an engineering tool. This is obtained by means of the parallel processing technique, together with an efficient CFD solution of the flow with de Discrete Vortex Method and the solution of the structure with the Finite Element Method.

Neste trabalho, a resposta dinâmica de um riser marítimo devido à geração e desprendimento alternado de vórtices é investigada numericamente. O riser é dividido em seções bidimensionais ao longo de seu comprimento. O Método dos Vórtices Discretos é empregado para a determinação das forças hidrodinâmicas que agem nestas seções bidimensionais. As seções hidrodinâmicas são resolvidas independentemente, e o acoplamento entre as mesmas é feito através da solução da estrutura no domínio do tempo pelo Método dos Elementos Finitos. Os resultados numéricos são comparados com resultados obtidos experimentalmente. Processamento paralelo é empregado para melhorar a performance do método. As simulações são realizadas através de uma metodologia mestre-escravo, utilizando MPI Message Passing Interface para explorar o paralelismo. A escalabilidade do algoritmo é mostrada e discutida. Este trabalho representa o desenvolvimento de um simulador que permite, efetivamente, a análise dinâmica de um riser com características e dimensões representativas das condições reais encontradas em campo, a um custo computacional factível para seu uso como uma ferramenta de engenharia. Isto é obtido por meio da técnica de processamento paralelo, aliada à solução do escoamento através de um método eficiente de CFD Método dos Vórtices Discretos e à solução da estrutura através do Método dos Elementos Finitos. / In this work the dy6namic response of a marine riser due to vortex shedding is numerically investigated. The riser is divided in two-dimensional sections along the riser length. The Discrete Vortex Method is employed for the assessment of the hydrodynamic forces acting on these two-dimensional sections. The hydrodynamic sections are solved independently, and the coupling among the sections is taken into account by the solution of the structure in the time domain by the Finite Element Method. The numerical results are compared with results obtained experimentally. Parallel processing is employed to improve the performance of the method. The simulations are carried out through a master-slave approach using MPI Message Passing Interface to exploit the parallelism. Scalability of the algorithm is shown and discussed. This work represents the development of a simulator that effectively allows the dynamic analysis of a riser with representative characteristics and dimensions of real field conditions, with a feasible computational cost for its use as an engineering tool. This is obtained by means of the parallel processing technique, together with an efficient CFD solution of the flow with de Discrete Vortex Method and the solution of the structure with the Finite Element Method.

Made available in DSpace on 2014-12-17T14:58:12Z (GMT). No. of bitstreams: 1 ElsonANM.pdf: 918866 bytes, checksum: 54ffb16de332a5a9c76c3afd577bcea1 (MD5) Previous issue date: 2005-09-16 / Conselho Nacional de Desenvolvimento Cient?fico e Tecnol?gico / The study of aerodynamic loading variations has many engineering applications, including helicopter rotor blades, wind turbines and turbo machinery. This work uses a Vortex Method to make a lagrangian description of the a twodimensional airfoil/ incident wake vortex interaction. The flow is incompressible, newtonian, homogeneus and the Reynolds Number is 5x105 .The airfoil is a NACA 0018 placed a angle of attack of the 0? and 5?simulates with the Painel Method with a constant density vorticity panels and a generation poit is near the painel. The protector layer is created does not permit vortex inside the body. The vortex Lamb convection is realized with the Euler Method (first order) and Adans-Bashforth (second order). The Random Walk Method is used to simulate the diffusion. The circular wake has 366 vortex all over positive or negative vorticity located at different heights with respect to the airfoil chord. The Lift was calculated based in the algorithm created by Ricci (2002). This simulation uses a ready algorithm vatidated with single body does not have a incident wake. The results are compared with a experimental work The comparasion concludes that the experimental results has a good agrement with this papper / Este estudo simula numericamente, atrav?s do M?todo de V?rtices, a intera??o entre uma nuvem de circula??o incidente e um aerof?lio NACA 0018 posicionado com ?ngulos de ataque de 0? e 5?. O n?mero de Reynolds utilizado ? de 5x105 e o escoamento ? considerado incompress?vel. Devido a complexidade do fen?meno de geometria tridimensional, apenas uma se??o do aerof?lio ? simulada transformando o trabalho em bidimensional. Utiliza-se o M?todo dos Pain?is com vorticidades distribu?das para simular o corpo e o ponto de gera??o de vorticidade ? deslocado da superf?cie do painel. Uma camada protetora ? criada para evitar a entrada de v?rtices no corpo. A convec??o dos v?rtices de Lamb presentes ? realizada atrav?s dos esquemas de Euler (1?ordem) e de Adams-Bashforth (2?ordem). A difus?o utiliza o M?todo do Avan?o Rand?mico. As cargas de sustenta??o s?o calculadas atrav?s de um algoritmo desenvolvido por Ricci (2002). A nuvem utilizada tem 366 v?rtices com circula??o de +0.001375692 ou -0.001375692 formando um c?rculo de raio 0,2. A simula??o utilizou um algoritmo existente que foi modificado para o problema em quest?o. Validado atrav?s da simula??o do perfil sem a presen?a da nuvem. Os resultados encontrados mostram uma boa concord?ncia com os encontrados em trabalhos experimentais

The thesis presents a two pronged approach for predicting aerodynamics of air- foils/wings in the vicinity of the ground. The first approach is effectively a model for ground effect studies, employing an inexpensive Discrete Vortex Method for the 2D pre- dictions and the well known Numerical lifting line theory for the 3D predictions. The second one pertains to the dynamic ground effect analysis which employs the state of the art moving mesh methodology based time accurate CFD. In that sense, the thesis deals with two ends of spectrum in the ground effect analysis; one, a model to be used in the concept design phase and the other an advanced CFD tool for analysis. The proposed model for ground effect studies is based on the well known Discrete Vortex Method (DVM). An important aspect of this method is that it employs what is referred to as the Generalized Kutta Joukowski Theorem (GKJ), meant for interaction problems with multiple vortices, for predicting the lift (and drag) within a potential flow framework. After ascertaining the correctness of using the GKJ theorem for lift prediction for airfoils in ground effect, a modified DVM is presented as a model for ground effect predictions. As per this model, knowing the free stream lift and drag (either from an ex- periment or from a RANS computation) the aerodynamics of the section in ground effect can be predicted. The model is effectively built by constraining the DVM to produce the reference lift/drag in the free stream. The accuracy of the model, particularly for the more relevant high lift sections used during take-off and landing, is systematically estab- lished for a number of test cases. Knowing the sectional ground effect, the extension to 3D analysis is very simple and this is achieved through the well known Numerical Lifting Line theory. The efficacy of the proposed method for the 3D applications is demonstrated using a high lift wing in ground effect. It is worth noting that the proposed model predicts the lift and drag very accurately, practically at no computational cost as compared to modern RANS based CFD tools requiring over 40 or 50 million volumes at a high computational cost and intense human intervention for generating the grids for every ground clearance. The other aspect of the thesis pertains to what is referred to as the Dynamic Ground Effect. Normally the CFD computations mimic the ground effect experiments in simulat- ing the ground effect. These simulations do not maintain geometric similarity with the actual landing or take-off sequence of the aircrafts and this can only be achieved when the simulations are dynamic. Dynamics is also important in case of combat aircrafts (particularly their naval versions) with an aggressive landing and take-off. The dynamic ground effect simulations also provides a framework for simulating varied gust conditions. This dynamic simulation of the ground effect is accomplished using a novel sinking grid methodology, which allows the grids to sink in the ground as the aircraft approaches the ground along the glide path. These simulations make use of the state of the art, time accurate moving grid methods and therefore can be computationally expensive. Never- theless, the utility of such computations in terms of their ability to produce continuous data has been highlighted in the thesis. In that sense, these dynamic computations will be cheaper as compared to the static simulations to produce data at the same level of resolution.