Computational Fluid Dynamics Investigation Research Articles

Efficiency Improvement of Darrieus Wind Turbine Using Oscillating Gurney Flap

In this work, a new model of Darrieus wind turbines with an oscillating gurney flap (OGF) is proposed. A detailed 2D computational fluid dynamics (CFD) investigation is carried out using ANSYS-Fluent 22.0 to assess the turbine performance. The OGF can alter its position between the upper and lower blade surfaces during the turbine rotation. Equations related to the combined motion are implemented through a user-defined function (UDF). The proposed model is validated where a good coincidence is achieved. The overset dynamic mesh method is used. It was found that a judicious synchronization of OGF and turbine blades creates beneficial vortex interactions, which correct the pressure distribution and lead to an overall improvement in the lift force. The magnitude of the improvement is highly dependent on the OGF length and the phase motion φ. The average torque coefficient Cm for the controlled case increased by more than 19% in comparison with the nominal case.

Computational Fluid Dynamics Investigation of the Spacing of the Aerodynamic Characteristics for Multiple Wingsails on Ships

Wind energy, as an inexhaustible energy source, has become a focal point in the development of new energy for ships. Sail-assisted technology, which leverages wind power, has been successfully applied to ship propulsion. The propulsion performance of sail-assisted ships is affected by the interference characteristics among multiple wingsails. To investigate interference characteristics, an arrangement scheme involving two-element wingsails and considering the relative wind direction angle was established. To obtain the inter-stage interference characteristics of wingsails, the Reynolds average N-S equation was used in the numerical simulation conducted under steady operating conditions. The results indicate that, at the relative wind angles of 30°, 90°, and 120°, the minimum horizontal spacing in a single row arrangement scheme is 1.5c. However, at relative wind angles of 90° and 120°, inter-stage interference may induce stall conditions in the wingsails. In a double-row arrangement scheme, the wake of the upstream wingsail interferes with the flow of the downstream sail at relative wind angles of 90°. An optimal propulsion performance is achieved with a horizontal spacing of 4c and a longitudinal spacing of 10c. Moreover, the interference performance of the two-element wingsails can be enhanced through a horizontal offset arrangement. This study provides a reference for the arrangement of wingsails on ships.

CFD investigation on start-up transient of HPR1000 secondary side passive residual heat removal heat exchanger

In advanced pressurized water reactor HPR1000, the secondary side passive residual heat removal heat exchanger (SPRS HX) is a crucial component. It removes the decay heat using boiling and convection heat transfer in the event of accidents. The robust functioning of the SPRS HX significantly influences overall reactor safety. To investigate the thermal hydraulic attributes of SPRS HX, computational fluid dynamics (CFD) is used in this study. The tube region is simulated using the porous media method (PMM). To model the two-phase flow in the tank, the drift flux model (DFM) is employed. The assessment of heat transfer rate from the primary to the secondary sides is conducted by empirical correlations. The governing equations are addressed through the utilization of FLUENT. Specifically, the momentum and energy equations are augmented with additional source terms that account for flow resistance and heat transfer in the tube region. These modifications are enacted through User Defined Functions (UDF). Three-dimensional profiles of the void fraction, fluid temperature and velocity are acquired, allowing for a comprehensive analysis of the heat transfer attributes of tube region.

Enhancing CO2 capture efficiency: Computational fluid dynamics investigation of gas-liquid vortex reactor configurations for process intensification

This study employs non-thermal Computational Fluid Dynamics (CFD) simulations to explore the efficacy of a gas–liquid vortex reactor (GLVR) for intensifying CO2 capture. The investigation concentrates on the multiphase flow and mass transfer behavior in diverse GLVR experimental units. Employing a multiphase Euler-Euler CFD model integrated with a Population Balance Model (CFD-PBM), a mass transfer model, and reaction kinetics, allows us to accurately simulate the reactive absorption of CO2 into aqueous Monoethanolamine (MEA). Experiments are conducted using a 30 wt% MEA solution for CO2 absorption, serving the purpose of model validation. This comprehensive approach enables a simulation of complex dynamics within the GLVR, emphasizing bubble breakage, coalescence, and reactive mass transfer processes. Examining bubble size distribution, pressure drop, CO2 absorption efficiency, and energy input systematically across various reactor geometries and operational conditions, our findings demonstrate that an optimized GLVR configuration significantly enhances CO2 absorption compared to the original design. Furthermore, the optimized GLVR outperforms state-of-the-art process intensification equipment in terms of CO2 absorption rate per unit reactor volume and energy efficiency.

