Computational Fluid Dynamics Code Research Articles

Cavitating flow through a Venturi using CFD: effects of inlet and outlet pressures

Hydrodynamic cavitation is a complex physical phenomenon that appeared in hydraulic systems (pumps, turbines, valves, Venturi tubes …etc.) when the fluid pressure decreases below the saturated vapor pressure. The works carried out in this study aimed to get a better understanding of the cavitating flow phenomenon. For this, we have numerically studied a cavitating bubbly flow through a venturi nozzle. The cavitation model is selected and solved using commercial code computational fluid dynamics (CFD). The numerical results are compared to the previous experimental and numerical data.The obtained results showed that the effect of the inlet pressure (10, 7, 5 and 2 bars) and outlet pressures (6.5, 6, 5, 4, 3.5, 1.5, and 1 bars) of the Venturi on pressure, velocity of the fluid flow and the vapor fraction formation. We found that the inlet and outlet pressures of the venturi are strongly affected on the evolution of the pressure, velocity and vapor fraction formation in the cavitating flow. The results presented in this study provide insight and useful guidance to the designer of hydrodynamic cavitation devices, identifying key geometrical and physical parameters of cavitation onset and intensification that can be manipulated to achieve the desired level of cavitation flow.

Numerical investigation on the core thermal hydraulic behavior of pool-type sodium-cooled fast reactor (SFR)

A thorough understanding of the reactor core thermal hydraulic behavior is essential for the design and safety analysis of Sodium-cooled Fast Reactors (SFR). Due to the application of hexagonal subassembly, the core thermal hydraulic behavior is significantly affected by the flow field within the subassemblies, the inter-wrapper region and the hot pool. Analysis of the core thermal hydraulic behavior requires a model coupling the three regions mentioned above, which has been identified as one of the thermal hydraulic challenges in SFR. In the present study, a 3D model that covers the three regions was developed for the core of the China Experimental Fast Reactor (CEFR) with the Computational Fluid Dynamic (CFD) code, Fluent. The inter-wrapper region and the hot pool were modeled in detail, while the subassemblies were modeled with a special porous medium model. The core thermal hydraulics behavior under steady state was studied, more specifically, information for the flow field distribution at the core outlet, the inter-wrapper flow and the duct wall temperature distribution was obtained. Under steady state, liquid sodium in the inter-wrapper region is supplied by the inner region of the hot pool. And it enters the inter-wrapper region from the core outer region and returns back to the hot pool inner region from the core central region. The inter-wrapper flow is cooled by non-fuel subassemblies and heated up by fuel subassemblies. For non-fuel subassembly, the ratio of the total heat transfer rate between the inter-wrapper flow and the subassemblies to the heat generated within subassemblies could reaches 96%; for fuel subassemblies, the maximum ratio of the total heat transfer rate between the inter-wrapper flow and the subassemblies to the heat generated within subassemblies is 2.45%. Significant temperature gradients have been observed on the duct wall, with maximum values of 156.69 K/m in the vertical direction and 2,196.00 K/m in the circumferential direction. The largest temperature gradient appears on the duct of subassemblies adjacent to the transition region of fuel subassemblies and non-fuel subassemblies.

CFD Simulations of a Loss of Heat Sink Experiment in the McMaster Nuclear Reactor

The McMaster nuclear reactor (MNR) is an important research facility that not only provides researchers with neutrons for fundamental science, but also supplies the medical industry with isotopes used for cancer treatment. To ensure the safety and performance of the MNR, modeling of the thermal hydraulics during nominal and accidental conditions is required. For such a task, system codes are customarily used. While system codes can assess the safety aspects of complex thermal–hydraulic systems, the question arises whether such systems can be modeled appropriately in a three-dimensional manner, such as computation fluid dynamics (CFD), while capturing the thermal mixing of the coolant. Modeling the thermal hydraulics of nuclear research reactors using CFD is relatively new. Therefore, validating such methods against experimental data is of utmost importance. The validation of CFD is the main focus of this work. More specifically, the influence of the loss of heat sink on the pool temperature is assessed using the CFD code Ansys CFX. The validation basis was provided by an experiment performed at the MNR in March 2023. In this experiment, the secondary heat removal system was intentionally shut down, and the pool temperature was measured in a few positions. The results obtained by modeling the loss of heat sink in the MNR using Ansys CFX agree well with the experiment.

