---

The present paper investigated the capability of various non-linear k–ε models for predicting flow field and pollutant dispersion around a cubical model building with a stack vent located on its roof center within the turbulent boundary layer. One quadratic model proposed by Nisizima and Yoshizawa, and two cubic models, proposed by Lien et al. and Ehrhard and Moussiopoulos were examined by comparing their simulation results with the wind tunnel data and standard k–ε model. All the computations were performed by using the self-developed object-oriented C++ programming in OpenFOAM CFD package, which contains applications and utilities for finite volume solvers. The standard k–ε model provided inadequate results for the flow field, because it could not reproduce the basic flow structures, such as reverse flow on the roof. By contrast, the non-linear models were able to predict anisotropic stresses and correctly showed the dominant stress over the roof to be the streamwise Reynolds stress. The non-linear models were able to predict the concentration field better than the SKE model due to inclusion of the quadratic and cubic terms. Among the RANS models, the Ehrhard model showed the best agreement with the experimental data. It was shown that concentrations predicted by all turbulence models were less diffusive than those of the experiment, although the non-linear k–ε models have reduced this difference.

---

In this article, two-phase slug flow is simulated numerically in a horizontal duct with rectangular cross-section using Volume Of Fluid (VOF) method. Conservation equations of mass, momentum and advection equation are solved in open source OpenFOAM code accompanying k-ω SST turbulence equations. Simulation is conducted based on the experimental results in the duct with rectangular cross-section. The results shows, due to Kelvin-Helmholtz (K-H) instability criteria slug initiation forms in the air-water interface during three dimensional turbulence modeling. Water level was increased slightly at interface in both numerical simulation and experiment. This level increase satisfies the K-H instability to generate a slug at interface. During slug initiation, the pressure behind slug is increased significantly. Big pressure gradient at the beginning of the slug in compare to the end of it causes the slug length to be increased as propagate along the duct. The numerical simulation of present research is capable of predicting the slug length accurately in accordance with experiment; however, the slug position with 22% inaccuracy was obtained. Comparison of the results with the numerical and experimental results of other researchers confirms higher accuracy of flow prediction in the present work.

---

In present study, volume of fluid method in OpenFOAM open source CFD package will be extended to consider phase change phenomena due to condensation process. Both phases (liquid – vapor) are incompressible and immiscible. Vapor phase is assumed in saturated temperature. Interface between two phases are tracked with color function volume of fluid (CF-VOF) method. ُSurface Tension is taken accounted by Continuous Surface Force (CSF) model and mass transfer occurs along interface is considered by Lee mass transfer model. Pressure-Velocity coupling will be solved with PISO algorithm in the collocated grid. This solver is validated with Stefan problem. In one dimensional Stefan problem, the desistance of interface motion from cold wall is compared by the analytical solution. Then condensate laminar liquid film flow over vertical plate is simulated in the presence of gravity. Numerical result shows calculated film thickness from numerical simulation is thinner than analytical solution. Also, it shows Nusselt number is a function of vapor specific heat which neglected in existing correlations, therefore analytical solution and experimental correlation should modified to consider this effect on the Nusselt Number.

---

The fThe flow field around the axisymmtric stream lined bodis which forms the main body of the airplaines and submarines has been the subject of several researches. Turning maneuvers of submarines result in cross flow separation that generates large hydrodynamic forces. The separation of a simple axisymmetric body is very complex in nature. Understanding these vortical flows is paramount to improving vehicle performance and design. A suitable way to reduce the effects of this separated flow is to use vortex generators. The main goal of the present study is to investigate the flow field around a Suboff standard underwater model employing the vortex generator by using the oil flow visualization method and CFD method (OpenFOAM code) at 0° ≤ α ≤ 30° angles of attack. The novelty of the this study is the application of oil flow visualizing method and CFD simulation which can help us to precisely study the structure of three-dimensional vortical flow field. The results show that Vortex Generators placed along the submarine do indeed significantly reduce cross flow separation, size of vortices and drag forces.

