• Nuclear Engineering and Technology
• /
• 제56권1호
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• 대한의용생체공학회:의공학회지
• /
• 제28권6호
• /
• pp.744-749
• /
• 2007

Kinematic and kinetic analysis has been performed for Soft Golf swings utilizing realistic three dimensional computer simulations based on three dimensional motion tracking data. Soft Golf is a newly developed recreational sport in South Korea aimed to become a safe and easy-to-learn sport for all ages. The advantage of Soft Golf stems from lighter weight of the club and much larger area of the sweet spot. This paper tries to look into kinematic and kinetic aspects of soft golf swings compared to regular golf swing and find the advantages of lighter Soft Golf clubs. For this purpose, swing motions of older aged participants were captured and kinematic analysis was performed for various kinematic parameters such as club head velocity, joint angular velocity, and joint range of motions as a pilot study. Kinetic analysis was performed by applying kinematic data to computer simulation models constructed from anthropometric database and the measurements from the participants. The simulations were solved using multi-body dynamics solver. Firstly, the kinematic parameters such as joint angles were obtained by solving inverse dynamics problem based on motion tracking data. Secondly, the kinetic parameters such as joint torques were obtained by solving control dynamics problem of making joint torque to follow pre-defined joint angle data. The results showed that mechanical loadings to major joints were reduced with lighter Soft Golf club.

• International Journal of Aeronautical and Space Sciences
• /
• 제2권2호
• /
• pp.1-18
• /
• 2001

A numerical technique for simulating the separation dynamics of strap-on boosters jettisoned in the dense atmosphere is presented. Six degree of freedom rigid body equations of motion are integrated into the three-dimensional unsteady Navier-Stokes solution procedure to determine the dynamic motions of strap-ons. An automated Chimera overlaid grid technique is introduced to achieve maximum efficiency for multi-body dynamic motion and a domain division technique is implemented in order to reduce the computational cost required to find interpolation points in the Chimera grids. The flow solver is validated by comparing the computed results around the Titan IV launch vehicle with experimental data. The complete analysis process is then applied to the. H-II launch vehicle, the central rocket in japans space program, the CZ-3C launch vehicle developed in China and the KSR-III, a three-stage sounding rocket being developed in Korea. From the analyses, separation trajectories of strap-on boosters are predicted and aerodynamic characteristics around the vehicles at every time interval are examined. In addition, separation-impulse devices generally introduced for safe separation of strap-ons are properly modeled in the present paper and the jettisoning force requirements are examined quantitatively.

• Nuclear Engineering and Technology
• /
• 제55권9호
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• 항공우주시스템공학회지
• /
• 제7권3호
• /
• pp.1-6
• /
• 2013

Numerical optimization of two dimensional high lift configuration was performed with flow solver and optimization method based on RSM(Response Surface Model). Navier-Stokes solver with Spalart-Allmaras turbulence model was selected for the simulation of highly complex and separated flows on the flap. For the simultaneous optimization of both flap shape and setting (gap/overlap), 10 design variables (eight variables for flap shape variation and two variables for flap setting) were chosen. In order to generate the response surface model, 128 experimental points were selected for 10 design variables. The objective function considering maximum lift coefficient, lift to drag ratio and lift coefficient at specific angle of attack was selected to reduce flow separation on the flap surface. The present method was applied to two dimensional fowler flap in landing configuration. After applying the present method, it was shown that the optimized high lift configuration had less flow separation on the flap surface and lift to drag ratio was suppressed over entire angle of attack range.

• 한국전산유체공학회지
• /
• 제11권1호
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• 한국전산유체공학회지
• /
• 제19권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.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2009년도 추계학술대회 논문집
• /
• pp.600-605
• /
• 2009

A general purpose viscous flow solver Ansys CFX is used to study a Savonius type wave energy converter in a 3D numerical viscous wave tank. This paper presents the results of a computational fluid dynamics (CFD) analysis of the effect of blade configuration on the performance of 3 bladed Savonius rotors for wave energy extraction. A piston-type wave generator was incorporated in the computational domain to generate the desired incident waves. A complete OWC system with a 3-bladed Savonius rotor was modeled in a three dimensional numerical wave tank and the hydrodynamic conversion efficiency was estimated. The flow over the rotors is assumed to be two-dimensional (2D), viscous, turbulent and unsteady. The CFX code is used with a solver of the coupled conservation equations of mass, momentum and energy, with an implicit time scheme and with the adoption of the hexahedral mesh and the moving mesh techniques in areas of moving surfaces. Turbulence is modeled with the k.e model. Simulations were carried out simultaneously for the rotor angle and the helical twist. The results indicate that the developed models are suitable to analyze the water flows both in the chamber and in the turbine. For the turbine, the numerical results of torque were compared for all the cases.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2009년도 춘계학술대회 논문집 특별세미나,특별/일반세션
• /
• pp.225-230
• /
• 2009

From the view point of railway vehicle dynamics, the interaction between wheel and rail have an huge effect on the behavior of the vehicle. This phenomenon is an unique motion, only for railway vehicles. Furthermore, close investigation of the backgrounds of the interaction is the key to estimate the dynamic behavior of the vehicle, successfully. To evaluate the model including flexible bodies such as car body and catenary system of the next generation express train, it is necessary to develop proper dynamic solver including a wheel rail contact module. In this study, wheel-rail contact module is developed using the general purpose dynamic solver. First of all, the procedure for calculation of the wheel-rail contact force has been established. Generally, yaw angle of the wheelset is ignored. Sets of information are summarized as tables and splined for further uses. With this information, normal force and creep coefficient can be extracted and used for FASTSIM algorithm, which has been shown good reliability over years. Normal force and longitudinal, lateral force at the contact surface are also calculated. Those data are verified by commercial railway simulation program 'VAMPIRE'. This procedure and program can offer a basic process for estimation of the dynamic behavior and wear of the wheel-rail system, even while running on the curved rail. Finally, multi-dimensional inspection tool will be developed including the prediction of the derailment.

• 한국전산구조공학회논문집
• /
• 제36권1호
• /
• pp.67-74
• /
• 2023

A meshless technique using the geometric conservation least-squares method (GC-LSM) was devised to discretize the governing equation of linear elasticity. Although the finite-element method is widely used for structural analysis, a meshless method was developed because of its advantages in a moving grid system. This work is the preliminary phase for developing a fully meshless-based fluid-structure interaction solver. In this study, Cauchy's momentum equation was discretized in strong form using GC-LSM for the structural domain, and the Newmark beta method was used for time integration. The solver was validated in 1D, 2D, and 3D benchmarking problems. Static and dynamic results were obtained. The results are more accurate than those of analytic solutions.