Coupled computational fluid dynamics and computational thermodynamics simulations for fission product retention and release: A molten salt fast reactor application

This study presents a computational capability for fission product retention and release in two-phase, multi-species systems representing Molten Salt Reactors (MSR) with coupled thermal-hydraulics and fuel coolant chemical behaviours. This is demonstrated through four simulated cases centred on the proposed Molten Salt Fast Reactor (MSFR). This is achieved by two-way coupling the Computational Fluid Dynamics (CFD) code OpenFOAM and the Computational Thermodynamics (CT) code Thermochimica, using the Joint Research Centre Molten Salt Database (JRCMSD). Local chemical equilibrium is assumed, implying that chemical kinetics are predominantly governed by mass transport. Four simulations address normal operating conditions, exploring: (i) dilution of fission products injected within the molten salt coolant, (ii) molten salt coolant evaporation rate, (iii) release of radioactive gaseous species, (iv) shifts in the UF4/UF3 ratio, and (v) comparison of vapour pressures of gaseous species. The influence of temperature-dependent viscosity on retaining fission products, compared to consistent values, is also discussed. The feasibility of integrating CFD with Thermochimica showed promising results, broadening insights into multiphysics systems and setting the stage for its application in more intricate scenarios.

The application of the state-of-the-art turbulence models in the OpenFOAM CFD codes and the validation for different flow control cases in the continuous casting process

Computational fluid dynamics (CFD) is an established approach to study flow behaviour. Despite advances in mathematical numerical modelling in recent decades, accurately predicting the turbulent flow behaviour of steel using CFD remains a challenge. In the continuous casting process, two main flow control systems are employed to transfer liquid steel from the tundish to the mould, namely stopper control and slide gate control. In additional to these two main flow control systems, the casting nozzle geometries, casting rate, and casting mould dimensions vary considerably, and this complexity makes it difficult to precisely simulate turbulent flow in the casting channel. In this paper, examples are presented showing how different turbulence models influence the simulation results, using the latest version of the open source CFD simulation software OpenFOAM. Prior to the CFD modelling, RHI Magnesita's water modelling facilities in Leoben (Austria) were used to validate the state-of-the-art turbulence models selected for a particular casting nozzle geometry and flow control system.

Improving efficiency of road transport through the use of hydrodiodes in hydraulic systems taking into account cavitation phenomena

Introduction. One of the promising ways to reduce dynamic loads during operation of hydraulic crane manipulator installations is the use of hydrodiodes. In hydraulic systems operating at high speeds of the working fluid, cavitation phenomena may occur in the hydrodiode, which are accompanied by increased noise and vibration, and may lead to the destruction of the hydrodiode, which is not acceptable. The paper compares the results of calculating the flow of liquid in the flow part of a vortex hydrodiode, taking into account cavitation phenomena and without cavitation phenomena, with the results of research tests. The analysis of the effect of cavitation on the working processes in a vortex hydrodiode for crane manipulator installations is presented.Materials and method. Computational fluid dynamics (CFD) models using the FLUENT CFD code to study the working processes occurring in the working chamber of a vortex hydrodiode were developed. The commercial CFD code ANSYS FLUENT to simulate the flow of liquid in the flow part of a vortex hydrodiode was used.Results. The paper verifies the results of a numerical experiment with the results of research tests. A quantitative and qualitative analysis of the effect of cavitation on the working processes of a vortex hydrodiode has been carried out.Discussion and conclusion. It has been found that the values of pressure and diode in calculations taking into account cavitation and without cavitation practically do not differ, thus, the effect of cavitation at Reynolds numbers Re<30000 does not significantly affect the quantitative values of the parameters of the vortex hydrodiode and the cavitation calculation module can not be used. However, at higher values of the Reynolds numbers, cavitation appears in the working cavity of the vortex hydrodiode and the calculated values of the parameters of the hydrodiode without cavitation and taking cavitation into account differ significantly. Therefore, when calculating high-speed flows, it is necessary to use the cavitation calculation module. The analysis of the effect of cavitation on the working processes of a vortex hydrodiode showed that in the forward direction of the flow, cavitation does not significantly affect the parameters of the hydrodiode and the place of its formation, the upper inlet region of th tangential chamber. In the opposite direction of flow, cavitation has a significant effect on the pattern of fluid flow in the vortex hydrodiode. Cavitation covers almost the entire inlet volume of the radial tube and partially captures the central part of the vortex chamber

