• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Korean Journal of Mathematics
• /
• v.13 no.2
• /
• pp.161-176
• /
• 2005

We consider a model of population dynamics whose mortality function is unbounded. We note that the regularity of the solution depends on the growth rate of the mortality near the maximum age. We propose Gauss-Legendre methods along the characteristics to approximate the solution when the solution is smooth enough. It is proven that the scheme is convergent at fourth-order rate in the maximum norm. We also propose discontinuous Galerkin finite element methods to approximate the solution which is not smooth enough. The stability of the method is discussed. Several numerical examples are presented.

• Transactions of Materials Processing
• /
• v.9 no.4
• /
• pp.395-403
• /
• 2000

Formation of Shear Band under the adiabatic condition is widely observed In the engineering materials during rapidly forming process lot a thermally rate-dependent material. The shear band stems from evolution of a narrow region in which an intensive plastic flow occurs. The shear band often plays a role of a precursor of the ductile fracture during a forming process. The objective of this study is to investigate the localization behavior using numerical method. In this work, the implicit finite difference scheme is employed due to the ease of convergence and the numerical stability It is noted that physical and mechanical properties of materials determine how the shear band is formed and then localized. Material properties can be characterized with inertia number dissipation number and diffusion number. It is observed that the dimensionless numbers effect on localization. Using a parametric study, comparison was made between CRS-1018 steel with WHA (tungsten heavy alloy). The deformation behavior of material in this study include an isotropic hardening as well as thermal softening. Moreover, this study suggests that a kinematic hardening constitutive relation be required to predict a more accurate strain level at a shear band.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.9
• /
• pp.3583-3589
• /
• 2010

The semi-implicit scheme proposed by Backhaus is introduced to solve two-dimensional shallow water equation. This mothod is applied for the numerical model solving surface elevation and velocity field of Geum River estuary. For the verification of the method, numerical solutions by this model are compared with ones by Heap's well known explicit model. Solutions of two models resemble each other. The time-step chosen for the semi-implicit scheme turned out to be 3 to 6 times longer than explicit model depending on the stringent CFL criterion. The computation time could be reduced at least 50%. It was proved that this scheme is easy to handle dry banks which can be seen in Geum River estuary and numerical stability is obtained for long time computation.

• Structural Engineering and Mechanics
• /
• v.85 no.4
• /
• pp.485-500
• /
• 2023

The non-linear static analysis of reinforced concrete (RC) structures using the three-dimensional (3D) finite element method is a time-consuming and challenging task. Moreover, this type of analysis encounters numerical problems such as the lack of convergence of results in the stages of growth and propagation of cracks in the structure. The time integration analysis along with the mass scaling (MS) technique is usually used to overcome these limitations. Despite the use of this method in the 3D finite element analysis of RC structures, a comprehensive study has not been conducted so far to assess the effects of the MS method on the accuracy of results. This study aims to evaluate the accuracy of the MS method in the non-linear quasi-static finite element analysis of RC structures. To this aim, different types of RC structures were simulated using the finite element approach based on the implicit time integration method and the mass scaling technique. The influences of effective parameters of the MS method (i.e., the allowable values of increase in the mass of the RC structure, the relationship between the duration of the applied load and fundamental vibration period of the RC structure, and the pattern of applied loads) on the accuracy of the simulated results were investigated. The accuracy of numerical simulation results has been evaluated through comparison with existing experimental data. The results of this study show that the achievement of accurate structural responses in the implicit time integration analyses using the MS method involves the appropriate selection of the effective parameters of the MS method.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.163-170
• /
• 2000

The distinct element method has been effectively applied to the analysis of stability and behavior of jointed rock masses. In this paper the modelling characteristics of different types of distinct element model were investigated. Arch tunnel examples were chosen to compare the calculation results of two computer codes, NURBM and CBLOCK, where the former is based on implicit algorithm, and the other on explicit one. CBLOCK calculations show that joint properties are very important parameters in the stability analysis and that the joint stiffness ratio associated with joint configuration could be used as an indicator, whereas NURBM differ from that. Some other disagreements were also identified.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.410-417
• /
• 2000

The distinct element method has been effectively applied to the analysis of stability and behavior of jointed rock masses. In this paper the modelling characteristics of different types of distinct element model were investigated. Arch tunnel examples were chosen to compare the calculation results of two computer codes, NURBM and CBLOCK, where the former is based on implicit algorithm, and the other on explicit one. CBLOCK calculations show that joint properties are very important parameters in the stability analysis and that the joint stiffness ratio associated with joint configuration could be used as an indicator, whereas NURBM differ from that. Some other disagreements were also identified.

• Proceedings of the IEEK Conference
• /
• 2002.07a
• /
• pp.137-140
• /
• 2002

In this paper, a new algorithm for two-dimensional (2-D) finite-difference time-domain (FDTD) method is presented. By this new method, the maximum time step size can be increased over the Courant-Friedrich-Levy (CFL) condition restraint. This new algorithm is adapted from an Alternating-Direction Implicit FDTD (ADI-FDTD) method. However, unlike the ADI-FDTD algorithm. the alternation is performed with respect to the blocks of fields rather than with respect to each respective coordinate direction. Moreover. this method can be efficiently simulated with parallel computation. and it is more efficient than the conventional FDTD method in terms of CPU time. Numerical formulations are shown and simulation results are presented to demonstrate the effectiveness and efficiency of our proposed method.

• 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.