• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.458-467
• /
• 2016

In this paper, a strongly coupled method for investigating the interaction between a cable and tow-fish is presented. The nodal position finite element method was utilized to analyze the nonlinear cable dynamics, and 6DOF equations of motion were employed to describe the 3D rigid body motion of the tow-fish. Combining cable and tow-fish systems into a single formulation allowed the two nonlinear systems to be strongly coupled into a unified nonlinear system. This strongly coupled system was numerically integrated in the time domain using a predictor/multi-corrector Newmark algorithm. To demonstrate the validity, efficacy, and applicability of the current approach, two different scenarios (virtual and sea trial) were simulated, and the simulation results were validated using the physical plausibility and the sea trial test.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.388-399
• /
• 2016

A motion analysis of an underwater towed cable is a complex task due to its nonlinear nature of the problem. The major source of the nonlinearity of the underwater cable analysis is that the motion of the cable involves large rigid-body motion. This large rigid-body motion makes difficult to use standard displacement-based finite element method. In this paper, the authors apply recently developed nodal position-based finite element method which can deal with the geometric nonlinearity due to the large rigid-body motion. In order to enhance the stability of the large-scale nonlinear cable motion simulation, an efficient time-integration scheme is proposed, namely predictor/multi-corrector Newmark scheme. Three different predictors are introduced, and the best predictor in terms of stability and robustness for impulsive cable motion analysis is proposed. As a result, the nonlinear motion of underwater cable is predicted in a very efficient manner compared to the classical finite element of finite difference methods. The efficacy of the method is demonstrated with several test cases, involving static and dynamic motion of a single cable element, and also under water towed cable composed of multiple cable elements.

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.54-63
• /
• 2002

In recent years, mechanical systems such as high speed vehicles and railway trains moving on elastic beam structures have become a very important issue to consider. In this paper, a general approach, which can predict the dynamic behavior of a constrained mechanical system moving on a flexible beam structure, is proposed. Various supporting conditions for the foundation support are considered for the elastic beam structure. The elastic structure is assumed to be a non-uniform and linear Bernoulli-Euler beam with a proportional damping effect. Combined differential-algebraic equation of motion is derived using the multi-body dynamics theory and the finite element method. The proposed equations of motion can be solved numerically using the generalized coordinate partitioning method and predictor-corrector algorithm, which is an implicit multi-step integration method.

• International Journal of Concrete Structures and Materials
• /
• v.2 no.2
• /
• pp.123-136
• /
• 2008

A new plastic-damage constitutive model applicable to lightweight concrete (LWC) and normal weight concrete (NWC) is proposed in this paper based on both continuum damage mechanics and plasticity theories. Two damage variables are used to represent tensile and compressive damage independently. The effective stress is computed in the Drucker-Prager multi-surface plasticity framework. The stress is then computed by multiplication of the damaged part and the effective part. The proposed model is coded as a user material subroutine and incorporated in a finite element analysis software. The constitutive integration algorithm is implemented by adopting the operator split involving elastic predictor, plastic corrector and damage corrector. The numerical study shows that the algorithm is efficient and robust in the finite element analysis. Experimental investigation is conducted to verify the proposed model involving both static and dynamic tests. The very good agreement between the numerical results and experimental results demonstrates the capability of the proposed model to capture the behaviors of LWC and NWC structures for static and impact loading.

• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2699-2708
• /
• 2020

Helical-coil steam generator (HCSG) technology is a major design candidate for small modular reactors due to its compactness and capability to produce superheated steam with high generation efficiency. In this paper, we investigate the feasibility of the passively autonomous power maneuvering by coupling the 3-D transient multi-physics of a soluble-boron-free (SBF) core with a time-dependent HCSG model. The predictor corrector quasi-static method was used to reduce the cost of the transient 3-D neutronic solution. In the numerical system simulations, the feedwater flow rate to the secondary of the HCSGs is adjusted to extract the demanded power from the primary loop. This varies the coolant temperature at the inlet of the SBF core, which governs the passively autonomous power maneuvering due to the strongly negative coolant reactivity feedback. Here, we simulate a 100-50-100 load-follow operation with a 5%/minute power ramping speed to investigate the feasibility of the passively autonomous load-follow in a 450 MW_th SBF PWR. In addition, the passively autonomous frequency control operation is investigated. The various system models are coupled, and they are solved by an in-house Fortran-95 code. The results of this work demonstrate constant steam temperature in the secondary side and limited variation of the primary coolant temperature. Meanwhile, the variations of the core axial shape index and the core power peaking are sufficiently small.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.145-151
• /
• 1999

An Interactive grid generation program(KGRID) with graphical user interface(GUI) has been improved. KGRID works on the UNLX environment and GUI has been implemented with OSF/Motif and X Toolkit and the graphics language is Open GL for visualization of the 3D objects. It supports more convenient user environment to generate 2D and 3D multi-block structured grid systems. It provides various useful field grid generation methods, which are the algebraic methods, the elliptic partial differential equations method and the predictor-corrector method. It also supports 3D surface grid generation with NURBS(Non-Uniform Rational B-Spline) and various stretching functions to control grid points distribution on curves and surfaces. And some menus are added to perform flexible management, for the objects. We generated surface and field grid system about full aircraft configuration using KGRID. The performance and stability of the KGRID is verified through the generation of the grid system about a complex shape.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.11
• /
• pp.165-175
• /
• 2000

In recent years, it becomes a very important issue to consider the mechanical systems such as high-speed vehicles and railway trains moving on elastic beam structures. In this paper, a general approach, which can predict the dynamic behavior of constrained mechanical system and elastic beam structure, is proposed. Also, various supporting conditions of a foundation support are considered for the elastic beam structures. The elastic structure is assumed to be a nonuniform and linear Bernoulli-Euler beam with proportional damping effect. Combined Differential-Algebraic Equations of motion are derived using multibody dynamics theory and Finite Element Method. The proposed equations of motion can be solved numerically using generalizd coordinate partitioning method and Predictor-Corrector algorithm, which is an implicit multi-step integration method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.5
• /
• pp.361-366
• /
• 2017

As the recent development of computing architecture and application software technology, real world simulation, which is the ultimate destination of computer simulation, is emerging as a practical issue in several research sectors. In this paper, metal plate motion in a square shock tube for small time interval was calculated using a supercomputing-based fluid-structure-combustion multi-physics simulation tool called Illinois Rocstar, developed in a US national R amp; D program at the University of Illinois. Afterwards, the simulation results were compared with those from experiments. The coupled solvers for unsteady compressible fluid dynamics and for structural analysis were based on the finite volume structured grid system and the large deformation linear elastic model, respectively. In addition, a strong correlation between calculation and experiment was shown, probably because of the predictor-corrector time-integration scheme framework. In the future, additional validation studies and code improvements for higher accuracy will be conducted to obtain a reliable open-source software research tool.