A hybrid Cartesian grid and gridless method for compressible flows

A hybrid Cartesian grid and gridless method is presented to compute unsteady compressible flows for complex geometries. In this method, a Cartesian grid is used as baseline mesh to cover the computational domain, while the boundary surfaces are addressed using a gridless method. This hybrid method combines the efficiency of a Cartesian grid method and the flexibility of a gridless method for the complex geometries. The developed method is used to compute a number of test cases to validate the accuracy and efficiency of the method. The numerical results obtained indicate that the use of this hybrid method leads to a significant improvement in performance over its unstructured grid counterpart for the time-accurate solution of the compressible Euler equations. An overall speed-up factor of about eight and a saving in storage requirements about one order of magnitude for a typical three-dimensional problem in comparison with the unstructured grid method are demonstrated

A hybrid building‐block and gridless method for compressible flows

A hybrid building‐block Cartesian grid and gridless method is presented to compute unsteady compressible flows for complex geometries. In this method, a Cartesian mesh based on a building‐block grid is used as a baseline mesh to cover the computational domain, while the boundary surfaces are represented using a set of gridless points. This hybrid method combines the efficiency of a Cartesian grid method and the flexibility of a gridless method for the complex geometries. The developed method is used to compute a number of test cases to validate the accuracy and efficiency of the method. The numerical results obtained indicate that the use of this hybrid method leads to a significant improvement in performance over its unstructured grid counterpart for the time‐accurate solution of the compressible Euler equations. An overall speed‐up factor from six to more than one order of magnitude and a saving in storage requirements up to one order of magnitude for all test cases in comparison with the unstructured grid method are demonstrated

A parallel meshless dynamic cloud method on graphic processing units for unsteady compressible flows past moving boundaries

This paper presents an effort to implement a recently proposed meshless dynamic cloud method (Hong Wang et al., 2010) on modern high-performance graphic processing units (GPUs) with the compute unified device architecture (CUDA) programming model. Within the framework of the meshless method, clouds of points used as basic computational stencils are distributed in the whole flow domain. The spatial derivatives of the governing equations are discretised by the moving-least square scheme on every cloud of points. Roe’s approximate Riemann solver is adopted to compute the convective flux. A dual-time stepping approach, which iterates in physical and pseudo temporal spaces, is employed to obtain the time-accurate solution. Simulation of steady compressible flows over a fixed aerofoil is firstly carried out to verify the GPU implementation of the method. Then it is extended to compute unsteady flows past oscillatory aerofoils. Numerical outcomes are compared with experimental and/or other reference results to validate the method. Significant performance speedup of more than an order of magnitude is verified by the numerical results. Systematic analysis shows that GPU is more energy efficient than CPU for solving aerodynamic problems. This demonstrates the potential of the proposed method to solve fluid–structure interaction problems

Unstructured Grid Generation Techniques and Software

The Workshop on Unstructured Grid Generation Techniques and Software was conducted for NASA to assess its unstructured grid activities, improve the coordination among NASA centers, and promote technology transfer to industry. The proceedings represent contributions from Ames, Langley, and Lewis Research Centers, and the Johnson and Marshall Space Flight Centers. This report is a compilation of the presentations made at the workshop

ジユウ ヒョウメンリュウ カイセキ ノ タメ ノ アタラシイ ビブン エンザンシ モデル ニ ヨル カイリョウガタ リュウシホウ

京都大学0048新制・課程博士博士(工学)甲第14147号工博第2981号新制||工||1442(附属図書館)26453UT51-2008-N464京都大学大学院工学研究科都市環境工学専攻(主査)教授 後藤 仁志, 教授 細田 尚, 准教授 牛島 省学位規則第4条第1項該当Doctor of EngineeringKyoto UniversityDFA

A Multi-Solver Scheme for Viscous Flows Using Adaptive Cartesian Grids and Meshless Grid Communication

This work concerns the development of an adaptive multi-solver approach for CFD simulation of viscous flows. Curvilinear grids are used near solid bodies to capture boundary layers, and stuctured adaptive Cartesian grids are used away from the body to fill the majority of the computational domain. An edge-based meshless scheme is used in the interface region to connnect the near-body and off-body codes. We show that the combination of a body-fitted grid near the surface coupled with an adaptive Cartesian grid system away from the surface leads to a highly efficient scheme with sharp feature resolution. The use of a meshless flow solver to interface the body-fitted and Cartesian grid systems leads to seamless grid communication without many of the complexities inherent in traditional Chimera overset grid interpolation schemes. The hierarchical structure of the nested Cartesian grids may be exploited to achieve multigrid convergence for steady problems and for use in dual-time stepping algorithms for unsteady problems. Results of two-dimensional steady airfoil calculations are presented. I

Meshless Direct Numerical Simulation of Turbulent Incompressible Flows