Computational Fluid Dynamics Investigation of Hydrodynamic Forces and Moments Acting on Stern Rudder Plane Configurations of a Submarine

This study presents the predicted hydrodynamic characteristics of different rudder plane configurations on the stern of a full-scale submarine in deep water, which are obtained using the Reynolds-Averaged Navier–Stokes method in Ansys Fluent Solver. First, the results obtained for the X-rudder plane configuration are verified according to previous numerical and experimental results in order to assess the accuracy of the simulation procedure. The X-rudder plane, Y-rudder plane, and Cross-rudder plane configurations in deep water with deflection angles ranging from −21 degrees to +21 degrees are then simulated. Next, the hydrodynamic forces and moments of the Cross-plane, X-plane, and Y-plane rudder configurations obtained through simulation are analyzed using Taylor’s expansion to estimate the hydrodynamic coefficients. The obtained results demonstrate that the X-force of the X-plane rudder configuration is larger than the corresponding forces acting on the Cross-plane rudder and Y-plane rudder configurations. Meanwhile, the Y-force and Z-force of the X-plane rudder configuration are significantly greater than the corresponding forces of the left configurations. The same tendency can be seen in the moment of the X-plane rudder about the y- and z-axes. However, the roll moment induced by the Y-plane and Cross-plane rudder configurations is significantly larger than that under the X-plane rudder configuration.

Influence of baffle height on the in-cylinder mixing and performance characteristics of a gasoline direct injection engine - a CFD investigation

In modern GDI engines, in-cylinder flows play a crucial role in the performance and emission characteristics. Adding baffles to the cylinder head is one way to enhance these flows, as baffle effectiveness is linked to its height. This study aims to investigate the influence of baffle height on the in-cylinder flows and its impact on the performance and emission characteristics of a four-stroke, four-valve, spray-guided GDI engine using computational fluid dynamics (CFD) analysis. The maximum baffle height is determined based on the geometry of the engine. Throughout the analysis, the engine operates at a constant speed of 1000 rev/min., under part-load conditions. From the results, the baffle heights of 2 mm and 3 mm significantly improve mixture stratification compared to the base engine. Furthermore, the indicated mean effective pressure (IMEP) for the engine with 2 mm and 3 mm baffle height improves by about 3% and 4%, respectively, compared to the base engine. In addition, with the increase in the baffle height, there is a significant reduction in the HC emissions with a marginal increase in NOX emissions.

3D CFD simulation of room air temperature and velocity with ceiling fan and exhaust fan and the resulting thermal comfort

The integration of energy-efficient solutions remains an essential hurdle in the pursuit of sustainable and comfortable indoor settings. As building energy consumption continues to climb, there is an urgent need to implement methods that reduce environmental effects, and at the same time, maintain adequate thermal comfort for occupant well-being and productivity. Ceiling fans, a standard fixture in many indoor settings, have gained popularity for their ability to provide better air movement and temperature adjustment while using substantially less energy than standard HVAC systems. This work presents a 3D Computational Fluid Dynamics (CFD) investigation of thermal comfort that can be achieved with a ceiling fan and an exhaust fan. A 3D fan zone model available in ANSYS Fluent was used to model the ceiling fan, considering the ceiling fan as a solid disc. Numerical methodology has been validated by verifying the transient simulation results against the experimental data available in the literature. The work analyzes thermal comfort parameters - temperature and velocity - in a room with a typical ceiling fan or an exhaust fan available in India. The room temperature and velocity distribution are presented for different times of a day in summer for a warm and humid climatic zone in India. The results show that the thermal comfort obtained with the exhaust fan is better than that with the ceiling fan for the evening, night, and early morning.

Thermal Hydraulic Analysis of Turbulent Flow of Superheated Steam Through Multiple Configurations of 3D Angular Piping Bend in Power Plant Installations

