In the case of a severe accident, natural convection plays an important role in the atmosphere mixing of nuclear reactor containments. In the previous study, to simulate the natural convection in the accident scenario within a nuclear reactor containment, the steady thermal boundary conditions (BCs) were prescribed on either cooled or heated wall. The present study, therefore, aims at the transient 3D numerical simulations of natural convection of air around a cylindrical containment with unsteady thermal BCs at the vessel wall. For that purpose, the experiment series was done in the CIGMA facility at Japan Atomic Energy Agency (JAEA). The upper vessel or both the upper vessel and the middle jacket was cooled by subcooled water, while the lower vessel was thermally insulated. A 3D model was simulated with OpenFOAM®, applying the Unsteady Reynolds-averaged Navier–Stokes equations (URANS) model. Different turbulence models were studied, such as the standard k-ε, standard k-ω, k-ω Shear Stress Transport (SST) and, low-Reynolds-k-ε Launder-Sharma. The results of the four turbulence models were compared versus the results of experimental data. The k-ω SST showed a better prediction compare to other turbulence models. Also, the accuracy of the predicted temperature and pressure were improved when the heat conduction on the internal structure, i.e., flat bar, was considered in the simulation. Otherwise, the predictions on both temperature and pressure were underestimated compared with the experimental results. Hence, the conjugate heat transfer in the internal structure inside the containment vessel must be modeled accurately.

In this manuscript, an automated framework dedicated to design space exploration and design optimization studies is presented. The framework integrates a set of numerical simulation, computer-aided design, numerical optimization, and data analytics tools using scripting capabilities. The tools used are open-source and freeware, and can be deployed on any platform. The main feature of the proposed methodology is the use of a cloud-based parametrical computer-aided design application, which allows the user to change any parametric variable defined in the solid model. We demonstrate the capabilities and flexibility of the framework using computational fluid dynamics applications; however, the same workflow can be used with any numerical simulation tool (e.g., a structural solver or a spread-sheet) that is able to interact via a command line interface or using scripting languages. We conduct design space exploration and design optimization studies using quantitative and qualitative metrics, and to reduce the high computing times and computational resources intrinsic to these kinds of studies, concurrent simulations and surrogate-based optimization are used.

High-Rise buildings with their particular features can affects on surrounding environment and makes new microclimates. In the windy conditions, the spaces that are between building blocks changes to passages and affects on the wind velocity, intensity and it’s other parameters.The importance of this effect is different in each level of building height. The Pedestrian-Level is the lowest and one of important areas. Markets, playgrounds and pedestrian access had located in this area and any unwanted microclimate changes like high velocity and turbulence in this level can makes discomfort and dangerous condition for residents. So this research tries to consider the pedestrian- level wind comfort in some High-Rise building complexes arrangement that had located in Tehran district 22 with Computational Fluid Dynamics (CFD) modeling and reaching to a suitable arrangement pattern. It had collected the required data through field study and librarian databases and then compared them with standard guidelines and analyzed them by comparative comparison method. As a result a linear arrangement that placed crossover to wind direction for providing wind comfort and preventing wind danger is been suggested in this region.

In order to perceive the behavior presented by the multiphase chemical reactors, the ant colony optimization algorithm was combined with computational fluid dynamics (CFD) data. This intelligent algorithm creates a probabilistic technique for computing flow and it can predict various levels of three-dimensional bubble column reactor (BCR). This artificial ant algorithm is mimicking real ant behavior. This method can anticipate the flow characteristics in the reactor using almost 30 % of the whole data in the domain. Following discovering the suitable parameters, the method is used for predicting the points not being simulated with CFD, which represent mesh refinement of Ant colony method. In addition, it is possible to anticipate the bubble-column reactors in the absence of numerical results or training of exact values of evaluated data. The major benefits include reduced computational costs and time savings. The results show a great agreement between ant colony prediction and CFD outputs in different sections of the BCR. The combination of ant colony system and neural network framework can provide the smart structure to estimate biological and nature physics base phenomena. The ant colony optimization algorithm (ACO) framework based on ant behavior can solve all local mathematical answers throughout 3D bubble column reactor. The integration of all local answers can provide the overall solution in the reactor for different characteristics. This new overview of modelling can illustrate new sight into biological behavior in nature.

Accurate prediction of the self-propulsion performance is one of the most important factors for energy-efficient design of a ship. In general, the hydrodynamic performance of a full-scale ship could be achieved by model-scale simulation or towing tank test with extrapolations. With the development of CFD methods and computing power, directly predict ship performance with full-scale CFD is an important approach. In this article, a numerical study on the full-scale self-propulsion performance with propeller operating behind ship at model- and full-scale is presented. The study includes numerical simulations using RANS method with double-model and VOF model respectively and scale effect analysis based on overall performance, local flow fields and detailed vortex identification. Verification study on grid convergence is also performed for full-scale simulation with global and local mesh refinements. And a series of sea trail tests were performed to collect reliable data for the validation of CFD predictions. The analysis of scale effect on hull-propeller interaction shows that the difference on hull boundary layer and flow separation is the main source of scale effect on ship wake. And the results of the fluctuations of propeller thrust and torque along with circulation distribution and local flow field show that propeller’s loading is significantly higher for model-scale ship. It is suggested that the difference on vortex evolution and interaction is more pronounced and have larger effects on ship’s powering performance at model-scale than full-scale according to the simulation results. From the study on self-propulsion prediction, it could be concluded that the simplification on free surface treatment does not only affect the wave-making resistance for bare hull but also the propeller performance and propeller induced ship resistance which can produced up to 5% uncertainty to the power prediction. Roughness is another important factor in full-scale simulation because it has up to approximately 7% effect on the delivery power. As a result of validation study, the numerical simulations of full-scale ship self-propulsion shows good agreement with the sea trail data especially for cases that have considered both roughness and free surface effects. This result will largely enhance our confidence to apply full-scale simulation in the prediction of ship’s self-propulsion performance in the future ship designs.

