• International Journal of Aeronautical and Space Sciences
• /
• v.9 no.2
• /
• pp.79-86
• /
• 2008

High performance direct-iterative hybrid linear solver for large scale finite element problem is developed. Direct solution method is robust but difficult to parallelize, whereas iterative solution method is opposite for direct method. Therefore, combining two solution methods is desired to get both high performance parallel efficiency and numerical robustness for large scale structural analysis problems. Hybrid method mentioned in this paper is based on FETI-DP (Finite Element Tearing and Interconnecting-Dual Primal method) which has good parallel scalability and efficiency. It is suitable for fourth and second order finite element elliptic problems including structural analysis problems. We are using the hybrid concept of theses two solution method categories, combining the multifrontal solver into FETI-DP based iterative solver. Hybrid solver is implemented for our general structural analysis code, IPSAP.

• Journal of Intelligence and Information Systems
• /
• v.10 no.2
• /
• pp.165-187
• /
• 2004

Constraint Programming in AS and Optimization in OR started and have grown in different backgrounds to solve common decision-making problems in real world. This paper tries to integrate results from those different fields by suggesting a hybrid solver as an integration framework. Starting with an integrating modeling language, a way to implement a hybrid solver will be discussed using Warren's abstract machine and an finite domain constraint programming solver structures. This paper will also propose some issues rising when implementing the hybrid solver and provide their solutions.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.2
• /
• pp.44-53
• /
• 2007

A three-dimensional viscous flow solver has been developed for the prediction of the aerodynamic performance of hovering helicopter rotor blades using unstructured hybrid meshes. The flow solver utilized a vertex-centered finite-volume scheme that is based on the Roe's flux-difference splitting with an implicit Jacobi/Gauss-Seidel time integration. The eddy viscosity are estimated by the Spalart- Allmaras one-equation turbulence model. Calculations were performed at three operating conditions with varying tip Mach number and collective pitch setting for the Caradonna-Tung rotor in hover. Additional computations are made for the UH-60A rotor in hover. Reasonable agreements were obtained between the present results and the experiment in both blade loading and overall rotor performance. It was demonstrated that the present vertex-centered flow solver is an efficient and accurate tool for the assessment of rotor performance in hover.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.115-126
• /
• 1999

A two-dimensional Navier-Stokes solver using hybrid meshes is parallelized with a domain decompostion method. The focus of this paper is placed on minimizing the amount of effort in parallelizing the serial version of the solver, and this is achieved by adding an additional layer of cells to each decomposed domain. Most subroutines of the serial solver are used without modification, and the information exchange between neighboring domains is achieved using MPI(Message Passing Interface) library. Load balancing among the processors and scheduling of the message passing are implemented to reduce the overhead of parallelization, and the speed-up achieved by parallelization is measured on the transonic invisicd and turbulent flow problems. The parallelization efficiencies of the explicit Runge-Kutta scheme and the implicit point-SGS scheme are compared and the effects of various factors on the results are also studied.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.27-34
• /
• 2006

In this paper, a recently proposed parallel hybrid particle-continuum (DSMC-NS) scheme employing 3D unstructured grid for solving steady-state gas flows involving continuum and rarefied regions is described [1]. Substitution of a density-based NS solver to a pressure-based one that greatly enhances the capability of the proposed hybrid scheme and several practical experiences of implementation learned from the development and verifications are highlighted. At the end, we present some simulation results of a realistic RCS nozzle plume, which is considered very challenging using either a continuum or particle solver alone, to demonstrate the capability of the proposed hybrid DSMC-NS method.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.4
• /
• pp.320-329
• /
• 2016

An implicit method of solving the six degree-of-freedom rigid body motion equations based on the second order Adams-Bashforth-Moulten method was utilised as an improvement over the leapfrog scheme by making modifications to the rigid body motion solver libraries directly. The implementation will depend on predictor-corrector steps still residing within the hybrid Pressure Implicit with Splitting of Operators - Semi-Implicit Method for Pressure Linked Equations (PIMPLE) outer corrector loops to ensure strong coupling between fluid and motion. Aitken's under-relaxation is also introduced in this study to optimise the convergence rate and stability of the coupled solver. The resulting coupled solver ran on a free floating object tutorial test case when converged matches the original solver. It further allows a varying 70%-80% reduction in simulation times compared using a fixed under-relaxation to achieve the required stability.

• Ocean Systems Engineering
• /
• v.7 no.3
• /
• pp.195-223
• /
• 2017

This paper presents 3D numerical simulations of a Free Standing Hybrid Riser under Vortex Induced Vibration, with prescribed motion on the top to replace the motion of the buoyancy can. The model is calculated using a fully implicit discretization scheme. The flow field around the riser is computed by solving the Navier-Stokes equations numerically. The fluid domain is discretized using the overset grid approach. Grid points in near-wall regions of riser are of high resolution, while far field flow is in relatively coarse grid. Fluid-structure interaction is accomplished by communication between fluid solver and riser motion solver. Simulation is based on previous experimental data. Two cases are studied with different current speeds, where the motion of the buoyancy can is approximated to a 'banana' shape. A fully three-dimensional CFD approach for VIV simulation for a top side moving Riser has been presented. This paper also presents a simulation of a riser connected to a platform under harmonic regular waves.