Abstract One of the most vulnerable parts in high-pressure and temperature pipeline networks in steam turbine power plant installations is pipe bends, due to intense fluctuation of velocity and pressure of the superheated steam flowing through the piping bend. On the other hand, because of unique chemical and physical properties, the transportation of superheated steam can adversely effect on the integrity of the pipe bends in the piping network. Therefore, the potential risk of pipeline failure becomes more probable at bends. Therefore, to ensure the survivability of the piping network, thorough investigations of the superheated steam flowing through pipe bends are of great significance in better understanding its flow behavior. However, to predict flow behavior at piping bend, intensive research by direct experiment at high temperature and pressure might not be possible to some extent. In that case, computer codes and models can help us analyze flow behavior of superheated steam at pipe bends. The current work is a computational fluid dynamics investigation, where numerical examination is carried out by commercial code ansys fluent for thermal hydraulic analysis of flow behavior of superheated steam through multiple configurations of 3D angular piping bend: 90 deg, 120 deg, and 135 deg. Validation of credibility of the computational approach is done against Sudo et al.'s experiment. With relatively good agreement between the experiment and numerical result using air flowing through 90 deg elbow, single-phase turbulent flow of superheated steam is examined by the application of realizable turbulence model (k–ɛ). The final assessment of the computed results has been done using qualitative and quantitative analyses. The study reveals that unlike the pressure profiles, the radial velocity profiles at the bend exit are noticeably different from those at the bend inlet, since they exhibit more rounded peaks and gentler velocity gradients. The study also shows that the radial pressure and velocity profiles do not display any sign of the existence of a pair of vortices.

Computational Fluid Dynamics (CFD) Validation and Investigation the Effect of Piston Bowl Geometries Performance on Port Fuel Injection-Homogeneous Charge Compression Ignition (PFI-HCCI) Engines

Computational Fluid Dynamics (CFD) Validation and Investigation the Effect of Piston Bowl Geometries Performance on Port Fuel Injection-Homogeneous Charge Compression Ignition (PFI-HCCI) Engines

Molding the acoustic cavity–analyzing the influence of toroidal vortex development on acoustic multi-bubble macrostructures under different ultrasonic horn tip diameters

The acoustic cavity structure typically experiences a sequence of transfigurations during its sinusoidal growth–collapse cycle. However, upon examining the cavity structure in aqueous bodies, it appears that the growth structure attained falls between two geometrical structures, namely, mushroom-like structure (MBS) and cone-like bubble structure (CBS), based on the actuated ultrasonic horn tip diameter. With the recurring observations of the emergence of proximal toroidal vortices, the present investigation conducts a numerical analysis exploring the vortex development under 3, 6, 13, 16, and 19 mm horn tips to establish a potential correlation between the vortex and the cavity structure. The study presents a computational fluid dynamic investigation to capture the nature of the vortex evolution, in terms of size and position, and its respective cavitation development. The first indicator of potential correlation was the equivalency of the vortex expansion–contraction frequency and the cavity's sub-harmonic frequency. It has been found that the cavity structure is molded into MBS by the presence of a symmetric locomotive vortex structure that extends up to 1.5 times the horn tip diameter. Meanwhile, CBS is observed to take shape in the presence of an eccentric locomotive vortex that attains a size within 0.2–0.6 times the horn tip diameter. The significance of the vortex size and position is also observed in the cavity's collapse, as the vortex appears to govern the ability of the cavity impinging jet to initialize the collapse phase.

CFD investigation on a novel pen partition-attached jet air supply for mitigating heat stress among lactating sows

To effectively mitigate heat stress among lactating sows, this study proposed a novel pen partition-attached jet air supply (PJAS) and compared its performance with commercial ceiling jet air supply (CJAS) by using computational fluid dynamics (CFD). Each PJAS unit, comprising a circular distribution duct with 5 slot orifices at each side, was installed at the pen partition between two adjacent pens to direct air towards lactating sows to enhance convective heat removal. The cooling performances of air supplier were primarily evaluated by assessing the convective heat flux from lying sows. The effect of jet angles induced by PJAS and CJAS on convective heat flux were investigated as well as the consequence of reducing the inlet dimension of CJAS. The results revealed that reducing the inlet dimension of CJAS doubled the speed of supplied air. However, it also limited the amount of air impinging on sow due to the narrowing of the inlet air jet. Consequently, the combined effect only marginally increased the convective heat flux when the inlet area of CJAS was reduced. By placing air supplier closer to lactating sows, the air jet induced by PJAS impinged on sows’ bodies directly and therefore significantly increased the convective heat removal. Compared to CJAS, PJAS required approximately 45% less ventilation rate to induce the same convective heat flux. When the ventilation rate was identical, PJAS could increase the convective heat flux by up to 27.6%. These results highlight PJAS as a promising energy-efficient localized cooling device for mitigating heat stress among lactating sows. To determine the optimal design and control strategy of PJAS based on industrial requirements, further research is still needed.

CFD investigation of chalcopyrite flotation prediction coupled with flotation kinetic model