Based on K-ε Standard Wall turbulence model (2-Equation) and Navier-Stokes (N-S) equations defined for incompressible fluids, fluid flow behaviour around hyperloop pods in an evacuated tube was simulated using ANSYS fluent solver assuming steady state and two dimensional conditions. In this research, to develop the case studies, using combination of different head and tail shape profile, four kind of hyperloop pods were developed with the aid of SolidWorks. These four pods have been investigated for their aerodynamic behaviour as four different case scenarios. The results of simulation depicts that an atmospheric pressure of 100 Pa with blockage ratio of 0.36 in tube provides the best possible aerodynamic behaviour for the designed hyperloop pod models. This research finds that overall aerodynamic behaviour of hyperloop pods can be varied by changing the head and tail shape profile of pods and a particular combination of head and tail shape profile can provide optimally best aerodynamic capabilities. Thus, this research paper provides a novel method of obtaining best aerodynamic capabilities in hyperloop pods by designing head profile optimally in combination with tail profile. This outcome will provide major contribution towards the development of hyperloop pods in future with better aerodynamic behaviour resulting in lesser electrical energy required to propel the hyperloop pods in evacuated tube.

This paper analyses the health and performance of 12-stage axial compressor of the TV3-117VM/VMA turboshaft operated in a desert environment. The results of the dimensional control of 4,800 worn blades are analysed to model the wear process. Operational experience and numerical simulations are used to assess the effectiveness of an Inlet Particle Separator. Numerical modal analysis is performed to generate the Campbell diagram of worn blades and identify resonant blade vibration which can lead to high cycle fatigue (HCF). It is shown that the gradual loss of the stall margin over time determines the serviceability limits of compressor blades. Recommendations setting out go / no-go criteria are made to maintenance and repair organisations.

Efficiency improvement is the new challenge in all fields of design. In this scenario the reduction of power losses is becoming more and more a main concern also in the design of power transmissions. Appropriate models to predict power losses are therefore from the earliest stages of the design phase. The aim of the project is to carry on lubrication simulations of several variants of a cylindrical-roller-bearing to understand the lubricant distribution and the related churning power losses. Several strategies to reduce the computational effort have been used. Among them the sectorial symmetry and three innovative meshing strategies (purely analytical with and without interfaces and analytical/subtractive) that have been implemented in the OpenFOAM® environment. The results of the different approaches were compared among them and with experimental observations showing good agreement and reasonable savings in terms of computational effort.

Among the different applicable irrigants for root canal disinfection, sodium hypochlorite 5.25% is one of the most attractive ones. The quality of root canal disinfection is dependent on some factors such as the employed approach, type of flow rate of irrigant and the size of needle. The majority of studies in the field of root canal disinfection are experimentally carried out. In the current article, Computation Fluid Dynamic (CFD) is used for modeling the antimicrobial liquid flow in the root canal and evaluate the effects of needle size and flow rate. Two needles, G28 and G30, are used for irrigation in three volumetric rates of flow including 0.10 mL⁄s , 0.20 mL⁄s and 0.30 mL⁄s. The results of numerical simulations revealed the improved quality of root canal disinfection by augmentation in the rate of flow and decrease in the inner diameter of the needle. According to the outcomes of the modeling, the highest average wall shear stress obtained in the case of using G28 needle and 30 mL⁄s flow rate, which was approximately 10.21 Pa.

In this study, the field measurement-based validation of a novel sediment transport calculation method is presented. River sections with complex bed topography and inhomogeneous bed material composition highlight the need for an improved sediment transport calculation method. The complexity of the morphodynamic features can result in the simultaneous appearance of the gravel and finer sand dominated sediment transport (e.g. parallel bed armoring and siltation) at different regions within a shorter river reach. For the improvement purpose of sediment transport calculation in such complex river beds, a novel sediment transport method was elaborated. The base concept of it is the combined use of two already existing empirical sediment transport models. The method was already validated against laboratory measurements. The major goal of this study is the verification of the novel method with a real river case study. The combining of the two sediment transport models is based on the implementation of a recently presented classification method of the locally dominant sediment transport nature (gravel or sand transport dominates). The results are compared with measured bed change maps. The verification clearly refers to the meaningful improvement in the sediment transport calculation by the novel manner in case of spatially varying bed content.

Intelligent algorithms are recently used in the optimization process in chemical engineering and application of multiphase flows such as bubbling flow. This overview of modeling can be a great replacement with complex numerical methods or very time-consuming and disruptive measurement experimental process. In this study, we develop the adaptive network-based fuzzy inference system (ANFIS) method for mapping inputs and outputs together and understand the behavior of the fluid flow from other output parameters of the bubble column reactor. Neural cells can fully learn the process in their memory and after the training stage, the fuzzy structure predicts the multiphase flow data. Four inputs such as x coordinate, y coordinate, z coordinate, and air superficial velocity and one output such as pressure gradient are considered in the learning process of the ANFIS method. During the learning process, the different number of the membership function, type of membership functions and the number of inputs are examined to achieve the intelligent algorithm with high accuracy. The results show that as the number of inputs increases the accuracy of the ANFIS method rises up to almost for all cases, while the increment in the number of rules has a effect on the intelligence of artificial algorithm. This finding shows that the density of neural objects or higher input parameters enables the moded for better understanding. We also proposed a new evaluation of data in the bubble column reactor by mapping inputs and outputs and shuffle all parameters together to understand the behaviour of the multiphase flow as a function of either inputs or outputs. This new process of mapping inputs and outputs data provides a framework to fully understand the flow in the fluid domain in a short time of fuzzy structure calculation.

The open source CFD code OpenFOAM has emerged as one of the most popular alternative to commercial CFD solvers. The recent version of OpenFOAM supports overlapping grids, so called Overset mesh. In this type of mesh, one or more sub-collection of control volumes (or cells) are allowed to overlap with other set of cells. This allows for great flexibility in modelling Fluid Structure Interaction (FSI) problems. Fluid mesh can be standard structured mesh of a rectangular/cuboidal/or cylindrical domain. The structure mesh can be generated separately by creating multiple layers around a given body. The overset mesh functionality allows for overlapping or immersing this structured mesh inside the fluid mesh and a relevant FSI or moving dynamics problems can be solved. The idea of overset mesh has been around since eighties but its support in OpenFOAM is very recent. Most CFD codes which support overset mesh have either been in-house or commercial CFD codes. OpenFOAM's support for overset is the first major open source code resource available for CFD problems. The aim of this research is to solve a classical benchmark airfoil problem using OpenFOAM overset mesh and compare the numerical results with experimental result. We report here flow simulations results around NACA 0018. The result obtained from overset mesh compares convincingly well with experimental results. Computations have been carried out for Reynolds numbers in the range of 10⁵ with angle of attack ranging from α= 5 degree to α= 30 degree with an interval of 5 degree. Turbulence is incorporated using k - ε turbulence model. This study helps developing confidence in using OpenFOAM overset support for more complicated flows and moving dynamics. This report is complemented with a brief description of finite volume discretization using overset mesh.

