• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• Ocean Systems Engineering
• /
• v.5 no.4
• /
• pp.245-259
• /
• 2015

The numerical simulation of wave slamming on a 3D platform deck was investigated using a coupled Level-Set and Volume-of-Fluid (CLSVOF) method for overset grid system incorporated into the Finite-Analytic Navier-Stokes (FANS) method. The predicted slamming impact forces were compared with the corresponding experimental data. The comparisons showed that the CLSVOF method is capable of accurately predicting the slamming impact and capturing the violent free surface flow including wave slamming, wave inundation and wave recession. Moreover, the capability of the present CLSVOF method for overset grid system is a prominent feature to handle the prediction of wave slamming on offshore structure.

• Computers and Concrete
• /
• v.10 no.5
• /
• pp.489-505
• /
• 2012

This paper presents a simplified grid beam model for simulating the nonlinear response of reinforced concrete flat-plate structures. The beam elements are defined with nonlinear behavior for bending moment and torsion. The flexural stiffness and torsional strength of the beam elements are defined based on experimental data to implicitly account for slab two-way bending effects. A failure criterion that considers the interaction between the punching strength and slab flexural behavior is incorporated in the model. The effects of bond-slip of slab reinforcement on connection stiffness are examined. The proposed grid beam model is validated by simulating large-scale tests of slab-column connections subjected to concentric gravity loading and unbalanced moment. This study also determines the critical parameters for a hysteretic model used to simulate flat-plates subjected to cyclic lateral loading.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.653-661
• /
• 1994

In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 supersonic/transonic compressor cascade flow. H-type grid was chosen for its simplicity in applying cyclic tridiagonal matrix algorithm along with conventional slip/no-slip boundary conditions. The thin-layer algebraic model of Baldwin-Lomax was employed for the calculation of turbulent flows. The test case inlet Mach No. was 1.612 and inlet/exit pressure ratio($P_2/P_1$) was 2.15. The results were compared with experimental results from current method were compared well in suction surface with the experiments and other computational results; however, not well in pressure surface. It might be due to the complex flowfields such as shock/boundary layer interaction, turbulence, and flow separation, etc. In the future, a proper turbulence modelling and adaptive grid system will be studied to improve the solution quality.

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

An iterative time marching procedure for solving incompressible internal flow has been applied to the flow around a rotor-stator. This procedure solves three-dimensional incompressible Reynolds-averaged Navier-Stokes equation on a moving, time-deforming, non-orthogonal body-fitted grid using second-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. To handle rotationg geometry, the multiblock technique is applied and the overall flow domain is subdivided into two blocks. In each block, a grid is generated and flowfield is solved independently of the other blocks. The boundary data for each block is provided by the neighboring blocks using bilinear interpolation technique.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.293-299
• /
• 2001

Since the previous cut-and-try design algorithm require much cost and time, it has recently been concerned the automatic design technique using the CFD and optimum design algorithm. In this study, the Navier-Stokes equations is solved to consider the more detail viscous flow informations of cascade interaction and O-H multiblock grid system is generated to impose an accurate boundary condition. The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. To validate present procedure, the time averaged aerodynamic loads are compared with experiment and good agreement obtained. Once the N-S equations have been solved, the computed aerodynamic loads may be used to computed the sensitivities of the aerodynamic objective function. The Modified Method of feasible Direction(MMFD) is usef to compute the

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.04a
• /
• pp.658-659
• /
• 2009

In this study, fluid-induced vibration (FIV) analyses have been conducted for tall building structure. In order to investigate the aeroelastic responses of tall building due to wind load, advanced computational analysis system based n computational fluid dynamics(CFD) and computational structural dynamics (CSD) has been developed. Fluid domains are modeled using the computational grid system with local grid deforming technique. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of tall structure for fluid-structure interaction (FSI) problems. Detailed aeroelastic responses and results are presented to show the physical phenomenon of the tall building.

• Steel and Composite Structures
• /
• v.31 no.1
• /
• pp.43-51
• /
• 2019

This paper is concerned with the numerical calculation of mixed-mode stress intensity factors (SIFs) of 2-D isotropic functionally graded materials (FGMs) by the natural element method (more exactly, Petrov-Galerkin NEM). The spatial variation of elastic modulus in non-homogeneous FGMs is reflected into the modified interaction integral ${\tilde{M}}^{(1,2)}$. The local NEM grid near the crack tip is refined, and the directly approximated strain and stress fields by PG-NEM are enhanced and smoothened by the patch recovery technique. Two numerical examples with the exponentially varying elastic modulus are taken to illustrate the proposed method. The mixed-mode SIFs are parametrically computed with respect to the exponent index in the elastic modulus and external loading and the crack angle and compared with the other reported results. It has been justified from the numerical results that the present method successfully and accurately calculates the mixed-mode stress intensity factors of 2-D non-homogeneous functionally graded materials.

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.373-377
• /
• 2023

The Monte Carlo (MC) method has become an indispensable part of the nuclear radiation research field. Several widely used and well-known MC packages were developed for simulation of radiation transport and interaction with matter. All these MC packages require users to prepare an input script. The input script can become lengthy for complex models. The process of preparing these input scripts is time-consuming and error-prone. In the present work, we have developed an open-source GUI computer program for modelling radiation transport and interaction in multi-segmented slab phantoms using grid-based system for the widely used PHITS MC package. The developed tools would be useful for future users of PHITS MC package and particularly inexperienced users. The present program is distributed under GPL license and all users can freely download, modify and redistribute the program without any restrictions.