Search published articles
Showing 55 results for Numerical Model
Volume 0, Issue 0 (8-2024)
Abstract
So far, the performance of masonry walls against in-plane lateral loads such as earthquake loads has been extensively studied, but less attention has been paid to out-of-plane loads such as explosions. Due to their large surface area, walls endure significant forces during explosions, leading to extensive damage and potentially causing severe financial losses and casualties. Given the increase in terrorist and sabotage attacks, reinforcing these structures seems necessary. In recent years, fiber-reinforced polymers (FRP) have been widely and effectively used in the reinforcement and performance improvement of these structures. Their light weight, high stiffness and high strength, and corrosion resistance are among the properties that have attracted researchers to use these materials. Finite element modeling not only provides a basis for better understanding the behavior of masonry walls but also is very useful in predicting the behavior of these members after reinforcement, especially in the absence of experimental results. In this study, using numerical modeling in ABAQUS software, the behavior of masonry walls reinforced with FRP strips against a blast with an explosive charge equivalent to 150 kg of TNT (the weight of explosive likely to be carried in a sabotage attempt via a vehicle) at a distance of 5 meters was investigated. Lagrangian equations were used to model the mechanical behavior of the structure, and the solver used in this research is an explicit solver to account for the time factor in the software’s integration process. The total time considered for the entire explosion process is 1 second, and the explosive load was applied to the studied structure using the Conwep method. The type of fibers, width, thickness, area, and angle of the FRP strips were important and influential parameters that were examined for the efficiency of this reinforcement method. The modeling results indicate that this reinforcement technique is highly effective in strengthening masonry walls against explosions, as it has reduced the deflection of the wall by at least 70% and its energy by up to 90%. It can also be inferred that an arrangement for reinforcing masonry walls with FRP strips is suitable if it covers the areas prone to damage, which in masonry walls are the mortar joints between the bricks. Therefore, the horizontal arrangement shows better performance compared to the vertical and diagonal arrangements. Similarly, reinforcing 100% of the wall area performs much better than reinforcing 50% and 25% of the wall area, but it is not economically acceptable. In general, similar to structural elements, non-structural elements can also exhibit plastic behavior in critical situations, preventing the collapse of these elements due to the absorbed energy. Therefore, walls with higher plastic energy show better behavior against explosive loads. Additionally, based on the hysterical displacement and kinetic energy diagrams of the wall, it can be seen that the behavior of reinforced walls is oscillatory, while the behavior of unreinforced walls is noticeably pulsating. Finally, the optimal arrangement of FRP strips proposed for reinforcing masonry walls against explosions in this study is the use of CFRP strips horizontally, with a thickness of 1 mm, a width of 24 cm, and covering 50% of the surface area; This configuration successfully decreased the deflection of the wall from 63.1 cm to 7.7 cm and damped approximately 13% of the blast wave energy.
Volume 0, Issue 0 (8-2024)
Abstract
The significant wave height is a critical parameter in the design and analysis of marine structures, as well as in their operational use. Consequently, predicting this parameter greatly contributes to improving the design and analysis of marine structures. Various modeling approaches for wave characteristics include numerical, empirical, and artificial intelligence models. This study employs the SWAN model, which is a third-generation model for the simulation and estimation of wave characteristics. Furthermore, soft computing models, including individual and hybrid artificial intelligence models such as Adaptive Neuro-Fuzzy Inference System (ANFIS), Support Vector Machine (SVM), and Emotional Artificial Neural Networks (EANN), have been utilized for wave height prediction, using data from the Amirabad buoy for validation purposes. In this research, the model inputs consist of wind speed, while the outputs are the wave heights. The analysis of the different models was carried out using statistical metrics, including bias, root mean square error, coefficient of variation, and coefficient of determination. The evaluation of the models using these statistics indicates an acceptable agreement between the significant wave heights estimated by the SWAN model and the buoy data. Additionally, each of the three artificial intelligence models mentioned demonstrates a relatively accurate capability in predicting wave height. A comparison of the results from the artificial intelligence models revealed that the Support Vector Machine model exhibited higher accuracy than the others. The Support Vector Machine model serves as an alternative method to the SWAN model or other numerical techniques, enhancing modeling outcomes when wave height data is unavailable or lacks the necessary statistical quality.
Volume 8, Issue 3 (7-2006)
Abstract
The partial differential equations for water flow and solute transport in a two-dimensional saturated domain are rendered discrete using the finite difference technique; the resulting system of algebraic equations is solved using a dynamic programming (DP) method. The advantage of the DP algorithm is that the problem is converted from solving an algebraic system of order NC(NL-1) NC(NL-1) into one of solving a difference equa-tion of order NCNC over NL-1 steps and involving NL-1 matrix inversions of order NCNC. The accuracy and precision of the solutions are shown by comparing the results with an analytical solution and calculation of mass the balance. In addition, the perform-ance of the DP model was compared with the results of the MOC model developed by US Geological Survey. In all cases, the DP model showed good results with sufficient accu-racy.
Akbar Rashidi Ebrahim Hesari, Navid Bazregar, Hossein Farjami,
Volume 9, Issue 1 (1-2020)
Volume 9, Issue 1 (1-2020)
Abstract
Aims: In the present study, a three-dimensional numerical hydrodynamic model with the capability to simulate the diffusion of chemical pollutants released in marine basins was developed and used as a case study to simulate the diffusion of phosphate released by fish culture cages located in the Sisangan marine basin.
Materials & Methods: The equations of the model including momentum equations, continuity of mass equation, temperature, salinity, representative of vertical velocity, the tendency of bottom pressure equations and an extra three-dimensional advection-diffusion equation for simulation of pollutant’s diffusion rewritten in the earth’s spherical coordinates with a vertical Sigma coordinate were solved using finite difference method. To provide the open boundary conditions, the model was used for simulation of wind driven currents in the Caspian Sea from 20 October 2018 to 20 May 2019. For the application of wind field and real geometric condition, it was used the time series of wind fields supplied by ECMWF reanalysis dataset and GEBCO bathymetry with 0.125 degrees resolution and 15 seconds of geographical resolution, respectively.
Findings: Considering the concentration of 17ppb for the phosphate as a concentration of pollutant in the source of the pollution in the location of the fish culture cages, wind induced currents and the diffusion of the phosphate were simulated for 8 months in both horizontal and vertical directions. The results of the simulations were demonstrated and analyzed within the framework of the horizontal surface current, distribution of the phosphate’s concentration in both horizontal and a vertical latitude-depth cross section.
Conclusion: The phosphate’s diffusion is affected by the wind induced currents and after 8 months, it could be extended to the distance of 11, 8.5, 9.5, and 5.7 kilometers far from the cages in eastern, western, southern, and northern directions, respectively. With the generation of vertical velocity and the turbulence effects in the upper layers, phosphate might be diffused in the vertical direction up to 400m depth, as well.
Volume 11, Issue 2 (4-2009)
Abstract
Mathematical simulation of flow toward drains is an important and indispensable stage in drainage design and management. Many related models have been developed, but most of them simulate the saturated flow toward drains without a due consideration of the unsaturated zone. In this study, the two dimensional differential equation governing saturated and unsaturated flow in porous media is numerically solved and water table variations between drains predicted. By introducing and linking a proper optimization model to the numerical one, saturated and unsaturated soil hydrodynamic parameters were estimated within the inverse problem technique context. Data for calibration and verification were provided through a conduction of laboratory experimentation. Other laboratory data were also employed for the proposed model evaluation. The results indicated that in addition to a prediction of the water table variations between drains, the inverse problem model can be employed to estimate the unsaturated soil hydrodynamic parameters with a high degree of precision.
Volume 11, Issue 4 (9-2011)
Abstract
In the current study fume extraction systems are studied numerically and the effect of various parameters as fresh air inlet gap size, fume temperature and composition as well as the dust size is investigated. To aim this goal a precise 3D model of the entire system and proper grid is generated and system of governing equations for the reactive turbulent two-phase flow is solved using a Finite-Volume based code. The results confirm that although increasing the gap size may lead to a reduction in CO volume fraction, but an increase in products temperature is inevitable. Besides, it is shown that, despite the high efficiency of settling chamber in removing the large size particles (greater than 45 micron), it has a poor efficiency in eliminating the smaller size dust particles.
Volume 11, Issue 4 (10-2009)
Abstract
The purpose of this paper is to present a 2D depth-averaged model for simulating and examining unsteady flow patterns in open channel bends. In particular, this paper proposes a 2D depth-averaged model that takes into account the influence of the secondary flow phenomenon through calculation of the dispersion stresses. The dispersion terms which arose from the integration of the product of the discrepancy between the mean and the actual vertical velocity distribution were included in the momentum equations in order to take into account the effect of the secondary current. This model used a time-splitting method for solving advection, diffusion and other momentum equation terms. The proposed model uses an orthogonal curvilinear coordinate system efficiently and accurately to simulate the flow field with irregular boundaries; it also used a finite volume projection method approach for solving the governing equation in a staggered grid. Two sets of experimental data were used to demonstrate the model's capabilities. The comparison of the simulated water surface elevation with the measurements shows good agreement and indicates that inclusion of the dispersion terms improved the simulation results.
Volume 13, Issue 5 (12-2013)
Abstract
Study on the physics of sediment particle movement in micro scale is essential for better understanding sediment transport phenomenon and estimating the rate of sediment transport in rivers and marine environment. Sediment particles basically transport in two modes of bed and suspended load. Bed load takes place through sliding, rolling and saltation. Many parameters influence on this process, which their effects are not fully understood. In this research the influence of the affecting parameters on movement of sediment particles in saltation under unidirectional steady flow are investigated. First, a numerical model is developed to simulate the particle motion in bed load saltation. Then the influencing parameters such as particle shape and its position between other particles, upon the jump length and average velocity of the particles are studied. The result of the study improves our understanding and results in better estimation of sediment transport rate for engineering application.
Volume 13, Issue 14 (3-2014)
Abstract
In this study, a two-dimensional, axisymmetric, computational Algorithm has been developed to simulate the plasma flowfield in a MPD thruster for the purpose of determining the flow behavior and electromagnetic characteristics distribution. The solution employs Roe’s flux vector difference method in combination with Powell’s characteristics-splitting scheme. To ensure the stable high-accuracy solution, new modification of MUSCL technique so called OMUSCL2 method is used. According to being supersonic strong gasdynamic expansion near the electrodes tip, HHT entropy correction is employed. Further improvements to the physical model, such as the inclusion of relevant classical transport properties, a real equation of state, multi-level equilibrium ionization models, anomalous transport, and multi-temperature effects, that are essential for the realistic simulation MPD flows, are implemented. Numerical results of a lab-scale thruster are presented, whereby comparison with experimental data shows good agreement between the predicted and measured enclosed current and electric potential.
Volume 14, Issue 1 (6-2014)
Abstract
Study on the physics of sediment particle movement at grain scale is essential for better understanding sediment transport phenomenon and estimating the rate of sediment transport in rivers and marine environment. Sediment particles basically transport in two modes of bed and suspended load. Bed load takes place through sliding, rolling and saltation, from which the latter is dominant. Many parameters influence on saltation phenomenon, which their effects are not fully understood. These influencing parameters make the saltation a stochastic phenomenon. In the present article the influence of the affecting parameters on movement of sediment particles at saltation mode of transport under unidirectional steady flow are investigated. A numerical model is developed to simulate the particle motion in bed load saltation with considering the main contributor forces. Then the influencing parameters that effect on the jump length and average velocity of the particles are studied. Among them are the initial condition, the particle position between other particles and the shape of particles. The influence of the velocity profile on the jump length and average velocity of the particles are also studied. In summary, the change in the initial condition including the initial velocity and angle produces less than 10% variation on the particle jump length and velocity. On the other hand the position of the grain between the other particles is considerably influential with 40% change in the jump length and average velocity. The particle shape is most important parameter in term of the influence on the jump length and average velocity; there is a 50% difference between the jump length of spherical particles and flake-shape particles, for average velocity it is about 10%. The result of the study improves our understanding of particle motion at grain scale and ultimately results in the better estimation of sediment transport rate.
Volume 14, Issue 2 (5-2014)
Abstract
As the coefficient of performance and the cooling power of adsorption chillers are low, the irreversibility calculation can identify the sources which limit the increase of performance parameters and effectively be used in association with current performance improvement techniques. Adopting the numerical modeling and calculating the temporal distribution of temperature in adsorber elements, this study measures the exergy destruction in different parts and processes of the adsorbent bed. The results show the maximum exergy destruction rate in isosteric phases, yet the total exergy destruction is low due to the short phase times. The highest total exergy loss is observed in isobaric heating phase due to the high irreversibility of desorption process and also long phase duration. Furthermore the effects of fin height and fin spacing on the exergy destruction of adsorbent bed are investigated. The results show that increasing fin height and fin spacing increase the total exergy destruction; however the dependency of fin spacing on exergy destruction is relatively low.
Volume 14, Issue 9 (12-2014)
Abstract
Due to the extreme increase in computational power over the recent years, numerical methods have gained the most proportion in analyzing composite structures and components because of the consideration complicated failure mechanisms such as delamination, fiber buckling and fiber breakage, matrix cracking, debonding ribs of skin and a combination of mentioned failure mechanisms. However exact three - dimensional modeling damages caused by impact phenomena is still a challenge. In present numerical work, the most advanced modeling techniques have been used to predict the behavior of composite structure under high velocity impact. The ribs and layers have been modeled using solid elements and a user defined material model with modified puck and Hashin (3D) failure criteria was implemented. Because these failure criteria do not exist in Commercial version of the Abaqus software, we have used Fortran software for writing these criteria so this capability was added to the software. Figures of velocity variations and force variations of projectile, damaged area, different mechanisms of fracture were reported as results and commented upon. In this study, The numerical results have been validated with experimental data and show very good agreement.
Volume 15, Issue 4 (12-2015)
Abstract
Anchors play a special role in geotechnical structures such as excavations. The anchor section in soil is generally divided into five zones including reinforcement element, grout, grout and surrounding soil mixture, shear zone and soil media. The main objective of the present research is to determine the soil-anchor interaction parameters for numerical modeling of anchored wall using FLAC2D software. Basically, the injection area determining is the main challenge in the anchor force nomination. According to the proposed method, the diameter of the injected area is determined based on the injection pressure, grout volume, porosity and shear zone thickness. It is shown that the diameter of the injected area is approximately increased by 40% relatively to the drilling diameter. The diameter of the injected area in rock media, however, is equal to the drilling diameter. The other parameters are determined using equalization of rock media formulas for soil media. In order to ensure the validity of the proposed method, the pull-out test is numerically simulated in FLAC2D software. The numerical results have been then verified with anchor tension results in an excavation project. The results indicate that ultimate load of anchor calculated from the numerical model is comparable with equations proposed by many researches. Also, there is a negligible difference between the displacement obtained in numerical simulation and pull-out test results. This method is therefore can be used in numerical modeling of anchored wall in soil media with high precision. Anchors play a special role in geotechnical structures such as excavations. The anchor section in soil is generally divided into five zones including reinforcement element, grout, grout and surrounding soil mixture, shear zone and soil media. The main objective of the present research is to determine the soil-anchor interaction parameters for numerical modeling of anchored wall using FLAC2D software. Basically, the injection area determining is the main challenge in the anchor force nomination. According to the proposed method, the diameter of the injected area is determined based on the injection pressure, grout volume, porosity and shear zone thickness. It is shown that the diameter of the injected area is approximately increased by 40% relatively to the drilling diameter. The diameter of the injected area in rock media, however, is equal to the drilling diameter. The other parameters are determined using equalization of rock media formulas for soil media. In order to ensure the validity of the proposed method, the pull-out test is numerically simulated in FLAC2D software. The numerical results have been then verified with anchor tension results in an excavation project. The results indicate that ultimate load of anchor calculated from the numerical model is comparable with equations proposed by many researches. Also, there is a negligible difference between the displacement obtained in numerical simulation and pull-out test results. This method is therefore can be used in numerical modeling of anchored wall in soil media with high precision.
Volume 16, Issue 2 (6-2016)
Abstract
Scour around pier in the flow is an Inevitable issue. Estimation of scour depth and understanding the flow field around pier would help us to design with safer factor. The most important factor of scour around pier is changing of streamelines that leads to a system so called local scour. It consist of two vortices: horseshoe vortices and wake vortex. The obstacle creates downflow jet in front of pier that collide with bed sediments and carry them to the downstream, making the horseshoe vortice. The wake vortex is caused by splitting the streamelines and formation of low pressure flow field region and absorption of flow in rear of pier and pick up the bsd sdiments in this district. In this study we used the numerical model SSIIM as a CFD model to simulate flow and scour pattern Simultaneously around a group piers. This model can be used in hydraulic and environment engineering and has the ability of sediment transport calculation in bed transient movement with temporal dependent as the most important advantage in compare with the other CFD models. The verification of this model was implemented by data and results reported for side by side piers examinations as one ofe the group categorize. In this model we considered the k-ε and k-ω seperately as a turbulence model to solve the eddy viscousity of 3D Navier-Stokes flow equations and use their outputs as inputs of sediment transition equations, we used Power-Law scheme as one of the descritization method of First-Order upstreame scheme to solve the flow and sediment equations on the grids. The pressure term of Navier-Stokes equations in cells was calculated by SIMPLE algorithm which is the First-Order upstreame scheme too. Also by changing the G⁄D distant ratio on the other simulation runs, we generated the diagrams with comparative situation with experimental diagrams. Results in the last time of simulation showed there is much more value of horizontal and vertical velocity between the piers than the other sides. It was 57% of final maximum scour depth in first hour of calculation. Similarly to velocity, The final scour patterns showed there is more scour depth counters between the piers. In details it was deriven that the scour depth pattern was symmetric in early time of calculation, but with time passing it appeared more in the region of between the piers. Although Numerical results show the SSIIM model have calculated the erosion depth in front of piers with high accuracy resulted from good calculation of downflow, comparisons between model results and data show the scour depth pattern that the model calculated the wake vortices behind the piers and Interference the horseshoe vortex between the piers with overestimate value and there are deeper countors of scour depth than experiment diagram. Also the RMS index of scour depth has been calculated in the grid and it represented the values of 0.0353 for k-ε model and 0.0899 for k-ω model. Therefore, the k-ε turbulence model resulted better scour depth pattern calculated in compare with k-ω turbulece model.
Volume 16, Issue 8 (10-2016)
Abstract
Large amount of diesel engine waste heats make researchers design systems that utilize the engine waste heat to provide the cooling demand of the heavy-duty vehicles and improve the engine efficiency. Considerable advantages of adsorption cooling system lead to be nominated for this purpose. Coolant and exhaust gases are the main sources of waste heats of diesel engines and using each of them to drive the adsorption cooling system requires its own equipment and working pair. In this paper, a detailed numerical model has been developed and to examine the performance of the cooling system driven by the coolant waste heat with working pair of silica gel-water and also driven by exhaust waste heat with zeolite13x-water working pair. An identical absorbent bed and ambient conditions have been employed to compare the performance of both systems to identify the more appropriate system. The results show that exhaust driven adsorption cooling system has more capability to meet the vehicle cooling demand. Moreover, the performance of the both adsorption cooling systems were examined under variable ambient condition. Results indicate that increase in ambient temperature leads to almost a linear performance drop in both systems that is more considerable in the coolant- driven adsorption system.
Volume 16, Issue 12 (2-2017)
Abstract
In the present study, for the first time, adsorbent bed of SWS-1L/water adsorption chiller with rectangular and trapezoidal finned flat-tube heat exchanger with has been simulated three dimensionally based on the distributed parameters model and finite volume method. Effects of some important parameters on the chiller performance such as bed averaged pressure, temperature and uptake variations with cycle time have been examined for better understanding of bed dynamic behavior. Also, a comparative study between two different configurations of adsorbent bed including rectangular and trapezoidal fins has been conducted based on identical adsorbent mass. For this purpose, bed temperature, uptake and pressure distributions as well as the vapor flow patterns at the end of heating cycle phases and also effects of fin height and spacing on the system performance have been studied. In this investigation at fixed bed length of 20mm, fin height and spacing variations have been examined in the range of 8-20mm and 3-12mm, respectively. Results indicated that the system performance with rectangular and trapezoidal adsorbent beds are almost similar except for those conditions which fin spacing is 3mm and fin height are 14, 20mm. For the mentioned dimensions, the specific cooling power (SCP) of rectangular beds are almost 5% and 17% (for fin heights of 14 and 20mm, respectively) better than those of trapezoidal beds. Maximum and minimum SCP of adsorption chiller with flat-tube heat exchanger were obtained about 882 and 163W/kg for the smallest and the largest bed geometry and operating conditions considered in this study.
Volume 17, Issue 1 (3-2017)
Abstract
In this paper, flexural behavior of composite sandwich beams under four point bending loading has been studied experimentally and numerically. The skins and the core of the sandwich composite beam have been made of woven glass/epoxy composites and polyvinylchloride foam with 70 kg/m3 density, respectively. The experiments were performed on the beams with different lengths and two different types of layup sequence for the skins as 0/90 and ±45. Failure was initiated in the beams due to indentation of the foam and extended to the face sheet failure under the loading roller. Numerical simulation of the sandwich beam has been performed using ABAQUS commercial software to verify experimental results. During the numerical simulations, the nonlinear material models were employed for shear stress-strain behavior modeling of the foam and the face sheets. In addition, due to the large deformation during bending test geometrical nonlinearity assumption was used in FE analysis. Failure initiation was predicted in the face sheets using modified Hashin criteria. Nonlinear stress analysis and failure predictions in the face sheets and the foam were conducted using USDFLD subroutine in ABAQUS software. Also crushable foam model was employed to simulate the plastic behavior of the foam core. The load-displacement curves and failure mechanisms predicted by the numerical simulations illustrated good correlation with the experimental data.
Volume 17, Issue 3 (5-2017)
Abstract
In the present study, adsorbent bed of an adsorption chiller with finned flat-tube heat exchanger has been simulated three dimensionally based on the heat and mass transfer model with finite volume method. To examine the inter-particle mass transfer resistance effects on the system performance parameters, two different configurations of adsorbent bed including rectangular and trapezoidal fins with identical length and adsorbent mass have been considered and the effects of bed length on the system performance for different fin height and fin pitch have been studied. Moreover, effects of bed length for different particle diameters and also heating source temperatures have been investigated. Results indicated that increasing of bed length (or in the other words increasing of inter-particle resistance) increases and decreases cycle time and specific cooling power, respectively, yet the coefficient of performance is not influenced. Also, increasing bed length reduces the difference between specific cooling power of rectangular and trapezoidal beds if there is any. Moreover it is clear that optimum particles size increase with bed length increase. Finally, it is shown that effect of higher heating fluid temperature on specific cooling power improvement for beds with smaller length is more significant than those with longer length.
Volume 17, Issue 3 (9-2017)
Abstract
In this paper, the dam break phenomena has been simulated in curved rivers using 3D numerical model, Flow-3D. It utilizes the finite volume scheme for structured meshes was used for solving the unsteady Reynolds-averaged Navier-Stokes equations in conjunction with the RNG k-ε closure model. In the utilized software, the Fractional Area/Volume Obstacle Representation (FAVOR) method is used to inspect the geometry in the finite volume mesh. FAVOR appoints the obstacles in a calculation cell with a factional value between 0 to 1 as obstacle fills in the cell. Fluid surface shape is illustrated by volume-of-fluid (VOF) function F(x,y,z,t). With the VOF method, grid cells are classified as empty, full, or partially filled with fluid. Cells are allocated in the fluid fraction varying from zero to one, depending on fluid quantity. The pressure and velocity are coupled implicitly by using the time-advanced pressures and time-advanced velocities in the momentum and continuity equations, respectively. FLOW3D solves these semi-implicit equations iteratively using relaxation techniques. In this paper the GMRES technique has been used as pressure implicit solver. A flux surface is a diagnostic feature in FLOW-3D for computing fluid flow rates. It can be used to obtain time-dependent information about the flow in different parts of the domain. A typical flux surface is a 100% porous baffle with no flow losses, so it does not affect the flow in any way. This feature gives the opportunity to determine the flood hydrograph at various stations downstream of the dam. Effects of curve angle and radious of curvature on the flood wave propagation and unsteady flow features along the curved reach, downstream of the dam has been investigated. Results showed that at the initial instants of the dam break in the straight channel, due to the effects of the dynamic wave, flood hydrographs at the dam location and at a distance downstream of the dam have local peak values, while in the curved chnnel cases, the flood wave becomes unstable immediately after the dam break and the local peak occures just at the dam section. The curved reach decelerate the flood wave propagation compared to the straight channel. Effect of channel curvature on the movement of the flood wave along the inner bank is higher than the outer bank and also the centerline of the curved channel. By decreasing the central radious of the bend, slope of the rising limb of the hydrograph and also the peak discharge, attenuates. Furthermore, the peak discharge time reduces. Unlike to effects of the curvature of the bend, increasing the bend angle does not affect the peak discharge. Changing the bend curvature and curve angle has no effect on the falling limb of the flood hydrograph at various stations downstream of the dam.
Volume 17, Issue 4 (6-2017)
Abstract
Accurate investigation of physical phenomena is one of the important challenges in engineering fields. The present study is about a wet tank which entrance of water is investigated in three cases. When the water wave moves into a tank, complex flow regimes are created. This complexity is mainly associated with different flow mechanisms during the entrance of water and propagation of waves at the bottom bed that should be modelled by means of Navier-Stokes equations with free-surface capability and in 3D phase. Due to complexity and time consuming of Navier-Stokes equations modelling, Shallow water equations are used with the assumption of hydrostatic pressure. First case is about efflux over a wet bed. Second, water influx from the middle top is investigated and then influx from top edges is modelled. A dimensionless number is introduced for each case based on water velocity, gap length and drop height which shows acceptable domain for appropriate compatibility between results. Finally, results of numerical modelling are compared with Navier-Stokes solutions which are obtained from STAR-CD software. Results show admissible compatibility with each other based on observations and inspections.