• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.390-403
• /
• 2017

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

• Journal of international Conference on Electrical Machines and Systems
• /
• v.1 no.1
• /
• pp.135-142
• /
• 2012

In this paper, optimal selection of the unified power flow controller (UPFC) damping controller parameters in order to improve the power system dynamic response and its stability based on two modified intelligent algorithms have been proposed. These algorithms are based on a modified intelligent particle swarm optimization (PSO) and continuous genetic algorithm (GA). After extraction of UPFC dynamic model, intelligent PSO and genetic algorithms are used to select the effective feedback signal of the damping controller; then, to compare the performance of the proposed UPFC controller in damping the critical modes of a single-machine infinite-bus (SMIB) power system, the simulation results are presented. The comparison shows the good performance of both presented PSO and genetic algorithms in an optimal selection of UPFC damping controller parameters and damping oscillations.

• Membrane and Water Treatment
• /
• v.9 no.2
• /
• pp.129-136
• /
• 2018

In this contribution, the control of multivariable reverse osmosis (RO) desalination plant using proportional-integral-derivative (PID) controllers is presented. First, feed-forward compensators are designed using simplified decoupling method and then the PID controllers are tuned for flux (flow-rate) and conductivity (salinity). The tuning of PID controllers is accomplished by minimization of the integral of squared error (ISE). The ISEs are minimized using a recently proposed algorithm named as teacher-learner-based-optimization (TLBO). TLBO algorithm is used due to being simple and being free from algorithm-specific parameters. A comparative analysis is carried out to prove the supremacy of TLBO algorithm over other state-of-art algorithms like particle swarm optimization (PSO), artificial bee colony (ABC) and differential evolution (DE). The simulation results and comparisons show that the purposed method performs better in terms of performance and can successfully be applied for tuning of PID controllers for RO desalination plants.

• Resources Recycling
• /
• v.31 no.6
• /
• pp.36-43
• /
• 2022

The general screen used to separate the particle size of recycled aggregate has restrictions when dealing with moisturized materials because of the blockage phenomenon. Therefore, in this study, to improve the separation efficiency of the conventional screen, the excellence of additional vibrating device based on spring was decided by a simulation experiment based on the discrete element method (DEM). The motion characteristic was investigated by analyzing the displacement, amplitude, and strain angle based on the spring design. Further, the particle motion was simulated by applying spring motion. The material flow and separation efficiency of the screen applied spring were confirmed as 9.2 kg/s and 97 %, respectively. Consequently, the improvement in the screen applied with blockage prevention spring was confirmed by comparing with the conventional screen.

• Structural Engineering and Mechanics
• /
• v.86 no.5
• /
• pp.607-619
• /
• 2023

This study, it was tried to evaluate the asphalt behavior under tensile loading conditions through indirect Brazilian and direct tensile tests, experimentally and numerically. This paper is important from two points of view. The first one, a new test method was developed for the determination of the direct tensile strength of asphalt and its difference was obtained from the indirect test method. The second one, the effects of particle size and loading rate have been cleared on the tensile fracture mechanism. The experimental direct tensile strength of the asphalt specimens was measured in the laboratory using the compression-to-tensile load converting (CTLC) device. Some special types of asphalt specimens were prepared in the form of slabs with a central hole. The CTLC device is then equipped with this specimen and placed in the universal testing machine. Then, the direct tensile strength of asphalt specimens with different sizes of ingredients can be measured at different loading rates in the laboratory. The particle flow code (PFC) was used to numerically simulate the direct tensile strength test of asphalt samples. This numerical modeling technique is based on the versatile discrete element method (DEM). Three different particle diameters were chosen and were tested under three different loading rates. The results show that when the loading rate was 0.016 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis till coalescence to the model boundary. When the loading rate was 0.032 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis. The branching occurs in these cracks. This shows that the crack propagation is under quasi-static conditions. When the loading rate was 0.064 mm/sec, mixed tensile and shear cracks were initiated below the loading walls and branching occurred in these cracks. This shows that the crack propagation is under dynamic conditions. The loading rate increases and the tensile strength increases. Because all defects mobilized under a low loading rate and this led to decreasing the tensile strength. The experimental results for the direct tensile strengths of asphalt specimens of different ingredients were in good accordance with their corresponding results approximated by DEM software.

• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.347-360
• /
• 1998

Numerical simulations for fluid flow and solute transport in a fracture rock masses are performed by using a transient flow model, which is based on the three-dimensional stochastic and discrete fracture network model (DFN model) and is coupled hydraulic model with mechanical model. In the numerical simulations of the solute transport, we used to the particle following algorithm which is similar to an advective biased random walk. The purpose of this study is to predict the response of the tracer test between two deep bore holes (GPK1 and GPK2) implanted at Soultz sous Foret in France, in the context of the geothermal researches.l The data sets used are obtained from in situcirculating experiments during 1995. As the result of the transport simulation, the mean transit time for the non reactive particles is about 5 days between two bore holes.

