Implementations of numerical simulations for solving systems of partial differential equations are often not verified and are falsely assumed to work correctly. As a result, the implementations are prone to coding errors that could degrade the accuracy of the solution. In order to ensure that a code is written correctly, rigorous verification of all parts of the code is necessary. Code verification is the task of ascertaining whether a numerical algorithm is solving the governing equations of the problem correctly. If an exact solution existed for the governing equations then verification would be easier but these solutions are rare because of the non-linearity of common Computational Fluid Dynamics (CFD) problems. In the absence of exact solutions, grid refinement studies are the most commonly used methods to verify codes using simulations on a sequence of grids but even these studies have limitations. The Method of Manufactured Solutions (MMS) is a novel and a recently developed technique that verifies the observed order-ofuracy of the implementation of a numerical algorithm. The method is more general and overcomes many of the limitations of the method of exact solutions and grid refinement studies. The central idea is to modify the governing equations and the boundary conditions by adding forcing functions or source terms in order to drive the discrete solution to a prescribed or ``manufactured' solution chosen a priori. A grid convergence study is performed subsequently to determine the observed orders. Two methods of accuracy assessment are presented here - solution accuracy analysis and residual error analysis. The method based on the error in the spatial residual is computationally less expensive and proved to be a valuable debugging tool. In the present work, the Method of Manufactured Solutions (MMS) is implemented on a compressible flow solver that solves the two-dimensional Euler equations on structured grids and an incompressible code that solves the two-dimensional Navier-Stokes equations on unstructured meshes. Exponential functions are used to ``manufacture' steady solutions to the governing equations. Solution and residual error analyses are presented. The influence of grid non-uniformity on the numerical accuracy is studied.

A detailed code verification study of an unstructured finite volume Computational Fluid Dynamics (CFD) code is performed. The Method of Manufactured Solutions is used to generate exact solutions for the Euler and Navier-Stokes equations to verify the correctness of the code through order of accuracy testing. The verification testing is performed on different mesh types which include triangular and quadrilateral elements in 2D and tetrahedral, prismatic, and hexahedral elements in 3D. The requirements of systematic mesh refinement are discussed, particularly in regards to unstructured meshes. Different code options verified include the baseline steady state governing equations, transport models, turbulence models, boundary conditions and unsteady flows. Coding mistakes, algorithm inconsistencies, and mesh quality sensitivities uncovered during the code verification are presented. In recent years, there has been significant work on the development of algorithms for the compressible Navier-Stokes equations on unstructured grids. One of the challenging tasks during the development of these algorithms is the formulation of consistent and accurate diffusion operators. The robustness and accuracy of diffusion operators depends on mesh quality. A survey of diffusion operators for compressible CFD solvers is conducted to understand different formulation procedures for diffusion fluxes. A patch-wise version of the Method of Manufactured Solutions is used to test the accuracy of selected diffusion operators. This testing of diffusion operators is limited to cell-centered finite volume methods which are formally second order accurate. These diffusion operators are tested and compared on different 2D mesh topologies to study the effect of mesh quality (stretching, aspect ratio, skewness, and curvature) on their numerical accuracy. Quantities examined include the numerical approximation errors and order of accuracy associated with face gradient reconstruction. From the analysis, defects in some of the numerical formulations are identified along with some robust and accurate diffusion operators. / Ph. D.

