• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.380-388
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation of condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new numerical scheme to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the cupid code.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.311-317
• /
• 2011

The present paper deals with the continuous works of extending the multi-dimensional limiting process (MLP) for compressible flows, which has been quite successful in finite volume methods, into discontinuous Galerkin (DG) methods. From the series of the previous, it was observed that the MLP shows several superior characteristics, such as an efficient controlling of multi-dimensional oscillations and accurate capturing of both discontinuous and continuous flow features. Mathematically, fundamental mechanism of oscillation-control in multiple dimensions has been established by satisfaction of the maximum principle. The MLP limiting strategy is extended into DG framework, which takes advantage of higher-order reconstruction within compact stencil, to capture detailed flow structures very accurately. At the present, it is observed that the proposed approach yields outstanding performances in resolving non-compressive as well as compressive flaw features. In the presentation, further numerical analyses and results are going to be presented to validate that the newly developed DG-MLP methods provide quite desirable performances in controlling numerical oscillations as well as capturing key flow features.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.37 no.7C
• /
• pp.601-609
• /
• 2012

In this paper, a limited-view CT image reconstruction method was studied to reduce the scan times and the X-ray dose for the patients. To reduce streak artifacts which is caused by insufficient number of views, we introduce a sinogram interpolation method based on image matching. Image matching is achieved using the characteristics of the neighboring views including intensity, gradient and distance between the pixels. Interpolation is performed using the image matching results.. A numerical phantom and Al-acryl phantom were used for evaluating the effectiveness of the proposed interpolation method. The results showed that streak artifacts were reduced in the reconstructed images while the details of the images were preserved. Moreover, maximum 5% improvements in terms of PSNR were observed.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.4
• /
• pp.263-273
• /
• 2011

The diffuse interface model of Saurel et al. (2008) is used for the computation of compressible cavitating flows around underwater missiles. Such systems use gas injection and natural cavitation to reduce drag effects. Consequently material interfaces appear separating liquid and gas. These interfaces may have a really complex dynamics such that only a few formulations are able to predict their evolution. Contrarily to front tracking or interface reconstruction method the interfaces are computed as diffused numerical zones, that are captured in a routinely manner, as is done usually with gas dynamics solvers for shocks and contact discontinuity. With the present approach, a single set of partial differential equations is solved everywhere, with a single numerical scheme. This leads to very efficient solvers. The algorithm derived in Saurel et al. (2009) is used to compute cavitation pockets around solid bodies. It is first validated against experiments done in cavitation tunnel at CNU. Then it is used to compute flows around high speed underwater systems (Shkval-like missile). Performance data are then computed showing method ability to predict forces acting on the system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.123-130
• /
• 2005

A digital in-line holographic particle image velocimetry (HPIV) which can be applied to measure three-dimensional velocity fields of turbulent flows was developed. There are three different implementation methods of HPIV: traditional film-based HPIV, intermediate HPIV and digital HPIV. The traditional film-based HPIV and intermediate HPIV method is rather troublesome to do experiments and takes long calculation time, compared with the digital HPIV, Configuration of the digital in-line HPIV is simple and the data processing routine is similar to conventional 2D PIV methods. The digital HPIV velocity field measurement consists of four steps: recording, numerical reconstruction, particle extraction and velocity extraction. In the velocity extraction process, we improved PTV algorithm to extract the displacement of particle each placed in 3D space. The developed digital in-line HPIV system was applied to a vertical jet flow. The 3D velocity vectors measured by the digital HPIV method in the near field are in a good agreement with 2D PIV results.

• Korean Journal of Optics and Photonics
• /
• v.14 no.3
• /
• pp.244-248
• /
• 2003

When confocal lenses are applied to a digital holography system, the interference pattern of a larger object size can be recorded on CCD. The angle of incoming light to the CCD can be reduced by the ratio of the two focal lengths of confocal lenses. The recordable spatial frequency is limited by the unit cell size of the CCD. Therefore the spatial frequency of interference on the CCD is lowered by reduction of the incoming light angle. By using confocal lenses , another merit can be achieved that the area of the zero order diffraction is reduced by the square of the ratio of two focal lengths at the numerical reconstruction.

• Journal of the Korean Society of Visualization
• /
• v.19 no.3
• /
• pp.115-122
• /
• 2021

Oral cancer surgery typically consists of resection of lesion, neck dissection and reconstruction, and it has an impact on the position of hyoid bone. Therefore, morphological change of airway can occur since the geometric parameter of airway is correlated with the hyoid bone. Airflow is affected by geometry of the airway. In this study, flow characteristics were compared between pre- and post-surgery models by both particle image velocimetry (PIV) and numerical simulation. 3D model of upper airway was reconstructed based on CT data. Velocity is accelerated by the reduced channel area, and vortex and recirculation region are observed in pre- and post-surgery models. For the post-surgery model, high pressure distribution is developed by significantly decreased hydraulic diameter, and the longitudinal flow stream is also interrupted.

• Structural Engineering and Mechanics
• /
• v.88 no.4
• /
• pp.301-309
• /
• 2023

In this paper, a new modified conjugate gradient (MCG) method is presented which is based on a new gradient regularizer, and this method is used to identify the dynamic load on airfoil structure without and with considering random structure parameters. First of all, the newly proposed algorithm is proved to be efficient and convergent through the rigorous mathematics theory and the numerical results of determinate dynamic load identification. Secondly, using the perturbation method, we transform uncertain inverse problem about force reconstruction into determinate load identification problem. Lastly, the statistical characteristics of identified load are evaluated by statistical methods. Especially, this newly proposed approach has successfully solved determinate and uncertain inverse problems about dynamic load identification. Numerical simulations validate that the newly developed method in this paper is feasible and stable in solving load identification problems without and with considering random structure parameters. Additionally, it also shows that most of the observation error of the proposed algorithm in solving dynamic load identification of deterministic and random structure is respectively within 11.13%, 20%.

• Journal of IKEEE
• /
• v.18 no.2
• /
• pp.272-281
• /
• 2014

In electrical impedance tomography (EIT), modified Newton Raphson (mNR) method is widely used inverse algorithm for static image reconstruction due to its convergence speed and estimation accuracy. The unknown conductivity distribution is estimated iteratively by minimizing a cost functional such that the residual error namely the difference in measured and calculated voltages is reduced. Although, mNR method has good estimation performance, EIT inverse problem still suffers from ill-conditioned and ill-posedness nature. To mitigate the ill-posedness, generally, regularization methods are adopted. The inverse solution is highly dependent on the choice of regularization parameter. In most cases, the regularization parameter has a constant value and is chosen based on experience or trail and error approach. In situations, when the internal distribution changes or with high measurement noise, the solution does not get converged with the use of constant regularization parameter. Therefore, in this paper, in order to improve the image reconstruction performance, we propose a new scheme to determine the regularization parameter. The regularization parameter is computed based on residual error and updated every iteration. The proposed scheme is tested with numerical simulations and laboratory phantom experiments. The results show an improved reconstruction performance when using the proposed regularization scheme as compared to constant regularization scheme.