• 대한조선학회논문집
• /
• 제43권5호
• /
• pp.538-549
• /
• 2006

A code is developed to simulate a viscous flow field around a deformable body using the hybrid Cartesian/immersed boundary method. In this method, the immersed boundary(IB) nodes are defined near the body boundary then velocities at the IB nodes are reconstructed based on the interpolation along the normal direction to the body surface. A new method is suggested to define the IB nodes so that a closed fluid domain is guaranteed by a set of IB nodes and the method is applicable to a zero-thickness body such as a sail. To validate the developed code, the vorticity fields are compared with other recent calculations where a cylinder orbits and moves into its own wake. It is shown the code can handle a sharp trailing edge at Reynolds number of $10^5$ under moderate requirements on girds. Finally the developed code is applied to simulate the vortex shedding behind a deforming foil with flapping tail like a fish. It is shown that the acceleration of fluids near the flapping tail contributes to the generation of the thrust for propulsion.

• 한국전산유체공학회지
• /
• 제21권2호
• /
• pp.81-89
• /
• 2016

Numerical simulations are carried out for the fluid flow and particle transport around two nearby circular cylinders in a side-by-side arrangement. The present study aims to understand the effects of the particle Stokes number and the spacing between two cylinders on particle dispersion and deposition characteristics. Simulations are based on an Eulerian-Lagrangian approach where the motion of particles is calculated by a Lagrangian approach based on one-way coupling. Results show that the flow structure is very different depending on the cylinder spacing, eventually affecting the overall pattern of particle dispersion significantly. It is also found that particles with smaller Stokes number tend to be distributed more uniformly in the wake of two cylinders, being located even inside the vortex cores. Meanwhile, particle deposition is analyzed in terms of the deposition efficiency and deposition location. The deposition efficiency of particles strongly depends on the Stokes number, whereas it is slightly affected by the cylinder spacing. The deposition location gets wider as the Stokes number increases, and it becomes asymmetric about the center of each cylinder as the cylinders get close.

• 한국항공우주학회지
• /
• 제46권7호
• /
• pp.567-574
• /
• 2018

스월은 연료의 후퇴율 증가뿐 아니라 연소 압력의 진동을 감소시키며, 하이브리드 로켓의 연소안정화에 기여하는 것으로 알려져 있다. 따라서 스월 산화제 분사에 의한 주연소실 내부의 유동 구조의 변화를 이해하고 연소안정화의 물리적 과정을 실험적으로 연구하였다. 결과에 의하면, 스월은 주연소실의 유동 구조를 변화시켜 후연소실의 500Hz 대역 p'과 q'의 발생을 억제할 뿐 아니라 두 진동의 위상차를 변화하여 상호결합(coupling)에 영향을 주고 있음을 확인했다. 또한 후연소실 화염 가시화를 통하여 스월에 의한 선회운동량의 증가로 와류발생과 흘림 등이 변화하여 연소안정화에 기여하는 것으로 분석됐다.

• 대한기계학회논문집B
• /
• 제30권12호
• /
• pp.1218-1227
• /
• 2006

In this paper a parametric study using an immersed boundary method has been carried out to investigate the effects of stable density stratification on the wakes past two identical three-dimensional hills aligned in tandem. The Reynolds number based on the uniform inlet velocity and twice the hill height was fixed at Re=300 while the Froude number based on the inlet velocity and the hill height was retained at Fr=0.2. Neutral flow without density stratification was also computed for comparison. Under a strong stratification, vertical motion of fluid particles over the three-dimensional hills is suppressed and the wake structures behind the hills become planar. Depending on the distance between the two hills, the flow pattern of each wake is significantly affected by the stratification. There is a critical hill distance at which flow characteristics drastically change. Qualitative and quantitative features of the wake interaction are reported.

• Wind and Structures
• /
• 제9권2호
• /
• pp.95-108
• /
• 2006

Wind turbines are highly complex structures for numerical flow simulation. They normally comprise of a turbine mounted on a tower thus the movement of the turbine blades and the blade/tower interaction must be captured. In addition the ground effect should also be included. There are many more important features of wind turbines and it is difficult to include all of them. A simplified set of features is chosen here for both the turbine and the tower to show how the method can begin to identify the main points connected with wind turbine wake generation and tip vortex tower interaction. An approach to modelling the rotating blades of a turbine is proposed here. The model uses point forces based on blade element theory to model the blades and takes into account their time dependent motion. This means that local instantaneous velocities can be used as a basis for the blade element theory. The model is incorporated into a large eddy simulation code and, although many important features are left out of the model, the velocity/power performance relation is generally of the correct order of magnitude. Suggested improvements to the method are discussed.