Traditional winglets are designed as fixed devices attached at the tips of the wings. The primary purpose of the winglets is to reduce the lift-induced drag, therefore improving aircraft performance and fuel efficiency. However, because winglets are fixed surfaces, they cannot be used to control lift-induced drag reductions or to obtain the largest lift-induced drag reductions at different flight conditions (take-off, climb, cruise, loitering, descent, approach, landing, and so on). In this work, we propose the use of variable cant angle winglets which could potentially allow aircraft to get the best all-around performance (in terms of lift-induced drag reduction), at different flight phases. By using computational fluid dynamics, we study the influence of the winglet cant angle and sweep angle on the performance of a benchmark wing at Mach numbers of 0.3 and 0.8395. The results obtained demonstrate that by adjusting the cant angle, the aerodynamic performance can be improved at different flight conditions.

In this study, the multi-relaxation-time lattice Boltzmann method is applied to investigate the oscillatory instability of lid-driven flows in two-dimensional semi-elliptical cavities with different vertical to horizontal aspect ratios K in the range of 1.0--3.0. The program implemented in this study is parallelized using CUDA (compute unified device architecture), a parallel computing platform, and computations are carried out on NVIDIA Tesla K40c GPU. To carry out precise calculations, the CUDA algorithm is extensively investigated, and its parallel efficiency indicates that the maximum speedup is 47.6 times faster. Furthermore, the steady--oscillatory Reynolds numbers are predicted by implementing the CUDA-based programs. The amplitude coefficient is defined to quantify the time-dependent oscillation of the velocity magnitude at the monitoring point. The simulation results indicate that the transition Reynolds numbers correlate negatively with the aspect ratio of the semi-elliptical cavity, and are smaller than those of the rectangular cavity at the same aspect ratio. In addition, the detailed vortex structures of the semi-elliptical cavity within a single period are also investigated when the Reynolds number is larger than the steady--oscillatory value to determine the effects of periodic oscillation of the velocity magnitude.

A common reactor type in the chemical and process industry is the fixed-bed reactor. Accurate modeling can be achieved with particle-resolved Computational Fluid Dynamic (CFD) simulations. However, the underlying bed morphology plays a paramount role. Synthetic bed-generation methods are much more flexible and faster than image-based approaches. In this study, we look critically at the two different bed generation methods: Discrete Element Method (DEM) (in the commercial software STAR-CCM+) and the rigid-body model (in the open-source software Blender). The two approaches are compared in terms of synthetically generated beds with experimental data of overall and radial porosity, particle orientation, as well as radial velocities. Both models show accurate agreement for the porosity. However, only Blender shows similar particle orientation than the experimental results. The main drawback of the DEM is the long calculation time and the shape approximation with composite particles.

The combination of artificial intelligence algorithms and numerical methods has recently become popular in the prediction of macroscopic and microscopic hydrodynamics parameters of bubble column reactors. The multi inputs and outputs machine learning can cover small phase interactions or large fluid behavior in industrial domains. This numerical combination can develop the smart multiphase bubble column reactor with the ability of low-cost computational time. It can also decrease case studies for the optimization process when big data is appropriately used during learning. There are still many model parameters that need to be optimized for a very accurate artificial algorithm, including data processing and initialization, the combination of inputs and outputs, number of inputs and model tuning parameters. For this study, we aim to train four inputs big data during learning process by an adaptive neuro-fuzzy inference system or adaptive-network-based fuzzy inference system (ANFIS) method, and we consider the superficial gas velocity as one of the input variables, while for the first time, one of the computational fluid dynamics (CFD) outputs named gas velocity is used as an output of the artificial algorithm. The results show that the increasing number of input variables improves the intelligence of the ANFIS method up to , and the number of rules during learning process has a significant effect on the accuracy of this type of modeling. The results also show that propper selection of model parameters results in more accuracy in prediction of the flow characteristics in the column structure.

A ducted photovoltaic façade (DPV) unit Studied using experimental Prototype and simulated in a full scale Computational Fluid Dynamics CFD Model. The Study comes in two parts; This is Part I with the title detailed above and Part II titled ‘A Ducted Photovoltaic Façade Unit with Buoyancy Cooling: Part II CFD Simulation’.. The process adopted in the experimental study is replicated in the simulation Part. The aim was to optimize the duct width behind the solar cells to allow for maximum buoyancy-driven cooling for the cells during operation. Duct widths from 5 to 50 cm were tested in a Proto-type. A duct width of 45 cm had the maximum calculated heat removed from the duct; however, the lowest cell-operating temperature was reported for duct width of 50 cm. It was found that the DT between ducts' inlets and outlets range from 5.47 °C to 12.32 °C for duct widths of 5–50 cm, respectively. The ducted system enhanced module efficiency by 12.69% by reducing PV temperature by 27 °C from 100°C to 73 °C. The maximum calculated heat recovered from the ducted PV system is 422 W. This is 47.98% from the incident radiation in the test. Total summation of heat recovered and power enhanced by the ducted system is 61.67%.

In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation field was applied for the solution of the problems of local strong wind areas in railway fields, and the mechanism of wind generation was discussed. At the same time, the affectivity of the application of computational fluid dynamic approach to railway field was discussed. The problem of local wind that occurs on the railway line in winter was taken up in this research. A computational simulation for the prediction of wind conditions by LES was implemented and it was clarified that the local strong wind area is mainly caused by separated flows originating from the small‐scale terrain positioned at its upstream (at approximately 180.0 m above sea level). Meanwhile, the effects of the size of calculation area and spatial grid resolution on the result of calculation and the effect of atmospheric stability were also discussed. It was clarified that when the air flow characteristic of the separated flow originating from the small‐scale terrain (at altitude of approximately 180.0 m) targeted in this research is reproduced at high accuracy by computational simulation of wind conditions, approximately 10.0 m of spatial resolution of computational grid in horizontal direction is required. As a result of the computational simulation of wind conditions of stably stratified flow (Fr = 1.0), lee waves were excited at the downstream of the terrain over time. As a result, the reverse‐flow region lying behind the terrain that had been observed at a neutral time was inhibited. Consequently, local strong wind area was generated at the downstream of the terrain and the strong wind area passing through the observation mast was observed. By investigating the speed increasing rate of local strong wind area induced at the time of stable stratification, it was found that the wind was approximately 1.2 times stronger than what was generated at a neutral time.

