Neste trabalho apresenta-se um novo esquema prático tipo upwind de alta resolução, denominado EPUS (Eight-degree Polynomial Upwind Scheme), para resolver numericamente equações de conservação TVD e é implementado no contexto do método das diferenças finitas. O desempenho do esquema é investigado na resolução de sistemas hiperbólicos de leis de conservação e escoamentos incompressíveis complexos com superfícies livres. Os resultados numéricos mostraram boa concordãncia com outros resultados numéricos e dados experimentais existentes / Is this work a new practical high resolution upwinding scheme, called EPUS (Eight-degree Polynomial Upwind Scheme), for the numerical solution of transient convection-dominated conservation equations is present. The scheme is based on TVD stability criterion and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving hyperbolic systems of conservation laws and complex incompressible flows with free surfaces. The numerical results displayed good agreement with other existing numerical and experimental data

Neste trabalho apresenta-se um novo esquema prático tipo upwind de alta resolução, denominado EPUS (Eight-degree Polynomial Upwind Scheme), para resolver numericamente equações de conservação TVD e é implementado no contexto do método das diferenças finitas. O desempenho do esquema é investigado na resolução de sistemas hiperbólicos de leis de conservação e escoamentos incompressíveis complexos com superfícies livres. Os resultados numéricos mostraram boa concordãncia com outros resultados numéricos e dados experimentais existentes / Is this work a new practical high resolution upwinding scheme, called EPUS (Eight-degree Polynomial Upwind Scheme), for the numerical solution of transient convection-dominated conservation equations is present. The scheme is based on TVD stability criterion and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving hyperbolic systems of conservation laws and complex incompressible flows with free surfaces. The numerical results displayed good agreement with other existing numerical and experimental data

This thesis studies advection modeling for heat transfer and fluid flow problems using a new Non--Inverted Skew Upwind Scheme (called NISUS). Variants of the new scheme are formulated and developed with 8-noded hexahedral elements using the Finite Element Method (FEM)and rectangular elements based on a Finite Volume Method (FVM). A new method of mass weighting to predict convective fluxes of each scalar from the nodal point values is developed. Due to an explicit representation in terms of nodal variables, local inversion of the upwind coefficient matrix is not needed. Also, this thesis evaluates two variants of the new scheme (i.e., 3-node / 3-point and 4-node / 8-point formulations) within a 3--D FEM and a third variant within a 2--D FVM. The 3--D FEM variants are applied to a variety of test problems involving the transport of a scalar variable, while the 2--D FVM variant is applied to fluid flow problems including natural convection in an enclosure and micro--channel flow simulations. The promising performance of NISUS, as compared with exact and previous solutions, is demonstrated both in terms of accuracy and stability. Furthermore, a new data storage format called Compressed Banded Data (CBD) is developed for sparse banded matrices generated by the control volume finite element method (CVFEM). The platform of the new CBD structure permits dynamic switching between various solvers, without any procedural change in the implementation of existing simulation software. The performance of different Krylov techniques with an ILU(0) preconditioner is observed and compared in three test problems with a direct solver. / October 2006

This thesis studies advection modeling for heat transfer and fluid flow problems using a new Non--Inverted Skew Upwind Scheme (called NISUS). Variants of the new scheme are formulated and developed with 8-noded hexahedral elements using the Finite Element Method (FEM)and rectangular elements based on a Finite Volume Method (FVM). A new method of mass weighting to predict convective fluxes of each scalar from the nodal point values is developed. Due to an explicit representation in terms of nodal variables, local inversion of the upwind coefficient matrix is not needed. Also, this thesis evaluates two variants of the new scheme (i.e., 3-node / 3-point and 4-node / 8-point formulations) within a 3--D FEM and a third variant within a 2--D FVM. The 3--D FEM variants are applied to a variety of test problems involving the transport of a scalar variable, while the 2--D FVM variant is applied to fluid flow problems including natural convection in an enclosure and micro--channel flow simulations. The promising performance of NISUS, as compared with exact and previous solutions, is demonstrated both in terms of accuracy and stability. Furthermore, a new data storage format called Compressed Banded Data (CBD) is developed for sparse banded matrices generated by the control volume finite element method (CVFEM). The platform of the new CBD structure permits dynamic switching between various solvers, without any procedural change in the implementation of existing simulation software. The performance of different Krylov techniques with an ILU(0) preconditioner is observed and compared in three test problems with a direct solver.

