• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.53 no.8
• /
• pp.433-442
• /
• 2004

Various models and various boundary conditions have been suggested for fluid transport simulations of high density plasma discharges such as the inductively coupled plasma discharge. In this work, we compare the various models using one-dimensional simulations based on the FDM(finite difference method), the upwind scheme, the power-law scheme, and the dielectric relaxation scheme［l］ Comparing the exactness, the numerical stability and the efficiency of the various models. the most adoptable model is suggested.

• Communications of the Korean Mathematical Society
• /
• v.37 no.3
• /
• pp.939-955
• /
• 2022

In this paper, we study a first-order non-linear singularly perturbed Volterra integro-differential equation (SPVIDE). We discretize the problem by a uniform difference scheme on a Bakhvalov-Shishkin mesh. The scheme is constructed by the method of integral identities with exponential basis functions and integral terms are handled with interpolating quadrature rules with remainder terms. An effective quasi-linearization technique is employed for the algorithm. We establish the error estimates and demonstrate that the scheme on Bakhvalov-Shishkin mesh is O(N^-1) uniformly convergent, where N is the mesh parameter. The numerical results on a couple of examples are also provided to confirm the theoretical analysis.

• Nuclear Engineering and Technology
• /
• v.55 no.11
• /
• pp.4228-4237
• /
• 2023

The Improved Deterministic Truncation of Monte Carlo (iDTMC) is a powerful acceleration and variance reduction scheme in the Monte Carlo analysis. The concept of the iDTMC method and correlated sampling-based real variance estimation are briefly introduced. Moreover, the application of the iterative scheme to the correlated sampling is discussed. The iDTMC method is utilized in a 3-dimensional small modular reactor (SMR) model problem. The real variances of burnup-dependent criticality and power distribution are evaluated and compared with the ones obtained from 30 independent iDTMC calculations. The impact of the inactive cycles on the correlated sampling is also evaluated to investigate the consistency of the correlated sample scheme. In addition, numerical performances and sensitivity analysis on the real variance estimation are performed in view of the figure of merit of the iDTMC method. The numerical results show that the correlated sampling accurately estimates the real variances with high computational efficiencies.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.120-126
• /
• 1991

An expanded scheme of a direct numerical solution method for solving the Navier-Stokes equation considering the modified boundary conditions for gas lubrication with ultra low clearance at high .LAMBDA. region is presented. Many examples are calculated by this scheme and their results are compared to the previous solutions using P$^{2}$H$^{[-992]}$ . This scheme has the advantages of fast calculation time and stable convergence in high .LAMBDA. region, and gives very good results in the case of fluid film thickness discontinuity.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.17 no.2 s.105
• /
• pp.165-170
• /
• 2006

A two dimensional(2-D) finite-difference time-domain(FDTD) method based on a novel finite difference scheme is developed to eliminate the numerical dispersion errors. In this paper, numerical dispersion and stability analysis of the new scheme are given, which show that the proposed method is nearly dispersionless, and stable for a larger time step than the standard FDTD method.

• International Journal of Aeronautical and Space Sciences
• /
• v.11 no.4
• /
• pp.319-325
• /
• 2010

For decades, researchers have rigorously studied the characteristics of flow traveling around blunt objects in order to gain greater understanding of the flow around aircraft, vehicles or vessels. Many different types of flow exist, such as boundary layer flow, flow separation, laminar and turbulent flow, vortex and vortex shedding; such types are especially observed around circular cylinders. Vortex shedding around a circular cylinder exhibits a two-dimensional flow structure possessing a Reynolds number within the range of 47 and 180. As the Reynolds number increases, the Karman vortex changes into a three-dimensional flow structure. In this paper, a numerical analysis was performed examining the flow and aero-acoustic field characteristics around a circular cylinder using an optimized high-order compact scheme, which is a high order scheme. The analysis was conducted with a Reynolds number ranging between 300 and 1,000, which belongs to B-mode flow around a circular cylinder. For a B-mode Reynolds number, a proper spanwise length is analyzed in order to obtain the characteristics of three-dimensional flow. The numerical results of the Strouhal number as well as the lift and drag coefficients according to Reynolds numbers are coincident with the other experimental results. Basic research has been conducted studying the effects an unstable three-dimensional wake flow on an aero-acoustic field.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.12
• /
• pp.1624-1634
• /
• 1997

A numerical study on two-dimensional unsteady transonic cascade flow has been performed by adopting dual time stepping and the k-.omega. turbulence model. An explicit 4 stage Runge-Kutta scheme for the compressible Navier-Stokes equations and an implicit Gauss-Seidel iteration scheme for the k-.omega. turbulence model are proposed for fictitious time stepping. This mixed time stepping scheme ensures the stability of numerical computation and exhibits a good convergence property with less computation time. Typical steady-state convergence accelerating schemes such as local time stepping, residual smoothing and multigrid combined with dual time stepping shows good convergence properties. Numerical results are presented for unsteady laminar flow past a cylinder and turbulent shock buffeting problem for bicircular arc cascade flow is discussed.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.179-182
• /
• 2004

This paper deals with a new mathematical formulation of nonlinear wave profile based on Banach fixed point theorem. As application of the formulation and its solution procedure, some numerical solutions was presented in this paper and nonlinear equation was derived. Also we introduce a new operator for iteration and getting solution. A numerical study was accomplished with Stokes' first-order solution and iteration scheme, and then we can know the nonlinear characteristic of Stokes' high-order solution. That is, using only Stokes' first-oder(linear) velocity potential and an initial guess of wave profile, it is possible to realize the corresponding high-oder Stokian wave profile with tile new numerical scheme which is the method of iteration. We proved the mathematical convergence of tile proposed scheme. The nonlinear strategy of iterations has very fast convergence rate, that is, only about 6-10 iterations arc required to obtain a numerically converged solution.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.551-555
• /
• 2006

In this study, new dispersion-correction terms are added to a leap-frog finite difference scheme for the linear shallow-water equations with the purpose of considering dispersion effects of the linear Boussinesq equations for propagation of tsunamis. The numerical model developed in this study is tested to the problem that the initial free surface displacement is a Gaussian hump over a constant water depth, and the predicted numerical results are compared with analytical solutions. The results of the present numerical model are accurate in comparison with those of existing models.