Computational Fluid Dynamics (CFD) simulations, based on Reynolds Averaged Navier Stokes (RANS) models, are a useful tool for a wide range of coastal and offshore applications, providing a high fidelity representation of the underlying hydrodynamic processes. Generating input waves in the CFD simulation is performed by a numerical wavemaker (NWM), with a variety of different NWM methods existing for this task. While NWMs, based on impulse source methods, have been widely applied for wave generation in depth averaged, shallow water models, they have not seen the same level of adoption in the more general RANS based CFD simulations, due to difficulties in relating the required impulse source function to the resulting free surface elevation for non-shallow water cases. This paper presents an implementation of an impulse source wavemaker, which is able to self-calibrate the impulse source function to produce a desired wave series in deep or shallow water at a specific point in time and space. Example applications are presented, for a numerical wave tank (NWT), based on the opensource CFD software OpenFOAM, for wave packets in deep and shallow water, highlighting the correct calibration of phase and amplitude. Also, the suitability for cases requiring very low reflection from NWT boundaries is demonstrated. Possible issues in the use of the method are discussed and guidance for good application is given.

This manuscript discusses a novel method to map pressure results from one 3D surface shell mesh onto another. This technique is especially important when transferring results from one numerical analysis to another. This method works independent of the actual pressures, and only focuses on ensuring the surface areas consistently match. By utilizing this approach, the cumulative forces consistently match for all input pressures. This method is demonstrated to work for pressure profiles with precipitous changes in pressures, and with small quadrangular source elements being applied to a mix of large quadrangular and triangular target elements, and the forces at all pressure profiles match remarkably.

Trailing oil is the tail section of contamination. There are two main reasons for the formation of trailing oil, one is the effect of laminar flow boundary layer, the other is the outflow of the preceding batch remained in the dead-legs. In the batch transportation of refined oil, under the action of viscous force, the preceding batch forms laminar boundary layer near the pipe wall and stays on the pipe wall, resulting in the phenomenon of contamination trailing and formation of trailing oil. When oil passes through the valve chamber of the oil transportation station, dead-leg will be formed. Due to gravity and convection diffusion, preceding batch flowing from dead-legs will form trailing oil in the pipeline. The phenomenon of trailing oil exists in the process of batch transportation, which will have an effect on the quality of oil. In this paper, Reynolds time-averaged method is used to simulate turbulence.Computational Fluid Dynamics(CFD) software is used to simulate different flow rates and bypass lengths to obtain contamination-related experimental data.Matlab software is used to perform multi-nonlinear regression for the oil substitution time, the length of the bypass and the flow rate. The formula for calculating the length of the trailing oil produced by the dead-leg is obtained. The modified formula for calculating the length of the contamination is obtained by combining the existing formula for calculating the length of the contamination.

The Abdominal Aortic Aneurysm (AAA) is a local dilation of the abdominal aorta and it is a cause for serious concern because of the high mortality associated with its rupture. Consequently, the understanding of the phenomena related to the creation and the progression of an AAA is of crucial importance. In this work the complicated interaction between the blood flow and the AAA wall is numerically examined using a fully coupled Fluid-Structure Interaction (FSI) method. The study investigates the possible link between the dynamic behaviour of an AAA and the blood viscosity variations attributed to the haematocrit value, while it also incorporates the pulsatile blood flow, the non-Newtonian behaviour of blood and the hyperelasticity of the arterial wall. It was found that blood viscosity has no significant effect on von Mises stress magnitude and distribution, whereas there is a close relation between the haematocrit value and the Wall Shear Stress (WSS) magnitude in AAAs. This WSS variation can possibly alter the mechanical properties of the arterial wall and increase its growth rate or even its rupture possibility. The relationship between haematocrit and dynamic behaviour of an AAA can be helpful in designing a patient specific treatment.

This work studied the corrosion of welded pipes and how welding destroyed surface film of pipes. Surface reaction of a welded pipe is key to understanding phenomena and important factors during the corrosion. This paper presents experiment and CFD modeling approaches of a welded pipe corrosion under salt vapor condition. The pipes were welded at currents of 60 A,70 A and 80 A to observe the effect of welding current on corrosion. A welded pipe is a stainless-steel ASTM A312 grade 304L and period of experiment about 0-600 hours that they are tested in vertical and horizontal alignments. In CFD software, there is not direct model of corrosion but it can use surface reaction and create add-on species and chemical reaction technique for imitating the corrosion mechanism. The modeling approaches of corrosion have presented in 3-dimensional transient times in CFD simulation. Surface reactions were performed by Species Model which involve site species. Site species in Species Model took place at gas-solid interfaces and in this case are salt vapor and surface pipe. Chemical reaction rate on the surface controls lost weight of a welded pipe and the model can be validated with experiment. In conclusion, in period 0-600 hours error between CFD modeling and experiment have error trend decreased. The error at 600 hours is 6% both of vertical pipe and horizontal pipe test. The modeling approaches closely with the performed experiment and can be accepted. Moreover, the model is able to predict corrosion of a welded pipe of different sizes and their lost weight after 600 hours without experiment. Also the model can predict lifetime of pipe.

Into the frame of the French TANDEM project (Tsunamis in the Atlantic and the English ChaNnel: Definition of the Effects through numerical Modelling) Principia has been working on the development and qualification of two in-house CFD software’s: the 2D EOLE-SV (Saint Venant) model for simulation of large scale tsunami propagation from the source up to coastal scale and the 3D EOLE-NS (Navier-Stokes) model dedicated to tsunami coastal impact modelling. This paper presents a large range of test cases carried out into the frame of the project and dedicated to the validation of numerical codes in various tsunami wave conditions. The main aspects of phenomena such as wave generation, propagation and coastal impact are investigated on academic situations. A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.

Mixing is considered as a critical process parameter (CPP) during process development due to its significant influence on reaction selectivity and process safety. Nevertheless, mixing issues are difficult to identify and solve owing to their complexity and dependence on knowledge of kinetics and hydrodynamics. In this paper, we proposed an optimization methodology using Computational Fluid Dynamics (CFD) based compartmental modelling to improve mixing and reaction selectivity. More importantly, we have demonstrated that through the implementation of surrogate-based optimization, the proposed methodology can be used as a computationally non-intensive way for rapid process development of reaction unit operations. For illustration purpose, reaction selectivity of a process with Bourne competitive reaction network is discussed. Results demonstrate that we can improve reaction selectivity by dynamically controlling rates and locations of feeding in the reactor. The proposed methodology incorporates mechanistic understanding of the reaction kinetics together with an efficient optimization algorithm to determine the optimal process operation and thus can serve as a tool for quality-by-design (QbD) during product development stage.

