• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.2
• /
• pp.286-290
• /
• 2008

This paper summarizes the components of an explicit aeroelastic solver developed especially for the simulation of dynamic fracture events occurring during the flight of solid propellant rockets. The numerical method combines an explicit Arbitrary Lagrangian Eulerian (ALE) version of the Cohesive Volumetric Finite Element (CVFE) scheme, used to simulate the spontaneous motion of one or more cracks propagating dynamically through a domain with regressing boundaries, and an explicit unstructured finite volume Euler code to follow the flow field during the failure event. A key feature of the algorithm is the ability to adaptively repair and expand the fluid mesh to handle the large geometrical changes associated with grain deformation and crack motion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1142-1149
• /
• 2003

The present study describes a numerical analysis for simulation of the sloshing of flows with free-surface which contained in a rectangular tank moving in harmonic or pitching motion. The VOF function, representing the volume fraction of a cell occupied by the fluid, is calculated for each cells, which gives the location of the free-surface filling any some fraction of cells with fluid. The time-dependent changes of free-surface height are used for visualization subject to several conditions such as fluid height, horizontal acceleration, sinusoidal motion, and viscosity. The free-surface heights were used for comparing wall-force, which is caused by sloshing of flows. Damping effects by baffles were extensively investigated for various conditions in terms of baffle shape and position.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.210-215
• /
• 2003

A fine hot-wire is used both as a heating element and a temperature sensor in transient hot-wire method. The traditional sensor system is unnecessarily big so that it takes large fluid volume to measure the thermal conductivity. To dramatically reduce this fluid volume, a new sensor fabrication and a data processing method are proposed in this article. Contrast to the conventional and most popular two wire sensor, the new sensor system is made up of divided multiple long and short wires. Through validation experiments, it is found that the measured thermal conductivities of the glycerin are exactly same each other between the conventional and proposed new method. Also some technical considerations in arranging the multiple wires are briefly discussed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2009.04a
• /
• pp.581-586
• /
• 2009

This paper presents an optimal design of magnetorheological(MR) damper based on analytical methodology and finite element analysis. The proposed MR damper consists of MR valve and gas chamber. The MR valve is constrained in a specific volume and the optimization problem identifies geometric dimensions of the valve structure that maximize the pressure drop of the MR valve or damping force of the MR damper. In this work, the single-coil annular MR valve structure is considered. After describing the schematic configuration and operating principle of MR valve and damper, a quasi-static model is derived based on Bingham model of MR fluid. The magnetic circuit of the valve and damper is then analyzed by applying the Kirchoff’s law and magnetic flux conservation rule. Based on the quasi-static modeling and the magnetic circuit analysis, the optimization problem of the MR valve and damper is built. The optimal solution of the optimization problem of the MR valve structure constrained in a specific volume is then obtained and compared with the solution obtained from finite element method.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.11
• /
• pp.1110-1118
• /
• 2009

This paper presents an optimal design of magnetorheological(MR) damper based on analytical methodology and finite element analysis. The proposed MR damper consists of MR valve and gas chamber. The MR valve is constrained in a specific volume and the optimization problem identifies geometric dimensions of the valve structure that maximize the pressure drop of the MR valve or damping force of the MR damper. In this work, the single-coil annular MR valve structure is considered. After describing the schematic configuration and operating principle of MR valve and damper, a quasi-static model is derived based on Bingham model of MR fluid. The magnetic circuit of the valve and damper is then analyzed by applying the Kirchoff' s law and magnetic flux conservation rule. Based on the quasi-static modeling and the magnetic circuit analysis, the optimization problem of the MR valve and damper is built. The optimal solution of the optimization problem of the MR valve structure constrained in a specific volume is then obtained and compared with the solution obtained from finite element method.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.5
• /
• pp.50-60
• /
• 2016

A design methodology of the modeled test facility to conserve an injection performance of a passive safety injection system is proposed. This safety injection system is composed of a core makeup tank and a safety injection tank. Individual tanks are connected with pressure balance line on the top side and injection line on the bottom side. It is important to conserve the scaled initial injection flow rate and total injection time since this system can be operated by small gravity head without any active pumps. Differential pressure distribution of the injection line induced by the gravity head is determined by the vertical length and elevation of each tank. However, the total injection time is adjustable by the flow resistance coefficient of the injection line. The scaling methodology for the tank and flow resistance coefficient is suggested. A key point of this test facility design is a scaling analysis for the flow resistance coefficient. The scaling analysis proposed on this paper is based on the volume scaling law with the same vertical length to the prototype and can be extended to a model with a reduced vertical length. A set of passive injection test were performed for the tanks with the same volume and the different length. The test results on the initial flow rate and total injection time showed the almost same injection characteristics and they were in good agreement with the design values.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.95-98
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid} method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.405-414
• /
• 2015

Fluid transport is a key issue in the development of microfluidic systems. Recently, Myong (2014) has proposed a new concept for droplet transport without external power sources and numerically validated the results for a hypothetical 2D, initially having a hemicylindrical droplet. In this paper, the movement of an actual water droplet, initially having a 3D hemispherical shape, on horizontal hydrophilic/hydrophobic surfaces is simulated using a commercial computational fluid dynamics (CFD) package, Fluent, with VOF (volume of fluid) method. The results are compared with the 2D analysis of Myong (2014), and the transport mechanism for the actual water droplet is examined based on the numerical results of the time evolution of the droplet shape, as well as the total kinetic, gravitational, surface free and pressure energies inside the droplet.

• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.1-10
• /
• 1984

General time-volume averaged conservation equations and jump conditions for two-phase flows are derived here. The time-averaged equations for a single phase region in two-phase flow are obtained from local instant balance equations by a technique often used for single phase turbulent flow equations. The results obtained by integrating the time averaged equations over a flow volume are spatially averaged twice; first, they are averaged over a single phase region of the k-th phase and then averaged over the total volume of the k-th phase, in a flow volume. The mass, momentum, and energy conservation equations are obtained from the general time-volume averaged equations. The advantages of the present model are explained by comparing it with Ishii's model (1) and Banerjee's model (2). Finally, the assumptions and approximate terms of the equations of the THERMIT-6S are clarified.