Heat transfer augmentation in heat sink with fin-bars inserted in cooling channel

Adding obstacles and obstructing fluid flow in cooling channel to induce whirl velocity enhances thermal performance in heat sinks. An analysis is presented of conjugate heat transfer in a heat sink with transversely inserted bars in the flow channel. The numerical study seeks to maximise the dimensionless thermal conductance subject to constant total volume. The study contrasts a standard heat sink without a bar with cooling channel inserts with one, two, three, and 4 bars. To improve the surface area for heat transfer, the extended surfaces are placed longitudinally at various points in the cooling channel. A high-density heat flux is applied at the bottom wall of the heat sink and coolant with an inlet velocity vin between 4.26 and 7.26 m/s is pumped in a forced convective laminar condition to eliminate the heat. A computational fluid dynamic (CFD) code incorporated in Ansys fluent is used to solve the mathematical models. The numerical analysis indicates that the configuration featuring 3 bars inserted at the roof of the cooling channel exhibits the highest global thermal conductance, achieving a 53 % increase over other configurations. All configurations with 3 bars inserted up, down, left and right, maintain a constant pump power of 3.5 W across a Reynolds number range of 857–1461. Bars positioned at the roof of the cooling channel experience the lowest outlet and inner wall temperatures, whereas those at the floor exhibit the highest temperatures. The optimised design meets the highest porosity limit of 0.8.

Propensity of 21700 cylindrical cells to thermal runaway under slight overcharging cycling – A numerical study

Slight overcharge of lithium-ion batteries (LIBs) could occur due to inadequate design of battery management system or unexpected malfunction of charger. In comparison with other abuse conditions, the evolution of thermal runaway (TR) under slight overcharging often requires a lengthy incubation period, necessitating relatively long periods for experiments. Numerical simulations with validated models offer a viable alternative to evaluate such risks in modules/packs and associated mitigation measures. Since a current-cutting device is usually installed inside the LIB module/packs, the current will be cut off once the overcharge reaches a certain level, therefore, slight overcharging is of more potential relevance. However, previous simulations mainly focused on TR process induced by continuous overcharging, the conditions associated with slight overcharging have been overlooked. As slight overcharging is difficult to detect, the generated heat could easily accumulate as the cycle process prevails, result in unwanted cell temperature increases with increasing propensity for TR. A three-dimensional (3-D) predictive tool for the transition of cylindrical 21700 cell from slight overcharging cycle to TR has been developed by implementing the published models for heat generation from various thermal decomposition reactions and ISC into in-house version of the opensource computational fluid dynamics (CFD) code OpenFOAM. The code is validated with newly conducted experiments involving both normal and slight overcharging cycles. The validated model has then been used to conduct further parametric studies to investigate the effects of the overcharge voltage, current-rates and convective cooling on the propensity and run up time to TR as well as the maximum cell temperature during TR for cells with different cathode materials under different ambient temperatures. This work forms a sound basis for evaluating TR propagation in modules/packs under slight overcharging/discharging cycles to aid the development of mitigative measures to improve LIB process safety in practical applications.

CFD Modeling of Aerosol Transport and Deposition Using a Drift-Flux Model

Aerosol transport and deposition are important processes in modeling of accident scenarios for a small modular reactor. An aerosol drift-flux model is attractive because it is computationally less expensive than Lagrangian particle tracking. It must be determined, however, how well it performs when implemented in a commercial computational fluid dynamics (CFD) code. This work presents results of modeling aerosol transport and deposition using a full Eulerian three-dimensional drift-flux model implemented in the commercial CFD code STAR-CCM+. The forces due to gravity and thermophoresis are included in the present drift-flux model along with Brownian motion and turbulent diffusion. The forces are added as a source term to a passive scalar transport equation. In addition, a drift velocity representing the forces is used in a built-in electrochemical species transport equation. The results of these two approaches are compared. An appropriate deposition velocity is used to calculate the aerosol concentration deposited on surfaces. The semiempirical relation proposed by Lai and Nazaroff (2000) is used to compute the deposition velocity due to gravitational settling, and the present results are compared with the experimental and numerical data obtained from the work of Chen et al. (2006). It was found that the concentration profile obtained from the present drift-flux model showed reasonable agreement with the literature data. A thermophoresis model showed good agreement when compared with the analytical solution of Nazaroff and Cass (1987). In addition to the particle concentration results, this work presents details of the drift-flux model implementation and the bulk flows. These extra details will enable comparisons by others developing similar models.