---

In this research two-phase slug flow regime in a T-junction branching divider is examined in two regular and irregular groups. Simulation is accomplished by OpenFOAM™ open source software. Simulation uses single fluid with volume of fluid (VOF) method to follow gas-liquid two-phase flow interface. Constant velocity boundary condition for inlet, constant pressure for outlet boundaries and no slip boundary condition are considered for fixed walls. Since slug flow regimes are one of the most complex two-phase flow regimes which its behavior could result in serious damages to the downward equipment's; the present research concentrates on the examination of slug flow behavior in the downstream of the T-junction. This study has concluded that using T junction eliminates flow fluctuation so the pressure and air velocity values decrease. Although the inlet of the vertical branch with cross section of 5×5 cm2 is not fully effective in decreasing upward slugs, but with increasing size of the inlet vertical side-branch from 5×5 cm2 to 10×5 cm2 and 20×5 cm2, pressure value of two-phase flow in the whole duct decreases. The consequences are the slug flow decreases in downstream but the plug flow grows up which means the objectives of the research has been accomplished. To verify the numerical results, comparison was made with the well justified previous works. The agreement was encouraging.

---

Conjugate heat transfer is one of the most important aspects of energy conversion and plays an important role in the thermal efficiency and fuel consumption of chambers. In the present work, a two-dimensional model for reacting flow is presented to calculate transport equations of mass, momentum, energy and species. A new solver is developed for the open-source OpenFOAM software. This new solver is able to predict the conjugate heat transfer effects of reactions and transport processes in fluid and heat conduction in solid as well as radiation in surrounding surface. The coupled method is used and the continuity of temperature and heat flux on the fluid and solid interface is applied in order to analyze conjugate heat transfer through boundary conditions. Experimental data of honeycomb burner is used to validate the new solver. Numerical results are in a good agreement with experimental data. The results show that change of fluid inlet condition and geometry dimensions affect the interaction of conjugate heat transfer and location of released heat of combustion. The location of flame is moved toward outlet as the inlet velocity is increased and toward inlet as the equilibrium ratio is increased. Increasing the length and thickness of solid reduces the preheat area as well.

---

Applications of hollow spherical particles in industry and in thermal spraying process have been developed in recent years. Despite dense droplets, in hollow droplets, the volume changes of the gas play an important role in the dynamics of impact and the shape of the formed splats. In plasma thermal spraying, impact velocities of particles to the surface is in the range of 50 m/s-300 m/s, therefore, changes in pressure and volume of the trapped gas, is important. In this research, impact of hollow droplet on a flat surface and its solidification has been simulated. Volume of fluid model for compressible flows at real thermal spraying condition is used while the impact velocities in the range of 50 m/s-300 is considered. In a few moments after the impact of droplet on the surface, a pressure wave is formed in the air. This wave, increase the vorticity in vicinity of interface of two fluid, which has a great effect on shaping the formed splats. Simulations showed that shape of formed splats vary with velocities in the range of 50 m/s-300 m/s. In higher velocities, the surface of the formed splat is more porous.

---