The present studies deliver the computational investigations of a 10 MW turbine with a diameter of 205.8 m developed within the framework of the AVATAR (Advanced Aerodynamic Tools for Large Rotors) project. The simulations were carried out using two methods with different fidelity levels, namely the computational fluid dynamics (CFD) and blade element and momentum (BEM) approaches. For this purpose, a new BEM code namely B-GO was developed employing several correction terms and three different polar and spatial interpolation options. Several flow conditions were considered in the simulations, ranging from the design condition to the off-design condition where massive flow separation takes place, challenging the validity of the BEM approach. An excellent agreement is obtained between the BEM computations and the 3D CFD results for all blade regions, even when massive flow separation occurs on the blade inboard area. The results demonstrate that the selection of the polar data can influence the accuracy of the BEM results significantly, where the 3D polar datasets extracted from the CFD simulations are considered the best. The BEM prediction depends on the interpolation order and the blade segment discretization.

Simulation programs contain Computational Fluid Dynamics - CFD codes and are a useful tool used for gas flow through porous materials. Conducting numerical simulations allows for detailed analysis of hydrodynamic phenomena. The results of numerical modeling should always be verifiable based on experimental data. Only their compliance with the results of experimental tests is a determinant of the correctness of the applied method. As part of the work, experimental studies of hydrodynamics of gas flow through an isotropic porous material were carried out and numerical simulation for material of the same shape was used. In the CFD modeling Kolmogorov's hypothesis for the transport of kinetic energy of turbulence k and transport of dissipation rate of kinetic energy of turbulence ε was used.

Unmanned Aerial Vehicles (UAVs) have gained significant attention in recent times due to their suitability to a wide variety of civil, military and societal missions. Development of an unmanned amphibious vehicle integrating the features of a multi-rotor UAV and a hovercraft is focus of the present study. Components and subsystems of the amphibious vehicle are developed with due consideration on aerodynamic, structural and environmental aspects. Finite element analysis (FEA) on static thrust conditions and skirt pressure are performed to evaluate the strength of structure. For diverse wind conditions and angles of attack (AOA), computational fluid dynamic (CFD) analysis is carried out to assess the effect of drag and suitable design modification is suggested. A prototype is built with a 7 kg payload capacity and successfully tested for stable operations in flight and water-borne modes. Internet of Things (IoT) based water quality measurement is performed in a typical lake and water quality is measured using pH, dissolved oxygen (DO), turbidity and electrical conductivity (EC) sensors. The developed vehicle is expected to meet functional requirements of disaster missions catering to the water quality monitoring of large water bodies.

In this paper, an attempt has been made to analyze the effect of spool port/ groove geometry on the pressure drop and chamber pressures which effect the performance parameters of the flow distributor valve. The work mainly involves formulation of detailed mathematical model of the valve and compare them on the same platform. For mathematical modelling, Matlab has been used. The size of the orifices is considered same throughout the model for better comparison. Initially the construction and functioning of flow distributor valve along with working principles of hydrostatic motor (Rotary Piston) is shown. Next shown the analytical analysis of area change and pressure drops due to different geometry of the spool valve ports. After that the computational fluid dynamics (CFD) analysis has been shown. A complete mathematical model to describe such flow distributor valve is developed after having a comprehensive knowledge of orifice characteristics, flow interactions based on valve geometry. Equations of flow through different orifices (fixed and variable area) of the valve have been developed based on the relationships obtained earlier.

The aim of this study is to provide the scientists with a straightforward correlation that can be applied for predicting the Fanning friction factor and consequently the pressure drop during blood flow in small caliber vessels. Due to the small diameter of the conduit, the Reynolds numbers are low and thus the flow is laminar. The study has been conducted using CFD simulations validated with relevant experimental data acquired using an appropriate experimental set-up. The experiments concern pressure drop measurement during the flow of a blood analogue that follows the Casson model, i.e. an aqueous glycerol solution that contains a small amount of xanthan gum and exhibits similar behavior to blood, in a smooth, stainless steel microtube (L=5.6cm and D=400 μm). The interpretation of the resulting numerical data led to the proposal of a simplified model that incorporates the effect of the flow rate, the hematocrit value (35-55%) and the vessel diameter (300-1800 μm) and predicts with better than ±10% the Fanning friction factor and consequently the pressure drop during laminar blood flow in small caliber vessels.

Aerodynamic performance of a floating offshore wind turbine (FOWT) is significantly influenced by platform surging motions. Accurate prediction of the unsteady aerodynamic loads is imperative for determining the fatigue life, ultimate loads on key components such as FOWT rotor blades, gearbox and power converter. The current study examines the predictions of numerical codes by comparing with unsteady experimental results of a scaled floating wind turbine rotor. The influence of platform surge amplitude together with the tip speed ratio on the unsteady aerodynamic loading has been simulated through unsteady CFD. It is shown that the unsteady aerodynamic loads of FOWT are highly sensitive to the changes in frequency and amplitude of the platform motion. Also, the surging motion significantly influences the windmill operating state due to strong flow interaction between the rotating blades and generated blade-tip vortices. Almost in all frequencies and amplitudes, CFD, LR-BEM and LR-uBEM predictions of mean thrust shows a good correlation with experimental results.

Bioprocesses for the production of renewable energies and materials lack efficient separation processes for the utilized microorganisms such as algae and yeasts. Dissolved air flotation (DAF) and microflotation are promising approaches to overcome this problem. The efficiency of these processes depends on the ability of microorganisms to aggregate with microbubbles in the flotation tank. In this study, different new or adapted aggregation models for microbubbles and microorganisms are compared and investigated for their range of suitability to predict the separation efficiency of microorganisms from fermentation broths. The complexity of the heteroaggregation models range from an algebraic model to a 2D population balance model (PBM) including the formation of clusters containing several bubbles and microorganisms. The effect of bubble and cell size distributions on the flotation efficiency is considered by applying PBMs, as well. To determine the impact of the model assumptions, the modeling approaches are compared and classified for their range of applicability. Evaluating computational fluid dynamics (CFD) of a DAF system shows the heterogeneity of the fluid dynamics in the flotation tank. Since analysis of the streamlines of the tank show negligible backmixing, the proposed aggregation models are coupled to the CFD data by applying a Lagrangian approach.