Coupled modeling and simulation of tank self-pressurization and thermal stratification

Numerical analysis of thermo-fluid dynamics for cryogenic propellant storage primarily consists of nodal and computational fluid dynamics (CFD) modeling approaches. While nodal approaches prioritize faster computation over fidelity, CFD approaches promise accuracy and fidelity at the expense of significant computational resources. This paper presents an efficient coupling methodology developed to facilitate the simulation of a long-term self-pressurization process of a cryogenic propellant tank as well as associated thermal stratification phenomenon under both normal and reduced gravity conditions. Key highlight of the methodology is the development of a coupling scheme based on a domain decomposition approach that separates the tank into two computational regions at the liquid-vapor interface. SINDA/FLUINT, the nodal code, is utilized to simulate the liquid region, while ANSYS Fluent, the CFD code, handles the vapor region. A data exchange algorithm was developed to compute required boundary conditions for each region using the local thermodynamic properties of assumed infinitely thin interface. The coupling between the nodal and CFD codes was demonstrated by simulating a self-pressurization process in a small size tank using hydrogen as the working fluid. A sensitivity study on the grid sizing and coupling time step was conducted to determine appropriate spatial and ensure a divergent-free explicit data exchange time step size. Using the coupling approach, two numerical case studies were performed to study tank self-pressurization by considering two initial hydrogen fill levels. The effectiveness and efficiency of the coupling methodology was assessed by comparing the temperature and pressure evolution results of the coupled simulation with those obtained solely from the CFD simulation. Results showed that the coupling scheme and data exchange algorithm were successfully implemented, and the coupled simulation agreed well with the CFD simulations. In addition, a stratification number was used to characterize the transient dynamic development of the thermal stratification. Lastly, the computational costs associated with both approaches were compared. Significant reductions in simulation time were achieved for both cases.

Analyses of Jet Buoyant Flow in a Multicompartment Containment Using an Open-Source Solver

Hydrogen mixing and distribution in a nuclear power plant containment during a postulated severe accident is a phenomenon with high safety relevance for current and advanced light water reactors. The Organisation for Economic Co-operation and Development (OECD)/Nuclear Energy Agency (NEA) SESAR Thermal-Hydraulics (SETH) project was carried out to generate experimental data on basic containment phenomena with the required instrumentation to assess computational fluid dynamics (CFD) and lumped parameter codes. In this paper, we analyze the PANDA SETH Test 16 of the OECD/NEA SETH project using containmentFOAM, which is a tailored open-source CFD code package for containment analysis developed at Forschungszentrum Jülich GmbH based on OpenFOAM.® SETH Test 16 addressed steam-air mixture transport driven by jets and plumes in the multicompartment PANDA facility (two vessels and an interconnecting pipe) at the Paul Scherrer Institut in Switzerland. Steam as a wall plume was injected at 0.040 kg/s and 140°C in vessel 1, which was initially filled with dry air at 108°C. The steam-air mixture was vented out from the top of vessel 2. The system pressure was constant at 1.3 bar for the whole test duration (4000 s). To simulate SETH Test 16, the k-ω shear stress transport turbulence model was adopted with the generalized gradient diffusive hypothesis buoyancy turbulence model. The simulation represented the experimental phenomenology reasonably well. The trends for the calculated temperature distribution were similar to those of the experiments and only yielded a slight discrepancy of about 2.7%. In addition, the distribution of the steam molar fraction was well predicted above the injection tube in vessel 1. However, below the injection tube in vessel 1, some discrepancies between the simulations and the experimental results were observed. These may have been caused by the higher turbulent mass diffusivity in the physical model and the challenge of the mesh resolution in the middle region of vessel 1.

A Comparative Study of Different CFD Codes for Fluidized Beds

