• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• Journal of Environmental Impact Assessment
• /
• v.18 no.6
• /
• pp.331-346
• /
• 2009

The serious problem facing two-dimensional finite element hydraulic model is the treatment of wet and dry areas. This situation is encountered in most practical river and coastal engineering problems, such as flood propagation, dam break analysis and so on. Especially, dry areas result in mathematical complications and require special treatment. The objective of this study is to investigate the wet and dry parameters that have direct relevance to model performance in situations where inundation of initially dry areas occurs. Several numerical simulations were carried out, which examined the performance of the marsh porosity method of RMA-2 model to investigate for application of parameters. Experimental channel with partly dry side slopes, straight channel with irregular geometry and Han river were performed for tests. As a result of this study, effectively applied marsh porosity method provide a reliable results for flow distribution of wet and dry area, it could be further developed to basis for extending to water quality and sediment transport analysis.

• Journal of the Korean Society of Visualization
• /
• v.17 no.3
• /
• pp.12-18
• /
• 2019

Direct numerical simulation of a fully developed turbulent plane Couette-Poiseuille flow with a two-dimensional (2-D) rod-roughened wall is performed to investigate the impacts of the surface roughness. It is shown that the logarithmic region in the mean velocity profile over the rough wall Couette-Poiseuille flow is significantly shortened by the surface roughness compared to that over a turbulent Couette-Poiseuille flow with smooth wall. The Reynolds shear stress over the rough wall Couette-Poiseuille flow is decreased compared to that for a smooth case in the outer layer. These results are attributed to weakened turbulence activity or roll-cell mode over the rough wall Couette-Poiseuille flow near the channel centerline due to suppressed development of u'-structure on the top wall, as documented through spanwise energy spectra of the streamwise velocity fluctuations. Inspection of congregation motion near the bottom wall and time evolution of u'-structure reveal weakened co-supporting cycle for the rough wall case.

• Nuclear Engineering and Technology
• /
• v.30 no.5
• /
• pp.452-461
• /
• 1998

The pin power distribution is determined by analyzing the rotational gamma scanning data for 36 element fuel bundle of HANARO. A fission monitor of Nb$^{95}$ is chosen by considering the criteria of the half-life, fission yield, emitting ${\gamma}$-ray energy and probability. The ${\gamma}$-ray spectra were measured in Korea Atomic Energy Research Institute(KAERI) by using a HPGe detector and by rotating the fuel bundle at steps of 10$^{\circ}$. The counting rates of Nb$^{95}$ 766 keV ${\gamma}$-rays are determined by analyzing the full absorption peak in the spectra. A 36$\times$36 response matrix is obtained from calculating the contribution of each rod at every scanning angle by assuming 2-dimensional and parallel beam approximations for the measuring geometry. In terms of the measured counting rates and the calculated response matrix, an inverse problem is set up for the unknown distribution of activity concentrations of pins. To select a suitable solving method, the performances of three direct methods and the iterative least-square method are tested by solving simulation examples. The final solution is obtained by using the iterative least-square method that shows a good stability. The influences of detection error, step size of rotation and the collimator width are discussed on the accuracy of the numerical solution. Hence an improvement in the accuracy of the solution is proposed by reducing the collimator width of the scanning arrangement.

• Journal of the Korean Society of Combustion
• /
• v.9 no.2
• /
• pp.10-17
• /
• 2004

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman#s detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. The results show that fuel-side and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.39 no.1
• /
• pp.42-48
• /
• 2002

In this paper, a decision method using multiple threshold values for PN code acquisition in Direct Sequence Code Division Multiple Access (DS-CDMA) systems is described. We apply this technique to the conventional double dwell serial search algorithm and analyze it in terms of mean code acquisition time. For the analysis, we present mathematical model of proposed algorithm and also perform the simulation under IMT-2000 channel models. Numerical results show that our proposed scheme outperforms the conventional one by 0.2 - 0.5 sec with respect to the mean code acquisition time because multiple threshold values mitigate the possible decline in search performance caused by the use of a single threshold.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.19 no.5
• /
• pp.833-841
• /
• 1994

In this paper we develop an approximation formalism for the queue length distribution of general queueing models. Our formalism is based on two steps of analytic approximation employing both the lower and upper bound techniques. It is favorable to a fast numerical calcuation for the queue length distribution of a superposition of a superposition of arbitary type traffic sources. In the application. M+ N D /M/1 is considered. The calculated result for queue length distribution measured by arriving or leaving customers show a good agreement with the direct simulation of the system. Especially, we demonstrate that our formula for M/M/1 is equivalent to the exact solution, while that D/M/1 is simplified in an analytic form.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2004.04a
• /
• pp.67-70
• /
• 2004

An investigation was performed to study the impact damage of the laminated composite plates caused by a low- velocity foreign object with multi-scale modeling based on the concepts of Direct Numerical Simulation (DNS)[4]. In the micro-scale part, we discretize the composite plates through separate modeling of fiber and matrix for the local microscopic analysis. A micro-scalemodel was developed for predicting the initiation of the damage and the extent of the final damage as a function of material properties, laminate configuration and the impactor's mass, etc. Anda macro-scale model was developed for description of global dynamic behavior. The connection betweenmicroscopic and macroscopic is implemented by the tied interface constraints of LS-DYNA contact card. A transient dynamic finite element analysis was adopted for calculating the contact force history and the stresses and strains inside the composites during impact resulting from a point-nose impactor. The low-velocity impact events such as contact force, deformation, etc. are simulated in the macroscopic sense and the impact damages, fiber-breakage, matrix cracking and delamination etc. are examined in the microscopic sense.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1407-1416
• /
• 2004

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$/$N_2$-Air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this computation. The results showed that an initially flat stagnation plane, on which an axial velocity was zero, was deformed into a complex-shaped plane, and an initial stagnation point was moved far away from a vortex head when the counterflow field was perturbed by the vortex. It was noted that the movement of stagnation point could alter the species transport mechanism to the flame surface. It was also identified that the altered species transport mechanism affected the distributions of the mixture fraction and the scalar dissipation rate.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.31 no.7A
• /
• pp.682-688
• /
• 2006

In this paper, a soft-decision method for initial PN code acquisition in pilot-channel aided Direct Sequence Code Division Multiple Access (DS-CDMA) systems is proposed in order to improve the acquisition performance. We apply this technique to the conventional hybrid search algorithm and analyze it in terms of mean code acquisition time. For the analysis, we present mathematical model of proposed algorithm and also perform the simulation under IMT-2000 channel models. Numerical results show that our proposed scheme outperforms the conventional one by 0.2 - 0.4 sec with respect to the mean code acquisition time because the soft decision technique can mitigate the possible decline in search performance caused by the use of a hard-decision technique.