• Coupled systems mechanics
• /
• 제1권4호
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Structural Engineering and Mechanics
• /
• 제75권4호
• /
• pp.465-475
• /
• 2020

The analytical approach for stability and response of iced conductor under uniform wind or turbulent wind is presented in this study. A nonlinear dynamic model is established to describe the motion of iced conductor galloping. In the case of uniform wind, the stability condition is derived by analyzing the eigenvalue associated with linearized matrix; The first order and second order approximation of galloping amplitude are obtained using multi-scale method. However, real wind has random characteristics essentially. To accurately evaluate the performance of the galloping iced conductor, turbulence wind should be described by random processes. In the case of turbulence wind, the Lyapunov exponent is conducted to judge the stability condition; The probability density of displacement is obtained by using the path integral method to predict galloping amplitude. An example is proposed to verify the effectiveness of the previous methods. It is shown that the fluctuating component of wind has little influence on the stability of iced conductor, but it can increase galloping amplitude. The analytical results on stability and response are also verified by numerical time stepping method.

• 대한전기학회논문지:전기기기및에너지변환시스템부문B
• /
• 제49권2호
• /
• pp.93-101
• /
• 2000

This paper presents a novel method for the transient of eddy currents in the dampers of a super-conducting synchronous generator(SCG). The method proposes a 2-D corrected model which takes into account the influence of leakage fluxes of the field winding ends by increasing the effective air gap in order to consider the high precision of the analysis for the conventional 2-D model. The electromagnetic fields for the corrected model are analyzed by the time-stepping finite element method, thus the eddy currents in the dampers and electro-motive forces(EMF) in the stator windings are calculated. As the results, it is proved the presented method is comparatively accurate by comparing measured phase EMF values and the simulation ones, where about 6.4% error at the maximum value of EMF is occurred between them.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• Ocean Systems Engineering
• /
• 제4권2호
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• 대한전기학회논문지:전기기기및에너지변환시스템부문B
• /
• 제48권9호
• /
• pp.467-475
• /
• 1999

This paper presents an analysis method of a forward converter, using both the finite element method considering the external circuit and a state variables equation. The converter operates at 50kHz and its one period is divided into two modes for the simplicity of the analysis. In the first mode, the switching transistor turns on and an input power is transferred into the load by the electromagnetic conversion action of a ferrite transformer. In the second mode, the switching transistor turns off and the stored energy in an inductor is delivered to the load, and the transformer core is demagnetized by the reset winding current. In this paper, time-stepping finite element method taking into account the on-state electrical circuit of the converter in used to analyze both the electrical circuit and electromagnetic field of the magnetic device during the first mode and the demagnetization period of the transformer core. Then a state variables equation for the circuit which the inductor current flows is constituted and solved during the second mode. As a result, the simulation results have been good agreement with the results obtained form experiment.

• 한국항공우주학회지
• /
• 제38권7호
• /
• pp.637-646
• /
• 2010

포텐셜 기저 패널법과 시간전진법을 통합한 시간영역 패널법을 이용하여 융합익기 형상 초소형 무인기 FM07에 대한 정상/비정상 공력해석을 수행하였다. 융합익기 형상 무인기인 FM07의 공력해석에는 초기 설계형상(Case I)과 가로세로비를 증가시키고 무게중심을 후방으로 이동시킨 개선형상(Case II)이 사용되었다. 정상 공력해석을 통해 FM07 무인기의 개선형상이 초기형상에 비하여 보다 큰 양항비와 높은 피칭 안정성을 나타냄을 확인하였다. FM07의 비정상 급가속(발사체 이륙단계)을 나타낸 공력해석에서는 초기 급격한 증가를 보이는 공력계수들이 수초 후에 안정화되어 정상상태에 근사한 값을 나타내었다. FM07의 순항시 발생할 수 있는 피치진동운동에 대한 해석을 수행하여 진동에 따른 공력계수의 이력 현상을 확인하였으며, 개선형상이 보다 큰 민감도를 갖는 것을 확인하였다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2003년도 제21회 추계학술대회 논문집
• /
• pp.157-163
• /
• 2003

A systematic numerical experiment has been conducted to study droplet gasification in high pressure environments with pressure oscillations. The general frame of previous rigorous model[1] is retained but tailored for flash equilibrium calculation of vapor-liquid interfacial thermodynamics. Time-dependent conservation equations of mass, momentum, energy, and species concentrations are formulated in axisymmetric coordinate system for both the droplet interior and ambient gases. In addition, a unified property evaluation scheme based on the fundamental equation of state and empirical methods are used to find fluid thermophysical properties over the entire thermodynamic domain of interest. The governing equations with appropriate physical boundary conditions are numerically time integrated using an implicit finite-difference method with a dual time-stepping technique. A series of calculation have been carried out to investigate the gasification of an isolated n-pentane droplet in a nitrogen gas environment over a wide range of ambient pressures and frequencies. Results show that the mean pressures and frequencies of the ambient gas have strong influences on the characteristics of the droplet gasification. The amplitude of the response increases with increasing pressure, and the magnitude of the vaporization response increases with the frequency.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• 제5B권3호
• /
• pp.253-257
• /
• 2005

By using a coupled circuit, time-stepping, two-dimensional finite element method (2-D FEM), the performance of a stand-alone permanent-magnet synchronous generator (PMSG) with inset rotor can be computed without involving the classical two-axis model. The effects of interpolar air gap length and armature resistance on the load characteristics are investigated. It is shown that the interpolar flux density, and hence the amount of voltage compensation, is affected by magnetic saturation. Validity of the coupled circuit and field analysis is confirmed by experiments on a prototype generator. The machine exhibits an approximately level load characteristic when it is supplying an isolated unity-power-factor load.