Over the last few decades, ocean research and exploration have made underwater mechanical systems a necessity. Underwater vehicles provide a new kind of marine platforms that could represent a great necessity in many areas of oceanographic research. Until now, the underwater vehicles come in a verity of shapes, sizes and means of propulsion. Depending on these characteristics, the type and mission of the vehicle are also determined. The underwater robots are used for different inspection and intervention missions in e.g. the oil and gas industry, ocean science research. Due to multiple applications to which the vehicle can participate, it can be successfully used and to determine methods of re-use of marine energy. Environmental mapping provides accurate information about the main areas of interest of the energy, as well as the exploitation possibilities of that. Most of the time, biomimetic robots were inspired their senso structure, from different kind of animals, such as insects, fish and birds. Nowadays, the concept of a underwater robotic vehicles capable to move independently, autonomously or remotely, has a great potential and a large application. This is the reason that the last studies have been directed on biomimetic robots. The fish and other underwater animals have evolved superior swimming capabilities in many ways and represent a starting point to explain the fluid-mechanical principles. Furthermore, the underwater animals develop and achieve extraordinary propulsion efficiencies, acceleration and maneuverability. They can also achieve high speed under water. Implanting and creating a vehicle through a biomimetic approach reduces the energy used to maneuver the vehicle as it can automatically correct its position and displacement. The paper presents an examination of the state of biomimetic robotic fishes, underlining the reason why bio-inspiration can help us in the underwater locomotion technology.

Computational fluid dynamics (CFD) has, in the last decade, being an essential problem solving tool in industries such as pharmaceutical, pulp, petrochemical as well as Oil and Gas processing. The use of CFD for mixer design is unpopular in many countries in Africa. Therefore, this study investigates the characteristics of Brine-Surfactant mixing in a horizontal pipeline using CFD. The CFD is conducted by AnsysFluent software (licensed). A T-junction pipe is created and meshed with unstructured tetrahedral elements using design modeler. Discretization is done by Finite Volume Method (FVM), cell-centered scheme with Second-Order Upwind Scheme. The pressure term is introduced into the continuity equation by SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm. The kinetic epsilon model, also known as k-e model is administered to define the properties of the fluid and geometry, such as velocity, species mole fraction, pipe diameter. The boundary conditions is selected based on filed data. The rate of fluid flow in the primary region is 650 bbl/day at 400psi in 4in diameter pipeline, which is 100m long. The numerical simulation was based upon the governing equations such s continuity, Navier’s stoke, energy as well as specie transport equations. The findings shows that higher concentration results in increased mixing time, while 2% conc. of surfactant reaches homogeneity in 20 minutes at 72 meters of the pipe length. The result validated with field detail and the empirical result from literature, and is consistent. This study provide insight on industrial mixer design, chemical injection system, as well as gas pipeline design and optimization, especially in multiphase scale transport.

At the Atsumi Wind Farm in Aichi Prefecture, Japan, damage to wind turbines occurred frequently due to terrain-induced turbulence. In the present study, numerical analyses of terrain-induced turbulence were conducted by reproducing the topography in the vicinity of the wind turbine sites in high resolution and using RIAM-COMPACT natural terrain version, which is based on large eddy simulation (LES). The results of the diagnoses indicated that, in the case of south-easterly wind, terrain-induced turbulence is generated at a small terrain feature located upstream of Wind Turbine (WT) #2, which serves as the origin of the turbulence. At the Atsumi Wind Farm, a combination of the series of wind diagnoses and on-site operation experience led to a decision to adopt an "automatic shutdown program" for WTs #1 and #2. Here, "automatic shutdown program" refers to the automatic suspension of wind turbine operation upon the wind speed and direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. The adoption of the "automatic shutdown program" has successfully resulted in a large reduction in the number of occurrences of wind turbine damage, thus, creating major positive economic effects. 1) a reduction in the repair costs by 9.322 million yen per year per wind turbine, 2) an increase in the availability factor by 8.05%, and 3) an increase in the capacity factor by 1.7%.

This paper proposes a procedure for predicting the actual wind speed for flow over complex terrain with CFD. It converts a time-series of wind speed data acquired from field observations into a time-series of actual scalar wind speed by using non-dimensional wind speed parameters which are determined beforehand with the use of CFD output. The accuracy and reproducibility of the prediction procedure were examined by simulating the flow with CFD with the use of high resolution (5 m) surface elevation data for the Noma Wind Park in Kagoshima Prefecture, Japan. The errors of the predicted average monthly wind speeds relative to the observed values were less than approximately 20%.

We have developed an unsteady and non-linear wind synopsis simulator called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain) in order to simulate the airflow on a microscale, i.e., a few tens of km or less. In RIAM-COMPACT, the large-eddy simulation (LES) has been adopted for turbulence modeling. LES is a technique in which the structures of relatively large eddies are directly simulated and smaller eddies are modeled using a sub-grid scale model. In the present study, we have conducted the numerical wind diagnoses for Taikoyama Wind Farm nacelle separation accident in Japan. The simulation results suggest that all six wind turbines at the Taikoyama Wind Farm are subject to significant influence from separated flow (terrain-induced turbulence) which is generated due to the topographic irregularities in the vicinity of the wind turbines. A proposal has been also made on reconstruction of the wind farm.

In this paper, we studied and analyzed the behavior of velocity profiles and the distribution of pressures in a pipe connected to an accessory. In this case, it was a 90º elbow fitting in which water flows as an incompressible flow. This study was developed numerically through the fluent program of ANSYS. It is of practical interest in engineering because it is based on results regarding the internal structure of a model of incompressible flows. In the same way, it provides a basis for the proper design and implementation of velocity profiles within an elbow.In the agricultural industry, suitable systems can be implemented for crop irrigation in this way with the results obtained to improve the design in the fluid distribution system.We used a feasible κ-ε turbulence model because this model causes less turbulence than the standard κ-ε model. A mesh independence study was developed, in which it was observed that combining different densities of meshes for the same geometry causes some variables, so a fine meshing was used to carry out the study so that the results are as close to reality as possible.

