• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1834-1845
• /
• 2021

Liquid lead-bismuth cooled fast reactor is one of the most promising reactor types among the fourth-generation nuclear energy systems. The flow and heat transfer characteristics of lead-bismuth eutectic (LBE) are completely different from ordinary fluids due to its special thermal properties, causing that the traditional Reynolds analogy is no longer recommended and appropriate. More accurate turbulence flow and heat transfer model for the liquid metal lead-bismuth should be developed and applied in CFD simulation. In this paper, a specific CFD solver for simulating the flow and heat transfer of liquid lead-bismuth based on the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was developed based on the open source platform OpenFOAM. Then the advantage of proposed model was demonstrated and validated against a set of experimental data. Finally, the simulation of LBE turbulent flow and heat transfer in a 7-pin wire-wrapped rod bundle with the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was carried out. The influence of wire on the flow and heat transfer characteristics and the three-dimensional distribution of key thermal hydraulic parameters such as temperature, cross-flow velocity and Nusselt number were studied and presented. Compared with the traditional SED model with a constant Pr_t = 1.5 or 2.0, the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is more accurate on predicting the turbulence flow and heat transfer of liquid lead-bismuth. The average relative error of the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is reduced by 11.1% at most under the simulation conditions in this paper. This work is meaningful for the thermal hydraulic analysis and structure design of fuel assembly in the liquid lead-bismuth cooled fast reactor.

• Korean Chemical Engineering Research
• /
• v.57 no.1
• /
• pp.51-57
• /
• 2019

A plant layout problem has a large impact on the overall construction cost of a plant. When determining a plant layout, various constraints associating with safety, environment, sufficient maintenance area, passages for workers, etc have to be considered together. In general plant layout problems, the main goal is to minimize the length of piping connecting equipments as satisfying various constraints. Since the process may suffer from the heat and friction loss, the piping length between equipments should be shorter. This problem can be represented by the mathematical formulation and the optimal solutions can be investigated by an optimization solver. General researches have overlooked many constraints such as maintenance spaces and safety distances between equipments. And, previous researches have tested benchmark processes. What the lack of general researches is that there is no realistic comparison. In this study, the plant layout of a real industrial C3MR (Propane precooling Mixed Refrigerant) process is studied. A MILP (Mixed Integer Linear Programming) including various constraints is developed. To avoid the violation of constraints, penalty functions are introduced. However, conventional optimization solvers handling the derivatives of an objective functions can not solve this problem due to the complexities of equations. Therefore, the PSO (Particle Swarm Optimization), which investigate an optimal solutions without differential equations, is selected to solve this problem. The results show that a proposed method contributes to saving the capital expenditures.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.931-941
• /
• 2020

There is a growing interest this paper for ocean sensing where autonomous vehicles can play an essential role in assisting engineers, researchers, and scientists with environmental monitoring and collecting oceanographic data. This study was conducted to develop a rigid sail for the autonomous sailing drone. Our study aims to numerically analyze the aerodynamic characteristics of curvy twin sail and compare it with wing sail. Because racing regulations limit the sail shape, only the two-dimensional geometry (2D) was open for an optimization. Therefore, the first objective was to identify the aerodynamic performance of such curvy twin sails. The secondary objective was to estimate the effect of the sail's spacing and shapes. A viscous Navier-Stokes flow solver was used for the numerical aerodynamic analysis. The 2D aerodynamic investigation is a preliminary evaluation. The results indicated that the curvy twin sail designs have improved lift, drag, and driving force coefficient compared to the wing sails. The spacing between the port and starboard sails of curvy twin sail was an important parameter. The spacing is 0.035 L, 0.07 L, and 0.14 L shows the lift coefficient reduction because of dramatically stall effect, while flow separation is improved with spacing is 0.21 L, 0.28 L, and 0.35 L. Significantly, the spacing 0.28 L shows the maximum high pressure at the lower area and the small low pressure area at leading edges. Therefore, the highest lift was generated.

• Journal of Korea Water Resources Association
• /
• v.54 no.11
• /
• pp.985-998
• /
• 2021

In general, the upstream dam changes downstream flow regime dramatically, i.e., from natural flow regime to hydropeaking flows. This study investigates the impact of the natural flow pattern on downstream fish habitat in a regulated river in Korea using the physical habitat simulation. The study area is a 13.4 km long reach of the Geum-gang River, located downstream from the Yongdam Dam, Korea. A field monitoring revealed that three fish species are dominant, namely Zacco platypus, Coreoleuciscus splendidus, and Opsariichthys bidens, and they account for 70% of the total fish community. Specially, Opsariichthys bidens is an indigenous species in the Geum-gang River. The three fish species are selected as target fish species for the physical habitat simulation. The Nays2D model, a 2D shallow water equation solver, and the HSI (Habitat Suitability Index) model are used for hydraulic and habitat simulations, respectively. To assess the impact of the natural flow pattern, this study uses the annual natural flow regime and hydropeaking flows from the dam. It is found that the natural flow regime increases significantly the Composite Suitability Index (CSI) in the study reach. Then, using the Building Block Approach (BBA), the scenarios for the modifying dam operations are presented in the study reach. Both Scenario 1 and scenario 2 are proposed by using the hydrological method considering both magnitude and duration of the inflow and averaging the inflow over each month, respectively. It is revealed that the natural flow regime embodied in scenario 1 and scenario 2 increases the Weighted Usable Area (WUA) significantly, compared to the hydropeaking flows. In conclusion, the modifying the dam operations by restoring to the natural flow pattern is advantageous to fish community.

