• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.695-696
• /
• 2002

Three-dimensional blood flow in the sac of the KTAH(Korean total artificial heart) is simulated using fluid-structure interaction model. The aim of this study is to delineate the three-dimensional unsteady-blood flow in the sac of KTAH. Incompressible viscous flow is assumed for blood using the assumption of Newtonian fluid. The numerical method employed in this study is the finite element software called ADINA. Fluid-structure interaction model between blood and sac is utilized to represent the deformation of the sac by the rigid moving actuator. Three-dimensional geometry of cactus type KTAH is chosen for numerical model with prescribed pressure boundary condition on the sac surface. Blood flow is generated by the motion of moving actuator and strongly interacts with the solid material surrounding blood. High shear stress is observed mainly near the inlet and outlet of the sac.

• Interaction and multiscale mechanics
• /
• v.2 no.3
• /
• pp.263-276
• /
• 2009

This paper is intended to investigate interaction response of a train running over a suspension bridge undergoing support settlements. The suspension bridge is modeled as a single-span suspended beam with hinged ends and the train as successive moving oscillators with identical properties. To conduct this dynamic problem with non-homogeneous boundary conditions, this study first divides the total response of the suspended beam into two parts: the static and dynamic responses. Then, the coupled equations of motion for the suspended beam carrying multiple moving oscillators are transformed into a set of nonlinearly coupled generalized equations by Galerkin's method, and solved using the Newmark method with an incremental-iterative procedure including the three phases: predictor, corrector, and equilibrium-checking. Numerical investigations demonstrate that the present iterative technique is available in dealing with the dynamic interaction problem of vehicle/bridge coupling system and that the differential movements of bridge supports will significantly affect the dynamic response of the running vehicles but insignificant influence on the bridge response.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.04a
• /
• pp.5-29
• /
• 1995

A three-dimensional flow field induced by two trains passing by each other inside a tunnel is studied based on the numerical simulation of the three-dimensional compressible Euler/Navier-Stokes equations formulated in the finite difference approximation. Domain decomposition method with the FSA(fortified solution algorithm) interface scheme is used to treat this moving-body problem. The computed resluts show basic characteristic of the flow field created when two trains passing by each other. History of the pressure distributions and the aerodynamic forces acting on the trains are mailnly discussed. The results indicate that the phenomenon is complicated due to the interaction of the flow induced by two trains. Strong side force occurs between the two trains when the front portion of the opposite train passes by. It fluctuates rapidly and maximum suction force occurs when two trains are aligned side by side. The results also indicate the effectiveness of the present numerical method for moving boundary problems.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1626-1633
• /
• 2006

In the present paper, in the context of the meshless interpolation of a moving least squares (MLS) type, a novel method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The support domains for the shape functions in the MLS approximation are defined from the primary nodes, and the secondary nodes use the same support domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and moved without an additional mesh. Several numerical examples are presented to illustrate the effectiveness of the present method.

• The Pure and Applied Mathematics
• /
• v.6 no.2
• /
• pp.113-120
• /
• 1999

This paper studies the problem of an infinite confined aquifer which at time t ＜ 0 is assumed motionless. At time t = 0 crude oil seeps into the aquifer, thereby contaminating the valuable drinking water. Since the crude oil and water are im-miscible, the problem is posed as a one-dimensional two-phase unsteady moving boundary problem. A similarity solution is developed in which the moving front parameter is obtained by Newton-Ralphson iteration. A numerical scheme, involving the front tracking method, is devised employing the fourth order Runge-Kutta method. Comparison of the exact and numerical schemes shows an error of only 3%. Thus the developed numerical scheme is quite accurate in tackling more realistic problems where exact solutions are not possible.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.5-10
• /
• 2003

In order to reduce the excessive numerical dissipation, the new spatial discretization scheme is introduced. The present method in this paper has the formula that has an additional procedure of defining transferred properties at a cell-interface, based on AUSMPW+. The newly defined transferred property could eliminate numerical dissipation effectively in non-flow aligned grid system. In addition, the present method guarantees the monotonic characteristic in capturing a discontinuity. Through a stationary or moving contact discontinuity and a stationary or moving shock discontinuity, a vortex discontinuity and shock wave/ boundary layer interaction, it is verified that the accuracy of the present method is improved.

• Structural Engineering and Mechanics
• /
• v.25 no.4
• /
• pp.467-480
• /
• 2007

Two methods were developed for analyzing problems with two adjacent sub-domains modeled by different kinds of elements in finite element method. Each sub-domain can be defined independently without the consideration of equivalent division with common nodes used for the interface. These two methods employ an individual interface to accomplish the compatibility. The MLSA method uses the moving least square approximation which is the basic formulation for Element Free Galerkin Method to formulate the interface. The displacement field assumed by this method does not pass through nodes on the common boundary. Therefore, nodes can be chosen freely for this method. The results show that the MLSA method has better approximation than traditional methods.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2410-2414
• /
• 2008

Human heart valves diseased by congenital heart defects, rheumatic fever, bacterial infection, cancer may cause stenosis or insufficiency in the valves. Treatment may be with medication but often involves valve repair or replacement (insertion of an artificial heart valve). Bileaflet mechanical heart valves (BMHVs) are widely implanted to replace the diseased heart valves, but still suffer from complications such as hemolysis, platelet activation, tissue overgrowth and device failure. These complications are closely related to both flow characteristics through the valves and leaflet dynamics. In this study, the physiological flow interacting with the moving leaflets in a bileaflet mechanical heart valve (BMHV) is simulated with a strongly coupled implicit fluid-structure interaction (FSI) method which is newly organized based on the Arbitrary-Lagrangian-Eulerian (ALE) approach and the dynamic mesh method (remeshing) in FLUENT. The simulated results are in good agreement with previous experimental studies. This study shows the applicability of the present FSI model to the complicated physics interacting between fluid flow and moving boundary.

• Structural Engineering and Mechanics
• /
• v.35 no.4
• /
• pp.511-533
• /
• 2010

As a truly meshfree method, meshless equilibrium on line method (MELM), for 2D elasticity problems is presented. In MELM, the problem domain is represented by a set of distributed nodes, and equilibrium is satisfied on lines for any node within this domain. In contrary to conventional meshfree methods, test domains are lines in this method, and all integrals can be easily evaluated over straight lines along x and y directions. Proposed weak formulation has the same concept as the equilibrium on line method which was previously used by the authors for enforcement of the Neumann boundary conditions in the strong-form meshless methods. In this paper, the idea of the equilibrium on line method is developed to use as the weak forms of the governing equations at inner nodes of the problem domain. The moving least squares (MLS) approximation is used to interpolate solution variables in this paper. Numerical studies have shown that this method is simple to implement, while leading to accurate results.

• Steel and Composite Structures
• /
• v.47 no.2
• /
• pp.297-313
• /
• 2023

This study aims to extend the application of the spectral element method (SEM) to wave propagation and free vibration analysis of functionally graded (FG) plates integrated with thin piezoelectric layers, plates with tapered thickness and structure on elastic foundations. Also, the dynamic response of the smart FG plate under impact and moving loads is presented. In this paper, the dynamic stiffness matrix of the smart rectangular FG plate is determined by using the exact dynamic shape functions based on Mindlin plate assumptions. The low computational time and results' independence with the number of elements are two significant features of the SEM. Also, to prove the accuracy and efficiency of the SEM, results are compared with Abaqus simulations and those reported in references. Furthermore, the effects of boundary conditions, power-law index, piezoelectric layers thickness, and type of loading on the results are studied.