Fluidized beds are pivotal in the process industry and chemical engineering, with Computational Fluid Dynamics (CFD) playing a crucial role in their design and optimization. Challenges in CFD modeling stem from the scarcity or inconsistency of experimental data for validation, along with the uncertainties introduced by numerous parameters and assumptions across different CFD codes. This study navigates these complexities by comparing simulation results from the open-source MFIX and OpenFOAM, and the commercial ANSYS FLUENT, against experimental data. Utilizing a Eulerian–Eulerian framework and the kinetic theory of granular flow (KTGF), the investigation focuses on solid-phase properties through the classical drag laws of Gidaspow and Syamlal–O’Brien across varied parameters. Findings indicate that ANSYS Fluent, MFiX, and OpenFOAM can achieve reasonable agreement with experimental benchmarks, each showcasing distinct strengths and weaknesses. The study also emphasizes that both the Syamlal–O’Brien and Gidaspow drag models exhibit reasonable agreement with experimental benchmarks across the examined CFD codes, suggesting a moderated sensitivity to the choice of drag model. Moreover, analyses were also carried out for 2D and 3D simulations, revealing that the dimensional approach impacts the predictive accuracy to a certain extent, with both models adapting well to the complexities of each simulation environment. The study highlights the significant influence of restitution coefficients on bed expansion due to their effect on particle–particle collisions, with a value of 0.9 deemed optimal for balancing simulation accuracy and computational efficiency. Conversely, the specularity coefficient, impacting particle–wall interactions, exhibits a more subtle effect on bed dynamics. This finding emphasizes the critical role of carefully choosing these coefficients to effectively simulate the nuanced behaviors of fluidized beds.

CFD Investigations on Heavy Liquid Metal Alternative Target Design for the SPS Beam Dump Facility

This study introduces numerical advancements in an alternative design for the Super Proton Synchrotron (SPS) Beam Dump Facility (BDF) at the European Laboratory for Particle Physics (CERN). The design envisions a high-power operation target made of flowing liquid lead. The proposed BDF is a versatile facility for both beam-dump-like and fixed-target experiments. The target behavior is studied, assuming a proton beam with a momentum of 400 GeV/c, a pulse frequency of 1/7.2 Hz, and an average beam power of 355 kW. Using various Computational Fluid Dynamics (CFD) codes, we evaluate the behavior of liquid lead and predict the thermal stress on the target vessel induced by the pulsed heat source generated by the charged particle beam. The comparison increases the reliability of the results, investigating the dependencies on the CFD modeling approach. The beam is a volumetric heat source with data from the beam-lead interaction simulations provided by the European Laboratory for Particle Physics and obtained with a Monte Carlo code. Velocity field and stress profiles can enhance the design of the lead loop and verify its viability and safety when operated with a liquid metal target.

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The system thermal-hydraulic (STH) code SAM has been coupled to the computational fluid dynamics (CFD) code Simcenter STARCCM+ utilizing a hybrid domain overlapping coupling method in space and an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH standalone models to capture, such as three-dimensional (3D) mixing and thermal stratification. The coupling method’s numerical stability was previously verified against multiple test cases including two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately one-dimensional behavior that is well predicted by a STH code. First, a STAR-CCM+ CFD model of the 3D test section is validated against experimental temperature data. Then a SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states and includes SAM standalone results for comparison. Last, the coupled model is validated against two TALL-3D transients: one exhibiting flow reversal in the test section and one exhibiting nonlinear, limit cycle oscillations (LCOs). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this increase is not predicted by the SAM standalone model. The coupled model also better predicts initial flow recovery and subsequent oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCOs. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. Nevertheless, the coupled model does reproduce the general oscillatory feedback observed. This is a significant achievement because the present coupling only uses an explicit coupling in time, as opposed to a tighter, semi-implicit coupling. In comparison, STH/CFD coupling efforts previously performed by other authors required a semi-implicit coupling to obtain similar results to this work. The strengths of the novel coupling scheme validated in the present paper lie in the superior convergence characteristics and the much simpler, nonintrusive implementation.

High-resolution and reduced-order multi-physics simulation on thermal–hydraulic and thermoelectric characteristics of the thermionic space reactor core