Code verification is the process of ensuring, to the degree possible, that there are no algorithm deficiencies and coding mistakes (bugs) in a scientific computing simulation. In this work, techniques are presented for performing code verification of boundary conditions commonly used in compressible and incompressible Computational Fluid Dynamics (CFD) codes. Using a compressible CFD code, this study assesses the subsonic inflow (isentropic and fixed-mass), subsonic outflow, supersonic outflow, no-slip wall (adiabatic and isothermal), and inviscid slip-wall. The use of simplified curved surfaces is proposed for easier generation of manufactured solutions during the verification of certain boundary conditions involving many constraints. To perform rigorous code verification, general grids with mixed cell types at the verified boundary are used. A novel approach is introduced to determine manufactured solutions for boundary condition verification when the velocity-field is constrained to be divergence-free during the simulation in an incompressible CFD code. Order of accuracy testing using the Method of Manufactured Solutions (MMS) is employed here for code verification of the major components of an open-source, multiphase flow code - MFIX. The presence of two-phase governing equations and a modified SIMPLE-based algorithm requiring divergence-free flows makes the selection of manufactured solutions more involved than for single-phase, compressible flows. Code verification is performed here on 2D and 3D, uniform and stretched meshes for incompressible, steady and unsteady, single-phase and two-phase flows using the two-fluid model of MFIX. In a CFD simulation, truncation error (TE) is the difference between the continuous governing equation and its discrete approximation. Since TE can be shown to be the local source term for the discretization error, TE is proposed as the criterion for determining which regions of the computational mesh should be refined/coarsened. For mesh modification, an error equidistribution strategy to perform r-refinement (i.e., mesh node relocation) is employed. This technique is applied to 1D and 2D inviscid flow problems where the exact (i.e., analytic) solution is available. For mesh adaptation based upon TE, about an order of magnitude improvement in discretization error levels is observed when compared with the uniform mesh. / Ph. D.

This work examines the Method of Nearby Problems as a way to generate analytical exact solutions to problems governed by partial differential equations (PDEs). The method involves generating a numerical solution to the original problem of interest, curve fitting the solution, and generating source terms by operating the governing PDEs upon the curve fit. Adding these source terms to the right-hand-side of the governing PDEs defines the nearby problem. In addition to its use for generating exact solutions the MNP can be extended for use as an error estimator. The nearby problem can be solved numerically on the same grid as the original problem. The nearby problem discretization error is calculated as the difference between its numerical solution and exact solution (curve fit). This is an estimate of the discretization error in the original problem of interest. The accuracy of the curve fits is quite important to this work. A method of curve fitting that takes local least squares fits and combines them together with weighting functions is used. This results in a piecewise fit with continuity at interface boundaries. A one-dimensional Burgers' equation case shows this to be a better approach then global curve fits. Six two-dimensional cases are investigated including solutions to the time-varying Burgers' equation and to the 2D steady Euler equations. The results show that the Method of Nearby Problems can be used to create realistic, analytical exact solutions to problems governed by PDEs. The resulting discretization error estimates are also shown to be reasonable for several cases examined. / Master of Science

No description available.

Thesis (PhD)--Stellenbosch University, 2013. / ENGLISH ABSTRACT: In this work the Method of Manufactured Solutions (MMS) is introduced for the code veri cation of full-wave frequency dependent electromagnetic computational software. At rst the method is sketched in the context of the veri cation and validation process and the need for proper code veri cation is highlighted. Subsequently, the MMS is investigated in its natural context: the Finite Element Method, speci cally for the E- eld Vector Wave Equation. The usefulness of the method to detect error in a computational code is demonstrated. The selection of Manufactured Solutions is discussed and it is demonstrated how it can be used to nd the probable cause of bugs. Mutation testing is introduced and used to show the ability to detect errors present in code. The MMS is nally applied in a novel manner to a Method of Moments (MoM) code. The challenges of numerical integration associated with the application of the operator is discussed and correct integration is successfully demonstrated. Subsequently the MMS is demonstrated to be successfully applied to the MoM and mutation testing is used to demonstrate the practical e cacy of the method. The application of the MMS to the MoM is the main contribution of this work. / AFRIKAANSE OPSOMMING: Die Metode van Vervaardigde Oplossings (MVO) word hier bekend gestel vir die veri kasie van numeriese volgolf frekwensie-afhanklike elektromagnetise kode. Die metode word eerstens in die bre e konteks van algemene veri kasie en validasie geplaas en gevolglik word die noodsaaklikheid van kode veri kasie beklemtoon. Daarna, word die toets-metode in die konteks van die Eindige Element Metode vir die E-veld vektorgolf vergelyking bestudeer. Die MVO is oorspronklik ontwikkel in die di erentiaalvergelyking omgewing. Die bruikbaarheid van die metode vir elektromagnetiese simulasies word prakties gedemonstreer deur die opsporing van werklike foute. Die metode word ook verder ondersoek vir die oorsprong van foute. Mutasietoetsing word bekendgestel en word gebruik om die metode verder prakties te veri eer. Die MVO word laastens in 'n nuwe manier gebruik om 'n Moment Metode kode te veri eer. Die praktiese probleme betrokke by numeriese integrasie word ondersoek en die korrekte toepassing van die integraal operator word prakties gedemonstreer. Daarna, word die MVO in hierdie konteks gedemonstreer deur verskeie voorbeelde te ondersoek. Mutasietoetsing word weereens gebruik om na die e ektiewiteit van die MVO te kyk om 'n Moment Metode kode te toets. Die toepassing van die MVO op 'n Moment Metode kode is die hoof bydrae van hierdie werk.

