• Wind and Structures
• /
• 제20권3호
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Journal of electromagnetic engineering and science
• /
• 제3권1호
• /
• pp.39-44
• /
• 2003

An unconditionally stable compact 2D Alternating-Direction Implicit (ADI) FDTD method for calculating dispersion characteristics of waveguide structures is proposed. The numerical stability and numerical dispersion relation of the proposed method are also presented and discussed. Numerical wavelengths for the dominant and higher order modes in a hollow waveguide are obtained from numerical simulations and compared with those from the analytical dispersion relation. The numerical results show that the proposed scheme has the potential to successfully analyze a class of waveguides having locally fine geometry with reduced numerical costs.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2007년도 정기총회 및 학술발표대회
• /
• pp.181-185
• /
• 2007

A free-surface correction(FSC) method is presented to solve the 3-D shallow water equations. Using the mode splitting process, FSC method can simulate shallow water flows under the hydrostatic assumption. For the hydrostatic pressure calculation, the momentum equations are firstly discretized using a semi-implicit scheme over the vertical direction leading to the tri-diagonal matrix systems. A semi-implicit scheme has been adopted to reduce the numerical instability caused by relatively small vertical length scale compare to horizontal one. and, as the free surface correction step the final horizontal velocity fields are corrected after the final surface elevations are obtained. Finally, the vertical final velocity fields can be calculated from the continuity equation. The numerical model is applied to the calculation of the simulation of flow fields in a rectangular open channel with the tidal influence. The comparisons with the analytical solutions show overall good agreements between the numerical results and analytical solutions.

• Geomechanics and Engineering
• /
• 제25권4호
• /
• pp.267-273
• /
• 2021

In this work, we compare the analytical solutions with the numerical solutions for photothermal interactions in semiconductor medium containing cylindrical cavity. This paper is devoted to a study of the photothermal interactions in semiconductor medium in the context of the coupled photo-thermal model. The basic equations are formulated in the domain of Laplace transform and the eigenvalue scheme are used to get the analytical solutions. The numerical solution is obtained by using the implicit finite difference method (IFDM). A comparison between the analytical solution and the numerical solutions are obtained. It is found that the implicit finite difference method (IFDM) is applicable, simple and efficient for such problems.

• Structural Engineering and Mechanics
• /
• 제75권5호
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.

• 대한조선학회논문집
• /
• 제35권1호
• /
• pp.1-14
• /
• 1998

2차원 천수방정식에 기초한 조류계산 수치모델을 개발하였다. 복잡한 해안선의 입력을 단순화할 수 있는 직각격자를 사용하였으며, 방정식은 효율이 뛰어난 것으로 알려져 있는 인수분해음해법(Implicit Factorization Method)에 의하여 수치적으로 풀었다. 대류항은 상류차분법으로 처리하여 수치소산이 자동적으로 들어가도록 하였으며, 따라서 수치불안정에 의한 발산을 억제할 수 있었다. 모델의 검증을 위하여 선형화된 천수방정식을 풀고 계산결과를 해석해와 비교하였으며 잘 일치하고 있음을 확인하였다. 수치모델의 응용 예로서 인천항부근의 경인운하 개발예정인 수로에서의 조류 계산결과를 보였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.7-12
• /
• 1998

Optimized high-order compact(OHOC) schemes were proposed, which have high spatial order of truncation and resolution to simulate the aeroacoustic problems due to unsteady compressible flows. Generally, numerical schemes are categorized explicit or implicit by time-marching method. In this research, OHOC differences which were developed with explicit time-marching method is used to have implicit formulation and the implicit OHOC differences result in block hepta-diagonal matrix. This paper presents the comparisons between the explicit and implicit OHOC schemes with a second order accuracy of time in the 1-d linear wave convection problem, and between the explicit OHOC scheme of 4th-order accuracy in time and the implicit OHOC scheme of 1st-order accuracy in tine for the 1-d nonlinear wave convection problem. With these comparisons, the characteristics of implicit OHOC scheme are shown in the point of CFL number.

• Earthquakes and Structures
• /
• 제7권2호
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.290-297
• /
• 2009

가압 경수로의 주요 기기에서 발생할 수 있는 과도 2상 유동(Two-phase flow) 현상에 대한 해석을 수행하기 위해 원자로 기기 열수력 해석 코드를 개발 중에 있다. 개발 중인 기기 열수력 해석 코드는 지배 방정식으로 Two-phase, three-field model을 사용하고 있으며, 복잡한 기하학적 형상의 원자로 기기를 모사하기 위해 비정렬 격자계(Unstructured grid)를 활용하고 있다. 수치해석 기법으로는, 원자로 계통 해석코드 RELAP5가 사용 중이며 대부분의 원자로 내 2상 유동 조건에서 안정적이며 정확하다고 알려진 Semi-implicit 방법을 적용하였다. 그러나 기존의 Semi-implicit 방법은 1차원, 엇갈림격자(Staggered grid)에 대해 개발되었기 때문에, 이를 다차원, 비정렬, 비엇갈림 격자(Non-staggered grid)에 적용하기 위해 기존의 Semi-implicit 방법을 수정하였다. 본 논문에서는 Semi-implicit 방법의 대류항을 이차정확도를 갖도록 확장하였으며, 이차정확도에 의한 수치확산의 감소를 평가하기 위해 수행된 수치시험의 결과를 기술하였다. 이차정확도 및 일차정확도로 계산된 값을 해석해 또는 격자 수렴성 시험을 통해 평가해 본 결과, 이차정확도 계산시 수치 확산의 감소 확인하였다.