A numerical study of the flow over a NACA aerofoil is presented in this paper. The numerical simulations are achieved with the computer code CFX and the computational domain is created by the computer tool ANSYS ICEM CFD. The CFX code is based on the finite volume method to solve the equations of mass, momentum and energy. The purpose of this paper is to determine the pressure distribution, flow patterns and the forces acting on the airfoil. Effects of the attack angle and Reynolds number on the velocity and pressure distribution, on the lift and drag coefficients are also explored.

In our previous works we have proposed design equations that can predict with reasonable accuracy the transition point from homogeneous to heterogeneous regime as well as the gas holdup and the mean Sauter diameter at the homogeneous regime. The validity of the proposed correlations was checked with data obtained using different geometrical configurations and several Newtonian and non-Newtonian liquids as well as the addition of surfactants. However, in all the experiments the gas phase was atmospheric air. This work investigates the effect of gas phase properties by conducting experiments employing various gases (i.e., air, CO₂, He) that cover a wide range of physical property values. Experiments revealed that only the use of low-density gas (He) has a measurable effect on bubble column performance. More precisely, when the low-density gas (He) is employed, the transition point shifts to higher gas flow rates and the gas holdup decreases, a fact attributed to the lower momentum force exerted by the gas. In view of the new data, the proposed correlations have been slightly modified to include the effect of gas phase properties and it is found that they can predict the aforementioned quantities with an accuracy of ±15%. It has been also proved that CFD simulations are an accurate means for assessing the flow characteristics inside a bubble column.

The optimization methods are increasingly used to solve challenging problems of aeronautical engineering. Typically, the optimization methods are utilized in design of aircraft airframe or its structure. The presented study is focused on an improvement of aircraft-flight-control procedures through the numerical optimization approach. The optimization problems concern selected phases of flight of light gyroplane - a rotorcraft using an unpowered rotor in autorotation to develop lift and an engine-powered propeller to provide thrust. An original methodology of computational simulation of rotorcraft flight was developed and implemented. In this approach the aircraft-motion equations are solved step-by-step, simultaneously with the solution of the Unsteady Reynolds-Averaged Navier-Stokes equations, which is conducted to assess aerodynamic forces acting on the aircraft. As a numerical optimization method, the BFGS algorithm was adapted. The developed methodology was applied to optimize the flight-control procedures in selected stages of gyroplane flight in direct proximity of the ground, where properly conducted control of the aircraft is critical to ensure flight safety and performance. The results of conducted computational optimizations proved qualitative correctness of the developed methodology. The research results can be helpful in design of easy-to-control gyroplanes and also in the training of pilots of this type of rotorcraft.

In this work the efficiency of a new μ-mixer design is investigated. As in this type of devices the Reynolds number is low, mixing is diffusion dominated and it can be enhanced by creating secondary flows. In this study we propose the introduction of helical inserts into a straight tube to create swirling flow. The influence of the insert’s geometrical parameters (pitch and length of the propeller blades) and of the Reynolds number on the mixing efficiency and on the pressure drop are numerically investigated. The mixing efficiency of the device is assessed by calculating a number, i.e. the Index of Mixing Efficiency that quantifies the uniformity of concentration at the outlet of the device. The influence of the design parameters on the mixing efficiency is assessed by performing a series of “computational” experiments, in which the values of the parameter are selected using DOE methodology. Finally using the numerical data, appropriate design equations are formulated, which, for given values ​​of the design parameters, can estimate with reasonable accuracy both the mixing efficiency and the pressure drop of the proposed mixing device.

The Gurney flap (GF) is a miniature lift-enhancement device and is usually mounted at the trailing edge of an airfoil. The GF and has been successfully applied to isolated airfoils, multi-element airfoils, and aircraft wings as well as helicopter rotor due to its attractive features of simplicity, cost-effective and separation control. The GF also has aroused the attention of researchers in the turbomachinery industry. However, limited studies are currently available on the application of a GF to an axial fan.Hence, in this paper, we conduct a wind tunnel and computational fluid dynamic (CFD) investigation on an axial fan profiled with a NACA 65-(12)10 airfoil to evaluate the effect of the GF on the performance of the fan. We also present the detailed flow features of the fan with and without the GF after validating the simulation results with the experimental results. The experimental results show that as the GF is installed higher on the fan blade, it can produce a higher total pressure rise accompanied with a greater loss of efficiency. The installation of the GF also enlarges the work capacity of the fan. Detailed flow field analysis, including the surface pressure distribution, vorticity distribution at the trailing edge and streamline distribution of the fan, is carried out to understand the mechanisms of the effect of the GF on the performance of the fan.

This work deals with the numerical investigation of the delivery of potential therapeutic agents through dentinal discs (i.e. a cylindrical segment of the dentinal tissue) towards the dentin-pulp junction. The aim is to assess the main key features (i.e. molecular size, initial concentration, consumption rate, disc porosity and thickness) that affect the delivery of therapeutic substances to the dental pulp and consequently to define the necessary quantitative and qualitative issues related to a specific agent before its potential application in clinical practice. The CFD code used for the computational study is validated with relevant experimental data obtained using micro Laser Induced Fluorescence (μ-LIF) a non-intrusive optical measuring technique. As the phenomenon is diffusion dominated and strongly dependent on the molecular size, the time needed for the concentration of released molecules to attain a required value can be controlled by their initial concentration. Finally, a model is proposed which, given the maximum acceptable time for the drug concentration to attain a required value at the pulpal side of the tissue along with the aforementioned key design parameters, is able to estimate the initial concentration to be imposed and vice versa.

In this paper, a two-phase model of air shock wave induced by rock-fall was described. The model was made up of the uniform motion phase (velocity was close to 0 m·s^-1) and the acceleration movement phase. The uniform motion phase was determined by experience, meanwhile the acceleration movement phase was derived by the theoretical analysis. A series of experiments were performed to verify the two-phase model and obtained the law of the uniform motion phase. The acceleration movement phase was taking a larger portion when height of rock-fall was higher with the observations. Experimental results of different falling heights showed good agreements with theoretical analysis values. Computational fluid dynamics (CFD) numerical simulation had been carried out to study the variation velocity with different falling height. As a result of this, the two-phase model could accurately and convenient estimating the velocity of air shock wave induced by rock-fall. The two-phase model could provide a reference and basis for estimating the air shock waves' velocity and designing the protective measures.