O som gerado por escoamentos, também conhecido como aeroacústica, tem se tornado cada vez mais importante em áreas industriais diversas desde aviação comercial até aparelhos eletrodomésticos, afetando diretamente os requisitos necessários para o desenvolvimento de novos produtos. Um caso particular é o ruído gerado por válvulas de compressores herméticos de refrigeração, sendo o compressor a principal fonte de ruído em refrigeradores domésticos. O presente trabalho tem por objetivo iniciar o desenvolvimento de uma ferramenta confiável de simulação capaz de auxiliar engenheiros na predição de problemas de aeroacústica, especialmente um que possa no futuro ser utilizado para estudar o ruído gerado pelo escoamento em válvulas de compressores herméticos. Para isso, foi desenvolvido um código para simulação numérica direta de aeroacústica. Utilizou-se processamento paralelo com decomposição de domínio para usar Simulação Numérica Direta em um tempo factível; esquemas de discretização espaciais e temporais de alta ordem para minimizar ao máximo os fenômenos de dissipação e dispersão do escoamento e das ondas acústicas e uma série de tratamentos no domínio como filtragem e estiramento da malha como também condições de contorno características com o intuito de obter uma solução adequada para estudo de aeroacústica. Assim, são apresentadas todas as etapas desenvolvidas no equacionamento, implementação e verificação. A verificação foi realizada segundo um processo matemático formal (Método das Soluções Manufaturadas) com o qual obteve-se que a ordem de precisão dos cálculos era a mesma da ordem formal dos esquemas de discretização utilizados para todas as variáveis. Também obteve-se a mesma concordância para análise do divergente da velocidade, verificando o código para simulação numérica direta de aeroacústica. Posteriormente, foram realizadas simulações de escoamentos compressíveis cisalhantes e seus resultados comparados com dados apresentados em literatura. Também foram calculadas as taxas de amplificação de perturbações e comparadas com a Teoria de Estabilidade Linear. Novamente, foram obtidos resultados satisfatórios nessas etapas, mostrando que a implementação do código DNS está verificada. / Sound generated by flow, known as aeroacoustics, is becoming more important in several industrial areas from commercial aircraft to household appliances, affecting directly the requirements to the development of new products. A particular case is the noise generated by valves of refrigeration hermetic compressors, being the compressor the main noise source in household refrigerators. This work has the goal of initiate the development of a reliable tool able to help engineers to predict aeroacoustics problems, specially one that can be used in the future to study the noise generated by the flow in valves of hermetic compressors. To do so, it was developed a numerical code to perform direct numerical simulation of aeroacoustics. It was used parallel processing with domain decomposition to use Direct Numerical Simulation in a feasible time; high order temporal and spatial discretization schemes to minimize the most the dispersion and dissipation phenomena of the flow field and of the acoustics waves and a series of treatments in the domain as filtering and mesh stretching as well as characteristics boundary conditions aiming a proper solution to study aeroacoustics. Thus, here all the steps developed in the formulation, implementation and verification are presented . The verification was done according to a formal mathematical procedure (Method of Manufactured Solutions) with which was found that the precision order of the calculations was the same of the formal order of the used discretization schemes for all variables. The same agreement was also obtained to the analysis of the divergence of the velocity, verifying the code to direct numerical simulation of aeroacoustics. Posteriorly, it were simulated shear compressible flows and the results were compared to literature data. Also, it were calculated the amplification rates of the disturbances and compared to Linear Stability Theory. Once more, it was obtained satisfactory results in these steps, showing that the implementation of the DNS code is verified.