This thesis studies advection modeling for heat transfer and fluid flow problems using a new Non--Inverted Skew Upwind Scheme (called NISUS). Variants of the new scheme are formulated and developed with 8-noded hexahedral elements using the Finite Element Method (FEM)and rectangular elements based on a Finite Volume Method (FVM). A new method of mass weighting to predict convective fluxes of each scalar from the nodal point values is developed. Due to an explicit representation in terms of nodal variables, local inversion of the upwind coefficient matrix is not needed. Also, this thesis evaluates two variants of the new scheme (i.e., 3-node / 3-point and 4-node / 8-point formulations) within a 3--D FEM and a third variant within a 2--D FVM. The 3--D FEM variants are applied to a variety of test problems involving the transport of a scalar variable, while the 2--D FVM variant is applied to fluid flow problems including natural convection in an enclosure and micro--channel flow simulations. The promising performance of NISUS, as compared with exact and previous solutions, is demonstrated both in terms of accuracy and stability. Furthermore, a new data storage format called Compressed Banded Data (CBD) is developed for sparse banded matrices generated by the control volume finite element method (CVFEM). The platform of the new CBD structure permits dynamic switching between various solvers, without any procedural change in the implementation of existing simulation software. The performance of different Krylov techniques with an ILU(0) preconditioner is observed and compared in three test problems with a direct solver.