• Wind and Structures
• /
• 제20권3호
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Wind and Structures
• /
• 제13권2호
• /
• pp.169-189
• /
• 2010

This paper describes the use of coupled Computational Fluid Dynamics (CFD) and Rigid Body Dynamics (RBD) in modelling the aerodynamic behaviour of wind-borne plate type objects. Unsteady 2D and 3D Reynolds Averaged Navier-Stokes (RANS) CFD models are used to simulate the unsteady and non-uniform flow field surrounding static, forced rotating, auto-rotating and free-flying plates. The auto-rotation phenomenon itself is strongly influenced by vortex shedding, and the realisable k-epsilon turbulence modelling approach is used, with a second order implicit time advancement scheme and equal or higher order advection schemes for the flow variables. Sequentially coupling the CFD code with a RBD solver allows a more detailed modelling of the Fluid-Structure Interaction (FSI) behaviour of the plate and how this influences plate motion. The results are compared against wind tunnel experiments on auto-rotating plates and an existing 3D analytical model.

• 한국초전도학회:학술대회논문집
• /
• 한국초전도학회 1999년도 High Temperature Superconductivity Vol.IX
• /
• pp.261-261
• /
• 1999

We observed leveling-off of the resistance in granular In/InO$_x$ thin films in the zero-temperature limit. The temperature T$_b$ at which the leveling-off appears gets larger as the sheet resistance R$_n$ increases. This is consistent with the concept that the leveling-off of the resistance is due to the dissipation of the bosonic phase and that the dissipation is enhanced as the resistance increases. The magnetic field dependence of the saturated resistance R$_b$ at low temperatures fits the modified square-root cusp-like form R$_b$/R$_n$=α exp[-b(B/B$_c$-1)$^{-1/2}$] for the magnetic field in the range B$_c$$_f$ where B$_c$ is the onset magnetic field of the resistance leveling-off. α and b are constants of order 1. For B>B$_f$ tansport properties are described by the theory of the fermi insulator. From the results, we attribute the leveling-off to the dissipative quantum tunneling of vortices, which supports the models predicting the vortex-motion-induced insulating phase related with the concept like"dirty boson" [1]l and "hose metal" [2].

• International Journal of Aeronautical and Space Sciences
• /
• 제10권2호
• /
• pp.23-33
• /
• 2009

The inherent aeromechanical complexity of a rotor system necessitated the comprehensive analysis code for helicopter rotor system. In the present study, an aerodynamic analysis module has been developed as a part of rotorcraft comprehensive program. Aerodynamic analysis module is largely classified into airload calculation routine and inflow analysis routine. For airload calculation, quasi-steady analysis model is employed based on the blade element method with the correction of unsteady aerodynamic effects. In order to take unsteady effects - body motion effects and dynamic stall - into account, aerodynamic coefficients are corrected by considering Leishman-Beddoes's unsteady model. Various inflow models and vortex wake models are implemented in the aerodynamic module to consider wake induced inflow. Specifically, linear inflow, dynamic inflow, prescribed wake and free wake model are integrated into the present module. The aerodynamic characteristics of each method are compared and validated against available experimental data such as Elliot's induced inflow distribution and sectional normal force coefficients of AH-1G. In order to validate unsteady aerodynamic model, 2-D unsteady model for NACA0012 airfoil is validated against aerodynamic coefficients of McAlister's experimental data.

• 신재생에너지
• /
• 제3권3호
• /
• pp.45-53
• /
• 2007

This paper demonstrates a computational method in predicting aerodynamic noise generated from wind turbines. Low frequency noise due to displacement of fluid and leading fluctuation, according to the blade passing motion, is modelled on monopole and dipole sources. They are predicted by Farassat 1A equation. Airfoil self noise and turbulence ingestion noise are modelled upon quadrupole sources and are predicted by semi-empirical formulas composed on the groundwork of Brooks et al. and Lowson. Aerodynamic flow in the vicinity of the blade should be obtained first, while noise source modelling need them as numerical inputs. Vortex Lattice Method(VLM) is used to compute aerodynamic conditions near blade. In the use of program X-foil [M.Drela] boundary layer characteristics are calculated to obtain airfoil self noise. Wind turbine blades are divided into spanwise unit panels, and each panel is considered as an independent source. Retarded time is considered, not only in low frequency noise but also In turbulence ingestion noise and airfoil self noise prediction. Numerical modelling is validated with measurement from NREL [AOC15/50 Turbine) and ETSU [Markham's VS45] wind turbine noise measurements.