• Ocean Systems Engineering
• /
• 제8권2호
• /
• pp.223-246
• /
• 2018

In this paper the results of CFD simulations, that were carried out to study the impact pressures acting on a symmetric wedge during water entry under the influence of gravity, are presented. The simulations were done using a solver implementing finite volume discretization and using the VOF scheme to keep track of the free surface during water entry. The parameters such as pressure on impact, displacement, velocity, acceleration and net hydrodynamic forces, etc., which govern the water entry process are monitored during the initial stage of water entry. In addition, the results of the complete water entry process of wedges covering the initial stage where the impact pressure reaches its maximum as well as the late stage that covers the rebound process of the buoyant wedge are presented. The study was conducted for a few touchdown velocities to understand its influence on the water entry phenomenon. The simulation results are compared with the experimental measurements available in the literature with good accuracy. The various computational parameters (e.g., mesh size, time step, solver, etc.) that are necessary for accurate prediction of impact pressures, as well as the entry-exit trajectory, are discussed.

• 소음진동
• /
• 제8권4호
• /
• pp.632-642
• /
• 1998

A method based on finite element discretization is developed for optimizing the polarization profile of PVDF film to create the modal transducer for specific modes. Using this concept, one can design the modal transducer in two-dimensional structure having arbitrary geometry and boundary conditions. As a practical means for implementing this polarization profile without repoling the PVDF film the polarization profile is approximated by optimizing electrode patterns, lamination angles, and poling directions of the multi-layered PVDF transducer. This corresponds to the approximation of a continuous function using discrete values. The electrode pattern of each PVDF layer is optimized by deciding the electrode of each finite element to be used or not. Genetic algorithm, suitable for discrete problems, is used as an optimization scheme. For the optimization of each layers lamination angle, the continuous lamination angle is encoded into discrete value using binary 5 bit string. For the experimental demonstration, a modal sensor for first and second modes of cantilevered composite plate is designed using two layers of PVDF films. The actuator is designed based on the criterion of minimizing the system energy in the control modes under a given initial condition. Experimental results show that the signals from residual modes are successfully reduced using the optimized multi-layered PVDF sensor. Using discrete LQG control law, the modal peaks of first and second modes are reduced in the amount of 12 dB and 4 dB, resepctively.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 1992년도 추계학술대회논문집; 반도아카데미, 20 Nov. 1992
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• 대한토목학회논문집
• /
• 제12권4_1호
• /
• pp.1-8
• /
• 1992

면내거동과 휨거동을 받는 원형구멍을 갖는 유한평판에서 원형구멍 주위의 응력을 모델링하기 위해 p-version 유한요소법이 제시되었으며, 또한 동일한 문제로써 원형구멍으로부터 발생된 균열해석을 위해 균열확장법이 사용되었다. 적분형 르장드르함수에 기초한 p-version 유한요소법이 원형구멍 주위의 응력경사가 심한 기하형상을 모델링하는데 적합함을 보여준다. 한편, 원형 경계조건을 표현하는데 이산화오차를 피하기 위해 초유한사상기법이 사용되었다. 앞에서 제시된 방법을 통한 수치해석 결과는 Nisida, Howland, Newman 등의 실험 및 이론결과와 종래의 유한요소법에 의한 수치해석결과와 비교하여 우수한 값을 보여 주고 있다.

• 대한기계학회논문집B
• /
• 제25권10호
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.

• 에너지공학
• /
• 제21권3호
• /
• pp.228-236
• /
• 2012

큰 에디 모사(LES)는 복잡한 연소실 유동에서 RANS 모델과 비교해 난류 유동장에 대해 모델의 보편성과 더 정확한 결과를 제공하기 때문에 점차적으로 사용이 증대되고 있다. 내연 기관의 연소실 내 난류 유동장 해석을 큰 에디 모사를 사용하여 수행하였다. 이산화 방식, 초기 조건, 시간 간격과 SGS 모델과 같은 모델과 수치 인자에 따른 영향을 평가하였다. SGS 모델을 사용한 LES 모사는 실험치와 유사한 결과를 보여주었다.

• International Journal of Aeronautical and Space Sciences
• /
• 제2권2호
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2008년도 하계학술발표대회 논문집
• /
• pp.1366-1371
• /
• 2008

Four simulation models of plastic plate heat exchangers are designed and simulated. The flat plate type heat exchanger is designed as the reference model in order to evaluate how much thermal performance increases. The turbulence promoter type heat exchanger is fabricated with cylindrical-type vortex generators and rib-type turbulence promoters. The corrugate type is obtained from the conventional stainless steel compact heat exchangers, which are called the herringbone-type compact heat exchangers. The dimple type heat exchanger has a number of dimples on its surface. In this study, the flow and heat transfer characteristics of the plastic plate heat exchanger are investigated using numerical simulation and compared with experimental results. The flows are assumed as a three-dimensional, incompressible and turbulent model. The standard k-$\varepsilon$ model is used as the turbulent flow modeling, the SIMPLE algorithm is used to treat the coupling between pressure and velocity, and first order upwind scheme is used for discretization of momentum, turbulent and energy. The computational analysis and experimental results both show that the friction coefficient and Nu number is highest in the corrugate type.

• 대한조선학회논문집
• /
• 제42권6호
• /
• pp.583-592
• /
• 2005

The turbulent free surface flow around a self-propelled KRISO 138K LNG Carrier is numerically simulated using the finite volume based multi-block RANS code, WAVIS developed at HRISO. The realizable k-$\varepsilon$ turbulence model with a wail function is employed for the turbulence closure. The free surface is captured with the Level-Set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. In order to obtain an accurate free surface solution and stable convergence, the computations are executed with a proper fine grid refinement around the free surface and with an adoption of implicit discretization scheme for the Level-Set formulation. The computed velocity vectors at the several stations and wave patterns show a good agreement with the experimental results measured at the KRISO towing tank.

• International Journal of Ocean System Engineering
• /
• 제1권3호
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.