The present study is subjected to analytical, numerical, and experimental simulation of hydraulic characteristics of flow over the streamlined weirs. Numerical simulations were performed using an open source software namely OpenFoam. According to the objectives of the present study, to evaluate the results of numerical modeling, experimental investigation was conducted, studying different models of streamlined weirs, experimentally. The profiles of the experimental models as well as the simulated numerical models were designed based on the Joukowsky transform function. By analyzing the results of different turbulence models including standard k-ε model, realized k-ε model, RNG k-ε model, k-ω SST model and Reynolds stress LRR model, the k-ω SST model was chosen as the most accurate numerical turbulence model for the simulation of flow over the streamlined weirs. The results of the numerical simulations for different flow discharges and different geometrical characteristics, indicated that, increasing the flow discharge and the relative eccentricity in Joukowsky transform function, tends to increase the velocity and consequently decrease the pressure over the weir crest. Therefore, the lowest pressure and the most probable potential of cavitation belongs to the circular-crested weirs with λ = 1 and high flow discharges. Furthermore, the results show that the greatest bed shear stresses and the compressive forces occur at the downstream end of the circular-crested weirs, thus the downstream zone of the circular-crested weirs is responsible to large values of bed erosion. This is partly due to formation of shock waves, reduction of the flow depth and enhanced velocity of flow downstream of the circular-crested weirs. Furthermore, the lowest bed shear stresses occur at the upstream end of the circular-crested weirs. Therefore, potential of sedimentation upstream of the circular-crested weirs increases. Accordingly, by employing streamlined weirs with λ< 1, and an appropriate curvature, unfavorable flow conditions would be improved, leading to a more safe and economic hydraulic structure. The present study is subjected to analytical, numerical, and experimental simulation of hydraulic characteristics of flow over the streamlined weirs. Numerical simulations were performed using an open source software namely OpenFoam. According to the objectives of the present study, to evaluate the results of numerical modeling, experimental investigation was conducted, studying different models of streamlined weirs, experimentally. The profiles of the experimental models as well as the simulated numerical models were designed based on the Joukowsky transform function. By analyzing the results of different turbulence models including standard k-ε model, realized k-ε model, RNG k-ε model, k-ω SST model and Reynolds stress LRR model, the k-ω SST model was chosen as the most accurate numerical turbulence model for the simulation of flow over the streamlined weirs. The results of the numerical simulations for different flow discharges and different geometrical characteristics, indicated that, increasing the flow discharge and the relative eccentricity in Joukowsky transform function, tends to increase the velocity and consequently decrease the pressure over the weir crest. Therefore, the lowest pressure and the most probable potential of cavitation belongs to the circular-crested weirs with λ = 1 and high flow discharges. Furthermore, the results show that the greatest bed shear stresses and the compressive forces occur at the downstream end of the circular-crested weirs, thus the downstream zone of the circular-crested weirs is responsible to large values of bed erosion. This is partly due to formation of shock waves, reduction of the flow depth and enhanced velocity of flow downstream of the circular-crested weirs. Furthermore, the lowest bed shear stresses occur at the upstream end of the circular-crested weirs. Therefore, potential of sedimentation upstream of the circular-crested weirs increases. Accordingly, by employing streamlined weirs with λ< 1, and an appropriate curvature, unfavorable flow conditions would be improved, leading to a more safe and economic hydraulic structure.

---

In this article, rigid wedge water entry problem under different conditions are evaluated using numerical scheme. It continues to be one of the fundamental issues raised in the hydrodynamics studies and known as a reference for the study of slamming phenomena. The exact calculation of the pressure caused by the slamming phenomenon can be used to analyze the appropriate structural analysis of the ships. In the current study, important variables such as speed and fluid pressure is investigated using computational fluid dynamics method based on the open source OpenFOAM code by numerical solution of the governing equations of tow phase fluid. In order to verify the simulation results obtained from this research, he values of the maximum pressure and t he location and exact time of its occurrence and also pressure coefficient distribution at the impact region have been compared by experimental results of other studies. These investigations have been utilized at different impact velocities and angles. By comparing the numerical results and experimental values, an error was found in the range of 2 to 9%. In addition, variables affecting the pressure applied to the wedge such as water entry velocity and different deadrise angles have been studied.

---

Accurate modeling of the sub-grid scales (SGS) is crucial in determining the accuracy of the large eddy simulations (LES) in turbulent flow analysis. In recent years, new branches of the sub-grid scales models called gradient-based models were developed in computing the sub-grid scales stresses and heat fluxes and used in large eddy simulations. In this work, the modulated gradient model (MGM) equations were implemented in the OpenFOAM package, and pimpleFoam solver was modified to improve the solution accuracy. The modulated gradient model is based on the Taylor-series expansion of the sub-grid scales stress and employs the local equilibrium hypothesis to evaluate the sub-grid scales kinetic energy. To assess the accuracy of the modulated gradient model as well as the improved pimpleFoam solver, turbulent channel flow at a frictional Reynolds number of 395 was simulated via the OpenFOAM package and results were compared with the direct numerical simulation (DNS) data as well as the numerical solution of the Smagorinsky, Dynamic Smagorinsky, Deardorff models. The results show that modulated gradient model evaluates first and second order turbulence parameters with a high-level of accuracy.