The simulation of a thermionic space reactor core necessitates a comprehensive consideration of the intricate interplay between fluid flow, heat transfer, and thermionic energy conversion, constituting a multi-physics coupling challenge. Given the absence of a thermionic conversion model in contemporary CFD (Computational Fluid Dynamics) codes, accurately simulating the thermal–hydraulic and thermo-electric behavior of thermionic reactors in existing CFD software is particularly challenging. To analyze the multi-physics coupling phenomena in thermionic reactors, this paper revises the conjugate heat transfer solver in the proprietary 3D (three-dimensional) finite-volume CFD code, mFVM (modular Finite-Volume Method Code), to enable bidirectional coupling with the thermionic conversion model. This enhancement aims to predict the high-resolution three-dimensional temperature field and thermoelectric characteristics of the thermionic reactor core. Furthermore, a reduced-order multi-physics model of the thermionic reactor core is established using RESYS (Reactor System Simulation Code). The numerical simulation outcomes for both the high-resolution and reduced-order multi-physics models, pertaining to electrically heated and fission-heated TFEs (Thermionic Fuel Elements), are validated against experimental data. The results indicate that the emitter electron cooling effects caused by thermionic emission reduced the maximum temperature of emitter by about 200 K for both electrically heated and fission-heated TFE. Therefore, a bidirectional coupling between the thermal–hydraulic model and the thermionic conversion model is essential for precisely calculating the temperature field within the TFE and the thermionic reactor core. Thereafter, both the high-resolution model and the reduced-order model are utilized to investigate the thermal–hydraulic and thermoelectric characteristics of the TOPAZ-II reactor core. The calculated electrical power of the TOPAZ-II reactor is 5.39 kWe using the mFVM model and 5.47 kWe using the RESYS model. Finally, an analysis of the thermoelectric conversion characteristics based on the simulation results of the reduced-order model reveals that the TOPAZ-II reactor system is capable of load-following only when the load resistance exceeds 0.06 Ω under nominal full power operation conditions.

Optimizing Coupled Fluid-Structure Simulations for Nuclear-Relevant Geometries

Abstract The numerical simulation of fluid-structure interactions (FSI) has gained interest to study flow-induced vibrations. Nevertheless, the high computational resources required by such simulations can represent a significant limitation for their application to industrial configurations. Therefore, simplified modeling approaches, when physically applicable, can represent an interesting compromise. This can be the case of slender structures (tubes, rods) often encountered in nuclear power plants. In this paper, an Euler–Bernoulli beam finite element model is implemented inside the computational fluid dynamics (CFD) code code_Saturne. With the goal of finding CFD methods less expensive than large eddy simulations (LES), unsteady Reynolds Navier–Stokes (URANS) and hybrid URANS/LES approaches are considered. The resulting fluid-structure model is able to calculate the vibration response of cantilever beams under a fluid flow, avoiding the necessity of CFD-finite element method (FEM) code coupling. The first part of the paper describes the model and its implementation: it allows to perform 2-way explicit fluid-structure coupling, using the Arbitrary Lagrangian-Eulerian approach to account for the structure deformations. Validation test cases are presented in the second part: first, the model is validated in terms of frequency, added mass, and damping for a cylinder vibrating in static air and water; then, the model is validated toward the vortex-induced resonance and lock-in mechanisms for a cylinder subjected to water cross-flow. The model is then applied to a real experimental configuration of two in-line cylinders in water cross-flow: the calculated vibrations are found to be in good agreement with the experimental measurements.

Research of VDT scheme for porous media thermal-hydraulics analysis of plate-type fuel assembly

Thermal-hydraulics analysis of nuclear reactor core is crucial for the safe operation of nuclear power plants. Porous media models in computational fluid dynamics (CFD) code are commonly used to simulate flow and heat transfer in nuclear reactor cores because they simplify geometry and reduce computational cost. However, their low spatial resolution limits detailed analysis of internal flow and temperature fields. This research develops a Virtual Domain Technology (VDT) scheme to improve the resolution of porous media models without increasing mesh density. VDT sets localized resistance coefficients to create virtual solid and fluid domains within the porous core representation. The flow and thermal effects of these virtual domains mimic real structural effects. VDT was verified on a cylinder flow case and showed excellent agreement with CFD for velocity and pressure fields. Application to a plate-type research reactor core demonstrated VDT’s ability to capture detailed temperature profiles and flow traces consistent with CFD, using only 6% as many mesh cells. By improving spatial resolution, VDT enhances porous media-based modeling of reactor cores to provide more accurate thermal-hydraulics analysis, while maintaining computational efficiency. This will aid optimization and safety analysis for advanced reactor designs.

