• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.141-148
• /
• 2007

We present an immersed boundary (IB) method for 3D simulation of flappingflags in a uniform flow. The proposed formulation is manipulated on the basis of an efficient Navier-Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A direct numerical method is developed to calculate the flag motion, with the elastic force treated implicitly. The fluid motion defined on an Eulerian grid and the flag motion defined on a Lagrangian grid are independently solved and the mass of flag is handled in a natural way. An additional momentum forcing is formulated from the flag motion equation in a way similar with the direct-forcing IB formulation and acts as the interaction force between the flag and ambient fluid. A series of numerical tests are performed and the present results are compared qualitatively and quantitatively with previous studies. The instantaneous flag motion is analyzed under different conditions and surrounding vortical structures are identified. The effects of physical parameters on the flapping frequency are studied.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.379-388
• /
• 2008

A hybrid Cartesian/immersed boundary code is expanded to simulate flow field around a three-dimensional body which undergoes large dynamic deformation. Immersed boundary nodes are automatically distributed based on the edges crossing triangles on body boundary. Velocity vectors are reconstructed at those immersed boundary nodes along local normal lines to the boundary. The reconstruction of pressure is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other results on laminar flow over a sphere. The code is applied to simulate flow around a foil which is attached to a body of revolution and rotates under periodic deformation. The periodic variation of the tip vortex is observed and the effects of the deformation on hydrodynamic force acting on the body are investigated.

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

In this study, it is performed numerical simulation on multi-phase flow by means of CIP-CSI2 scheme. It is applied In a two-phase free surface flow problem at a high density ratio equivalent to that of an air-water system, for examining the computational capability. The method that is being developed and improved is a CIP(Constrained Interpolation Profile) and CSL2(Conservative Semi-Lagrangian) based Cartesian Grid Method.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.6
• /
• pp.459-479
• /
• 2014

This paper considers experimental and numerical studies on added resistance in waves. As the numerical methods, three different methods, strip method, Rankine panel method and Cartesian-grid method, are applied. The computational results of vertical motion response and added resistance are compared with the experimental data of Series 60($C_B=0.8$) hull, S175 containership and KVLCC2 hull. To investigate the influence of above-still water hull form, a Rankine panel method is extended to two nonlinear methods: weakly-nonlinear and weak-scatterer approaches. As nonlinear computational models, three ships are considered: original KVLCC2 hull, 'Ax-bow' and 'Leadge-bow' hulls. Two of the three models are modified hull forms of original KVLCC2 hull, aiming the reduction of added resistance. The nonlinear computational results are compared with linear results, and the improvement of computational result is discussed. As experimental approach, a series of towing-tank experiment for ship motions and added resistance on the three models (original KVLCC2 hull, 'Ax-bow' and 'Leadge-bow') are carried out. For the original KVLCC2 hull, uncertainty analysis in the measurement of vertical motion response and added resistance is performed in three waves conditions: ${\lambda}/L=0.5$, 1.1, 2.0. From the experimental results, the effects of hull form on added resistance are discussed.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.35-43
• /
• 2008

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.200-208
• /
• 2007

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.2
• /
• pp.188-195
• /
• 2010

Green water load is an important parameter to be considered in designing a modern ship or offshore structures like FPSO and FSRU. In this research, a numerical simulation method for green water phenomenon is introduced. The Navier-Stokes equations and the continuity equation are used as governing equations. The equations are calculated using Finite Difference Method(FDM) in rectangular staggered grid system. To increase the numerical accuracy near the body, the Cartesian cut cell method is employed. The nonlinear free-surface during green water incident is defined by Marker-density method. The green waters on a box in regular waves are simulated. The simulation results are compared with other experimental and computational results for verification. To check the applicability to moving ship, the green water of the ship which is towed by uniform force in regular wave, is simulated. The ship is set free to heave and to surge.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.45-50
• /
• 2000

We study a conceptual numerical model on shock-vortex interaction setting an impulsive shock in a compressible vertex. Navier-Stokes equations are solved for the investigation of interactive structure and acoustic wave propagation. The rotationally symmetric vortex enforces two compression-expansion pairs resultantly forming a quadrupolar shape. These compressive and expansive waves cylindrically propagate to the far field and turn to acoustic waves. Using a fine uniform Cartesian grid system and a TVD-high resolution method, the flow data irl: precisely obtained to extend our interest to the sound source.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.25 no.12
• /
• pp.1559-1565
• /
• 1988

In computed tomography, the errors asociated with interpolation in the reconstruction process degrade the reconstructed image and may cause divergence unless a large number of rays is used. A new data acquisition scheme without interpolation is developed in this paper. Samples (projection data ) are taken in phase with samples of the Cartesian grid to eliminated errors associated with interpolation process.