• Journal of the Society of Disaster Information
• /
• v.18 no.3
• /
• pp.542-549
• /
• 2022

Purpose: Pipes are widely used as applied devices in many industrial fields such as machinery, electronics, electricity, and plants, and are also widely used in safety-related fields such as firefighting and chemistry. With the diversification of products, the importance of technology in the piping field is also increasing. In particular, when changing the existing copper pipe to stainless steel, it is necessary to evaluate safety and flow characteristics through structural analysis or flow analysis. Method: This study investigated the structural stability of the 6.35 and 15.88 socket models, which are integrated insert type connectors developed by a company, using FEM. For structural analysis, HyperMesh as pre-processor, HYPER VIEW as post-processor, and LS-DYNA as solver were used. Result: In the case of 6.35 socket, the maximum stresses at hook, holder and hinge were 95.02MPa, 19.59MPa and 44.01MPa, respectively, and in case of 15.88 socket, it was 127.7 MPa, 40.09MPa and 45.23MPa, respectively. Conclusion: Comparing the stress distribution of the two socket models, the stress in the 15.88 socket, which has a relatively large outer diameter, appears to be large overall, but it is significantly lower than the yield stress of each material, indicating that there is no problem in structural safety in both models.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.261-269
• /
• 2021

Without any validation of the incompressible assumption, most of previous studies on cavitation flow and its noise have utilized numerical methods based on the incompressible Reynolds Average Navier-Stokes (RANS) equations because of advantage of its efficiency. In this study, to investigate the effects of the flow compressibility on the Tip Vortex Cavitation (TVC) flow and noise, both the incompressible and compressible simulations are performed to simulate the TVC flow, and the Ffowcs Williams and Hawkings (FW-H) integral equation is utilized to predict the TVC noise. The DARPA Suboff submarine body with an underwater propeller of a skew angle of 17 degree is targeted to account for the effects of upstream disturbance. The computation domain is set to be same as the test-section of the large cavitation tunnel in Korea Research Institute of Ships and Ocean Engineering to compare the prediction results with the measured ones. To predict the TVC accurately, the Delayed Detached Eddy Simulation (DDES) technique is used in combination with the adaptive grid techniques. The acoustic spectrum obtained using the compressible flow solver shows closer agreement with the measured one.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.289-301
• /
• 2009

In this work, wave run-up heights and resultant wave forces on a vertical revetment due to tsunami (solitary wave) are investigated numerically using a numerical wave tank model called CADMAS-SURF (CDIT, 2001. Research and Development of Numerical Wave Channel (CADMAS-SURF). CDIT library, No. 12, Japan.), which is based on a 2-D Navier-Stokes solver, coupled to a volume of fluid (VOF) method. The third order approximate solution (Fenton, 1972. A ninth-order solution for the solitary wave. J. of Fluid Mech., Vol. 53, No.2, pp.257-271) is used to generate solitary waves and implemented in original CADMAS-SURF code. Numerical results of the wave profiles and forces are in good agreements with available experimental data. Using the numerical results, the regression curves determined from the least-square analysis are proposed, which can be used to determine the maximum wave run-up height and force on a vertical revetment due to tsunami. In addition, the capability of CADMAS-SURF is demonstrated for tsunami wave forces acting on an onshore structure using various configuration computations including the variations of the crown heights of the vertical wall and the position of the onshore structure. Based on the numerical results such as water level, velocity field and wave force, the direct effects of tsunami on an onshore structure are discussed.

• Composites Research
• /
• v.36 no.5
• /
• pp.303-309
• /
• 2023

In this paper, a Progressive Damage and Failure Analysis (PDFA) modeling method was developed using ABAQUS/EXPLICIT to predict in-plane damage and delamination for Open-Hole Compression (OHC) testing. The proposed PDFA model was constructed based on Hashin criteria and cohesive behavior. The strength and stiffness of OHC specimens with three types of stacking sequences [(45/-45/0₂)₃]_s , [(45/0/-45/90)₃]_s and [45/-45/0/45/-45/90/(45/-45)₂]_s were compared to comprehensively evaluate the validity of the Finite Element(FE) model of PDFA. The strength and stiffness of the OHC specimens were predicted relatively well, with less than a percentage error 10.0 %. For the numerical simulation case for each layup, the damage initiation/evolution of OHC specimens were evaluated for delamination and tension/compression matrix damage before and after failure.