• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.

• 한국연소학회:학술대회논문집
• /
• 2003.05a
• /
• pp.57-62
• /
• 2003

A general procedure of obtaining reliable one-step kinetics model for hydrocarbon mixture from the fully detailed chemistry is described iin this study. One-step theoretical formulation of the induction parameter model IPM uses a theoretical reconstruction of the induction time database obtained from a detailed kinetics library. Non-dimensional induction time calculations is compared with that of detailed kinetics. The IPM was latter implemented to fluid dynamics code and applied for the numerical simulation of detonation wave propagation. The numerical results including the numerical smoked-foil record show the all the details of the detonation wave propagation characteristics at the cost around 1/100 of the detailed kinetics calculation.

• Journal of the Korean Society of Safety
• /
• v.31 no.4
• /
• pp.143-149
• /
• 2016

In this study, a series of numerical simulation of dam break flow was conducted using EFDC model, and input conditions including cell size, time step, and turbulent eddy viscosity were considered to analyze parameter sensitivity. In case of coarse mesh layout, the propagated length of the shock wave front was ${\Delta}_x$ longer than that of other mesh layouts, and the velocity results showed jagged edge, which can be cured by applying fine grid mesh. Turbulent eddy viscosity influenced magnitude of the maximum velocity passing through gate up to 20% and the cell Peclet number less than 2.0 ensured no numerical oscillations.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.453-458
• /
• 2004

This study proposes a new two-level condensation scheme for the construction of a reduced system. In the first step, the candidate area is selected for the construction of the reduced system by energy estimation in element-level. In the second step, primary degrees of freedom are selected by sequential elimination from the candidate degrees of freedom linked to the selected elements. Numerical examples demonstrate that the proposed method saves the computational cost effectively and provides a reduced system which predicts the eigenvalues accurately. Moreover, the well-constructed reduced system can present the reliable behavior of the structure under arbitrary dynamic loads comparing to that of global system. Time integration in a reduced system can save the computing time remarkably. Through a few numerical examples, the efficiency and reliability of the proposed scheme are verified.

• Journal of Mechanical Science and Technology
• /
• v.20 no.3
• /
• pp.376-381
• /
• 2006

In this paper, the quaternion, the dual Euler, and the direction cosine methods are numerically compared using a non-aerodynamic 6 degree-of-freedom rigid model at all-attitude angles of an aircraft. The dual Euler method turns out to be superior to the others in the applications because it shows better numerical accuracy, stability, and robustness in integration step sizes. The dual Euler method is affordably less efficient than the quaternion method in terms of computational cost. Numerical accuracy and stability, which allow larger integration step sizes, are more critical in modern real-time applications than computational efficiency because of today's increased computational power. If the quaternion method is required because of constraints in computation time, then a suppression mechanism should be provided for algebraic constraint errors which will eventually add computational burden.

• Journal of Korea Water Resources Association
• /
• v.44 no.5
• /
• pp.363-376
• /
• 2011

1-D hydrodynamic numerical models have been most widely used in the field of flood analysis. The model's input data are upstream/downstream boundaries, roughness coefficients, cross-sections, and so on, and computational distance step and time step are the most important factors in order to guarantee the computational accuracy, stability, and efficiency. In this study, a theoretical explanation is presented for the basis of the previous empirical selection criteria of cross-section's location; also, the estimation technique of computationally variable distance step is proposed to reflect the properties of flow at every computational time step. Combining this technique with 1-D unsteady numerical model, it was applied to two events of Teton dam failure flood and the Han River flood. The numerical experimental results demonstrate that the accuracy and stability is increased when used more interpolated cross-sections and show that the proposed technique of computationally variable distance step has the same order of accuracy with smaller numbers of cross-section than previous empirical selection criteria. The practical use of this technique will be possible to analyze the river floods with high efficiency as well as accuracy and stability.

• Structural Engineering and Mechanics
• /
• v.81 no.5
• /
• pp.591-606
• /
• 2022

When the friction pendulum system and shear keys work together to resist the ground motion, which inclined inputs (non 45°) to the bridge structure, the shear keys in XY direction will be sheared asynchronously, endowed the friction pendulum system with a violent curvilinear motion on the sliding surface during earthquakes. In view of this situation, firstly, this paper abandons the equivalent linearization model of friction and constructs a Spring-Coulomb friction plane isolation system with XY shear keys, and then makes a detailed mechanical analysis of the movement process of friction pendulum system, next, this paper establishes the mathematical model of structural time history response calculation by using the step-by-step integration method, finally, it compiles the corresponding computer program to realize the numerical calculation. The results show that the calculation method in this paper takes advantage of the characteristic that the friction force is always µmg, and creatively uses the "circle making method" to express the change process of the friction force and resultant force of the friction pendulum system in any calculation time step, which can effectively solve the temporal nonlinear action of the plane friction; Compared with the response obtained by the calculation method in this paper, the peak values of acceleration response and displacement response calculated by the unidirectional calculation model, which used in the traditional research of the friction pendulum system, are smaller, so the unidirectional calculation model is not safe.

• Nuclear Engineering and Technology
• /
• v.55 no.7
• /
• pp.2454-2465
• /
• 2023

Film boiling may lead to burnout of the heating element. Even though burnout does not occur, the heating element is subject to deformation because it is not sufficiently strong to withstand external loads. In particular, the ballooning and rupture of a tube under film boiling are important phenomena in the field of nuclear reactor safety. If the tube-type cladding of nuclear fuel ruptures owing to high internal pressure and thermal load, radioactive materials inside the cladding are released to the coolant. Therefore, predicting the ballooning and rupture is important. This study presents numerical simulations to predict the ballooning behavior and rupture time of a horizontal tube at high internal pressure under saturated film boiling. To do so, a multi-step coupled simulation of conjugated film boiling heat transfer and ballooning using creep model is adopted. The numerical methods and models are validated against experimental values. Two different nonuniform heat flux distributions and four different internal pressures are considered. The three-step simulation is enough to obtain a convergent result. However, the single-step simulation also successfully predicts the rupture time. This is because the film boiling heat transfer characteristics are slightly affected by the tube geometry related to creep ballooning.