• Journal of Soil and Groundwater Environment
• /
• v.14 no.6
• /
• pp.1-18
• /
• 2009

A series of three-dimensional numerical simulations using a generalized multidimensional hydrodynamic dispersion numerical model is performed to simulate effectively and to evaluate quantitatively impacts of urbanization on density-dependent groundwater flow and salt transport in a coastal aquifer system, Suyeong-Gu, Busan, Korea. A series of steady-state numerical simulations of groundwater flow and salt transport before urbanization with material properties of geologic formations, which are established by numerical modeling calibrations considering all the urbanization factors, is performed first without considering all the urbanization factors. A series of transient-state numerical simulations of groundwater flow and salt transport after urbanization is then performed considering the urbanization factors individually and all together. Finally, the results of both numerical simulations are compared with each other and analyzed. The results of the numerical simulations show that density-dependent groundwater flow, salt transport, and seawater intrusion in the coastal aquifer system are intensively and extensively impacted by the urbanization factors. Especially, these urbanization factors result in the changes of the total groundwater volume and salt mass in the coastal aquifer system. However, such impacts of each urbanization factor are not spatially uniform but locally different.

• Wind and Structures
• /
• v.37 no.4
• /
• pp.255-274
• /
• 2023

The aim of this paper is to enhance the accuracy of predicting time-averaged external surface pressures on low-rise buildings by utilizing Computational Fluid Dynamics (CFD) simulations. To achieve this, benchmark studies of the Silsoe cube and the Texas Tech University (TTU) experimental building are employed for comparison with simulation results. The paper is structured into three main sections. In the initial part, an appropriate domain size is selected based on the precision of mean pressure coefficients on the windward face of the cube, utilizing Reynolds Averaged Navier-Stokes (RANS) turbulence models. Subsequently, recommendations regarding the optimal computational domain size for an isolated building are provided based on revised findings. Moving on to the second part, the Silsoe cube model is examined within a horizontally homogeneous computational domain using more accurate turbulence models, such as Large Eddy Simulation (LES) and hybrid RANS-LES models. For computational efficiency, transient simulation settings are employed, building upon previous studies by the authors at the Windstorm Impact, Science, and Engineering (WISE) Lab, Louisiana State University (LSU). An optimal meshing strategy is determined for LES based on a grid convergence study. Three hybrid RANS-LES cases are investigated to achieve desired enhancements in the distribution of mean pressure coefficients on the Silsoe cube. In the final part, a 1:10 scale model of the TTU building is studied, incorporating the insights gained from the second part. The generated flow characteristics, including vertical profiles of mean velocity, turbulence intensity, and velocity spectra (small and large eddies), exhibit good agreement with full-scale (TTU) measurements. The results indicate promising roof pressures achieved through the careful consideration of meshing strategy, time step, domain size, inflow turbulence, near-wall treatment, and turbulence models. Moreover, this paper demonstrates an improvement in mean roof pressures compared to other state-of-the-art studies, thus highlighting the significance of CFD simulations in building aerodynamics.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.5
• /
• pp.1-9
• /
• 2013

Variable area nozzle, where both throat and exit area vary, is required for optimal expansion and optimal nozzle shape upon operation of after-burner. Steady-state and transient analyses are carried out for each condition with and without afterburner operation and as a function of the location of the nozzle flap. Effects of that nozzle displacement on flow and thrust characteristics are analyzed from numerical results. With variable area nozzle adopted, the combustion field is variable in time, leading to periodically variable thrust. For off-design conditions, flow separation shows up due to over expansion at the flap tips and shock wave does in the nozzle due to under expansion. The undesirable phenomena can be solved by control of variable area nozzle.

• Proceedings of the KIEE Conference
• /
• 2005.11b
• /
• pp.46-48
• /
• 2005

This paper discusses a probabilistic method for power system security assessment. The security analysis relates to the ability of the electric power systems to survive sudden disturbances such as electric short circuits or unanticipated loss of system elements. It consists of both steady state and dynamic security analyses, which are not two separate issues but should be considered together. In steady state security analysis including voltage security analysis, the analysis checks that the system is operated within security limits by OPF (optimal power flow) after the transition to a new operating point. Until now, many utilities have difficulty in including dynamic aspects due to computational capabilities. On the other hand. dynamic security analysis is required to ensure that the transition may lead to an acceptable operating condition. Transient stability, which is the ability of power systems to maintain synchronism when subjected to a large disturbance. is a principal component in dynamic security analysis. Usually any loss of synchronism may cause additional outages and make the present steady state analysis of the post-contingency condition inadequate for unstable cases. This is the reason for the need of dynamic studies in power systems. Probabilistic criterion can be used to recognize the probabilistic nature of system components while considering system security. In this approach. we do not have to assign any predetermined margin of safety. A comprehensive conceptual framework for probabilistic static and dynamic assessment is presented in this paper. The simulation results of the Western System Coordinating Council (WSCC) system compare an analytical method with Monte-Carlo simulation (MCS).

• Nuclear Engineering and Technology
• /
• v.36 no.6
• /
• pp.528-539
• /
• 2004