---

In the present study, the effects of several global chemical kinetics in 3-dimensional numerical simulation of methane combustion in a horizontal combustion chamber which has lifted flame by a set of open source code OpenFOAM, is compared. The purpose of this comparison is to study the effects of 1, 2 and 4 step global kinetics on velocity, temperature and species distribution. In this simulation, conservation and state equations are solved simultaneously. Partial differential equations are discreted by finite volume method. The effects of turbulence by standard k-e, radiation by P1 model and turbulence-combustion interaction by PaSR are modeled. The results of numerical simulations have been validated by a cylindrical combustion chamber experimental data. The results show that the kinetics have considerable differences in results of velocity, temperature and species in the final third of the chamber where the flame is located, and differently predict locations of the flame. According to these results, 4-step mechanisms were more accurate than the 2-step type. Between 4 step mechanisms, JL is more accurate than Kim in overall; However, its calculation time is higher than the Kim. Single step kinetics were not able to keep the lifted flame.Towards the experimental results, 2-step model predicts the flame in downward and Kim mechanism estimates the flame in the upward.

---

In the present research, numerical simulation of the characteristic chart and steady-state Wakefield flow around a marine propeller is conducted. Solutions were performed using the open-source OpenFOAM software and the steady incompressible simple-Foam solver. The gradients were calculated using the linear Gauss algorithm, and the pressure equation was solved with the multi-grid method. In this research, characteristic chart simulation of the propeller was carried out for the entire operational conditions and the effect of using Realizable-k-ε and k-ε-v^2-f turbulence models on the results was investigated. The results were found to be in good agreement in all conditions except for the near bollard region. In this region, the propeller inlet angle of attack severely increased, and the two equation model predicted the thrust coefficient with 24% error, while implementing the four equation model significantly developed the results and decreased the error to 5%. The wake region parameters were also investigated in the numerical simulations at different longitudinal and radial cross sections behind the propeller which showed good agreement compared with the available experimental data. Wake region investigation showed that the flow behavior in downstream cross sections is similar to the corresponding upstream section with smaller variation ranges and for the swirling flow behind the propeller, the maximum and minimum angular position of the wake components rotates. The obtained results also show that the wake axial velocity component deviation is extremely large at the blade tip.

---

“Minimum-dissipation sub-grid models” are simple alternatives to the Smagorinsky-type approaches to imposing sub-grid scales (SGS)' effects in the large-eddy simulation (LES) approach. Recently, a new model in this family called “anisotropic minimum-dissipation (AMD)” model is represented. AMD is classified as a static type eddy-viscosity sub-grid scale model. The model is more cost effective than the dynamic Smagorinsky model, furthermore; it is not only able to consider the effect of various directions in computing sub-grid stress but also capable of operating for transitional flows from laminar to turbulent. In this study, this sub-grid model has been implemented in the open source package OpenFOAM and its performance is evaluated in the prediction of the flow field inside a channel with a pressure driven air flow. The accuracy of the model has been investigated at different Reynolds numbers including transient and fully turbulent flows and compared with the dynamic Smagorinsky model as well as direct numerical simulation (DNS) solutions. Results reveal that this sub-grid model is quite accurate over a broad range of Reynolds numbers once calculating velocity profiles as well as first and second-order turbulent quantities.

---

