• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.414-422
• /
• 2021

The aim of this work considers a second order point reactor kinetics model based on the P1 approximation of transport theory, called in this work as P1 point reactor model. The P1 point reactor model implicitly considers the time derivative of the neutron source which has not been thus considered previously. The inverse method to calculate the reactivity in nuclear reactors -chosen because its high accuracy- Matrix Formulation. The numerical results shown that the Matrix Formulation for the reactivity estimation constitutes a method with insignificant calculation errors.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.1
• /
• pp.267-285
• /
• 2014

An upstream scheme based on the pseudostress-velocity mixed formulation is studied to solve convection-dominated Oseen equations. Lagrange multipliers are introduced to treat the trace-free constraint and the lowest order Raviart-Thomas finite element space on rectangular mesh is used. Error analysis for several quantities of interest is given. Particularly, first-order convergence in $L^2$ norm for the velocity is proved. Finally, numerical experiments for various cases are presented to show the efficiency of this method.

• Structural Engineering and Mechanics
• /
• v.21 no.5
• /
• pp.571-590
• /
• 2005

An 8 node assumed stress hexahedral element with rotational degrees of freedom is proposed for static and free vibration analyses. The element formulation is based directly on an 8-node element. This direct formulation requires fewer computations than a similar element that is derived from an internal 20-node element in which the midside degrees of freedom are eliminated by expressing them in terms of displacements and rotations at corner nodes. The formulation is based on Hellinger-Reissner variational principle. Numerical examples are presented to show the validity and efficiency of the present element for static and free vibration analysis.

• Journal of Korean Society for Quality Management
• /
• v.21 no.1
• /
• pp.121-135
• /
• 1993

The derived expressions and computations of the system reliability in the consecutive k out of n failure structure with sink-source pole are discussed and a new simplified formulation to compute the reliability of consecutive k out of $n^*$ failure system is presented, and the relations with system reliability when the number of first consecutive falied components from source is greater than k are presented in this paper. The new simplified formulation is illustrated in the numerical experiments.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.10
• /
• pp.55-65
• /
• 2014

Recent researches on debris flow is focused on understanding its movement mechanism and building a numerical simulator to predict its behavior. However, previous simulators emulating fluid-like debris flow have limitations in numerical stability, geometric modeling and application of various boundary conditions. In this study, depth integration is applied to continuity equation and force equilibrium for debris flow. Thickness of sediment, and average velocities in x and y flow direction are chosen for main variables in the analysis, which improve numerical stability in the area with zero thickness. Petrov-Galerkin formulation uses a discontinuous test function of the weighted matrix from DG scheme. Presented mechanical constitutive model combines fluid and granular behaviors for debris flow. Effects on slope angle, inducing debris height, and bottom friction resistance are investigated for a simple slope. Numerical results also show the effect of embankment at the bottom of the slope. Developed numerical simulator can assess various risk factors for the expected area of debris flow, and facilitate embankment design in order to minimize damage.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Structural Engineering and Mechanics
• /
• v.68 no.1
• /
• pp.1-15
• /
• 2018

In ocean industry, free surface type ART (Anti Roll tank) system has been widely used to suppress the roll motion of floating structures. In those, various obstacles have been devised to obtain the sufficient damping and to enhance the controllability of freely rushing water inside the tank. Most of previous researches have paid on the development of simple mathematical formula for coupled ship-ARTs analysis although other numerical and experimental approaches exist. Little attention has been focused on the use of 3D panel method for preliminary design of free surface type ART despite its advantages in computational time and general capacity for hydrodynamic damping estimation. This study aims at developing a potential theory based hydrodynamic code for the analysis of floating structure with baffled ARTs. The sloshing in baffled tanks is modeled through the linear potential theory with FE discretization and it coupled with hydrodynamic equations of floating structures discretized by BEM and FEM, resulting in direct coupled FE-BE formulation. The general capacity of proposed formulation is emphasized through the coupled hydrodynamic analysis of floating structure and sloshing inside baffled ARTs. In addition, the numerical methods for natural sloshing frequency tuning and estimation of hydrodynamic damping ratio of liquid sloshing in baffled tanks undergoing wave exiting loads are developed through the proposed formulation. In numerical examples, effects of natural frequency tuning and baffle ratios on the maximum and significant roll motions are investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.

• Journal of Korean Society of Steel Construction
• /
• v.10 no.2 s.35
• /
• pp.317-326
• /
• 1998

This paper describes a reduced finite element formulation for nonlinear transient heat transfer analysis based on Lanczos Algorithm. In the proposed reduced formulation all material nonlinearities of irradiation boundary element are included using the pseudo force method and numerical time integration of the reduced formulation is conducted by Galerkin method. The results of numerical examples demonstrate the applicability and the accuracy of the proposed method for the nonlinear transient heat transfer analysis.

• Smart Structures and Systems
• /
• v.29 no.2
• /
• pp.279-286
• /
• 2022

In this article, a novel propagation formulation of Rayleigh waves in a compressible isotropic half-space with impedance boundary condition is proposed by embedding the normal stress. In a two-dimensional case, it is assumed that a design boundary is free of normal traction and a shear traction depends on linearly a normal component of displacements multiplied by frequencies. Therefore, impedance boundary conditions affect the normal stress, where the impedance parameters correspond to dimensions of stresses over velocity. On the other hand, vanished impedance values are traction-free boundary conditions. The main purpose of this article is to present theoretically the existence and uniqueness of a Rayleigh wave formulation relying on secular equation's mathematical analyses. Its velocity varies along with the impedance parameters. Moreover, numerical experiments with different values for the velocity of Rayleigh waves are carried out. The present Rayleigh waves study is a fundamental step in analyzing the cause and effect of physical states such as building or structure damages resulting from natural dynamics. The results of the study generate a basic design formulation theory to test the effects of Rayleigh waves affecting structures when an earthquake occurs. The presence and uniqueness of the proposed formulation is verified by mutual comparisons of several numerical examples.