It is essential in commercial reactors that the safety limits imposed on the fuel pellets and fuel clad barriers, such as the linear power density (LPD) and the departure from nucleate boiling ratio (DNBR), are not violated during reactor operations. In order to accurately monitor the safety limits of current reactor states, a detailed three-dimensional (3D) core power distribution should be estimated from the in-core detector signals. In this paper, we propose a calculation methodology for detailed 3D core power distribution, using in-core detector signals and core monitoring constants such as the 3D Coupling Coefficients (3DCC), node power fraction, and pin-to-node factors. Also, the calculation method for several core safety parameters is introduced. The core monitoring constants for the real core state are promptly provided by the core design code and on-line MASTER (Multi-purpose Analyzer for Static and Transient Effects of Reactors), coupled with the core monitoring program. through the plant computer, core state variables, which include reactor thermal power, control rod bank position, boron concentration, inlet moderator temperature, and flow rate, are supplied as input data for MASTER. MASTER performs the core calculation based on the neutron balance equation and generates several core monitoring constants corresponding to the real core state in addition to the expected core power distribution. The accuracy of the developed method is verified through a comparison with the current CECOR method. Because in all the verification calculation cases the proposed method shows a more conservative value than the best estimated value and a less conservative one than the current CECOR and COLSS methods, it is also confirmed that this method secures a greater operating margin through the simulation of the YGN-3 Cycle-1 core from the viewpoint of the power peaking factor for the LPD and the pseudo hot pin axial power distribution for the DNBR calculation.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.327-330
• /
• 2006

This paper represents the numerical analysis on effects of radius ratio in a concentric annulus with a rotating inner cylinder. The numerical model consisted of two cylinder which inner cylinder is rotating and outer cylinder is fix, and the axial direction is used the cyclic condition because of the length for axial direction is assumed infinite. The diameter of inner cylinder is assumed 86.8mm, the numerical parameters are angular velocity and radius ratio. Also, the whole walls of numerical model have no-slip and the working fluid is used water at $20^{\circ}C$. The numerical analysis is assumed the transient state to observe the flow variations by time and the 3-D cylindrical coordinate system. The calculation grid adopted a non-constant grid for dense arrangement near the wall side of cylinder, the standard $k-{\omega}$ high Reynolds number model to consider the effect of turbulence flow and wall, the fully implicit method for time term and the quick scheme for momentum equation. The numerical method is compared with the experimental results by Wereley and Lueptow, and the results are very good agreement. As the results, TVF isn't appeared when Re is small because of the initial flow instability is disappear by effect of the centrifugal force and viscosity. The vortex size is from 0.8 to 1.1 for TVF at various $\eta$, and the traveling distance for wavy vortex have the critical traveling distance for each case.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.941-951
• /
• 2016

The aim of this work is to simulate the thermohydraulic consequences of a main steam line break and to compare the obtained results with Rossendorf Coolant Mixing Model (ROCOM) 1.1 experimental results. The objective is to utilize data from steady-state mixing experiments and computational fluid dynamics (CFD) calculations to determine the flow distribution and the effect of thermal mixing phenomena in the primary loops for the improvement of normal operation conditions and structural integrity assessment of pressurized water reactors. The numerical model of ROCOM was developed using the FLUENT code. The positions of the inlet and outlet boundary conditions and the distribution of detailed velocity/turbulence parameters were determined by preliminary calculations. The temperature fields of transient calculation were averaged in time and compared with time-averaged experimental data. The perforated barrel under the core inlet homogenizes the flow, and therefore, a uniform temperature distribution is formed in the pressure vessel bottom. The calculated and measured values of lowest temperature were equal. The inlet temperature is an essential parameter for safety assessment. The calculation predicts precisely the experimental results at the core inlet central region. CFD results showed a good agreement (both qualitatively and quantitatively) with experimental results.

• Journal of Hydrogen and New Energy
• /
• v.23 no.4
• /
• pp.364-372
• /
• 2012

Transient mixing states of two different fuel oils, dimethylformamide (DMF) oil and JetA1 oil, were investigated by using a color image processing and a neural network. A tank ($D{\times}H$, $310{\times}370mm$) was filled with JetA1 oil. The DMF oil was filled at a top tank, and was mixed with the JetA1 oil in the tank mixing tank via a sudden opening which was performed by nitrogen gas with 1.9 bar. An impeller was rotated with 700 rpm for mixing enhancements of the two fuel oils. To visualize the mixing state of the DMF oil with the JetA1 oil, the DMF oil was coated with Rhodamine B whose color was red. A LCD monitor was used for uniform illumination. The color changes of the DMF oil were captured by a camcoder and the images were transferred to a host computer for quantifying the information of color changes. The color images of two mixed oils were captured with the camcoder. The R, G, B color information of the captured images was used to quantify the concentration of the DMF oil. To quantify the concentration of the DMF oil in the JetA1 oil, a calibration of color-to-concentration was carried out before the main experiment was done. Transient mixing states of DMF oil with the JetA1 oil since after the sudden infiltration were quantified and characterized with the constructed visualization technique.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.16 no.1
• /
• pp.11-21
• /
• 2020

This study aims to develop and advance the evaluation technology for assessing PHWR safety. For this purpose, the complete loss of AC power or station blackout (SBO) was selected as a target accident scenario and the analysis model to evaluate the plant responses was envisioned into the MARS-KS input model. The model includes the main features of the primary heat transport system with a simplified model for the horizontal fuel channels, the secondary heat transport system including the shell side of steam generators, feedwater and main steam line, and moderator system. A steady state condition was achieved successfully by running the present model to check out the stable convergence of the key parameters. Subsequently, through the SBO transient analyses two cases with and without the coolant leakage via the PHTS pumps were simulated and the behaviors of the major parameters were compared. The sensitivity analysis on the amount of the coolant leakage by varying its flow area was also performed to investigate the effect on the system responses. It is expected that the results of the present study will contribute to upgrading the evaluation technology of the detailed thermal hydraulic analysis on the SBO transient of the operating PHWRs.