A presente tese se refere a três principais objetivos. Um objetivo foi desenvolver um código de simulação numérica direta para simulação de ondas de instabilidade em um escoamento de Poiseuille plano. O outro objetivo foi analisá-Io através do Método das Soluções Manufaturadas (MMS), e por fim, um terceiro objetivo foi estudar o regime não-linear da evolução de trens de ondas modulados em um escoamento de Poiseuille plano. O código resolve numericamente, com diferenças finitas de ordem de precisão alta e métodos pseudo espectrais, as equações de Navier-Stokes tri-dimensionais e incompressíveis numa formulação vorticidade-velocidade. O MMS é um método de verificação de código mais completo que os normalmente usados, por exemplo, comparação com teoria de estabilidade linear. O código usa diferenças finitas de ordem alta de precisão, mas com diferentes ordens em diferentes regiões do domínio. O MMS é pouco utilizado neste tipo de código. Concluiu-se que estes códigos em geral não operam na chamada faixa assintótica de erro. Na faixa de trabalho, a ordem do erro varia no domínio computacional de forma consistente com os métodos numéricos empregados. Isto permite, entre outras coisas, a otimização do esquema numérico. Após os testes de verifição, simulações numéricas dos trens de ondas foram realizadas. A análise dos resultados foi feita através das teorias de instabilidade primária e secundária e teoria fracamente não-linear. Foi estudado um trem de ondas modulado numa região do diagrama de instabilidade onde vários estudos para ondas mono cromáticas foram realizados. Os resultados sugeriram que nesta região o regime não-linear de transição dos trens de ondas modulados é governado pela instabilidade tipo-K. Com a redução da amplitude inicial de perturbação um cenário mais complexo, que pode estar envolvendo outros mecanismos, foi observado. Casos mais próximos ao primeiro ramo do diagrama também foram estudados. Nesta região a teoria linear prevê ondas tri-dimensionais sendo as mais instáveis, fato que poderia favorecer a ocorrência da chamada transição oblíqua. Os resultados indicaram que o trem de ondas modulado se divide em duas regiões que tendem-se afastar uma da outra. Apesar de estar associado com a instabilidade linear, este comportamento não havia sido antecipado na literatura. Neste cenário, para trens de ondas modulados isolados, o regime não-linear não pode ser correlacionado claramente com nenhum dos cenários clássicos de transição. Neste contexto transição oblíqua pode estar restrito a situações que envolvem a interação entre trens de ondas modulados. / The current thesis had three objectives. The first objective was to develop a code of direct numerical simulation (DNS) to simulation of waves of instability in a plane Poiseuille flow. The other main objective was to analyze it through the method of manufactured solutions (MMS). Finally, a third objective was to study the nonlinear regime of spanwise modulated wavetrains in a plane Poiseuille flow. Using high-order finite differences and pseudo-spectral methods, the DNS code solved the incompressible three-dimensional Navier-Stokes equations in a vorticity-velocity formulation. The MMS is a verification method more complete than the others more often used, for example, comparison with linear stability theory. The code used different high-order finite differences in different regions of the domain. The MMS has been little used for this type of code. It was concluded that these codes generally do not operate in the so-called errar asymptotic range. In the working range, the observed arder changes in the computational domain in a manner consistent with the numerical methods employed. This allows, among other things, optimization of the numerical scheme. After testing, numerical simulations of the wavetrains were performed. The analysis of the results was made based on the primary and secondary instability theories and weakly non-linear theory. A modulated wavetrain was studied in a region of the stability diagram where several studies for monochromatic waves were performed. The results suggested that in this region the non-linear regime of transition of the modulation wavetrain is governed by the K-type instability. With the reduction of the magnitude of the initial disturbance, a more complex scenario, which may involve other mechanisms, was observed. Cases near the first branch of the instability diagram were also studied. In this region, the linear theory predicts tri-dimensional waves are the most unstable, a fact that could lead to the so-called oblique transition. The results indicated that the modulated wavetrain divided into two regions, which tend move further from each other. Despite being associated with the linear instability, this behavior was not anticipated in the literature. In this scenario, for isolated modulated wavetrains, the non-linear system could not be clearly correlated with any of the classic scenarios of transition. In this context, oblique transition may perhaps to situations involving the interaction between modulated wavetrains.