A meshless direct pressure-velocity coupling procedure is presented to perform Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of turbulent incompressible flows in regular and irregular geometries. The proposed method is a combination of several efficient techniques found in different Computational Fluid Dynamic (CFD) procedures and it is a major improvement of the algorithm published in 2007 by this author. This new procedure has very low numerical diffusion and some preliminary calculations with 2D steady state flows show that viscous effects become negligible faster that ever predicted numerically. The fundamental idea of this proposal lays on several important inconsistencies found in three of the most popular techniques used in CFD, segregated procedures, streamline-vorticity formulation for 2D viscous flows and the fractional-step method, very popular in DNS/LES. The inconsistencies found become important in elliptic flows and they might lead to some wrong solutions if coarse grids are used. In all methods studied, the mathematical basement was found to be correct in most cases, but inconsistencies were found when writing the boundary conditions. In all methods analyzed, it was found that it is basically impossible to satisfy the exact set of boundary conditions and all formulations use a reduced set, valid for parabolic flows only. For example, for segregated methods, boundary condition of normal derivative for pressure zero is valid only in parabolic flows. Additionally, the complete proposal for mass balance correction is right exclusively for parabolic flows. In the streamline-vorticity formulation, the boundary conditions normally used for the streamline function, violates the no-slip condition for viscous flow. Finally, in the fractional-step method, the boundary condition for pseudo-velocity implies a zero normal derivative for pressure in the wall (correct in parabolic flows only) and, when the flows reaches steady state, the procedure does not guarantee mass balance. The proposed procedure is validated in two cases of 2D flow in steady state, backward-facing step and lid-driven cavity. Comparisons are performed with experiments and excellent agreement was obtained in the solutions that were free from numerical instabilities. A study on grid usage is done. It was found that if the discretized equations are written in terms of a local Reynolds number, a strong criterion can be developed to determine, in advance, the grid requirements for any fluid flow calculation. The 2D-DNS on parallel plates is presented to study the basic features present in the simulation of any turbulent flow. Calculations were performed on a short geometry, using a uniform and very fine grid to avoid any numerical instability. Inflow conditions were white noise and high frequency oscillations. Results suggest that, if no numerical instability is present, inflow conditions alone are not enough to sustain permanently the turbulent regime. Finally, the 2D-DNS on a backward-facing step is studied. Expansion ratios of 1.14 and 1.40 are used and calculations are performed in the transitional regime. Inflow conditions were white noise and high frequency oscillations. In general, good agreement is found on most variables when comparing with experimental data

The study on adaptive Cartesian grid methods for compressible flow and their applications

This research is mainly focused on the development of the adaptive Cartesian grid methods for compressibl  e flow. At first, the ghost cell method and its applications for inviscid compressible flow on adaptive tree Cartesian grid are developed. The proposed method is successfully used to evaluate various inviscid compressible flows around complex bodies. The mass conservation of the method is also studied by numerical analysis. The extension to three-dimensional flow is presented. Then, an h-adaptive Runge–Kutta discontinuous Galerkin (RKDG) method is presented in detail for the development of high accuracy numerical method under Cartesian grid. This method combined with the ghost cell immersed boundary method is also validated by well documented test problems involving both steady and unsteady compressible flows over complex bodies in a wide range of Mach numbers. In addition, in order to suppress the failure of preserving positivity of density or pressure, which may cause blow-ups of the high order numerical algorithms, a positivity-preserving limiter technique coupled with h-adaptive RKDG method is developed. Such a method has been successfully implemented to study flows with the large Mach number, strong shock/obstacle interactions and shock diffraction. The extension of the method to viscous flow under the adaptive Cartesian grid with hybrid overlapping bodyfitted grid is developed. The method is validated by benchmark problems and has been successfully implemented to study airfoil with ice accretion. Finally, based on an open source code, the detached eddy simulation (DES) is developed for massive separation flow, and it is used to perform the research on aerodynamic performance analysis over the wing with ice accretion

무격자 유동 해석을 위한 질점 생성 프로그램 개발

학위논문 (석사)-- 서울대학교 대학원 : 기계항공공학부, 2016. 2. 김규홍.This study aims to develop meshless point generation technique which can be applied to complicated geometry or moving boundary. Unlike the conventional finite volume method, meshless method requires only point system as a computational domain. Therefore Generation of computational domain is relatively easier than conventional FVM. In this study, meshless point generation technique is developed. For the validation, the results obtained from the meshless method were compared with the results obtained from the conventional FVM. Practical models, such as Space shuttle or Missile, were selected as a validation model. For the steady calculation, governing equation is Euler equation. And for the unsteady calculation, Euler equation in Arbitrary Lagrangian Eulerian form is selected as governing equation Meshless method developed by Huh is used for a flow solver. In both meshless method and FVM, AUSMPW+ was used for numerical flux scheme, minmod limiter was used for limiting process and LU-SGS was used for time integration. From the results, the robustness and the accuracy of meshless point generation technique are verified.Chapter 1. Introduction 1.1 Meshless method 1.2 Motivation Chapter 2. Point generation technique 2.1 Meshless point system 2.2 Near surface point system 2.3 Background point system 2.4 Local points cloud 2.5 Moving principle of points Chapter 3. Numerical Method 3.1 Governing Equation 3.2 Least square method 3.3 Spatial Discretization 3.4 Time Integration 3.5 Dual-time stepping for meshless method Chapter 4. Numerical Analysis 4.1 Steady problems 4.2 Unsteady problems Chapter 5. Conclusions Chapter 6. References 국문 초록Maste