In this paper, 3D investigation has been employed to study the wake instability of viscoelasic fluid flow behind unconfined sphere. For estimating the proper properties of the viscoelastic fluid in this study a non-linear Giesekus model is used as the constitutive equation of viscoelastic fluid. Numerical computations are carried out by solving the governing and the onstitutive equations of the viscoelasic fluid flow using the finite volume technique and OpenFOAM which is an open source code is used as the CFD solver. At first velocity field and flow streamlines of Newtonian fluid around the sphere for various Reynolds numbers have been plotted and by plotting the velocity magnitude and pressure at a point behind the sphere versus time, the value of Recr in which the flow become unstable has been reported. Furthermore, for validating the present numerical code, variation of drag coefficient around the sphere versus Reynolds number has been compared with previous investigations. In the following, the effect of Reynolds and Wisenberg number on fluid flow and instability of wake formation behind a sphere have been investigated at high values of Reynolds number for the first time. Results show that at high values of Reynolds number the effect of Wisenberg number has less effect in contrast with Reynolds number on flow instability behind the sphere.

---

---

The motion of the liquid free surface in a container (sloshing phenomenon) inserts a momentum on the container walls. This makes a great disorder in the movement of the carrier vehicle or inserts a large force and momentum on the container walls. The reason for this phenomenon is the establishment of destructive waves and hydrodynamic forces. The side effects of this phenomenon in various industries, such as ship industries carrying liquid fuels, liquid fuel rocket industries, fuel tanks or water tanks, increase the importance of predictions of the behaviors of this phenomenon. One way of controlling is to use baffles or plates in the transverse direction of the tank. In this study, the governing equations on this phenomenon have been solved using the OpenFOAM software. This software solves partial differential equations using the finite volume method, which by default considers geometry to be three dimensional. In order to solve the two-phase flow, a modified volume of the fluid model (VOF) is applied and the moving mesh model is used for the movement of the container body. In the VOF method, the phases are expressed as a fraction of one (volume fraction). To determine this parameter, based on the continuity equation, a differential equation is regulated and solved. For the turbulent flow model, a modified k-e model is used by considering the effects of free-surface flows. Also, an experimental model of a real moving liquid container has been used for validation of the predictions of the presented simulation. The results show that the experimental and numerical results are in good accordance. In addition, the results show that using vertical baffles up to 50% can reduce the fluctuations caused by this phenomenon.

---

Roll motion is one of the most important and dangerous motions of the ship and can even result in capsizing of the vessel. Therefore, its control has been always of interest for marine industry researchers. Among the various methods and equipment for controlling the roll motion, the use of free surface anti-roll tanks has been one of the most important methods and which used in many cases due to its simplicity in construction and design. The high efficiency of these tanks at all speeds and even without speed is another strength of these tanks. This study investigates the effect of the free surface anti-roll tank on the roll motion numerically and experimentally. In the numerical simulation, a CFD sloshing solver, based on the “Open source Field Operation And Manipulation”, known in short as Open-FOAM, and assuming 2D laminar flow conditions, is customized to calculate the sloshing loads from the tank. The predicted roll damping and moments due to the anti-roll tank are validated against experimental results. This simulator could be used as a sloshing simulator to couple with seakeeping solvers.

---

Numerical methods as one of the subcategories of theoretical models can predict the behavior of energetic materials with appropriate accuracy and away of experimental tests limitations. In this investigation, computational fluid dynamics tool has been used to predict the blast wave propagation with Consideration of geometrical obstacles. Two solvers (extendedSonicFoam and blastFoam) from the open source technology module, OpenFOAM  have been used for simulations and To enhance confirmation with reality, large eddy simulation method was employed for turbulence modeling. In addition to the ideal gas equation of state (EOS), the BKW EOS, which is a complete EOS with an explicit temperature dependence, have been used to correlate the various thermodynamic parameters. Several gauges were positioned to record the pressure-time signals and the experimental data reported in the resources were used for validation. It should be noted that the maximum error of simulations was 12.29% for different blast wave parameters. deviation from standard for ideal gas numerical results was greater than that of real gas assumption and blastFoam solver has been predicted maximum positive phase overpressure, arrival time and positive phase impulse, which are the important parameters of blast wave, with less error in comparison to extendedSonicFoam solver.