A presente tese se refere a três principais objetivos. Um objetivo foi desenvolver um código de simulação numérica direta para simulação de ondas de instabilidade em um escoamento de Poiseuille plano. O outro objetivo foi analisá-Io através do Método das Soluções Manufaturadas (MMS), e por fim, um terceiro objetivo foi estudar o regime não-linear da evolução de trens de ondas modulados em um escoamento de Poiseuille plano. O código resolve numericamente, com diferenças finitas de ordem de precisão alta e métodos pseudo espectrais, as equações de Navier-Stokes tri-dimensionais e incompressíveis numa formulação vorticidade-velocidade. O MMS é um método de verificação de código mais completo que os normalmente usados, por exemplo, comparação com teoria de estabilidade linear. O código usa diferenças finitas de ordem alta de precisão, mas com diferentes ordens em diferentes regiões do domínio. O MMS é pouco utilizado neste tipo de código. Concluiu-se que estes códigos em geral não operam na chamada faixa assintótica de erro. Na faixa de trabalho, a ordem do erro varia no domínio computacional de forma consistente com os métodos numéricos empregados. Isto permite, entre outras coisas, a otimização do esquema numérico. Após os testes de verifição, simulações numéricas dos trens de ondas foram realizadas. A análise dos resultados foi feita através das teorias de instabilidade primária e secundária e teoria fracamente não-linear. Foi estudado um trem de ondas modulado numa região do diagrama de instabilidade onde vários estudos para ondas mono cromáticas foram realizados. Os resultados sugeriram que nesta região o regime não-linear de transição dos trens de ondas modulados é governado pela instabilidade tipo-K. Com a redução da amplitude inicial de perturbação um cenário mais complexo, que pode estar envolvendo outros mecanismos, foi observado. Casos mais próximos ao primeiro ramo do diagrama também foram estudados. Nesta região a teoria linear prevê ondas tri-dimensionais sendo as mais instáveis, fato que poderia favorecer a ocorrência da chamada transição oblíqua. Os resultados indicaram que o trem de ondas modulado se divide em duas regiões que tendem-se afastar uma da outra. Apesar de estar associado com a instabilidade linear, este comportamento não havia sido antecipado na literatura. Neste cenário, para trens de ondas modulados isolados, o regime não-linear não pode ser correlacionado claramente com nenhum dos cenários clássicos de transição. Neste contexto transição oblíqua pode estar restrito a situações que envolvem a interação entre trens de ondas modulados. / The current thesis had three objectives. The first objective was to develop a code of direct numerical simulation (DNS) to simulation of waves of instability in a plane Poiseuille flow. The other main objective was to analyze it through the method of manufactured solutions (MMS). Finally, a third objective was to study the nonlinear regime of spanwise modulated wavetrains in a plane Poiseuille flow. Using high-order finite differences and pseudo-spectral methods, the DNS code solved the incompressible three-dimensional Navier-Stokes equations in a vorticity-velocity formulation. The MMS is a verification method more complete than the others more often used, for example, comparison with linear stability theory. The code used different high-order finite differences in different regions of the domain. The MMS has been little used for this type of code. It was concluded that these codes generally do not operate in the so-called errar asymptotic range. In the working range, the observed arder changes in the computational domain in a manner consistent with the numerical methods employed. This allows, among other things, optimization of the numerical scheme. After testing, numerical simulations of the wavetrains were performed. The analysis of the results was made based on the primary and secondary instability theories and weakly non-linear theory. A modulated wavetrain was studied in a region of the stability diagram where several studies for monochromatic waves were performed. The results suggested that in this region the non-linear regime of transition of the modulation wavetrain is governed by the K-type instability. With the reduction of the magnitude of the initial disturbance, a more complex scenario, which may involve other mechanisms, was observed. Cases near the first branch of the instability diagram were also studied. In this region, the linear theory predicts tri-dimensional waves are the most unstable, a fact that could lead to the so-called oblique transition. The results indicated that the modulated wavetrain divided into two regions, which tend move further from each other. Despite being associated with the linear instability, this behavior was not anticipated in the literature. In this scenario, for isolated modulated wavetrains, the non-linear system could not be clearly correlated with any of the classic scenarios of transition. In this context, oblique transition may perhaps to situations involving the interaction between modulated wavetrains.

