• Proceedings of the KIEE Conference
• /
• 1988.07a
• /
• pp.950-953
• /
• 1988

The equations of motion for linearly elastic bodies undergoing large displacement motion are derived. This produces a set of equations which are efficient to numerically integrate. The equations for the elastic bodies are formulated and simplified to provide as much efficiency as possible in their numerical solution. A futher efficiency is obtained through the use of floating reference frame. The equation are presented in two forms for numerical integration. 1) Explicit numerical integration 2) Implicit numerical integration. In this paper, there was used the numerical integration. The implicit numerical integration is extended to solved second order equation, futher reducing the numerical effort required. The formulation given is seen to be occulate and is expected to be efficient for many types of problems.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.173-195
• /
• 2021

A new implicit discontinuous Galerkin spectral element method (DGSEM) based on the first order hyperbolic system(FOHS) is presented for solving elliptic type partial different equations, such as the Poisson problems. By utilizing the idea of hyperbolic formulation of Nishikawa[1], the original Poisson equation was reformulated in the first-order hyperbolic system. Such hyperbolic system is solved implicitly by the collocation type DGSEM. The steady state solution in pseudo-time, which is the solution of the original Poisson problem, was obtained by the implicit solution of the global linear system. The optimal polynomial orders of 𝒪(𝒽^𝑝+1)) are obtained for both the solution and gradient variables from the test cases in 1D and 2D regular grids. Spectral accuracy of the solution and gradient variables are confirmed from all test cases of using the uniform grids in 2D.

• Journal of the Korean Mathematical Society
• /
• v.59 no.2
• /
• pp.279-298
• /
• 2022

In this paper, the existence and uniqueness for the global solution of neutral stochastic functional differential equation is investigated under the locally Lipschitz condition and the contractive condition. The implicit iterative methodology and the Lyapunov-Razumikhin theorem are used. The stability analysis for such equations is also applied. One numerical example is provided to illustrate the effectiveness of the theoretical results obtained.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.2
• /
• pp.196-202
• /
• 2012

In this study, numerical analysis of viscous flows is carried out based on the unstructured grid. There exist some difficulties in expressing and computing numerical derivatives on the unstructured grid due to lack of the structured characteristics. The general computer algorithms are developed to perform numerical derivatives easily and extended to be applicable to various geometries composed of hybrid meshes. And the optimal method of strongly implicit procedure is newly contrived to accelerate the rate of convergence in solving the pressure Poisson equation. To verify numerical schemes, the driven cavity problems of 2 and 3 dimension are simulated. The numerical results are compared with others and our numerical schemes are shown to be valid.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1624-1629
• /
• 2008

The governing equations in generalized curvilinear coordinates for a 3D pulsatile flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. The computational technique implements the implicit approximate factorization method of the Beam and Warming method (1978), which is the extension of the Alternate Direction Implicit (ADI) method. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Rindt & Steenhoven, 1991).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.8 s.227
• /
• pp.1140-1151
• /
• 2004

This paper reports on the finite elements analysis for die compaction process of cemented carbide tool parts. Experimental data were obtained under die compaction and triaxial compression with various loading conditions. The elastoplastic constitutive equations based on the yield function of Shima and Oyane were implemented into an explicit finite element program (ABAQUS/Explicit) and implicit finite element program (PMsolver/Compaction-3D) to simulate compaction response of cemented carbide powder during die compaction. For simulation of die compaction, the material parameters for Shima and Oyane model were obtained by uniaxial die compaction test. Explicit finite element results were compared with implicit results for cemented carbide powder.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.90-97
• /
• 1999

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.3
• /
• pp.138-155
• /
• 2022

To solve the initial value problem we present a new single-step implicit method based on the Euler method. We prove that the proposed method has convergence order 2. In practice, numerical results of the proposed method for some selected examples show an error tendency similar to the second-order Taylor method. It can also be found that this method is useful for stiff initial value problems, even when a small number of nodes are used. In addition, we extend the proposed method by using weighted averages with a parameter and show that its convergence order becomes 2 for the parameter near $\frac{1}{2}$. Moreover, it can be seen that the extended method with properly selected values of the parameter improves the approximation error more significantly.