CFD (computational fluid dynamics) simulation plays an important role in flotation prediction research and application. The flotation process of chalcopyrite in a three-dimensional mineralized tube was simulated coupled with flotation kinetic model. Regarding kinetic prediction model in the flotation process is explored at gas flow rate of 0.06 m3/h, slurry flow rate of 0.4 and 0.5 m3/h, including collision probability model, attachment probability model, detachment probability model and bubble slip velocity model. Flotation experiments of chalcopyrite is performed, and the simulation of flotation process is realized by adding UDF (user-defined function). Furthermore, the flotation rate constant is taken to verify the validity of the model. The prediction error of flotation rate constant is 19 % and 17 % between flotation experiments and simulations at slurry flow rate of 0.4 and 0.5 m3/h, respectively. According to the combinations of flotation kinetic model, it’s found that the prediction results of Yoon-Luttrell collision probability model, Yoon-Luttrell attachment probability model, Schulze detachment probability model and Kostoglou bubble slip velocity model are superior to other combinations. the prediction error of flotation rate constant is 8 % at the slurry flow rate of 0.4 m3/h.

Computational fluid dynamics investigation on the irrigation of a real root canal with a side-vented needle

BackgroundRoot canal therapy is one of the main treatments for root canal diseases, and effective irrigation is the key to successful treatment. Side-vented needle is one of the commonly used needle types in clinic. In the real root canal, due to the influence of the curvature of the root canal, the irrigation flow field in different needle directions shows obvious differences. At the same time, changes in root canal curvature and working depth will lead to changes in irrigation efficiency and the flow field. Both the mainstream of the irrigation flow and the shear stress near the wall changes significant. Consequently, either the replacement in the root canal or the removal efficiency of the smear layers is apparently modified.Materials and methodsIn this paper, the permanent root canal of the maxillary first molar prepared until 15/04 were scanned by micro-CT, and then imported into the software for 3D reconstruction. The key parameters of flushing efficiency of 30G side needle at different working depths of 4.75 mm, 5 mm, 5.25 mm and 5.5 mm were compared. Meanwhile, the simulated models with different curvatures of 0°, 5°, 10°, 20° and 30° based on the real root canal were reconstructed to investigate the curvature effect on the irrigation efficiency.ResultsThe results show that moderate working depth (such as 4.75 mm and 5.25 mm in present paper) helps to improve the replacement capacity of irrigation flow. At the same time, the apical pressure decreased as the working depth increased. The curvature of the root canal seriously affects the removal depth of the smear layers of the root canal. A root canal with a large curvature (especially 20° and 30°) can significantly improve the difficulty of irrigation.Conclusions(1) Moderate working depth helps to improve the displacement capacity, the ERD of the irrigation flow is generally improved at the working depths of 4.75 mm and 5.25 mm, and the apical pressure will decrease with the increase of working depth. (2) The large curvature of the root canal can significantly improve the difficulty of irrigation. The curvature of the root canal can severely influence the removal depth of the smear layer on the wall. It can be found both the span and the depth of the ESS for little curvatures (5° and 10°) root canals are higher than those for large curvatures (20° and 30°).

A computational fluid dynamics investigation of endothelial cell damage from glaucoma drainage devices

Glaucoma drainage devices (GDDs) are prosthetic-treatment devices for treating primary open-angle glaucoma. Despite their effectiveness in reducing intraocular pressures (IOP), endothelial cell damage (ECD) is a commonly known side-effect. There have been different hypotheses regarding the reasons for ECD with one being an induced increase in shear on the corneal wall. A computational fluid dynamics (CFD) model was used to investigate this hypothesis in silico. The Ahmed Glaucoma Valve (AGV) was selected as the subject of this study using an idealised 3D model of the anterior chamber with insertion angles and positions that are commonly used in clinical practice. It was found that a tube-cornea distance of 1.27 mm or greater does not result in a wall shear stress (WSS) above the limit where ECD could occur. Similarly, a tube-cornea angle of 45° or more was shown to be preferable. It was also found that the ECD region has an irregular shape, and the aqueous humour flow fluctuates at certain insertion angles and positions. This study shows that pathological amounts of WSS may occur as a result of certain GDD placements. Hence, it is imperative to consider the associated fluid force interactions when performing the GDD insertion procedure.

CFD Investigation of Methane Combustion with Excess Air in Can-Type Gas Turbine Combustor