O som gerado por escoamentos, também conhecido como aeroacústica, tem se tornado cada vez mais importante em áreas industriais diversas desde aviação comercial até aparelhos eletrodomésticos, afetando diretamente os requisitos necessários para o desenvolvimento de novos produtos. Um caso particular é o ruído gerado por válvulas de compressores herméticos de refrigeração, sendo o compressor a principal fonte de ruído em refrigeradores domésticos. O presente trabalho tem por objetivo iniciar o desenvolvimento de uma ferramenta confiável de simulação capaz de auxiliar engenheiros na predição de problemas de aeroacústica, especialmente um que possa no futuro ser utilizado para estudar o ruído gerado pelo escoamento em válvulas de compressores herméticos. Para isso, foi desenvolvido um código para simulação numérica direta de aeroacústica. Utilizou-se processamento paralelo com decomposição de domínio para usar Simulação Numérica Direta em um tempo factível; esquemas de discretização espaciais e temporais de alta ordem para minimizar ao máximo os fenômenos de dissipação e dispersão do escoamento e das ondas acústicas e uma série de tratamentos no domínio como filtragem e estiramento da malha como também condições de contorno características com o intuito de obter uma solução adequada para estudo de aeroacústica. Assim, são apresentadas todas as etapas desenvolvidas no equacionamento, implementação e verificação. A verificação foi realizada segundo um processo matemático formal (Método das Soluções Manufaturadas) com o qual obteve-se que a ordem de precisão dos cálculos era a mesma da ordem formal dos esquemas de discretização utilizados para todas as variáveis. Também obteve-se a mesma concordância para análise do divergente da velocidade, verificando o código para simulação numérica direta de aeroacústica. Posteriormente, foram realizadas simulações de escoamentos compressíveis cisalhantes e seus resultados comparados com dados apresentados em literatura. Também foram calculadas as taxas de amplificação de perturbações e comparadas com a Teoria de Estabilidade Linear. Novamente, foram obtidos resultados satisfatórios nessas etapas, mostrando que a implementação do código DNS está verificada. / Sound generated by flow, known as aeroacoustics, is becoming more important in several industrial areas from commercial aircraft to household appliances, affecting directly the requirements to the development of new products. A particular case is the noise generated by valves of refrigeration hermetic compressors, being the compressor the main noise source in household refrigerators. This work has the goal of initiate the development of a reliable tool able to help engineers to predict aeroacoustics problems, specially one that can be used in the future to study the noise generated by the flow in valves of hermetic compressors. To do so, it was developed a numerical code to perform direct numerical simulation of aeroacoustics. It was used parallel processing with domain decomposition to use Direct Numerical Simulation in a feasible time; high order temporal and spatial discretization schemes to minimize the most the dispersion and dissipation phenomena of the flow field and of the acoustics waves and a series of treatments in the domain as filtering and mesh stretching as well as characteristics boundary conditions aiming a proper solution to study aeroacoustics. Thus, here all the steps developed in the formulation, implementation and verification are presented . The verification was done according to a formal mathematical procedure (Method of Manufactured Solutions) with which was found that the precision order of the calculations was the same of the formal order of the used discretization schemes for all variables. The same agreement was also obtained to the analysis of the divergence of the velocity, verifying the code to direct numerical simulation of aeroacoustics. Posteriorly, it were simulated shear compressible flows and the results were compared to literature data. Also, it were calculated the amplification rates of the disturbances and compared to Linear Stability Theory. Once more, it was obtained satisfactory results in these steps, showing that the implementation of the DNS code is verified.