• Journal of applied mathematics & informatics
• /
• 제29권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.

• 전기학회논문지
• /
• 제61권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
• /
• 제14권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년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• 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.

• 수산해양기술연구
• /
• 제28권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
• /
• 제9권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• 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.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2005년도 제36회 하계학술대회 논문집 C
• /
• 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
• /
• 제87권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.

• 전산구조공학
• /
• 제6권2호
• /
• pp.103-114
• /
• 1993

유전적분형 물성방정식에 근거한 선형 점탄성문제의 효율적인 수치해석을 위해서 새로운 유한요소해법을 공식화하였다. 각 시간구간에서 변수변화를 선형적으로 가정하고 유전적분의 계산을 매우 효율적으로 처리하였다. 그 결과 과거의 해석법에 비하여 수치정확도 및 경제성에서 큰 향상을 얻었다.