This paper presents an investigation on the effect of excess air to combustion characteristics in a full-scale, gas turbine combustor, commonly used in power plant. The investigation was carried out using Computational Fluid Dynamics (CFD), and prior validation was made with the actual operation data of a power station, as well as the adiabatic flame temperature of methane. The mass fraction of CO<sub>2</sub>, O<sub>2</sub> and NO<sub>x</sub> emissions for blended methane with different percentages of excess air explosions was also investigated. The stoichiometric excess air varies from 0% to 30% with air-fuel mixture of 2.7 kg/s. The geometric model of the combustor is extracted from actual gas turbine combustor using 3D scanner and converted into CAD model for simulation. The Navier-Stokes equations were solved using commercial CFD code, ANSYS Fluent, with RNG k-ε chosen to close the turbulence. For reacting species, the species transport model is assumed. Results showed that the addition of excess air during combustor has little effect on the velocity and temperature distribution, both at the combustor interior, as well as at the exit. For the emission of CO<sub>2</sub> and O<sub>2</sub>, though there was no clear trend on the relations between the emission of these species and the excess air, the impact was quite significant. The production of NO<sub>x</sub> was also found to be independent on the excess air ratio, but instead, was a strong function of combustion and exhaust gas temperature.

Theoretical and CFD investigation on nonlinear dynamics of a cable under vortex-induced vibration with different aerodynamic shapes

This paper investigates the effect of the aerodynamic shape on nonlinear dynamics of a cable induced by vortex-induced vibration (VIV) with theoretical and computational fluid dynamics (CFD) method. Firstly, a two-dimensional (2-D) CFD model is established and calculated by using the shear stress transport (SST) k-ω model to supply the necessary drag and lift coefficients for perturbation analysis. Secondly, the planar transverse vibration equation of the cable is established and the wind load is described by the van der Pol wake oscillator. Then, the above two equations are discretized by the Galerkin method, and the modulation equations are obtained by the multiscale method. The results are verified by Runge-Kutta method. Thereafter, three different aerodynamic shapes of the cable are considered and five wind attack angles are defined to extensively investigate the influence on nonlinear dynamic behaviors. The results show that the variation of the cable cross-section significantly changes its aerodynamic performance, and the aerodynamic coefficient of the asymmetric structure is very sensitive to the wind attack angle. Arranging aerodynamic measures (such as soft tail) in the downwind of the cable can effectively reduce the aerodynamic coefficient and dynamic response of the structure.

Experimental and 3D computational fluid dynamics (CFD) investigation of multiple arrays of jet impingement on a flat surface

Experimental and 3D computational fluid dynamics (CFD) investigation of multiple arrays of jet impingement on a flat surface

Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms.

Computational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient-specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient-specific viscosity values derived from hematocrit measurements. The objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient-specific viscosity in future CFD investigations of cerebral aneurysms. CFD simulations of forty-one cerebral aneurysms were conducted using patient-specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient-specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t-tests for differences in the time-averaged maps of these hemodynamic parameters between standardized and patient-specific viscosity simulations were conducted for each aneurysm. Bland-Altman analysis was used to examine the cohort-wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t-test was used to compare the percent change in WSS, WSSG, and OSI between patient-specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t-tests, a Bonferroni-corrected significance level of 0.0167 was used. 63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient-specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study-wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after-the-fact correction for patient-specific viscosity is not feasible. Standardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient-specific viscosity. We recommend implementing hematocrit-based patient-specific blood viscosity values for all CFD simulations of cerebral aneurysms.

Computational Fluid Dynamics Investigation of the Impact of Wall Treatment on Fluidization in a Liquid-Solid Fluidized Bed

Computational Fluid Dynamics Investigation of the Impact of Wall Treatment on Fluidization in a Liquid-Solid Fluidized Bed

Experimental and Computational Fluid Dynamics Investigation on Tanning Process in a Rotating Drum

Mass transfer of chemicals greatly affects leather production efficiency and product quality. Leather shows different motions in a rotating drum during processing, which is strongly associated with chemicals’ mass transfer. However, how leather motions affect mass transfer remains unclear, which disfavors highly efficient leather manufacturing process. Here, different leather motion states were obtained by adjusting the drum rotation speed. Experimental results showed that the duration of leather rolling motion greatly increased by 41% when the rotation speed increased from 5 r/min to 20 r/min, and the uptake of the tanning agent was consequently improved, which indicated that the rolling motion is beneficial to mass transfer. Computational fluid dynamics simulation results showed that the mass transfer rate under rolling motion was higher than those under slipping, elevating and hanging motions, because the flow velocity and concentration gradient near the leather surface were higher under rolling motion. Accordingly, increasing the rolling motion enhanced the mass transfer in leather processing. This work identifies the leather motion beneficial for mass transfer and provides guidance on operating condition optimization and drum design for high-efficiency leather production.