• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.52-61
• /
• 2004

In the present study, a solution algorithm for the computation of unsteady flows on unstructured meshes is presented. Dual time stepping is incorporated to achieve the second-order temporal accuracy while reducing errors associated with linearization and factorization. This allows any time step size, which is suitable for considering physical phenomena of interest. The Gauss-Seidel scheme is used to solve the linear system of equations. A special treatment based on spring analogy is made to handle meshes with high aspect-ratio cells. The present method was validated by comparing the results with experimental data and those obtained from rigid motion.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.1
• /
• pp.66-78
• /
• 2007

A multi-point aerodynamic shape optimization technique has been developed for helicopter rotor blades in hover based on a continuous adjoint method on unstructured meshes. The Euler flow solver and the continuous adjoint sensitivity analysis were formulated on the rotating frame of reference. The 'objective function and the sensitivity were obtained as a weighted sum of the values at each design point. The blade section contour was modified by using the Hicks-Henne shape functions. The mesh movement due to the blade geometry change was achieved by using a spring analogy. In order to handle the repeated evaluation of the design cycle efficiently, the flow and adjoint solvers were parallelized based on a domain decomposition strategy. A solution-adaptive mesh refinement technique was adopted for the accurate capturing of the wake. Applications were made to the aerodynamic shape optimization of the Caradonna-Tung rotor blades and the UH-60 rotor blades in hover.

• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.33-48
• /
• 2009

In the present study, a fast grid deformation technique has been incorporated into the unsteady compressible and incompressible viscous flow solvers on unstructured hybrid meshes. An algebraic method based on the basis decomposition of normal edge vector was used for the deformation of viscous elements, and a ball-vertex spring analogy was adopted for inviscid elements among several spring analogy methods due to its robustness. The present method was validated by comparing the results obtained from the grid deformation and the rigid motion of entire grids. Fish swimming motion of an NACA0012 airfoil and flapping wing motion of a generic fighter were also simulated to demonstrate the robustness of the present grid deformation technique.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.47-50
• /
• 2008

The unsteady vortical flow over ONERA 70-degree delta wing was simulated using RANS and DES flow solvers on hybrid unstructured meshes. A study of time accuracy is accomplished to determine the effects of time step and number of iteration in pseudo-time stepping on numerical solutions. The grid size test is also performed to demonstrate that DES can be used to capture more accurately the unsteady vortical flow features over RANS simulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.25-30
• /
• 2003

The numerical simulations with unstructured mesh by cell-centered and vertex-centered approaches were performed for the quadrilateral and triangular meshes. For the 2-D incompressible supersonic vortex flow, the simulation results and the analytic solution were compared and the accuracy was assessed. The calculation efficiency was measured by the parameter defined by the consumed CPU time multiplied by absolute error, As a results, equilateral triangular mesh yielded the best accuracy and efficiency among the tested meshes.

• Structural Engineering and Mechanics
• /
• v.14 no.6
• /
• pp.625-647
• /
• 2002

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.

• Korean Journal of Computational Design and Engineering
• /
• v.4 no.2
• /
• pp.153-161
• /
• 1999

An atomatic mesh generation scheme with unstructured quadrilateral elements on trimmed surfaces has been developed. Trimmed surfaces are often encountered in modeling of structures with complex shapes such as aircrafts, automobile structures, pressure vessels and etc. For unstructured mesh generation with quadrilateral elements, a domain decomposition algorithm employing loop operators has been used. Mesh generation on trimmed surface is performed in three steps. First, trimmed surfaces with holes or cuts are transformed to th largest projection planes in which the meshes are constructed. The constructed meshes are transformed to the u-v parametric plane and then finally to the original 3D surfaces. Th exact locations of holes or cuts in projection planes are determined by the Newton-Raphson method. Sample meshes are constructed to demonstrate the effectiveness of the proposed algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.294-301
• /
• 2011

A two-dimensional hybrid flaw solver has been developed for the accurate and efficient simulation of steady and unsteady flaw fields. The flow solver was cast to accommodate two different topologies of computational meshes. Triangular meshes are adopted in the near-body region such that complex geometric configurations can be easily modeled, while adaptive Cartesian meshes are, utilized in the off-body region to resolve the flaw more accurately with less numerical dissipation by adopting a spatially high-order accurate scheme and solution-adaptive mesh refinement technique. A chimera mesh technique has been employed to link the two flow regimes adopting each mesh topology. Validations were made for the unsteady inviscid vol1ex convection am the unsteady turbulent flaws over an NACA0012 airfoil, and the results were compared with experimental and other computational results.