• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1129-1141
• /
• 2011

In this paper, a class of large time-stepping method based on the finite element discretization for the Cahn-Hilliard equation with the Neumann boundary conditions is developed. The equation is discretized by finite element method in space and semi-implicit schemes in time. For the first order fully discrete scheme, convergence property is investigated by using finite element analysis. Numerical experiment is presented, which demonstrates the effectiveness of the large time-stepping approaches.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.61 no.8
• /
• pp.1123-1129
• /
• 2012

This paper presents the fast steady state analysis using time-stepping finite element method for a high power three-phase transformer. The high power transformer spends huge computational cost of the time-stepping finite element method. It is because that the high power transformer requires a lot of time to reach steady state by its large inductance component. In order to reduce computational cost, in this paper, the adaptive time-step control algorithm combined with the embedded 2nd 4th singly diagonally implicit Runge-Kutta method and the analysis strategy using variation of the winding resistance are studied, and their numerical results are compared with those from the typical time-stepping finite element method.

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1269-1289
• /
• 2014

With the sub-stepping technique, the numerical analysis in real-time dynamic hybrid testing is split into the response analysis and signal generation tasks. Two target computers that operate in real-time may be assigned to implement these two tasks, respectively, for fully extending the simulation scale of the numerical substructure. In this case, the integration time-step of solving the dynamic response of the numerical substructure can be dozens of times bigger than the sampling time-step of the controller. The time delay between the real and desired feedback forces becomes more striking, which challenges the well-developed delay compensation methods in real-time dynamic hybrid testing. This paper focuses on displacement prediction and force correction for delay compensation in the real-time dynamic hybrid testing with a large integration time-step. A new displacement prediction scheme is proposed based on recently-developed explicit integration algorithms and compared with several commonly-used prediction procedures. The evaluation of its prediction accuracy is carried out theoretically, numerically and experimentally. Results indicate that the accuracy and effectiveness of the proposed prediction method are of significance.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.181-182
• /
• 2003

A preconditioned numerical method for gas-liquid to-phase flow is applied to solve cavitating flow. The present method employs a density based finite-difference method of dual time-stepping integration procedure and Roe's flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristics at low Mach number. By this method, two-dimensional internal flows through a venturi tuve and decelerating cascades are computed and discussed.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.28 no.1
• /
• pp.79-92
• /
• 1992

This is analyzed using the finite element method which is appling excellent isoparametric curve element in the aspect of large usages of dynamic responses in which is regarding geometric and material nonlinear of a large scale shell structure of an airplane, a submarine, a ship, and an ocean structure. The solution of dynamic equations is got by direct integration method using time-stepping procedure and regarding Central Difference Method of the both solutions. But because formal matrix factorization is not necessary in each time step and it does not take less time to compute relatively, this method must be regarded very few time steps on the condition. Axisymmatric shell problems are inspected using 8 node Isoparametric element in this paper. Partial axisymmatric spherical shell is used as a model to analyze axisymmatric nonlinear dynamic behavior regarding. Total Lagrangian formulation in geometric nonlinear behavior and elastio-viscoplastic in material nonlinear behavior.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.390-403
• /
• 2017

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1903-1908
• /
• 2004

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

• Proceedings of the KIEE Conference
• /
• 2005.07c
• /
• pp.2315-2317
• /
• 2005

In this paper, an iterative method to model the electromagnetic heating of electrically large lossy dielectrics is presented. Frequency domain finite element (FEM) solutions of the wave equation are determined for the lossy inhomogeneous dielectric as the material properties are change with temperature and time. The power absorbed from microwave losses is applied to a finite element time domain (FETD) calculation of the heat diffusion equation. Time steps appropriate for updating the piecewise material properties in the wave equation and the time stepping of the heat equation are presented. The effects of preheating and source frequency are investigated.

• Structural Engineering and Mechanics
• /
• v.87 no.4
• /
• pp.363-373
• /
• 2023

Since the first family of structure-dependent methods can simultaneously integrate unconditional stability and explicit formulation in addition to second order accuracy, it is very computationally efficient for solving inertial problems except for adopting auto time-stepping techniques due to no nonlinear iterations. However, an unusual stability property is first found herein since its unconditional stability interval is drastically different for zero and nonzero damping. In fact, instability might occur for solving a damped stiffness hardening system while an accurate result can be obtained for the corresponding undamped stiffness hardening system. A technique of using a stability factor is applied to overcome this difficulty. It can be applied to magnify an unconditional stability interval. After introducing this stability factor, the formulation of this family of structure-dependent methods is changed accordingly and thus its numerical properties must be re-evaluated. In summary, a large stability factor can result in a large unconditional stability interval but also lead to a large relative period error. As a consequence, a stability factor must be appropriately chosen to have a desired unconditional stability interval in addition to an acceptable period distortion.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.