Esta dissertação trata da solução numérica de problemas em dinâmica dos fluidos usando dois novos esquemas upwind de alta resolução, denominados FDPUS-C1 (Five-Degree Polynomial Upwind Scheme of \' C POT. 1\' Class) e SDPUS-C1 (Six-Degree Polynomial Upwind Scheme of \'C POT.1\' Class), para a discretização de termos convectivos lineares e não-lineares. Os esquemas são baseados nos critérios de estabilidade TVD (Total Variation Diminishing) e CBC (Convection Boundedness Criterion) e são implementados, nos contextos das metodologias de diferenças finitas e volumes finitos, no ambiente de simulação Freeflow (an integrated simulation system for Free surface Flow) para escoamentos imcompressíveis 2D, 2D-1/2 e 3D, ou no código bem conhecido CLAWPACK ( Conservation LAW PACKage) para problemaw compressíveis 1D e 2D. Vários testes computacionais são feitos com o objetivo de verificar e validar os métodos numéricos contra esquemas upwind populares. Os novos esqumas são então aplicados na resolução de uma gama ampla de problemas em CFD (Computational Fluids Dynamics), tais como propagação de ondas de choque e escoamentos incompressíveis envolvendo superfícies livres móveis. Em particular, os resultados numéricos para leis de conservação hiperbólicas 2D e equações de Navier-Stokes incompressíveis 2D, 2D-1/2 e 3D demosntram que esses novos esquemas convectivos tipo upwind polinomiais funcionam muito bem / This dissertation deals with the numerical solution of fluid dynamics problems using two new high resolution upwind schemes,. namely FDPUS-C1 and SDPUS-C1, for the discretization of the linear and non-linear convection terms. The Schemes are based on TVD and DBC stability criteria and are implemented in the context of the finite difference and finite volume methodologies, either into the Freeflow code for 2D, 2D-1/2 and 3D incompressible flows or in the well-known CLAWPACK code for 1D and 2D compressible flows. Several computational tests are performed to verify and validate the numerical methods against other popularly used upwind schemes. The new schemes are then applied to solve a wide range of problems in CFD, such as shock wave propagation and incompressible fluid flows involving moving free msurfaces. In particular, the numerical results for 2D hyperbolic conservation laws and 2D, 2D-1/2 and 3D incompressible Navier-Stokes eqautions show that new polynomial upwind convection schemes perform very well

In the quest for a high fidelity numerical scheme for CFD it is necessary to satisfy demands on accuracy, conservation, positivity and upwinding. Recently the requirement of rotational invariance has been added to this list. In the present work we are mainly interested in upwinding and rotational invariance of Least Squares Kinetic Upwind Method (LSKUM). The standard LSKUM achieves upwinding by stencil division along co-ordinate axes which is referred to as co-ordinate splitting method. This leads to symmetry breaking and rotational invariance is lost. Thus the numerical solution becomes co-ordinate frame dependent. To overcome this undesirable feature of existing numerical schemes, a new algorithm called KUMARI (Kinetic Upwind Method Avec Rotational Invariance, 'Avec' in French means 'with') has been developed. The interesting mathematical relation between directional derivative, Fourier series and divergence operator has been used effectively to achieve upwinding as well as rotational invariance and hence making the scheme truly or genuinely multidimensional upwind scheme. The KUMARI has been applied to the test case of standard 2D shock reflection problem, flow past airfoils, then to 2D blast wave problem and lastly to 2D Riemann problem (Lax's 3rd test case). The results show that either KUMARI is comparable to or in some cases better than the usual LSKUM.

Esta dissertação trata da solução numérica de problemas em dinâmica dos fluidos usando dois novos esquemas upwind de alta resolução, denominados FDPUS-C1 (Five-Degree Polynomial Upwind Scheme of \' C POT. 1\' Class) e SDPUS-C1 (Six-Degree Polynomial Upwind Scheme of \'C POT.1\' Class), para a discretização de termos convectivos lineares e não-lineares. Os esquemas são baseados nos critérios de estabilidade TVD (Total Variation Diminishing) e CBC (Convection Boundedness Criterion) e são implementados, nos contextos das metodologias de diferenças finitas e volumes finitos, no ambiente de simulação Freeflow (an integrated simulation system for Free surface Flow) para escoamentos imcompressíveis 2D, 2D-1/2 e 3D, ou no código bem conhecido CLAWPACK ( Conservation LAW PACKage) para problemaw compressíveis 1D e 2D. Vários testes computacionais são feitos com o objetivo de verificar e validar os métodos numéricos contra esquemas upwind populares. Os novos esqumas são então aplicados na resolução de uma gama ampla de problemas em CFD (Computational Fluids Dynamics), tais como propagação de ondas de choque e escoamentos incompressíveis envolvendo superfícies livres móveis. Em particular, os resultados numéricos para leis de conservação hiperbólicas 2D e equações de Navier-Stokes incompressíveis 2D, 2D-1/2 e 3D demosntram que esses novos esquemas convectivos tipo upwind polinomiais funcionam muito bem / This dissertation deals with the numerical solution of fluid dynamics problems using two new high resolution upwind schemes,. namely FDPUS-C1 and SDPUS-C1, for the discretization of the linear and non-linear convection terms. The Schemes are based on TVD and DBC stability criteria and are implemented in the context of the finite difference and finite volume methodologies, either into the Freeflow code for 2D, 2D-1/2 and 3D incompressible flows or in the well-known CLAWPACK code for 1D and 2D compressible flows. Several computational tests are performed to verify and validate the numerical methods against other popularly used upwind schemes. The new schemes are then applied to solve a wide range of problems in CFD, such as shock wave propagation and incompressible fluid flows involving moving free msurfaces. In particular, the numerical results for 2D hyperbolic conservation laws and 2D, 2D-1/2 and 3D incompressible Navier-Stokes eqautions show that new polynomial upwind convection schemes perform very well