• Coupled systems mechanics
• /
• v.1 no.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.

• Journal of computational fluids engineering
• /
• v.17 no.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
• /
• v.75 no.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.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.49 no.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.05a
• /
• 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
• /
• v.4 no.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.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.48 no.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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.7
• /
• pp.637-646
• /
• 2010

A time-domain panel method based on the potential flow theory and the time-stepping method is developed to predict the steady/unsteady aerodynamic characteristics of FM07, which is the BWB (Blended-wing body) type MAV. In the aerodynamic analyses, we used two types of the initial model(Case I) and the improved model(Case II), which is moved the gravity center toward the rear and has larger aspect ratio. In the steady aerodynamic analyses, it is revealed that improved model has higher lift to drag ratio(L/D) and more stable pitch characteristic than those of the initial model. In the unsteady aerodynamic analyses for sudden acceleration motion similar to the launch phase of MAV, it seemed that there is a rapid increase of the lift coefficient after the launch and unsteady results are good agreed compare with steady results in just a few times. In the analysis for pitch oscillation motion, which is occurred at the cruise condition of the FM07, it shows that unsteady aerodynamic coefficients looped around steady results and the improved model has more sensitive aerodynamic characteristics.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.10a
• /
• 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
• /
• v.5B no.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.