• Clean Technology
• /
• v.25 no.3
• /
• pp.206-222
• /
• 2019

Recently, various multiphase flows have been developed, and among them some models have been commercialized. However, most of them have been developed based on a pressure-based approach; therefore, various numerical difficulties were involved inherently. Accordingly, in order to overcome these numerical difficulties, a multiphase flow model, MultiPhaSe flow (MPS), following a fractional-flow based approach was developed. In this study, by combining a contaminant transport module describing an enhanced dissolution effect of a surfactant with MPS, a MultiPhaSe flow and TranSport (MPSTS) model was developed. The developed model was verified using the analytical solution of Clement. The MPSTS model can simulate the process of surfactant enhanced aquifer remediation including interphase mass transfer and contaminant transport in multiphase flow by using the coupled particle tracking method and Lagrangian-Eulerian method. In this study, a surfactant was used in a non aqueous phase liquid (NAPL) contaminated area, and the effect of hydro-geological heterogeneity in the layered media on remediation efficiency was studied using the developed model. According to the numerical simulation, when hydraulic conductivity in a lower layer is 10 times, 20 times, and 50 times larger than that in an upper layer, the concentration of dissolved diesel in the lower layer is much higher than that in the upper layer because the surfactant moves faster along the lower layer owing to preferential flow; thus, the surfactant enhances dissolution of residual non aqueous phase liquid in the lower layer.

• Journal of Korean Society of Water and Wastewater
• /
• v.27 no.3
• /
• pp.383-394
• /
• 2013

A sedimentation basin is used to remove suspended sediments which can cause abrasive and erosive wear on hydraulic turbines of hydropower plants. This sediment erosion not only decreases efficiency of the turbine but also increases maintenance costs. In this study, the three-dimensional numerical simulations were carried out on the overseas hydropower project. The simulations of flow and suspended sediment concentration were obtained using FLOW-3D computational fluid dynamics code. The simulations provide removal efficiency of a sedimentation basin based on particle sizes. The influence of baffles on the flow field and the removal efficiency of suspended sediments in the sedimentation basin has been investigated. This paper also provides the numerical simulations for sediment-induced density currents that may occur in the sedimentation basin. The simulation results indicate that the formation of density currents decreases the removal efficiency. When a baffle is installed in the sedimentation basin, the baffle provides intensive settling zones resulting in increasing the sediments settling. Thus the enhanced removal efficiency can be achieved by installing the baffle inside the sedimentation basin.

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.68-80
• /
• 1995

A series of parametric calculations have been performed in order to investigate the short-term and short-range plume and puff behavior of toxic gaseous and solid pollutant dispersion in an open atmosphere. The simulation is made by the use of the computer program developed by this laboratory, in which a control-volume based finite-difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling appeared In Wavier-Stokes equation. The Reynolds stresses are solved by the standard two-equation k-$\varepsilon$ model modified for buoyancy together with the RNG(Renormalization Group) k-$\varepsilon$ model. The major parameters considered in this calculation are pollutant gas density and temperature, the relative velocity of pollutants to that of the surrounding atmospheric air, and particulate size and density together with the height released. The flow field is typically characterized by the formation of a strong recirculation region for the case of the low density gases such as $CH_4$ and air due to the strong buoyancy, while the flow is simply declining pattern toward the downstream ground for the case of heavy molecule like the $CH_2C1_2$and $CCl_4$, even for the high temperature, $200^{\circ}C$. The effect of gas temperature and velocity on the flow field together with the particle trajectory are presented and discussed in detail. In general, the results are physically acceptable and consistent.

• Ocean Systems Engineering
• /
• v.10 no.4
• /
• pp.359-371
• /
• 2020

The performance of a multi-level overtopping wave energy converter (OWEC) has been numerically and experimentally investigated in a two-dimensional wave tank in order to study the effects of opening width of additional reservoirs. The device is a fixed OWEC consisting of an inclined ramp together with several reservoirs at different levels. A particle-based numerical simulation utilizing the Lattice Boltzmann Method (LBM) is used to simulate the flow behavior around the OWEC. Additionally, an experimental model is also built and tested in a small wave flume in order to validate the numerical results. A comparison in energy captured performance between single-level and multi-level devices has been proposed using the hydraulic efficiency. The enhancement of power capture performance is accomplished by increasing an overtopping flow rate captured by the extra reservoirs. However, a noticeably large opening of the extra reservoirs can result in a reduction in the power efficiency. The overtopping flow behavior into the reservoirs is also presented and discussed. Moreover, the results of hydrodynamic performance are compared with a similar study, of which a similar tendency is achieved. Nevertheless, the LBM simulations consume less computational time in both pre-processing and calculating phases.