• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.35-40
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convertive terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. In this paper k-$\varepsilon$ turbulence model with wall function is used to increase efficiency of computation times.

• Journal of the Institute of Electronics and Information Engineers
• /
• v.50 no.11
• /
• pp.165-170
• /
• 2013

When we build up a 3D model from the given 2D images, linear algorithms are often used to reduce computational cost or for initialization of nonlinear algorithms. However, contemporary linear algorithms have apparently linear structures, but virtually they are implemented using SVD. The SVD is also implemented using numerical analysis algorithms that need initialization. Moreover, solutions using SVD are more difficult to analyze than closed-form solutions. To avoid from such inconvenient numerical analysis algorithms of the contemporary methods and for convenient analysis of solutions, this paper proposes a convenient linear method that produces a closed-form solution.

• Structural Engineering and Mechanics
• /
• v.2 no.3
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• Journal of Soil and Groundwater Environment
• /
• v.24 no.3
• /
• pp.13-23
• /
• 2019

The purpose of this study is to suggest conceptual models based on finite numerical method that can be used to assess contaminant transport through subsurface and estimate exposed concentration at contaminated site. This study tested various assumptions of the numerical models for contaminant transport in unsaturated and saturated zones to simulate the pathways to the human exposal point. For this purpose, models for seven possible scenarios of contaminant transport were simulated using the numerical code MODFLOW and MT3D. The simulation results that showed different peak concentrations and travel times were compared. In conclusion, the potential utility of the numerical models in the site specific risk analysis suggested as well as future research ramifications.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.10
• /
• pp.55-62
• /
• 2012

In this research, the algorithm for simulating specific grain size distribution(GSD) with large diameter granular material was developed using the distinct element analysis program $PFC^{3D}$(Particle Flow Code). This modeling approach can generate the initial distinct elements without clump logic or cluster logic and prevent distinct element from escaping through the confining walls during the process. Finally the proposed distinct element model is used to simulate large triaxial compression test of the rockfill material and we compared the simulation output with lab test results. Simulation results of Assembly showed very well agreement with the GSD of the test sample and numerical modeling of granular material would be possible for various stress conditions using this application through the calibration.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.95-98
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid} method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.55 no.3
• /
• pp.123-130
• /
• 2006

This paper presents characteristics analysis of skewed Synchronous Reluctance Motor using equivalent magnetic circuit and compares with the result of Finite Element Method. Torque ripple must be reduced, because it is producing noise and vibration. There is many kinds of method to reduce torque ripple, but generally we apply skewing stator or rotor. The 2D Finite Element Method(FEM) or 3D FEM is used to analyze the motor, since skew influence the average torque in the motor. However, the FEM takes much time in spite of the advanced computer and numerical technique. This paper will analyze characteristics of skewed synchronous reluctance motor using equivalent magnetic circuit.

• International Journal of Control, Automation, and Systems
• /
• v.2 no.3
• /
• pp.319-332
• /
• 2004

This paper proposes some further improvements on N.K. Bose's 2D stability test for polynomials with real coefficients by revealing symmetric properties of the polynomials, resultants occurring in the test and by generalizing Sturm's method. The improved test can be fulfilled by a totally algebraic algorithm with a finite number of steps and the computational complexity is largely reduced as it involves only certain real variable polynomials with degrees not exceeding half of their previous complex variable counterparts. Nontrivial examples for 2D polynomials having both numerical and literal coefficients are also shown to illustrate the computational advantage of the proposed method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.4
• /
• pp.407-415
• /
• 2011

The stability condition and wideband characteristics of 3D ID-FDTD algorithm which has low dispersion error with isotropic dispersion are presented in this paper. 3D ID-FDTD method was proposed to improve the defect of the Yee FDTD such as the anisotropy and large dispersion error. The published paper calculated the stability condition of 3D ID-FDTD algorithm by using numerical method, however, it is thought that the examples were not sufficient to verify the stability condition. Thus, in this paper, various simulations are included in order to hold reliability under the conditions that the plane wave propagation is assumed with a single frequency and a wideband frequency. Also, the 3D ID-FDTD algorithm is compared to those that have the similar FDTD algorithm with ID-FDTD such as Forgy's method and non-standard FDTD method in a wideband. Finally, the radar cross section(RCS) for the large sphere with high dielectric constant is calculated.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.77-94
• /
• 2016

A forked tunnel, as a special complicated underground structure, is composed of big-arch tunnel, multi-arch tunnel, neighborhood tunnels and separate tunnels according to the different distances between two separate tunnels. Due to the complicated process of design and construction, surrounding jointed rock mass stability of the big-arch tunnel which belongs to the forked tunnel during excavation is a hot issue that needs special attentions. In this paper, an elasto-plastic damage constitutive model for jointed rock mass is proposed based on the coupling method considering elasto-plastic and damage theories, and the irreversible thermodynamics theory. Based on this elasto-plastic damage constitutive model, a three dimensional elasto-plastic damage finite element code (D-FEM) is implemented using Visual Fortran language, which can numerically simulate the whole excavation process of underground project and perform the structural stability of the surrounding rock mass. Comparing with a popular commercial computer code, three dimensional fast Lagrangian analysis of continua (FLAC3D), this D-FEM has advantages in terms of rapid computing process, element grouping function and providing more material models. After that, FLAC3D and D-FEM are simultaneously used to perform the structural stability analysis of the surrounding rock mass in the forked tunnel considering three different computing schemes. The final numerical results behave almost consistent using both FLAC3D and D-FEM. But from the point of numerically obtained damage softening areas, the numerical results obtained by D-FEM more closely approach the practical behaviors of in-situ surrounding rock mass.