The Effects of a Seagull Airfoil on the Aerodynamic Performance of a Small Wind Turbine

Birds’ flight characteristics such as gliding and dynamic soaring have inspired various optimizations and designs of wind turbines. The implementation of biological wing geometries such as the airfoil profile of seabirds has improved wind turbine performance. However, the field can still benefit from further investigation into the aerodynamic characteristics of an inspired design. Therefore, this study evaluated the effect of a seagull airfoil design on the aerodynamic performance of the National Renewable Energy Laboratory (NREL) Phase VI wind turbine. By replacing its S809 airfoil with the laser-scanned profile of the seagull airfoil, the aerodynamic behavior at key locations of the NREL Phase VI wind turbine blade was numerically simulated in a three-dimensional environment using the Ansys Fluent 2022 R1 computational fluid dynamics (CFD) code. The results were validated against the experimental data, and analysis of the torque outputs, pressure distributions, and velocity profiles that were generated by both the baseline and modified models demonstrated the ability of the seagull airfoil profile to modify regions of minimum and maximum local velocities to achieve highly favorable pressure differentials, significantly increasing the torque output of the NREL Phase VI wind turbine by 350, 539, 823, and 577 Nm at 10, 15, 20, and 25 m/s inlet velocities, respectively.

Numerical investigation of vortex induced vibrations on cylinders

The paper addresses Vortex Induced Vibrations (VIV) caused by tower vortex shedding through LES simulations of a wind tunnel experiment. A high-fidelity (HiFi), Fluid Structure Interaction (FSI) framework for aeroelastic analysis of 3D flexible cylinders is established, based on the open-source Computational Fluid Dynamics (CFD) code OpenFOAM and the coupling library preCICE. In the present study, the FSI framework is employed with the aim to replicate a low Reynolds (Re=26000) forced vibrations wind tunnel experiment of a circular cross-section rigid cylinder. The paper compares predictions of the aerodynamic damping of the oscillating cylinder for different oscillation frequencies within the lock-in range and amplitudes against experimental results. While the model consistently predicts the lock-in range of frequencies and the trends in the variation of the aerodynamic damping with the oscillation amplitude, it falls short of accurately capturing the maximum value of the aerodynamic damping at the onset of lock-in.

Enhancing performance in solar air channels: A numerical analysis of turbulent flow and heat transfer with novel shaped baffles

Solar air heaters are emerging as a promising solution for sustainable heating, as they harness abundant solar energy to provide heat for a variety of applications. However, their optimum efficiency remains a major challenge. Deflectors have significant potential to improve the performance of solar air heaters by redirecting airflow and optimizing heat transfer, making them more efficient and more suitable for large-scale adoption. This article examines the performance of a new baffle design aimed at improving heat transfer in the channel by increasing the exchange surface between fluid and solid surface. It consists of a “≠” -shaped baffle mounted on the top and bottom walls of the channel. An in-depth analysis of the spacing between these baffles was carried out to assess their impact on the thermal and aerodynamic properties of the solar air channel. Computational Fluid Dynamics (CFD) was used to perform the calculations, using the Finite Volume Method (FVM) in conjunction with the SIMPLE algorithm. All studies are carried out using the CFD code Fluent. Crucial steps in this research include studying the influence of these new baffles at different Reynolds numbers, ranging from 15000 to 35000, as well as varying the separation distance between the two baffles, chosen periodically (DBaffle/2, 3DBaffle/4, DBaffle, 5DBaffle/4 and 3DBaffle/2). Some of the numerical results were validated with the available experimental data and showed satisfactory agreement. The study identified an optimal case for the new baffles, by reducing the distance between the baffles, an improvement in system performance was observed, by increasing the flow strength and widening the recycling cells, which has an impact on dynamic behavior and heat transfer. The presence of this new baffle design, with minimal separation distance between them, demonstrated a significant thermal improvement, with a thermal improvement factor increasing by 44% compared to the simple channel without baffles, thus highlighting the effectiveness of the new solar air channel with “≠” -shaped baffles.