• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.203-206
• /
• 2004

In this study, a numerical analysis has been performed for a three-dimensional pulsatile blood flow associated with the elastic blood vessel and curved bileaflet for multiple cycles in terms of fluid-structure interaction. Here, blood has been assumed as a Newtonian, incompressible fluid. Pressure profiles have been used as boundary conditions at the ventricle and the aorta. From this analysis, the motion of the leaflet has been observed with fluttering phenomenon and rebound, and the flow fields of blood have been obtained with recirculation and regurgitation. The results can contribute to the development of design methodology for the curved bileaflet mechanical heart valve.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.591-596
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a planar two-dimensional channel, and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-section area, indicative of wave speed flow limitation or "choking" (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, flow rate decreases. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is ${\le}$ 2 the area throat is located near the downstream end; as wall taper is increased to ${\ge}$ 3 the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.

• International Journal of Aeronautical and Space Sciences
• /
• v.1 no.1
• /
• pp.21-28
• /
• 2000

The flowfield of hypersonic shock-shock interaction has been simulated using a two-dimensional Navier-Stokes code based on AUSMPW+ scheme. AUSMPW+ scheme is a new hybrid flux splitting scheme, which is improved by introducing pressure-based weight functions to eliminate the typical drawbacks of AUSM-type schemes, such as non-monotone pressure solutions. To study the real gas effects, three different gas models are taken into account in the present paper: perfect gas, equilibrium flow and non equilibrium flow. It has been investigated how each gas model influences on the peak surface loading, such as wall pressure and wall heat transfer, and unsteady structure of flowfield in the region of shock-shock interaction. With the results, the value of peak pressure is not sensitive to the real gas effects nor to the wall catalyticity. However, the value of peak heat transfer rates is affected by the real gas effects and the wall catalyticity. Also, the structure of the flowfield changes drastically in the presence of real gas effects.

• International Journal of Fluid Machinery and Systems
• /
• v.6 no.2
• /
• pp.87-93
• /
• 2013

In order to study the effect of the fluid-structure interaction (FSI) on the simulation results, the external characteristics and internal flow features of a diffuser pump were analyzed with a two-way flow solid coupling method. And the static and dynamic structure analysis of the blade was also caculated with the FEA method. The steady flow field is based on Reynolds Averaged N-S equations with standard $k-{\varepsilon}$ turbulent model, the unsteady flow field is based on the large eddy simulation, and the structure response is based on elastic transient structural dynamic equation. The results showed that the effect of FSI on the head prediction based on CFD really exists. At the same radius, the van mises stress on the nodes closed shroud and hub was larger than other nodes. A large deformation region existed near inlet side at the middle of blades. The strength of impeller satisfied the strength requirement with static stress analysis based on the fourth strength theory. The dynamic stress varied periodically with the impeller rotating. It was also found that the fundamental frequency of the dynamic stress is the rotating frequency and its harmonic frequency. The frequency of maximum stress amplitude at node 1626 was 7 times of the rotating frequency. The frequency of maximum stress amplitude at node 2328 was 14 times of the rotating frequency. No matter strength failure or fatigue failure, the root of blades near shroud is the key region to analyse.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.126-131
• /
• 2004

The interaction between the hypersonic free stream and the side jet flow at high altitudes is investigated by direct simulation Monte Carlo(DSMC) method. Since there is a great difference in density between the free stream and the side jet flow, the weighting factor technique which could control the number of simulation particles, is applied to calculate these two flows simultaneously. Chemical reactions are not considered in the calculation. For validation, the corner flow passing between a pair of plates that are perpendicularly attached is solved. The side jet flow is then injected into this comer flow and solution is found for the merged flow. Results are compared with the experiments. For a more realistic rocket model, the flow past a blunted cone cylinder shape is solved. The leeward or windward jet injection is merged with this flow. The effect on the rocket surface is observed at various flow angles. The lambda effect and the wake structure are found like low attitudes. High interaction between the free stream and the side jet flow is observed when the side jet is injected in the windward direction.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.7
• /
• pp.642-648
• /
• 2005

This study has been conducted to investigate the thermopneumatic and flow characteristics of diffuser/nozzle based thermopneumatic micropumps. In this study, a transient three-dimensional numerical analysis using FSI (Fluid-Structure Interaction) model has been employed to analyze the effects of the interaction between the membrane and two fluids (air and water) in the thermopneumtic micropump. The transient temperature and pressure in the cavity, the transient displacements of the membrane and the net flow rate of the micropump have been closely calculated for the frequency of 1 Hz. It has been found that the difference of the flow rates at the inlet and outlet is larger in the cooling period than in the heating period and that the duty ratio is very important in association with pump performance because the temperature in the cavity ascends drastically in the heating period and descends slowly in the cooling period. This study can be regarded as fundamental understandings for the design and analysis of thermopneumatic micropumps.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.151-158
• /
• 2009

To implement the insects' flapping flight for developing flapping MAVs(micro air vehicles), the unsteady flow characteristics of the insects' forward flight is investigated. In this paper, two-dimensional FSI(Fluid-Structure Interaction) simulations are conducted to examine realistic flow features of insects' flapping flight and to examine the flexibility effects of the insect's wing. The unsteady incompressible Navier-Stokes equations with an artificial compressibility method are implemented as the fluid module while the dynamic finite element equations using a direct integration method are employed as the solid module. In order to exchange physical information to each module, the common refinement method is employed as the data transfer method. Also, a simple and efficient dynamic grid deformation technique based on Delaunay graph mapping is used to deform computational grids. Compared to the earlier researches of two-dimensional rigid wing simulations, key physical phenomena and flow patterns such as vortex pairing and vortex staying can still be observed. For example, lift is mainly generated during downstroke motion by high effective angle of attack caused by translation and lagging motion. A large amount of thrust is generated abruptly at the end of upstroke motion. However, the quantitative aspect of flow field is somewhat different. A flexible wing generates more thrust but less lift than a rigid wing. This is because the net force acting on wing surface is split into two directions due to structural flexibility. As a consequence, thrust and propulsive efficiency was enhanced considerably compared to a rigid wing. From these numerical simulations, it is seen that the wing flexibility yields a significant impact on aerodynamic characteristics.

• Structural Engineering and Mechanics
• /
• v.37 no.4
• /
• pp.401-413
• /
• 2011

The aerodynamic stability of orthotropic tensioned membrane structures with rectangular plane is theoretically studied under the uniform ideal potential flow. The aerodynamic force acting on the membrane surface is determined by the potential flow theory in fluid mechanics and the thin airfoil theory in aerodynamics. Then, based on the large amplitude theory and the D'Alembert's principle, the interaction governing equation of wind-structure is established. Under the circumstances of single mode response, the Bubnov-Galerkin approximate method is applied to transform the complicated interaction equation into a system of second order nonlinear differential equation with constant coefficients. Through judging the stability of the system characteristic equation, the critical divergence instability wind velocity is determined. Finally, from different parametric analysis, we can conclude that it has positive significance to consider the characteristics of orthotropic and large amplitude for preventing the instability destruction of structures.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.607-610
• /
• 2002

As machines become smaller and faster multileaf foil bearings are used to overcome the problems with heat, friction and wear Systems with foil bearings do not need a separate system for lubrication. These bearings are self acting and are therefore green systems. Until now, there have been many studies on the structural and dynamical performances. Therefore the object of the present study is to predict the flow and structural characteristics by using the Fluid/structure interaction method. The increase in RPM led to the increase in pressure, temperature difference, maximum velocity, Mach number, shear stress and torque. In the case of 90,000 RPM effects such as choking led to a non-lineararity in the system. Also the effect of eccentricity ratio was observed and showed that eccentricity increased the maximum pressure and the density difference while decreasing the shear stress and torque.