Performing optimisation and scale-up studies of crystallisation systems requires accurate and computationally efficient mathematical models. The assumption of the ideal mixing conditions in batch reactors typically produce inaccurate results while the computational expense of CFD models is still prohibitively high. Therefore, in this work, a new intermediary approach is proposed that takes into account the non-ideal mixing conditions in the reactor and requires less computational resources than full CFD simulations. Starting with the Danckwerts concept of the intensity of segregation, an analogy between its application to chemical reactions and the kinetics of the crystallisation phenomena (such as nucleation and growth) has been made. As a result, the modified kinetics expressions have been derived which incorporate the effect of non-idealities present in stirred reactors. This way, based on the experimental measurements of the mixing time using the Laser Induced Fluorescence (LIF) technique, computationally more efficient mathematical models can be developed in two ways: (1) the accurate semi-empirical correlations are available for standard mixing configurations with the most often used types of impellers, (2) CFD simulations can be utilised for estimation of the mixing time; in this case it is necessary to simulate only the mixing process. The benefits offered by the LIF experimental technique have been demonstrated and some frequent problems in its application analysed. The mixing time results for configurations with and without baffles for three types of impellers and four different rotational speeds have been presented. The false shorter mixing times in the non-baffled configurations have been observed and this phenomena explained by the existence of two segregated zones in the reactor and confirmed by additional experiments. The precise measurements in these cases have been shown as difficult using the LIF technique, particularly for higher rpms. The experimental data has been compared to the preliminary simulation results obtained from the Smoothed Particle Hydrodynamics method and the standard k-ε turbulence model with the modest success. The shortcomings of the SPH model have been recognized and the directions for the future work discussed.

The second article in the series presents the application of the Smoothed Particle Hydrodynamics (SPH) method to modelling of batch crystallisation in stirred tanks. A methodology to integrate the population balance equations (PBE) in parallel and independently from the Navier-Stokes equations is demonstrated. The benefits of the proposed methodology in terms of computational requirements, accuracy and availability of the crystal size distribution are discussed. The specific formulation of the SPH equations where the resulting system of ordinary differential equations is solved using the weighted contributions rather than numerically by solving a linear system of equations allows for massive parallelisation and a very loose coupling of the population balance and the fluid dynamics. It has been demonstrated, that the population balance equations can be solved on a Shared Memory Architecture (SMA) system using the OpenMP interface while the fluid dynamics equations being computed independently on a General Purpose Graphics Processing Unit (GPGPU) using the NVidia CUDA technology. This way, a significant portion of the computational overhead due to the large number of additional transport equations resulting from the discretisation of the population balance was removed: the SPH simulation coupled with 200 population balance equations was only 40% slower compared to SPH-only simulation. Two methods for the solution of population balance equations that preserve full crystal size distribution were implemented: discretised population balance (DPB) and method of characteristics (MOCH). The DPB equations are solved using the high-resolution finite-volume method with flux limiter and the effect of a large number of different flux limiters have been investigated. Both methods were validated using the case studies from the literature where an analytical solution can be derived. The developed models were applied to a numerical solution of coupled computational fluid dynamics and population balance equations to model a batch crystallization process. The effect of the hydrodynamics on the local temperature/supersaturation and the resulting crystal size distribution was captured and compared to the ideal mixing case. The simulation results from the DPB and MOCH methods were compared in terms of computational requirements and accuracy and MOCH selected as computationally more efficient and accurate.

NACA 0015 and NACA 4415 aerofoil are most common four digits and broadly used aerodynamic shape. Both of the shapes are extensively used for various kind of applications including turbine blade, aircraft wing and so on. NACA 0015 is symmetrical and NACA 4415 is unsymmetrical in shape. Consequently, they have big one-of-a-kind in aerodynamic traits at the side of widespread differences of their utility and performance. Both of them undergo the same fluid principle while applied in any fluid medium giving dissimilar outcomes in aerodynamics behavior. On this work, experimental and numerical investigation of each NACA 0015 and NACA 4415 is done to decide their performance. For this purpose, aerofoil section is tested for a prevalence range attack of angle (AOA). The study addresses the performance of NACA 0015 and NACA 4415 and evaluates the dynamics of flow separation, lift, drag, pressure and velocity contour and so on. This additionally enables to layout new optimistic aerofoil, which is critical to enhance the efficiency and performance of an aircraft in terms of lift enhancement and drag reduction.

In this study, a basic study was performed to analyze the seasonal temperature status of a research room in the Global Environment Research Building where ceiling-embedded indoor units are installed to study the room temperature status of the building as well as to improve its thermal environment. In addition, a direction for improvement of the indoor thermal environment in the winter was proposed through a CFD (computational fluid dynamics) simulation and was proven by an additional experiment. Through the results of this study, it appeared that if the ceiling-embedded indoor unit was installed in the small indoor space without considering the thermal vulnerability of its perimeter boundary, the air temperature of the upper part was greatly different from that of the bottom part in the winter. Hence, in this study, as a means to improve it, convectors were installed to minimize the effect of the external thermal environment and angle-controllable air flowing fans were installed to mitigate the stratification distribution. With such result, it was intended to present the essential data for improvement of the thermal environment as well as conservation of heating energy in the winter for buildings by reviewing the use of the ceiling-embedded indoor unit in the future.

Industries, which are mainly responsible for high energy consumptions, need to invest in research projects in order to develop new managing systems for rational energy use and to tackle the devastating effects of climate change caused by human behavior. The study reported in this paper concerns the forging industry, where the production processes generally start with the heating of the steel in furnaces and continue with other processes, such as heat treatments and mechanical machining. One of the most critical operations, in terms of energy loss, is the opening of the furnace doors for the insertion and extraction operations. During this time, the temperature of the furnaces decreases by hundreds of degrees in a few minutes. Because the dispersed heat needs to be supplied again through the combustion of fuel, increasing the consumption of energy and the pollutant emissions, the evaluation of the amount of the lost energy is crucial for the development of operating or mechanical systems able to contain this dispersion. To perform this study, CFD simulation software was used. Results show that at the door opening, because of temperature and pressure differences between the furnace and the ambient, turbulences are generated. Results also show that the amount of energy lost for an opening of 10 minutes for radiation, convection and conduction is equal to 5606 MJ where convection is the main contributor with 5020 MJ. The model created, after being validated, has been applied to perform other simulations in order to improve the energy performance of the furnace. Results show that a reduction of the opening time of the door allows energy savings and limits pollutant emissions.