Drafting is a common technique to reduce the drag experienced by elite runners on races, leading to better running efficiency and energy savings. The tactic has been successfully used in previous races where world records were set, especially by Eliud Kipchoge in the 2018 and 2022 editions of the Berlin marathon. In the 2023 edition of the Chicago Marathon, Kenyan runner Kelvin Kiptum broke the Marathon record after a 2:00:35 finish. This feat becomes even more impressive by the lack of use of drafting formations, despite the availability of two pacers for the majority of the race. In this study, the drag faced by Kiptum and his pacers during the race is calculated by means of computational fluid dynamics (CFD). The performance of each runner is evaluated from an energetic standpoint, and the analysis is extended to include more efficient drafting formations. Our results indicate that, by properly using the advantages of drafting, Kiptum could have finished the race at a staggering 1:57:34, a full three minutes better than his own record and 215 seconds better than the previous world record.

The present is a study of the CFD simulations intended to simulate the emissions of pollutants that are generated after the combustion of proposed alternative aircraft fuels (Hydrogen, Ethanol and Methane) to compare with the emissions generated after the combustion of Kerosene and Benzene in a 2D cylindrical combustion chamber. Given that air traffic is a main contributor to not only 3% of man-made greenhouse effects but also of the generation of smog over heavy air traffic urban areas generating an impact on the air quality and the population of those areas.

This paper presents some results from a computational fluid dynamics (CFD) model of a multi-megawatt crosswind kite spinning on a circular path in a straight downwind configuration. The unsteady Reynolds averaged Navier-Stokes equations closed by the k−ω SST turbulence model are solved in the three-dimensional space using ANSYS Fluent. The flow behaviour is examined at the rotation plane, and the overall (or global) induction factor is obtained by getting the weighted average of induction factors on multiple annuli over the swept area. The wake flow behaviour is also discussed in some details using velocity and pressure contour plots. In addition to the CFD model, an analytical model for calculating the average flow velocity and radii of the annular wake downstream of the kite is developed. The model is formulated based on the widely-used Jensen’s model (Technical Report Risø-M; No. 2411, 1983), which was developed for conventional wind turbines, and thus has a simple form. Expressions for the dimensionless wake flow velocity and wake radii are obtained by assuming self-similarity of flow velocity and linear wake expansion. Comparisons are made between numerical results from the analytical model and those from the CFD simulation. The level of agreement was found to be reasonably good. Such computational and analytical models are indispensable for kite farm layout design and optimization, where aerodynamic interactions between kites should be considered.

The pulsating flow is one of the techniques which can enhance heat transfer, therefore leading to energy saving in tubular heat exchangers. This paper investigates the heat transfer and flow characteristics in a two-dimensional in-line tube bundle with the pulsating flow by a numerical method using the Ansys Fluent. Numerical simulation is performed for Reynolds number Re = 500 with different frequencies and amplitude of pulsation. Heat transfer enhancement was estimated from the central tube of the tube bundle. Pulsation velocity had an asymmetrical character with a reciprocating flow. The technique developed by the authors to obtain asymmetric pulsations was used. This technique allows simulating an asymmetric flow in heat exchangers equipped with a pulsation generation system. Increase in both the amplitude and the frequency of the pulsations has a significant effect on heat transfer enhancement. Heat transfer enhancement is mainly observed in the front and back of the cylinder. At a steady flow in these areas, heat transfer is minimal due to the weak circulation of the flow. The increase in heat transfer in the front and back of the cylinder is associated with increased velocity and additional flow mixing in these areas.

In the present study, a numerical model was developed to predict the flow pattern inside a broiler building. The model intends to predict the velocities fields inside the domain and ammonia (NH3) emitted or released by litter from poultry housing. The numerical model developed in Computational Fluid Dynamics (CFD) commercial code, intends to represent a commercial broiler building, and intends to simulate the 3D and heat transfer, in steady state flow. The evaporative cooling pads were also included in the model. The validation of the model was based in experimental measurements obtained in previous studies. The simulations were focused on Summer, Winter and also Mid-Season situation. The numerical results of NH3 concentration were compared with the experimental measurements, and a quite good agreement was verified. The numerical results allowed the characterization of: the inside flow pattern developed for the summer and winter situation; the inside NH3 distribution, and the velocity field distribution inside the broiler building. It was found that NH3 concentration increases along the tunnel, especially in low flow rate imposed from the exhaust fan. Also, it was verified that the low velocities inside domain are no sufficient to remove the gaseous pollutants.

Solid asperity interactions are common and inevitable under severe loading conditions for any lubricated contact. Heavy-duty machine components (gear, bearings etc.) generally operate under Mixed Lubrication, where uneven surface features contact each other when generated fluid pressure is not enough to support the external load. The Reynolds equation is commonly used to simulate smooth lubricated contacts numerically. In rough lubricated interfaces where opposite surface asperities make contact, the Reynolds equation alone cannot accurately predict pressure using the traditional numerical simulation method. In this paper, Lubrication-Contact interface conditions (LCICs) have been implemented and extended to solve the multiple asperity contact problem using the full-multigrid approach. The developed novel algorithm has significantly accelerated the solution process and improved the accuracy and efficiency of pressure calculation for fluid-solid sub-interactions that can occur under ML regions. The Finite Difference Method (FDM) results have been compared with Computational Fluid Dynamics (CFD) simulation to validate the newly developed model. Hence, the proposed optimized solution method will provide valuable insight to researchers and industry engineers interested in simulating the ML problem where the effect of the fluid-solid interface can be captured effectively to improve reliability in the calculation of the life expectancy of the lubricated parts.

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.

Craters are the most prominent geomorphological features on the surface of celestial bodies, which is playing a crucial role in studying the formation and evolution of celestial bodies, as well as in landing and planning for surface exploration. Currently, the main automatic crater detection models and datasets focus on the detection of large and medium craters. In this paper, we created 23 small lunar crater datasets for model training based on the Chang’E-2 (CE-2) DOM, DEM, Slope, and integrated data with 7 kinds of visualization stretching methods. And then, we proposed the YOLO-Crater model for Lunar and Martian small crater detection by replacing EioU and VariFocal loss to solve the crater samples imbalance problem and introducing a CBAM attention mechanism to mitigate interference from the complex extraterrestrial environment. The results show that the accuracy (P = 87.86%, R = 66.04%, and F1 = 75.41%) of the Lunar YOLO-Crater model based on the DOM-MMS (Maximum-Minimum Stretching) dataset is the highest and better than that of YOLOX model. And the Martian YOLO-Crater, trained by the Martian dataset from the 2022 GeoAI Martian Challenge, achieves good performance with P = 88.37%, R = 69.25%, and F1 = 77.65%. It indicates that the YOLO-Crater model has a strong transferability and generalization capability, which can be applied to detect small craters on the Moon and other celestial bodies.

The output of non-ferrous metal smelters is 450,000 tons (3.1% of world output), and copper ranks 10th in the world for health impact evaluation. According to the latest information from smelters, the overall average concentration in 2020 was 8.2 ppb each of SO2, NO2, O3, PM10, and PM2.5. 7.5 ppb, 30.8 ppb, 21.3 μg/m3, and 13.2 μg/m3, and the maximum concentrations were 390 ppb, 50 ppb, 167 ppb, 148 μg/m3, and 107 μg/m3. Art Corporation, A2C (Atmosphere to CFD) model was used to investigate the resident health follow-up management (cohort study) study (SAS Enterprise Guide) as the basis for the source of diffuse pollution in the area. The average concentration of PAHs (polynuclear aromatic hydrocarbons) among carcinogens in the air by survey site was the highest in Seokpo 4-ri (60.7 ng/m3), followed by Seokpo 1-ri (48.6 ng/m3) and Seokpo. Among the respiratory diseases, the disease (Acute upper respiratory tract disease(J00-J06), Other upper respiratory tract diseases (J32-J39), Acute lower respiratory tract infection (excluding pneumonia) (J20-J22), Chronic lower respiratory disease (persons) (excluding J40-J47, J45-J46), Asthma (J45-J46), Rhinitis (J30-J31), Respiratory disease (J00-J99), Cough (R05)) that has a correlation between the exposed area and the control area is as follows.

This paper investigates urban microclimate in relation to the typical major streets and mobility layout in Turin, i.e. boulevards with lateral minor roads, called “controviali”, an almost unique urban mobility layout. The aim is to examine the possibility of intervention in relation to the type and scope of mobility for microclimate mitigation and climate adaptation purposes. To this, this paper evaluates to which extent urban morphology and mobility features contribute to determining microclimatic conditions, to finally identify site-specific prioritised strategies to cope with its mitigation. This research is structured in two phases. The first phase is dedicated to urban-scale analysis of the streets with “controviali”, resulting in the creation of a database containing their properties and typology, to finally derive six clusters of boulevards. The second phase involves the selection of one representative case study per cluster and simulation of their microclimate conditions, performed by the software ENVI-met. The results are discussed in relation to the parameters proven in the prior literature review to have a major impact on the microclimate. The assessment of the results indicated which elements affect the climatic conditions the most, allowing for the proposal of different prioritised mitigation strategies in the conclusion.

This study aims to investigate the impact of body heat loss on the thermal and aerodynamic conditions in a whole-body cryotherapy chamber. The underlying hypothesis is that the heat generated by the human body alters the thermal and aerodynamic environment inside the cabin. A numerical study was conducted to test this hypothesis to analyze the thermodynamic exchanges between the human body and the cabin during a 3-minute whole-body cryotherapy session. The computational fluid dynamics (CFD) approach was used to study the unsteady heat transfer between the human body and the interior of the cryotherapy cabin. A thermal boundary condition, based on a mathematical model developed from experimental data, was applied to simulate skin cooling kinetics over time. The post-processing of the 3D results, including temperature, velocity fields, and thermal flux maps at the body surface, provided insight into the thermo-convective mechanisms involved in a whole-body cryotherapy session. The study found that body heat loss significantly affects the temperature fields inside the cabin, leading to global modifications of the aeraulic and thermal conditions. These findings suggest that cryotherapy protocols may need to be adjusted or the cabin set temperature optimized to enhance the therapeutic benefits.

The automotive industry has witnessed a significant push towards increased fuel efficiency and reduced carbon emissions. Among many factors influencing overall vehicle efficiency, aerodynamic drag is a critical consideration, accounting for a substantial portion of fuel consumption. Our investigation focused on the drag coefficients of three prominent 2023 passenger car models – the Toyota Corolla, Hyundai Elantra, and Kia Soul. Using ANSYS Fluent to perform computational fluid dynamics, we demonstrated an inverse correlation between drag force and gas mileage, underlining the significance of aerodynamic design in fuel efficiency. Notably, the Hyundai Elantra emerged as the most aerodynamically efficient of the three vehicles, followed closely by the Toyota Corolla. In contrast, the boxier-designed Kia Soul exhibited significantly higher drag coefficients and drag forces, rendering it the least efficient of the trio. These findings emphasize the importance of environmental factors in car design. Shaping aerodynamic profiles directly impacts fuel efficiency, stability, and overall performance.

Tissue engineered heart valves can grow, repair, and remodel after implantation, presenting a more favorable long-term solution compared to mechanical and porcine valves. Achieving functional engineered valve tissue requires the maturation of human cells seeded onto valve scaffolds under favorable growth conditions in bioreactors. The mechanical stress and strain on developing valve tissue caused by different pressure and flow conditions in bioreactors are currently unknown. The aim of this study is to quantify the wall shear stress (WSS) magnitude in heart valve prostheses under different valve geometries and bioreactor flow rates. To achieve this, the study used fluid-structure interaction simulations to obtain the valve's opening geometries during the systolic phase. These geometries were then used in computational fluid dynamics simulations with refined near-wall mesh elements and ranges of prescribed inlet flow rates. The data obtained included histograms and regression curves that characterized the distribution, peak, and median WSS for various flow rates and opening configurations. The study also found that the upper region of the valve near the commissures experienced higher WSS magnitudes than the rest of the valve.

This paper presents the results of simulations of particle-laden air-solid jet flow in long straight tubes using CFD-DEM, along with an analysis of coarse-graining. Although previous studies have used CFD-DEM for similar flows, these have typically been in a dilute flow regime where uncoupled simulations are accurate. However, fully coupled simulations can introduce issues, necessitating validation studies to ensure that all coupling parameters are effectively used and that the physics is accurately represented. This paper validates the simulations against two different experimental studies, with fluid Reynolds numbers between 10,000 and 40,000 and Stokes numbers between 5.6 and 50. Interestingly, the profiles of mean particle velocity exhibit fewer discrepancies as the Stokes number increases but more discrepancies for the root mean square velocity compared to the experiments. Particle number flux is consistent with the experiments after the nozzle exit. Coarse-graining is also applied to the same simulations, achieving relatively accurate results. However, as expected, scaling of contact collision frequencies, forces, and stresses cannot be achieved, meaning that coarse-graining may be useful for comparing designs or operating parameters on an industrial scale but falls short when measuring the total energy dissipation of one experiment.

As the need to minimize environmental impact continues to rise, it is essential to incorporate, advance, and adopt renewable energy sources and materials to attain climate neutrality. It is established that economic growth is built upon infrastructure where cement industry plays a crucial role. However, it is also known that this industry is actively looking for ways to transition towards low-carbon practices to encourage sustainable and environmentally conscious practices. To this end, the use of refused derived fuels to substitute fossil fuels is very appealing, as these present the potential to lower the clinker production costs and CO2 emissions. Bearing this in mind, the primary objective of this work is to gain insights of the combustion behavior in an industrial rotary kiln, studying real-life scenarios and assessing the potential of substituting alternative fuels for fossil fuels to reduce CO2 emissions. A 3D CFD turbulent combustion model was formulated in Ansys® considering a Pillard NOVAFLAM® burner, where refused derived and petcoke fuels were used and different secondary air mass flows were considered. From the obtained results it was possible to conclude that the outcome of the combustion process is greatly influenced by the fuel-to-air ratio. Increasing the secondary air mass flow promotes the occurrence of a complete and efficient combustion process, leading to an enhanced fuel conversion and decreased formation of pollutants as CO, soot, and unburnt hydrocarbons.

This article presents a detailed examination of the methodology and modeling tools utilized to analyze gas flow in pipelines, rooted in the fundamental principles of gas dynamics. The methodology integrates numerical simulations with modern neural network techniques, particularly focusing on the Physics-Informed Neural Network (PINN) method. This innovative approach combines artificial neural networks (ANNs) with physical equations, offering a more efficient and accurate way to model various complex processes and phenomena. The proposed mathematical model, based on the Euler equation, has been meticulously implemented using the Python language. Verification with analytical solutions ensures the accuracy and reliability of the computations. In the research, a comprehensive comparative analysis was conducted between results obtained using the PINN method and those from conventional Computational Fluid Dynamics (CFD) approaches. The analysis highlighted the advantages of the PINN method, which produced smoother pressure and velocity fluctuation profiles while reducing computation time, demonstrating its potential as a transformative modeling tool. The data derived from this study are of paramount importance for ensuring ongoing energy supply reliability and can also be used to create predictive models related to gas behavior in pipelines. The application of modeling techniques for gas flow simulation has the potential to revolutionize the integrity of our energy infrastructure and utilization of gas resources. However, it is crucial to emphasize that the effectiveness of such models relies on continuous monitoring and frequent updates to ensure alignment with real-world conditions. This research not only contributes to a deeper understanding of compressible gas flows but also underscores the crucial role of advanced modeling methodologies in the sustainable management of gas resources for both current and future generations.

Plate fin-tube heat exchangers are widely used in air conditioning and refrigeration systems and other industry fields. Various errors made in the manufacturing process can result in the formation of an air gap between the tube and fin. Several numerical simulations were carried out for a symmetric section of plate fin-tube heat exchanger to study the influence of air gap on heat transfer under periodic flow conditions. Different locations and sizes of an air gap spanning 1/2 circumference of the tube were considered for the range of airflow velocities. Velocity and temperature fields for cases with air gap were compared with ideal thermal contact cases. Blocking of heat flow by the gap leads to the reduction of heat transfer rate. Fin discontinuity in the front of the tube causes the smallest reduction of the heat transfer rate in comparison to the ideal tube-fin contact, especially for thin slits. The rear gap position is the worst in the smallest gap range. Therefore, reversing the flow direction can lead to up to a 15% heat transfer increase, if mainly the rear gaps are present. The introduction of a thin slit in the front of the tube leads to convective heat transfer enhancement, which should be further investigated.

The second-order factor effect of burner optical ports and edge inter-matrices (EIM) and the first-order factor of pressure on the soot formation process and behavior of premixed sooting flames in a high-pressure burner are numerically investigated here. Three-dimensional computational fluid dynamics (CFD) simulations of a premixed flame C₂H₄/air at p = 1.01 and 10 bar using a one-step chemistry approach are first performed to justify the satisfied predictability of the prospective axisymmetric two-dimensional (2D) and one-dimensional (1D) simulations. The justified 2D simulation approach shows the generation of an axial vorticity around the EIM and axial multi-vorticities due to the high expansion rate of burnt gases at the high pressure of 10 bar. This leads to the development of axial multi-sooting zones, which are manifested experimentally by visible luminous soot streaks, and to the boosting of soot formation conditions of a relatively low-temperature field, < 1800 K, and a high mixing rate of gases in combustion around and above the EIM location. Nevertheless, a tolerable effect on the centerline soot volume fraction (f_V) profile, f_V < 3%, is manifested only at high heights above the burner of the atmospheric sooting flame C₂H₄/air ϕ = 2.1, and early at the high pressure of 10 bar of this flame, f_V < 10%. Enhancing the combustion process reactivity by decreasing the rich equivalence ratio of the fuel/air mixture and/or rising the pressure results in the prior formation of soot precursors, which shifts the sooting zone upstream.

This paper attempts to develop a Dynamic Mode Decomposition (DMD) based Reduced Order Model (ROMs) which can quickly but accurately predict the forces and moments experienced by a road vehicle such that they be used by an on-board controller to determine the vehicle's trajectory. DMD can linearize a large dataset of high-dimensional measurements by decomposing them into low-dimensional coherent structures and associated time-dynamics. This ROM then also can be applied to predict the future state of the fluid flow. Existing literature on DMD is limited to low Reynolds number applications. This paper presents DMD analyses of the flow around an idealized road vehicle, called the Ahmed body, at a Reynolds number of 2.7E6. The high dimensional dataset used in this paper was collected from a computational fluid dynamics (CFD) simulation performed using Menter's shear stress transport (SST) turbulence model within the context of Improved Delayed Detached Eddy Simulations (IDDES). The DMD algorithm, as available in literature, was found to suffer nonphysical dampening of the medium-to-high frequency modes. Enhancements to the existing algorithm were explored, and a modified DMD approach is presented in this paper which includes: (a) a requirement of higher sampling rate to obtain higher resolution of data, and (b) a custom filtration process to remove spurious modes. The modified DMD algorithm thus developed was applied to the high Reynolds number, separation dominated flow past the idealized ground vehicle. The effectiveness of the modified algorithm was tested by comparing future predictions of force and moment coefficients as predicted by the DMD-based ROM to the reference CFD simulation data, and found to offer significant improvement.

Marine vessels operating on the Arctic Sea route are constantly prone to collisions and friction with ice. This study discusses the wear of the hull plate caused by the collision of ice against vessels operating in Arctic Sea routes. The abrasive wear of the hull due to ice impact was numerically assessed based on both the incident behavior of ice particles interacting with the flow around the hull and the wear loss of the hull surface caused by the contact force of ice particles. A multi-phase approach was adopted to account for the behavior of ice particles continuously affected by the fluid force around the hull. The fluid force acting on the ice floe was evaluated using computational fluid dynamics (CFD) and the dynamic motion of the drift ice was evaluated using the discrete element method (DEM). The motion of the floating ice particles was updated in real time by iteratively coupling the fluid force and the motion of the ice floe at each time step of the numerical simulation. The results of the wear simulation models were presented in terms of the shape change of the hull surface due to wear. At first, the wear was evaluated for cases in which only the surface paint of the hull was damaged. Thereafter, a computation model considering the shape change of the hull surface experiencing long-term friction of ice particles was introduced. Finally, the numerical procedures to predict the abrasive wear of the hull surface by ice impact were discussed.

Backgroumd: The feeling of nasal obstruction is one of the most complaints of the Ear Nose and Throats (ENT) patients. The reason of this complaint is variable, such as Nasal Septal Deviation (NSD). NSD is a significant obstacle to normal respiration and the volume of air entering the lungs. Understanding how the deviation causes the obstruction is so important, because it leads to better management of this disease. Rhinomanometry and acoustic rhinometry are methods that are used for objective evaluation of upper airway, but they are semi invasive and the accuracy are not guaranteed. With development of Computed Tomography (CT scan) and Computational Fluid Dynamics (CFD), numerical studies on airflow in different parts of the nasal cavity have been performed. Objectives: The purpose of this study was to compare nasal airflow patterns between patient with NSD and healthy nasal cavity using a computational fluid dynamics (CFD) Methods: A patient with NSD and a healthy human were selected. Using CFD, the maximal airflow velocity, nasal resistance, maximal Wall Shear Stress (WSS), and pressure drop in bilateral nostrils of both cases were assessed. Results: Pressure and velocity profiles are evaluated at a mass flow rate of 17.4 Liters Per Minute (LPM). The velocity profile shows an abnormal flow profile in the septal deviated nasal cavity and the peak velocity is observed at the middle nasal region, whereas this point is at the nasal valve level in a healthy nasal cavity. Wall shear stress (WSS) caused by NSD is lower than in normal nasal cavities. The pressure drop in a nasal cavity with deviated septum is higher than in a healthy nasal cavity. Based on the results of this investigation nasal resistance in NSD was higher than in normal, healthy nasal cavity. Conclusions: NSD plays an important role in nasal airflow regime and it results in the feeling of obstruction comparing to normal nasal cavity. Also there is a promising future for CFD to become a useful tool in understanding airflow physiology and pathology.

The envelope (façade) of a building is the first barrier between the exterior and the interior of the building and withstands the highest variation of temperature and solar radiation. Trombe walls are used to take advantage of that and contribute to the heating of interior air, helping the heating system. In this study, a new Trombe wall design is presented with the aim to contribute significantly to the ventilation of residential buildings. For this an exterior wall equipped with a proposed Trombe system was studied in a numerical simulation analysis. The proposed systems consist of two parts, the air cavities on the exterior of the wall, covered with silica glass, and another layer used for heat transfer, on the interior, comprised of a MPCLBs wall of 115mm. Two cities in Romania were considered to represent the coldest and hottest climate in Romania. Brașov representing the cold climate and Constanța representing the hot climate. In the investigation, both the presence and absence of solar radiation were taken into account. In total 4 cases were established. For all four cases monthly research was made using monthly means temperature, direct and diffuse solar radiation. The results are promising and illustrate that the system works the best during the transitional seasons of spring and autumn. The lower the outdoor temperature the higher the temperature can be increased. Overall, solar radiation accounted for an average of 2 °C increase. The temperature increase varied between 3.4 °C to 15.99 °C for Brașov and between 6.42 °C to 12.07 °C in Constanța.

An aneurysm is a vascular malformation, it can be classified according to its location (cerebral, aortic); and according its shape (saccular, fusiform, and mycotic). The study of blood flow interaction with the aneurysms has gained attention from physicists an engineers in recent years. Shear stresses, oscillatory shear index (OSI), gradient oscillatory number (GON), and residence time had been used as potential variables to describe the hemodynamics as well as the origin and evolution of aneurysm. However, the causes of aneurysms and the hemodynamics conditions that promote the growth of these structures are still under debate. The present work presents numerical simulation on three aneurysms: 2 aortic and 1 cerebral. The results shows that the blood rheology is not relevant for aortic aneurysms, however, for cerebral aneurysm could play a relevant role in the hemodynamics. The turbulence models evaluated show equivalent results in both cases. Lastly, a simulation considering the fluid-structure interaction (FSI) shows that this phenomenon is the dominant factor for aneurysms simulation.

Water resource pollution, particularly in river channels, presents a grave environmental challenge that necessitates a comprehensive and systematic approach encompassing assessment, forecasting, and effective management. This article provides a comprehensive exploration of the methodology and modeling tools employed to scrutinize the process of river channel pollution due to silting, rooted in the fundamental principles of hydrodynamics and pollutant transport dynamics. The study's methodology seamlessly integrates numerical simulations with state-of-the-art neural network techniques, with a specific focus on the Physics-informed neural networks (PINN) method. This innovative approach represents a groundbreaking fusion of artificial neural networks (ANN) and physical equations, offering a more efficient and precise means of modeling a wide array of complex processes and phenomena. The proposed mathematical model, grounded in the Euler equation, has been meticulously implemented using the Ansys Fluent software package, ensuring accuracy and reliability in the computations. In a pivotal phase of the research, a thorough comparative analysis was conducted between the results derived from the PINN method and those obtained through conventional numerical approaches using the Ansys Fluent software package. The outcomes of this analysis revealed the superior performance of the PINN method, characterized by the generation of smoother pressure fluctuation profiles and significantly reduced computation time, underscoring its potential as a transformative modeling tool. The calculated data originating from this study assumes paramount significance in the ongoing battle against river sedimentation. Beyond this immediate application, these findings also serve as a valuable resource for creating predictive materials pertaining to river channel silting, thereby empowering decision-makers and environmental stakeholders with essential information. The utilization of modeling techniques to address pollution concerns in river channels holds the potential to revolutionize risk management and safeguard the integrity of our vital water resources. However, it is imperative to underscore that the effectiveness of such models hinges on ongoing monitoring and frequent data updates, ensuring that they remain aligned with real-world conditions. This research not only contributes to the enhanced understanding and proactive management of river channel pollution due to silting but also underscores the pivotal role of advanced modeling methodologies in the preservation of our invaluable water resources for present and future generations.

Snow cover plays a critical role in the global energy and water cycles. Snow cover on the Tibetan Plateau (TP) provides vital water sources in western China and Himalayan regions in addition to its weather and climate significance. The massive high mountain topography of the TP is main conditions for the presence and persistence of snow cover on the plateau at the mid-low latitudes of the Northern Hemisphere (NH). However, how mountain topography controls snow cover distribution on the TP is largely remain unclear and the relationship is not well quantified. Here, the spatial distribution of snow cover and topographic controls on snow cover on the TP are examined based on snow cover frequency (SCF) derived from MODIS snow cover product (MOD10A2 v5) and digital elevation model (DEM). The results show that snow cover on the TP is spatially unevenly distributed and is characterized by rich snow and high SCF on the interior and surrounding high mountain ranges, and less snow and low SCF in inland basins and river valleys. Snow cover on the TP presents elevation dependence with the higher the altitude, the higher the SCF, the longer the snow cover duration and the more stable the intra-annual variation. Annual mean SCF below 3000 m above sea level (masl) is less than 4% and it reaches 77% above 6000 masl. The intra-annual snow cover variation below 4000 masl features a unimodal distribution, while above 4000 masl it presents a bimodal distribution. The mean minimum SCF below 6000 masl occurs in summer, while above 6000 masl it occurs in winter. Because of differences in solar radiation and water vapor sources, mean SCF generally increases with mountain slopes and it is the highest on the north-facing aspect, whereas the lowest is observed on the south-facing aspect.

The four-way reversing valve is a key component determining whether a heat pump air conditioning system can work. Studying the four-way reversing valve is of great significance for improving the performance of the heat pump air conditioning system. In this paper, the control equation of the sliding block movement process is constructed to simulate the reversing process of the four-way reversing valve in fluid mechanics, and a detailed study is carried out from the sliding block speed, sliding block acceleration, sliding block force, refrigerant leakage, and other aspects. Research has shown that optimizing the structure of the slider can greatly improve the working performance of the four-way reversing valve. Increasing the height of the slider not only improves the pressure relief ability of the four-way reversing valve and prevents the slider from being damaged by high-pressure gas impact, but also reduces the leakage of refrigerant during the directional process; When reducing the length of the slider, although it will improve the pressure relief capacity of the four-way reversing valve, it also increases the amount of refrigerant leakage. The research results provide a theoretical basis and improvement direction for designers to improve the performance of four-way reversing valves.

Wall Shear Stress (WSS) abnormalities in carotid artery blood flow are key hemodynamic factors leading to arterial endothelial damage, atherosclerotic plaque formation, and carotid artery stenosis. Precise simulation of carotid artery blood flow WSS in vitro play an important role in understanding the underlying hemodynamic mechanism of carotid artery stenosis. This study presents a design for a microfluidic cell culture chamber with a variable cross-section and a step structure, leveraging microfluidic technology and principles of fluid dynamics. We propose an optimization strategy for microchannel dimensions based on Computational Fluid Dynamics (CFD) simulations. This strategy aims to replicate the distinct WSS waveform features found at various locations within the carotid bifurcation as observed in vivo. Our simulations indicate that with a step height of 0.09 mm, a width of 4 mm for the initial segment, and 5 mm for the subsequent wider channel, the device can effectively model the characteristics of low oscillatory WSS seen at sites of carotid sinus stenosis. Furthermore, it accurately represents the high-magnitude pulsatile WSS waveforms in the more uniform arterial sections downstream of the carotid sinus. The study also reveals that vortex formation and variations in low oscillatory WSS within the stepped section of the microchannel correlate with the step's height and width, as well as with the dimensions of the wider channel section. An increase in step width leads to a decrease in both the vortex area and the magnitude of negative WSS oscillations. When the step height is less than 0.03 mm, vortex formation is inhibited, challenging the simulation of low oscillatory WSS patterns. The variable cross-section microfluidic channel developed in this study provides a platform for simulating carotid artery WSS waveforms and facilitates research into the interplay between low oscillatory WSS at sites of carotid sinus stenosis and arterial endothelial dysfunction.

Hydrogen storage in high pressure tanks can occur with different filling strategies. Many studies have been carried out on supplies with increasing pressure rates. However, hydrogen filling with constant mass inflow rate have not received the same attention. The present work aims to carry out numerical simulations to study hydrogen flow, in a type 3 tank of 70 MPa NWP and constant inflow rate condition. From the analysis of the results, it was possible to verify that there is a linear relationship between the increase in the total equilibrium temperature and the final hydrogen temperature. There is also a linear increasing relationship between Inlet and final temperatures when the equilibrium temperature is fixed. Considering fully adiabatic walls results in significant increases in temperature. The difference between total and static temperatures found in the Inlet tube decreases with the decrease in mass flow rate. The choice of a polymer for the tank lining led to higher temperatures when compared to an aluminum lining.

Machine learning (ML) based on neural network (NN) facilitates data-driven techniques for handling large amounts of data, either obtained through experiments or simulations at multiple spatio-temporal scales, thereby finding the hidden patterns underlying these data and promoting efficient research methods. The main purpose of this paper is to extend the capabilities of a new solver, called realFluidReactingNNFoam, under development at University of Perugia in OpenFOAM with NN algorithm for replacing complex real-fluid thermophysical property evaluations, using the approach of coupling OpenFOAM and Python-trained NN models. Currently, NN models are trained against data generated using the Peng-Robinson equation of state (PR-EoS) assuming mixture frozen temperature. The OpenFOAM solver, where needed, calls the NN models in each grid-cell with appropriate inputs, and the returned results are used and stored in suitable OpenFOAM data structures. Such inference for the thermophysical properties is achieved via the “Neural Network Inference in C made Easy” (NNICE) library, which proved to be very efficient and robust. The overall model is validated considering a liquid-rocket benchmark comprised of liquid-oxygen (LOX) and gaseous-hydrogen (GH2) streams. The model accounts for real-fluid thermodynamics and transport properties, making use of the PR-EoS and the Chung transport models. First, the development of a real-fluid model (RFM) with artificial neural network is described in detail. Then, the numerical results of the transcritical mixing layer (LOX/GH2) benchmark are presented and analyzed in terms of accuracy and computational ef-ficiency. Results of the overall implementation indicate that the combined OpenFOAM and ML approach provides a speed-up factor higher than seven, while preserving the original solver accuracy.

Unsmooth reversing is one of the most common faults in the four-way reversing valve of air conditioning system, and the airflow in the valve chamber and the shape of the slider of the reversing valve are the main factors in unsmooth reversing. In order to study the influence of the airflow in the chamber and the slider shape on the reversing process, the fluid flow in the valve chamber of a four-way reversing valve is obtained using computational fluid dynamics (CFD) simulations, and the influence of slider structure on slider thrust and resistance of the four-way reversing valve is analyzed. Four-way reversing valves with no-cutting, straight-cutting and arc-cutting sliders are employed to evaluate the reversing characteristics. A pressure measuring device for a four-way reversing valve chamber is designed, which is performed to test the chamber pressure of three different types of sliders of four-way reversing valves. The maximum error of experimental tests and simulations is within 5% of engineering tolerance, so as to verify the reliability of numerical simulation. The results indicate that the chamber pressure and resistance-thrust ratio with the cutting slider are reduced, which is beneficial to the reversing of the four-way reversing valve. Compared with the straight cutting structure, the arc cutting slider four-way reversing valve has a more stable reversing process and better comprehensive reversing performance. The results can provide a theoretical basis for engineers to improve the performance of four-way reversing valves.

This work presents a design and model of the thermo-hydrodynamic behavior applying CFD for pin-fin heat sink with square shape (rectangular fin profile), through a novel pin arrangement that improves heat removal in the system. The novelty consist in the insertion of a deflector row that will help to direct the flow, so it hits the rectangular cross-sectional fins the presents in the system. The heat sink model is placed in channel designed in which air flows as suitable medium, and its overall performance is investigated. In addition, a constant heat flux is applied to the bottom wall of the heat sink model, equal values and amounts that​ correspond to the heat fluxes generated by current electronic equipment and devices. The numerical results of the global thermal resistance, the pressure drop and the Nusselt number are reported. The results allow determining best arrangement or geometry configuration, location and functionality of the fin-deflectors were collocated for evaluate thermal-hydraulic efficiency of system.

Polymerase chain reaction (PCR) technique is one of the molecular methods in amplifying DNA for the detection of malaria. However, the limitations of PCR especially when used for routine clinical practice can hamper its sensitivity and specificity. With that, this study focuses on designing a microfluidic device that will mimic the function of a conventional genus-specific PCR based on 18S rRNA gene to detect malaria parasites (Plasmodium falciparum) at low grade parasitemia. The design was drawn and simulated using ANSYS 14.5 Computational Fluid Dynamics (CFD). The simulation shows that adding loops to the design increases its relative deviation but to a minimal extent as compared to having only a straight path design, which indicates that looping is acceptable in designing a microfluidic device to minimize chip length. Also, increasing the cross-sectional area of the fluid path decreases the efficiency of the design, thus, the design with a relatively smaller cross-sectional area is favored. And lastly, among the three materials utilized, the chip made of polypropylene is the most efficient with a relative deviation of 0.94 as compared to polycarbonate and polydimethylsiloxane which have relative deviations of 2.78 and 1.92, respectively.

Recanalization poses challenge in coil embolization for cerebral aneurysms. Establishing predictive models for post-embolization recanalization is important for clinical decision-making. However, conventional statistical and machine learning (ML) models may overlook critical parameters during the initial selection process. In this study, we automated the identification of significant hemodynamic parameters using a PointNet-based deep neural network (DNN), leveraging their inherent 3D features. Further feature analysis was conducted utilizing saliency mapping, an explainable artificial intelligence (XAI) technique. The study encompassed the analysis of velocity, pressure, and wall shear stress (WSS) in both pre- and post-coiling models derived from computational fluid dynamic (CFD) simulations for 59 aneurysms. Velocity was identified as the most pivotal parameter, supported by the lowest P-value from statistical analysis and the highest area under the curves (AUROC)/precision-recall curves (AUPRC) values from DNN model. Moreover, visual XAI analysis revealed robust injection flow zones with notable impingement points in pre-coiling models, as well as pronounced interplay between flow dynamics and the coiling plane were important 3D features in identifying the recanalized aneurysms. The combination of DNN and XAI was found to be an accurate and explainable approach not only at predicting post-embolization recanalization but also at discovering unknown features in the future.

This paper summaries the main findings from the full-scale CFD analyses conducted at SINTEF Ocean on the case of MV REGAL, which is one of the benchmark vessels studied in the ongoing joint industry project JoRes. The numerical approach is described in detail, and comparative results are presented regarding the propeller open water characteristics, ship towing resistance, and ship self-propulsion performance. The focus of numerical investigations is on the assessment of the existing simulation best practices applied to a ship scale case in a blind simulation exercise, and performance of thereof with different turbulence modelling methods. The results are compared directly with full-scale performance predictions based on model tests conducted at SINTEF Ocean and sea trails data obtained in the JoRes project.

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in e.g. the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. The paper consists of three distinct parts – an introduction to and review of the field of computational rhinology, a review of published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of the paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

Abstract：Objective: To compare the effects of fermentation on the chemical composition of Gastrodia Tuder Halimasch Powder, and to provide a basis for the processing, processing and clinical application of this medicinal material. Methods: UHPLC-Q-Orbitrap HRMS was used to identify compound structure by comparing excimer ions, characteristic fragment ions and reference substances, and referring to relevant literature and database search. The main differential components were screened by orthogonal partial least squares-discriminant analysis (OPLS-DA), and the differential components were quantitatively studied by high performance liquid chromatography, and the differences in the types and contents of chemical components before and after fermentation were compared. Results: A total of 77 compounds were identified or preliminarily derived after fermentation, including 21 terpenes, 18 organic acids, 9 flavonoids, 7 nucleosides, 6 amides, 4 amino acids, 2 pyrrolidone derivatives, 2 sterols, 2 alkaloids and 6 others.Through high performance liquid chromatography, it was found that the content of daidzein genistein and ergosterol was significantly increased after fermentation, and a new component ergosterol was produced. Compared with Gastrodia Tuder Halimasch Powder before fermentation, the content of genistin and other components was significantly reduced, and there may be a process of genistin conversion to genistein.Conclusion: The fermentation process will have a certain influence on the types and contents of the chemical components of Gastrodia Tuder Halimasch Powder, and ergosterol, genistein and .daidzein may be the main pharmacodynamic components, which can provide a useful reference for the clinical rational application, quality control and mechanism of action of Gastrodia Tuder Halimasch Powder.

(1) Background: This article presents a study that aims to provide a precise understanding of the temperature distribution within a Whole-Body Cryostimulation (WBC) chamber, whether it is empty or occupied by one or several individuals.; (2) Methods: The study employs a mixed numerical and experimental approach, utilizing simplified Computational Fluid Dynamics (CFD) simulations and an experimental analysis that employs thermocouples to determine the 3D thermal field in actual conditions.; (3) Results: The results reveal a non-negligible temperature difference between the setpoint and actual temperature in the middle of the cryochamber. Furthermore, it is shown that the presence of individuals inside the chamber results in both an average temperature rise and a more heterogeneous thermal behavior, which appear to be associated with the number of individuals present.; (4) Conclusions: The findings of this study emphasize the need for further research to establish temperature guidelines and standardize measurement methods for effective WBC treatment.

Computational Fluid Dynamics (CFD) is revolutionizing the oil and gas industry by offering powerful insights and innovative solutions. This paper explores the cutting-edge applications of CFD in enhancing and innovating within the sector. CFD simulations provide engineers with invaluable insights into fluid behavior, enabling optimized designs and informed decision-making. The software facilitates foresight by predicting system performance, optimizing drilling techniques, and maximizing reservoir management strategies. Moreover, CFD optimizes equipment design, increases operational efficiency, and reduces costs through virtual testing, eliminating the need for costly physical prototypes. The application of CFD aligns with the United Nations Sustainable Development Goals (SDGs), promoting energy efficiency, climate action, responsible production, and sustainable consumption. It addresses flow assurance challenges, enhances safety protocols, and mitigates environmental impacts. CFD is a powerful tool driving innovation, resilience, and sustainability within the oil and gas industry, ensuring a more efficient, safe, and environmentally conscious future. As advancements in CFD methodologies continue, the potential for further enhancements and breakthroughs within the sector becomes even more promising, driving the industry towards a sustainable and prosperous future.

Seepage erosion is one of the main reasons for the local collapse or instability of embankments. To investigate the characteristics and mechanism of seepage erosion for cohesionless soils, the model tests by using an independently developed seepage erosion device and the numerical simulations based on the Discrete Element Method-Computational Fluid Dynamics (DEM-CFD) coupling model were carried out. The results show that the seepage erosion process of the cohesionless soil can be characterized by four stages: stable seepage stage, fine particles upward migration stage, sand samples boiling stage, and erosion damage stage. The skeleton structure of soil sample under seepage flow is continually changed due to the loss of fine soil particles, which results in a significant decrease in the sample strength and may ultimately lead to the failure of the sample. The results of this study can provide references and bases for the design, construction, and long-term service of embankments or earth dams under the complex seepage condition, reducing the risk of seepage erosion.

Abstract: In this study, a computational fluid dynamics (CFD) simulation was conducted to validate the CFD model with experimental data from Zimparov et.al (2012), which investigates the heat transfer performance in a spirally corrugated tube that has a twisted tape inserted. The heat transfer was then compared to a simple corrugated tube without the twisted tape and to a smooth tube with no corrugations and no twisted tape. This simulation is compared with a previous experimental study conducted by (Zimparov et al. 2012), which focuses on tubes with a height-to-diameter ratio e/Di > 0.04 and small relative pitches of the twisted tape, H/Di. The largest improvement was noticed in the tube with e/Di= 0.057 and ridge pitch-to-height ratio p/e=6.77 with a twisted tape of H/Di= 4.7. This tube is labeled 5035. Tube 5035 was found to have the most significant enhancement, hence the focus of the CFD simulation will be on this tube, and the simulation will range for Reynolds numbers, 3.5 x 103 < Re < 5.0 x 104. The focus of this simulation is the evaluation of heat transfer and friction factors; a metric was used to obtain an empirical representation of the tube’s performance. Keywords: Spirally corrugated tube; Twisted tape; Heat transfer enhancement; Friction factor; Nusselt number; Computational fluid dynamics (CFD); Star CCM+

Tablets are the most used solid dosage form in medical treatments. The realistic modeling 1 of such specific object shapes in computer simulations has been demonstrated already by the author 2 for bi-convex round tablets but needs to be extended for a wider range of tablet shapes. Moreover, 3 the number of tablets must be large enough to simulate industrial drum coaters, which are typically 4 used for film coating processes. Whilst this long lasting process where tablets flow continuously, 5 tablets are loaded by frequent changes of individual contacts. Using a tailored simulation tool named 6 pharmAssist , written for GPUs (modern massively parallel processors in graphics cards and HPC 7 accelerators), we implemented the widely used oblong tablet shape. We are able to simulate flows 8 involving industrial scale number of such tablets. Data mining allows for statistical analyzis of 9 pair-wise contacts of different geometric features. Our simulation results serves as an observation 10 space, which cannot be addressed by experimental analyzis and will help to avoid degradation of 11 tablets like chipping or other defect phenomena. This paper reflects current development work in 12 order to extend the DEM with a new shape class, which comprises the modelling of compacted or 13 extruded objects, which have a smooth exterior, a flat band and sharp contour edges

Many farms utilize closed-type livestock systems to enhance productivity and facilitate effective environmental management. However, the confined nature of these closed spaces poses an increased risk of exposure to harmful gases and organic dust for both workers and livestock. Additionally, the introduction of outside air through ventilation systems can lead to temperature fluctuations within the breeding environment, resulting in potential productivity issues. This research paper employs computational fluid dynamics (CFD) to develop ventilation operation management plans that address both the working environment and the breeding environment simultaneously. The proposed plans are designed to be easily implemented in practical farm settings. The findings of this study, based on simulation analysis, indicate that while ventilation is effective in reducing harmful gases and improving the working environment, its efficiency decreases after the initial 3 minutes of operation. Furthermore, uncontrolled ventilation can cause sudden temperature changes, which may adversely affect the well-being of the livestock. However, when upgraded ventilation structures are implemented, significant improvements in the working environment (an average of 27.3% improvement) can be achieved while maintaining temperature stability for the livestock. These results highlight the importance of referring to the provided ventilation operation management table before commencing work, as it enables workers to improve the working environment while minimizing the potential impact of ventilation on the breeding environment.

OpenFOAM is a CFD software widely used in both industry and academia. The exaFOAM project aims at enhancing the HPC scalability of OpenFOAM, while identifying its current bottlenecks and proposing ways to overcome them. For the assessment of the software components and the code profiling during the code development, lightweight but significant benchmarks should be used. The answer was to develop microbenchmarks, with a small memory footprint and short runtime. The name microbenchmark does not mean that they have been prepared to be the smallest possible test cases, as they have been developed to fit in a compute node, which usually has dozens of compute cores. The microbenchmarks cover a broad band of applications: incompressible and compressible flow, combustion, viscoelastic flow and adjoint optimisation. All benchmarks are part of the OpenFOAM HPC Technical Committee repository and are fully accessible. The performance using HPC systems with Intel and AMD processors (x86_64 architecture) and Arm processors (aarch64 architecture) have been benchmarked. For the workloads in this study, the AMD processor seems particularly suited resulting in an overall shorter time-to-solution.

This comprehensive study focused on the standard taxi sign used in Ireland, its impact on drag production, fuel expenses, and CO2 emissions. Experimental analysis revealed that the conventional taxi sign significantly increased drag, especially when mounted on streamlined vehicles such as saloon cars, due to flow separation issues on the rear roof and rear windshield. Longitudinal reorientation of the sign offered a 14-fold reduction in drag increase compared to the traditional placement. It was found that positioning the sign in the middle of the roof offered the greatest fuel efficiency. Furthermore, the study estimated that implementing longitudinal repositioning on all Irish taxi signs could save drivers approximately 832 per year and reduce national CO2 emissions by a substantial 22,464 tonnes annually. Comparative analyses with international taxi signs demonstrated that the Irish sign had significantly larger drag contributions, emphasizing the need for improved aerodynamics. To address the inherent drag issue, the study explored novel appendable devices and proposed alternative taxi sign designs. Among the tested solutions, a magnet-mounted front ramp proved the most effective, reducing total drag by nearly 30%. Additionally, a motorized flip-up taxi sign design demonstrated a remarkable 40% reduction in drag. Finally, a newly proposed taxi sign design, featuring longitudinal positioning and pointed triangular front and rear faces, resulted in a minimal 4.3% increase in vehicle drag compared to the baseline car.

Abstract：Ballast pick-up velocity is a crucial parameter in the design of a slurry system in excavating devices, it influences transportation stability, safety, and economy of the system. In this study, a vertical pipe hydraulic experiment with different ballast sizes under water condition was conducted to investigate the ballast pick-up velocity in the shield slurry system. Thereafter, the coupling method of computational fluid dynamics and the discrete element method (CFD-DEM) was employed to establish a numerical model of the ballast pick-up velocity in the shield slurry system, and its applicability was verified by comparing the experiment results. Additionally, the ballast pick-up process under actual working conditions was studied, and the influence of ballast size and concentration on pick-up velocity was investigated. Overall, the primary findings of the study are as follows: 1) The numerical method can be applied to predict both the single ballast and ballast group pick-up velocity. 2) The pick-up process of the ballast group can be divided into four stages. 3) The conveying velocity of the ballast should be higher than the pick-up velocity of the single ballast to facilitate smooth discharge. 4) For a certain increase in ballast concentration, the pick-up velocity tends to be stable.

In recent years discrete element method (DEM) has gradually been applied to the traditional fluidization simulation of fine particles in micro fluidized bed (MFB). The application of DEM in simulating fast fluidization of fine particles in MFB has not yet received attention. This article presents a drag model that relies on the surrounding environment of particles, namely particle circumstance-dependent drag modle or PCDD model. Fast fluidization in MFB of fine particles is simulated using DEM based on the PCDD model. It is found that the gas-solid two-phase flow in MFB meets the common law of fast fluidization in circulating fluidized bed (CFB) to a certain extent, and also presents a special law that is significantly different from that in CFB. Simulations indicate that the local structure in MFB exhibits particle aggregation which is a natural property of fast fluidized, forming a structure where continuous dilute phase and dispersed concentrated phase coexist. The formation and fragmentation of particle clusters in different local regions of the bed have time synchronization. There exists strong effect of solid back-mixing in MFB, leading to relatively low outlet solid flux. The gas back-mixing effect is, however, not so distinct. The axial porosity shows a monotonically increasing distribution with the bed height, but does not strictly follow the single exponential distribution. The solid volume fraction at the bottom of the bed is significantly lower than the correlated value in CFB. The axial heterogeneous distribution of the cross-sectional average porosity in the lower half of the bed is also weakened. The radial porosity shows a distribution pattern of higher in the central region and lower in the sidewall region. Compared with the correlation results in CFB, the porosity near the central region in MFB is relatively low, while the porosity near the sidewall region is relatively high. The geometric size of the container, especially the radial size, is much smaller than that of a CFB. The relative area of contact between the MFB wall and particles is larger, and the wall friction factor becomes more significant. This may be the main reason for the time synchronization of particle agglomeration, severe particle back-mixing, low outlet solid flux, monotonically increasing axial porosity and weak core-annular structure of radial porosity.

The hepatic sinusoids are lined by thin endothelium containing transcellular pores dubbed fenestra. These fenestrations are completely open channels connecting the sinusoidal lumen to the underlying Space of Disse (SoD) and the hepatocytes of the liver parenchyma. The fenestrations are in the size range 0.05-0.35µm in diameter and cover between 5-15% of the sinusoidal endothelial surface area depending on their location along the sinusoids. Motivation/Research Question: So far, the narrowness of the sinusoids prevented the direct measurement of hemodynamic features such as sinusoidal pressure and flow velocity. A better understanding of these parameters might help in better understanding the physiology of the hepatic niche and possible implications in liver diseases. Few simulations of liver blood flow focus on the level of the individual sinusoid, fewer still include the transcellular pores (fenestrations) of the sinusoidal endothelium, and none have included i) a porosity gradient along the sinusoid wall modelled with through-all pores rather than a porus medium, ii) the presence of the Space of Disse, iii) lymphatic drainage. Model: Computed fluid dynamic (CFD) simulations performed using a numerical model with relevant anatomical characteristics (length, diameter, porous size and porosity). Boundary conditions included physiological/pathological values of inlet/outlet pressure and lymphatic outflow from the portal region of the SoD. Results: The pressure and flow velocity of the sinusoidal lumen was entirely dependent on the shape, i.e. constant versus divergent radius and not dependent on porosity in the sinusoidal wall. The velocity through the space of Disse (SoD) was affected by the addition of lymphatic drainage and increases in porosity in the pericentral area of the model. Variations in porosity also affected flow velocity through the fenestrations. Conclusions: The flow velocity in the SoD was modified by differences in porosity, while the flow velocity in the lumen of the sinusoids was unaffected, even by the absence of fenestrations. The overall shape of the vessel is the singular most important factor in the pressure flow behavior of the sinusoidal lumen. The flow rate over the hepatocytes is the flow rate through the SoD, and is modestly affected by the distribution of porosity along the sinusoid, and by the addition of a lymphatic drainage, this parameter would be of interest for modelling blood exchange with the hepatic parenchyma, especially if accounting for zonation.

Flooding, exacerbated by climate change, poses a significant threat to certain areas, increasing in frequency and severity. In response, the construction of supplementary dams has emerged as a reliable solution for flood management. This study employs a Remote Sensing (RS) approach integrated with GIS to assess the feasibility of constructing a supplementary dam near Linville, Brisbane, Australia, with the aim of mit-igating floods and preventing overtopping failure at the Wivenhoe Dam. Using QGIS software and a 25-meter resolution DEM from the Queensland Spatial Catalogue ‘QSpatial’ website, four potential dam sites were analysed, considering cross-sections, watershed characteristics, and water volume calculations. Systematic selection criteria were applied to identify the most optimal option based on dam wall di-mensions, volume-to-area, and volume-to-cost ratios. The selected option was further assessed against predefined criteria, yielding one optimal choice. The study provides insights into the feasibility and ef-fectiveness of supplementary dam construction for flood mitigation in the region, with recommendations for future research and implementation.

The present state of research in Computational Fluid Dynamics (CFD) is marked by an ongoing process of refining numerical methods and algorithms with the goal of achieving accurate modeling and analysis of fluid flow and heat transfer phenomena. Remarkable progress has been achieved in the domains of turbulence modeling, parallel computing, and mesh generation, resulting in heightened simulation precision when it comes to capturing complex flow behaviors. Nevertheless CFD faces a significant challenge due to the time and expertise needed for meticulous simulation setup and intricate numerical techniques. To surmount this challenge, we introduce paint2sim—an innovative mobile application designed to enable on-the-fly 2D fluid simulations using a device’s camera. Seamlessly integrated with OpenLB, a high-performance Lattice Boltzmann-based library, paint2sim offers accurate simulations. The application leverages the capabilities of the Lattice Boltzmann Method (LBM) to model fluid behaviors accurately. Through a symbiotic interaction with the open-source OpenCV library, paint2sim can scan and extract hand-drawn simulation domains, affording the capability for instant simulation and visualization. Notably, paint2sim can also be regarded as a digital twin, facilitating just-in-time representation and analysis of 2D fluid systems. The implications of this technology extend significantly to both fluid dynamics education and industrial applications, effectively lowering barriers and rendering fluid simulations more accessible. Encouragingly, the outcomes of simulations conducted with paint2sim showcase promising qualitative and quantitative results. Overall, paint2sim offers a groundbreaking approach to mobile 2D fluid simulations, providing users with just-in-time visualization and accurate results, while simultaneously serving as a digital twin for fluid systems.

Efficient implementation of clean energy technologies, including pioneering solutions like mobile solar PV systems mounted on trailers (MSPT), is pivotal for addressing energy challenges, especially in regions with constrained power grid access. This endeavour is vital for meeting escalating elec-tricity demands and aligning with the UN Sustainable Development Goal (SDG) aimed at ensuring dependable and sustainable energy provision in developing countries. This paper investigates the aerodynamic behaviour of a designed MSPT using numerical and experimental approaches. Spe-cifically, the study focuses on the effects of wind velocity and tilt angels on the drag and lift forces, as well as drag and lift coefficients on the panel used in the MSPT system. Additionally, the overall wind force acting on the entire MSPT (including 9 large solar PV panels) was examined, considering the combined effects of wind flow and the system's geometry. The simulations were conducted using ANSYS-Fluent software and experimental testing was performed utilizing scaled-down models to validate the accuracy of the simulation. The insights and the knowledge from this study can be utilized to optimize the design, installation, and operation of such systems, ultimately enhancing their performance, reliability, and paving ways for further improvements in mobile solar PV technology.

Labyrinth seals are commonly used in turbomachinery in order to control leakage flows. Flutter is one of the most dangerous potential issues for them, leading to High Cycle Fatigue (HCF) life considerations or even mechanical failure. This phenomenon depends on the interaction between aerodynamics and structural dynamics; mainly due to the very high uncertainties regarding the details of the fluid flow through the component, it is very hard to predict accurately. In 2014, as part of the E-Break research project funded by the European Union (EU), an experimental campaign regarding the flutter behaviour of labyrinth seals was conducted at Centro de Tecnologias Aeronauticas (CTA). During this campaign, three realistic seals were tested at different rotational speeds, and the pressure ratio where the flutter onset appeared was determined. The test was reproduced using a linearized uncoupled structural-fluid methodology of analysis based on Computational Fluid Dynamics (CFD) simulations, with results only in moderate agreement with experimental data. A procedure to adjust the CFD simulations to the steady flow measurements was developed. Once this method was applied, the matching between flutter predictions and the measured data improved, but still some discrepancies could be found. Finally, a set of simulations to retain the influence of the external cavities was run, what further improved the agreement with testing data.

This study presents an extensive collection of data on the aerodynamic behavior at a low Reynolds number and geometric coefficients for 2900 airfoils obtained through the class shape transformation (CST) method. By employing a verified OpenFOAM-based CFD simulation framework, lift and drag coefficients were determined at a Reynolds number of 10^5. Considering the limited availability of data on low Reynolds number airfoils, this dataset is invaluable for a wide range of applications, including unmanned aerial vehicles (UAVs) and wind turbines. Additionally, the study offers a method for automating CFD simulations that could be applied to obtain aerodynamic coefficients at higher Reynolds numbers. The breadth of this dataset also supports the enhancement and creation of machine learning (ML) models, further advancing research into the aerodynamics of airfoils and lifting surfaces.

The application of numerical simulation methods for studying vertical-axis wind turbines offers significant advantages over conventional horizontal-axis wind turbines regarding wind energy utilisation, start-up wind speed, and environmental impact. This paper aims to investigate and establish the aerodynamic design performance of asymmetrical (NACA 4412) airfoils used in the application of vertical-axis wind turbines. The computational fluid dynamics (CFD) simulation uses Ansys software based on the Reynolds Averaged Navier-Stokes (RANS) equations. The aerodynamic parameters (pressure and velocity) are evaluated and compared within a range of angles of attack (AoA) 5-12 degrees at a fixed Reynolds number, as is a usual condition for wind turbine application. The lift coefficient (Cl), drag coefficient (Cd), and Cl/Cd ratio are obtained within the range AoA for the asymmetrical airfoils, and the performances have been predicted using the potential flow panel and CFD methods. This study suggests that while Xfoil may not be as precise as computational fluid dynamics (CFD), its predicted lift and drag coefficients generally align with experimental results. These help us to understand the best air-foil selection in a wind turbine application and predict its optimum AoA at which it performs at its best.

This paper presents the second stage of a tandem fixed wing unmanned aerial vehicle (UAV) aerodynamic development. In the initial stage, the UAV was optimized analyzing its characteristics only in symmetrical flight conditions. Posted requirements were that both wings should produce relevant positive lift, the initial stall must occur on the front wing first, the center of pressure should be close to the center of gravity, while longitudinal static stability should be in the optimum range. Computational fluid dynamics (CFD) analyses were performed, where applied calculation model was derived from the authors’ previous successful projects. The eight version TW V8 has satisfied all longitudinal requirements. Lateral-directional CFD analyses of V8 showed that the ratio of the lateral and directional stability at the nominal cruising regime was optimal, but both lateral and directional static stabilities were too high. On further development versions the lower vertical tail was eliminated, negative dihedral was implemented on the front wing, and four inverted blended winglets were added. Version TW V14 has largely improved lateral and directional stability characteristics, while their optimum ratio at cruising regime was preserved. Longitudinal characteristics were also well preserved. Maximum lift coefficient and lift-to-drag ratio were increased, compared to the V8.

A micro energy harvesting device proposed in the literature, is numerically studied. It consists of two bluff bodies in a micro-channel and a flexible diaphragm at its upper wall. Vortex shedding behind bodies induces pressure fluctuation causing vibration of the diaphragm that converts me-chanical energy to electrical by means of a piezoelectric membrane. Research on enhancing vortex shedding is justified due to the low power output of the device. Amplitude and frequency of un-steady pressure fluctuation on the diaphragm center are numerically predicted. Vortex shedding severity is mainly assessed in terms of pressure amplitude. The CFD model set-up is described in detail and appropriate metrics to assess energy harvesting potential are defined. Several cases are simulated to study the effect of inlet Reynolds number and channel blockage ratio on the prospec-tive performance of the device. Furthermore, the critical blockage ratio leading to vortex shedding suppression is sought. Higher inlet velocity for a constant blockage ratio is found to enhance vor-tex shedding and pressure drop. Great blockage ratio values, but lower than the critical one, seem to provide great pressure amplitudes in the expense of moderate pressure drop. There is evidence that the field is fruitful for research and relevant directions are provided.

Over the last decades, an important number of regional and global Digital Elevation Models (DEMs) have been released publicly. As a consequence, researchers need to choose between several of these models to perform their studies and to use these DEMs as third-party data to compute derived products (e.g., for orthorectification). However, the comparison of DEMs is not trivial. For most quantitative comparisons, DEMs require to be expressed in the same Coordinate Reference System (CRS), sampled over the same grid (i.e., be at the same ground sampling distance with the same pixel-is-area or pixel-is-point convention), and with heights relative to the same Vertical Reference System (VRS). Thankfully, many open tools allow us to perform these transformations precisely and easily. Despite these rigorous transformations, local or global planimetric displacements may still be observed from one DEM to another. These displacements, or disparities, may lead to significant biases in comparisons of DEM elevations or derived products such as slope, aspect or curvature. Therefore, the control of DEM planimetric accuracy is certainly the most important task to perform at first before any comparison. This paper presents the disparity analysis method enhanced to achieve a sub-pixel accuracy, by interpolating the linear regression coefficients computed within an exploration window. This new method is significantly faster than oversampling the input data because it uses the correlation coefficients that have been already computed in the disparity analysis. To demonstrate the robustness of this algorithm, artificial displacements have been performed in a 11x11 grid with a 0.1-pixel step in both directions. The Copernicus DEM GLO-30 has been interpolated by a “classical” bicubic to achieve these 11x11 displacements. This validation method has been applied in four 10 km x 10 km DEMIX tiles showing different height distributions. Globally, this new sub-pixel accuracy method is robust. Artificial displacements have been retrieved with with typical errors (e_b) ranging from 12 to 20 % of the pixel size (worst case in Croatia).These errors in displacement retrievals are not constant in the 11x11 grid and the overall error E_b depends on the roughness (altimetry variations) encountered in the different tiles. The second idea of this paper is to assess the impact of the bicubic parameter (slope of the weight function at a distance d=1 of the interpolated point) on the accuracy of the displacement retrieval. By considering E_b as a quality indicator, tests have been performed in the four DEMIX tiles making the bicubic parameter vary between -1.5 and 0.0 by step of 0.1. For each DEMIX tile, the best bicubic (BBC) parameter b* is interpolated from the four E_b minimal values. This BBC parameter b* is low for flat areas (around -0.95) and higher in mountainous areas (around -0.75). The roughness indicator is the standard deviation of the slope norms computed from all the pixels of a tile. A logarithmic regression analysis performed between the roughness indicator and the BBC parameter b* computed in 67 DEMIX tiles shows a high correlation (r=0.717). The logarithmic regression formula estimating the BBC parameter from the roughness indicator is generic and may be applied to estimate the displacements between two different DEMs. This formula may also be used to set up a future Adaptative Best BiCubic (ABBC) that will estimate the local roughness in a sliding window to compute a local BBC.

Even though inhalation dosimetry is determined by three factors (i.e., breathing, aerosols, and the respiratory tract), the first two categories have been more widely studied than the last. Both breathing and aerosols are quantitative variables that can be easily changed, while respiratory airway morphologies are difficult to reconstruct, modify, and quantify. Despite several methods are available for model reconstruction and modification, developing an anatomically accurate airway model and morphing it to various physiological conditions remains labor-intensive and technically challenging. The objective of this study is to explore the feasibility of using an adjoint-CFD model to understand airway shape effects on vapor deposition and control vapor flux into the lung. A mouth-throat model was used, with the shape of the mouth and tongue being automatically varied via adjoint morphing, and the vapor transport being simulated using ANSYS Fluent coupled with a wall absorption model. Two chemicals with varying adsorption rates, Acetaldehyde and Benzene, were considered, which exhibited large differences in dosimetry sensitivity to airway shapes. For both chemicals, large shape deformations were explored to find the maximal possible morphing; both the maximal and intermediate deformations were simulated for morphology parametric studies. Results show that changing the mouth-tongue shape can alter the oral filtration by 3.2% for Acetaldehyde and 0.27% for Benzene under a given inhalation condition. This study demonstrates that the hybrid adjoint-CFD approach can be a practical and efficient method to investigate morphology-associated variability in inhalation dosimetry of vapors and nanomedicines.

In this paper, dual-rotor counter-rotating (CR) configurations of vertical axis wind turbines (VAWT) are briefly inspected and divided into three types. This investigation was focused on one of these types – the CR-VAWT with co-axial rotors, in which two equal rotors are placed on the same shaft, displaced from each other along it and rotate in opposite directions. For this CR-VAWT with three-blade H-Darrieus rotors, the properties of the design in terms of aerodynamics, mechanical transmission and electric generator, as well as control system are analyzed. A new direct-driven dual-rotor (DR) permanent magnet synchronous generator (PMSG) was proposed, in which two built-in low-power PM electric machines have been added. They perform two functions – start-up and overclocking of the rotors to the angular velocity at which the lifting force of the blades is generated and stabilizing the CR-VAWT work as wind gusts act on the two rotors. Detailed in this paper is the evaluation of aerodynamic performance of the CR-VAWT via 3D computational fluid dynamics (CFD) simulations. The evaluation was conducted using the CONVERGE CFD software with the inclusion of the actuator line model for the rotor aerodynamics, which significantly reduces the computational effort. Obtained results show that both rotors, while they rotate in opposite directions, had a positive impact on each other: the power coefficients of the upper and lower rotors in the CR-VAWT increased by 5.5% and 13.3% respectively compared to the single-rotor VAWT at the optimal distance between the rotors of 0.3 from a rotor height.

The result of the multidisciplinary collaboration of researchers from different areas of knowledge to validate a solar radiation model is presented. The MAPsol is a 3D local-scale adaptive solar radiation model that allows to estimate direct, diffuse and reflected irradiance for clear sky conditions. The model includes the adaptation of the mesh to complex orography and albedo, and considers the shadows cast by the terrain and buildings. The surface mesh generation is based on surface refinement, smoothing and parameterisation techniques and allows the generation of high quality adapted meshes with a reasonable number of elements. Another key aspect of the paper is the generation of a high-resolution digital elevation model (DEM). This high-resolution DEM is constructed from LiDAR data and its resolution if two times more accurate than the publicly available DEMs. The validation process uses direct and global solar irradiance data obtained from pyranometers at the University of Salamanca located in an urban area affected by systematic shading from nearby buildings. This work provides an efficient protocol for studying solar resources, with particular emphasis on areas of complex orography and dense buildings where shadows can potentially make solar energy production facilities less efficient.

This study proposes a customized industrial nozzle for generating impulsive air jets. It is installed on an automatic machine where wine caps are stacked and shot consecutively into a rotating cylinder by an air jet. This process is very noisy; hence, this study aimed to investigate possible geometric variations of the nozzle that can reduce the emitted noise. First, the nozzle was tested in a laboratory to measure the air consumption at different supply pressures. Subsequently, 3D models of the nozzle and its variations were created and used for computational fluid dynamics simulations. Different boundary conditions were set, first to validate the model and compare it with the experimental test results, and then to simulate the real working conditions and determine the geometry that is less noisy while maintaining the velocity peak. Among the various possibilities, shortening the final ducts of the nozzle appears to be the most promising solution. These modified nozzles could be easily added to current machines to provide immediate benefits, and this study represents a promising start for action on other machines where this type of device is present.

To establish simulated micro-gravity conditions, biological organisms are placed on a simulated microgravity system, which primarily rotate in two or three dimensions and rely on slow sedimentation with latency to detect gravity. Fluid dynamics arise as a consequence of rotation of the chambers with in which the cells are filled. To minimize the effects of centrifugal forces, it is necessary to position the samples as close to the center of rotation as possible. Strong shear forces and velocities within the vessel can have an adverse impact on cells, making this a critical topic in cell biology. During experimental investigation, this specific phenomenon is not apparent, but numerical modelling readily illustrates and is not widely analyzed in literature's. Furthermore, bulk flask are frequently employed on random positioning machines (RPM), not merely because they demand a substantial quantity of reagents. The flask presented in the study can be easily developed and fitted in conventional cell culture plates containing medium and reagents similar to those used in ground experiments, thus allowing it to remain close to the center of rotation. In the present study, a series of fluid dynamics simulations are performed to examine the impact of shear forces within the microvessel presented under varying conditions using a foreign medium condition, and these results are compared with existing vessels by observing the effects of shear forces and maximum velocity zones.

This paper investigates a numerical code coupling technique to tackle multi-physics and multi-scale simulations using state-of-the-art software packages that typically address some specific modeling domain. The coupling considers the in-house FEM code, FEMuS, and the FVM code OpenFOAM, by exploiting the MED library from the SALOME platform. The present approach is tested on a buoyancy-driven fluid flow within a square cavity, where the buoyancy force constitutes the coupling term. In uncoupled scenarios, momentum and temperature equations are solved both in FEM and FVM codes without data exchange. In the coupled setting, only the OpenFOAM velocity and the FEMuS temperature fields are solved separately and shared at each time step (or vice versa). The MED library handles the coupling with the addition of suitable ad-hoc data structures that perform the field transfer between codes. Different Rayleigh numbers are investigated, comparing the outcomes of coupled and uncoupled cases with reference literature results. Additionally, a boundary data transfer application is presented to extend the capabilities of the coupling algorithm to coupled applications with separate domains. In this problem, the two domains share interfaces and boundary values on specific fields as fluxes are exchanged between the two numerical codes.

This study explored the feasibility of fluoride removal from simulated semiconductor industry wastewater and its recovery as calcium fluoride using fluidized bed crystallization. The continuous reactor showed best performance (>90% fluoride removal and >95% crystallization efficiency) at a calcium to fluoride ratio of 0.6 within first 40 days of continuous operation. The resulting particle size increased by more than double during this time along with a 36% increase in the seed bed height indicating deposition of CaF2 onto the silica seed. The SEM-EDX analysis showed the size and shape of the crystals formed along with high amount of Ca-F ions presence. The purity of the CaF2 crystals were determined to be 91.1% though ICP-OES analysis. Following the continuous experiment, different process improvement strategies are explored. The addition of excess amount of calcium resulted in an additional 6% fluoride removal, however, compared to this single stage process, a two-stage approach was found to be better strategy to achieve a low effluent fluoride concentration. The fluoride removal reached 94% with this two-stage approach under optimum condition of 4 + 1 h HRT combinations and a [Ca2+]/[F-] ratio of 0.55 and 0.7 for the two reactors, respectively. CFD simulation showed the impact of inlet diameter, bottom-angle shape and width and height ratio of the reactor on mixing inside reactor and possibility of further improvement in reactor performance by optimizing the FBR configuration.

The production of anatase titanium dioxide particles plays a crucial role in the sulfate process used for manufacturing white pigment. A key factor in improving the quality of white pigments is enhancing the thermal mixing process within the precipitation tank. This improvement ensures the uniform dispersion of seed particles instead of their agglomeration, leading to the formation of particles with uniform sizes. The objective of this study was to enhance three-phase CFD simulations involving the mixing process of H₂SO₄ solution, steam as a gas phase, and solid seed particles. By analyzing the trajectories of the seed particles using CFD, the optimal injection position for the seed particles within the mixing process was determined. Subsequently, lab scale test and real field test were conducted based on the insights gained from the CFD simulations. The particle size distribution of two different types of seed inlets was analyzed and compared using Transmission Electron Microscopy (TEM). The findings of this study demonstrate that the developed multi-phase CFD simulation can be effectively utilized to enhance the precipitation process for the production of anatase titanium dioxide particles.

The increasing demand for higher power in rocket launchers is driving the development of higher power nozzles, achieved essentially by increasing the expansion ratio. However, this can lead to flow separation and resulting unsteady, asymmetrical forces, called side loads, which can limit the life of both the nozzle itself and other engine components. If this problem can be overcome, engine performance can be significantly increased. Therefore, various separation control methods have been proposed, but none have been fully implemented to date due to the uncertainties associated with modeling and predicting the flow phenomena. A numerical study of high area ratio rocket engine is presented to analyze the aeroelastic performance of the structure under the condition of flow separation. The finite volume method, the 2nd order upwind scheme and the extended wall function are used to analyze the flow inside the rocket nozzle and the transition of the Free Shock Separation (FSS) to Restricted Shock Separation (RSS) separation pattern ) to be discussed in detail. Since the location of the flow separation point depends heavily on the flow.

The performance of OpenFOAM in strong scaling tests on HPC systems with AMD EPYC processors exhibits a pronounced superlinear speedup. Simple test cases show superlinear speedups of over 300%, which significantly impacts the efficient use of computing resources. With the last generation of HPC architectures, a superlinear speedup of about 10% to 20% was well expected and accepted by CFD users. The measured superlinear speedup is much more pronounced and withstands the communication overhead to even larger scales. A detailed performance analysis of OpenFOAM follows, employing various High-Performance Computing (HPC) architectures, including AMD, ARM and Intel systems. The performance metric FVOPS (Finite VOlumes solved Per Second) is introduced to compare the performance of Computational Fluid Dynamics (CFD) applications when varying the grid size, as occurs in a strong scaling test. The achievable FVOPS depends on various factors, including the simulation type, boundary conditions, and especially the grid size of a use case. Analysing FVOPS on a single node level with varying grid size shows a significant difference in performance and cache utilization, which explains the large superlinear speedups seen in the strong scaling tests. Furthermore, FVOPS can be used as a simple benchmark to determine the optimal number of grid elements per rank to simulate a given use case at peak efficiency on a given platform, resulting in time, energy, and cost savings. The FVOPS metric also facilitates the direct comparison between different HPC architectures. The tests using AMD, ARM, and Intel processors show a peak in performance when employing around 10,000 grid elements per core. The presence of a large L3 cache on AMD processors is particularly advantageous, as indicated by L3 cache miss rates observed on AMD EPYC processors. Our results suggest that future HPC architectures with larger caches and higher memory bandwidth would benefit the CFD community.

The cavitation behavior of cryogenic liquids in the field of cryogenic engineering is distinct, showcasing fundamental differences from the cavitation behavior observed in ambient temperature water. Cavitation occurrence imposes severe limitations on the operational stability of cryogenic media transportation equipment, leading to safety concerns and increased manufacturing costs. Investigating the spatiotemporal characteristics of cavitation in cryogenic fluids holds significant value for both academia and engineering. The complexity of cavitation in cryogenic fluids primarily manifests through the coupled effects of thermodynamics, vortices, and bubbles during the cavitation process. In order to further explore the coupling mechanism, numerical simulations employing the Delayed Detached Eddy Simulation (DDES) turbulence model and Sauer-Schnerr cavitation model are conducted to validate the accuracy of the numerical calculations using the Hord experiment hydrofoil as a benchmark. Subsequently, the unsteady cavitation flow of liquid nitrogen around the NACA 0015 hydrofoil is simulated by the same numerical calculation method. The results indicate that the simulated results best coincide with experimental results when the bubble number density in the Sauer-Schnerr cavitation model reaches 108. The upstream development of the re-entrant jet under the driving force of inverse pressure gradient is the fundamental reason for the detachment of the primary cavitation zone. Under the same inflow cavitation number, the thermodynamic effect significantly inhibits the generation of bubbles and changes the range of the inverse pressure gradient, thus affecting the separation behavior of the primary cavitation zone.

As the demand for fossil fuels continues to rise and their availability becomes increasingly limited, it is crucial to explore innovative approaches to optimize fuel combustion in existing systems. Among various engines, diesel engines are heavily reliant on fluid motion within the engine cylinder, which plays a critical role in fuel-air mixing. This abstract focuses on the potential of dynamic mathematical model simulation for multi-compartment rotor compressed combustion engines, aiming to revolutionize power generation by enhancing fuel burning efficiency. By leveraging advanced simulation techniques, this research seeks to unlock new insights into fluid dynamics and optimize fuel-air mixing, leading to more efficient and sustainable energy utilization. Through this exploration, we aim to pave the way for groundbreaking advancements in the field of combustion engines, ultimately contributing to a greener and more sustainable future. In the realm of internal combustion engines, the fluid dynamics within the engine cylinder greatly influence the combustion processes and heat transfer. For compression ignition (C.I) engines, both the bulk gas motion and turbulence characteristics play vital roles in engine performance. This study focuses on enhancing the combustibility of the fuel-air mixture in C.I engines by modifying the engine design to induce turbulence through squish and tumble flows. The modification involves creating multiple compartments on the rotor crown, comprising three small chambers spaced 120º apart. By implementing this design, the aim is to optimize the fuel-air mixing and improve the overall combustion efficiency. Through comprehensive analysis and experimentation, we seek to revolutionize C.I engine technology and pave the way for more efficient and environmentally-friendly power generation. In this study, the performance of a compression ignition (C.I) engine with a multiple compartment rotor has been analyzed under motoring conditions using computational fluid dynamics (CFD) simulations with FLUENT software. The obtained results have been compared with those of a base C.I engine. The analysis focused on the effects of the modified engine design on the tumble ratio and squish velocity, which are important parameters for fuel-air mixing and combustion efficiency. The results indicated a significant improvement with the modified engine, with a 35% increase in tumble ratio and a 31% increase in squish velocity compared to the base engine. These findings suggest that the modified engine design has the potential to enhance fuel-air mixing and combustion performance, leading to improved overall engine efficiency. This research presents promising insights into the optimization of C.I engine design and contributes to the development of more efficient and environmentally-friendly power generation technologies..

In Active Flow Control (AFC) applications, to activate the boundary layer the use of pulsating flow has notorious energy advantages versus constant blowing/suction jet injections. For a given AFC application, five parameters, jet location and width, inclination angle, frequency of injection and momentum coefficient need to be tuned. Two main devices are nowadays capable of injecting pulsating flow with the sufficient momentum coefficient to delay the boundary layer separation, these are the Zero Net Mass Flow actuators (ZNMFA) and the Fluidic Oscillators (FO). In the present research, a novel FO configuration is analyzed for the first time at relatively high Reynolds numbers, fluid is considered as incompressible. After obtaining the typical linear correlation between the incoming Reynolds number and the outlet flow oscillating frequency, the effect of outlet width and mixing chamber wedge inclination angle dimensional modifications is addressed, modifications of the outlet width are observed to create large variations on the FO performance. The origin of the self-sustained oscillations is also analyzed in the present manuscript and greatly helps in clarifying the forces acting on the jet inside the mixing chamber. In fact we can conclude saying that the current FO configuration is pressure driven although the mass flow forces appear to be a much more relevant than in previous FO configurations studied.

This article presents a study on the aerodynamic drag of a generic dual-axle flatbed trailer and explores ways to reduce the drag using appendable drag-reducing devices. The primary sources of drag originated from the van and trailer’s rear, along with the trailer’s wheels. The most effective initial device for reducing drag was a full trailer underside cover which offered a 7% drag reduction. Additionally, ladder racks, dropsides, and rear gates were studied, and it was found that protruding ladder racks significantly increased drag. Rear gates added large amounts of drag and should be removed and stored when not needed. The study also explored novel mid-section devices that increased the van’s base pressure and reduced drag. An axle test revealed that drag for single, dual, and triple axle trailers was very similar in direct flow, but different in yawed flow. A drawbar length test showed a near-linear relationship between drawbar length and drag, manifesting as a 1.7% change in drag per 250mm change in drawbar length. Several novel modifications were made to the trailer, including fitting six unique appendable devices, which offered a total 7.3% drag reduction. A novel rear van device known as the Multi-Stage Converging Cavity was introduced which reduced drag by nearly 18%. When all the devices were used together, a total 25% drag reduction was observed for the van-trailer combination.

The discrete element method (DEM) coupled with computational dynamics (CFD) has been considered one of the most sensitive ways of studying fluidized beds. This paper presents gas-solid fluidized bed simulations by use of CFD-DEM. The two-phase interaction is calculated in a way of decoupling. Within the effective neighborhood of the particles, the kernel approximation method is used to calculate the local porosity, thereby calculating the drag force on the particles. At the grid scale, the momentum exchange coefficient is calculated using a two fluid model based method, thereby calculating the source term of the phase interaction in the fluid control equation. The simulated form, size and motion process of the big bubble are all consistent with experimental observations. The simulated solid volume fraction and relative pressure averaged in a horizontal plane at 45 mm above the bottom are in reasonable agreement with the experimental data. The fluctuation of the simulated bed height over time is similar to the experimental measurement results. The reproduction of these qualitative and quantitative results indicates that the proposed method effectively improves simulation performance and accuracy.

Typhoons bring great damages to transmission line systems located in coastal areas. Strong wind and extreme precipitation are the main sources of damaging effects. Transmission lines suffered from wind-driven rain exhibit more susceptibility to damage due to the coupled effect of wind and rain water. This paper presents an integrated numerical simulation framework based on mesoscale WRF model, multiphase CFD model and FEM model to analyze the motions of a transmission line subjected to coupled wind and rain loads during typhoon events. A full-scale transmission line in Zhoushan island is employed to demonstrate the effectiveness of the proposed framework by simulating typhoon evolution in terms of wind fields and rainfall, solving the coupled wind and rain fields around the conductor and predicting the dynamic responses of the transmission line during Super Typhoon Lekima (2019). The results show that the horizontal displacements of the transmission line under the joint actions of wind and rain increases approximately 17%~18% compared to those of wind loads only. It is important to consider the coupled effects of wind-driven rain on conductors in the design of transmission lines under typhoon conditions.

Vertical-axis wind turbines (VAWTs) are receiving more and more attention as they involve simple design, cope better with turbulence and are insensitive to wind direction which has a huge impact on the cost since the yaw mechanism is not needed. However, VAWTs still suffer from low conversion efficiency. As a results tremendous efforts are exerted to improve its efficiency which mainly focus on two methods regardless of whether the study is a CFD simulation or a field test or lab test experiment. An active approach which involves the modification of the rotor itself such as the blade design, angle, trailing and leading edges, inner blades, chord thickness, contra-rotating rotor, etc. While the second approach involves passive techniques where the flow is directed to optimally face the downwind rotor by mounting a guiding vane such as a diffuser or other shapes at the upwind position of the rotor. Among all techniques undertaken, the CRWT rotor technique seems to be the most effective with an output even comparable to that of Horizontal Axis Wind Turbines (HAWT), while the Savonius rotor received more attention compared to other designs of VAWTs. Apart from technological issues, it was also suggested that geographical issues such as proper site siting of the wind turbine rotor at a particular location where a uniform flow could be guaranteed is of paramount importance to ensure an effective conversion capacity of wind turbines.

The coseismic surface three-dimensional deformation fields were significant for quantified the geometric and kinematic characteristics of earthquake rupture faults. Traditional geodetic techniques were constrained by their intrinsic limitations, such as InSAR, which could only extract far-field deformation fields due to its incoherence, and GNSS, which could only acquire deformations at discrete points. The recently developing optical pixel correlation technique based on high resolution remote sensing images could acquire the near-field coseismic horizontal deformation. In this study, the vertical displacement field of the 2022 Mw 6.6 Menyuan earthquake was obtained by differencing the pre- and post-earthquake DEMs from GaoFen-7 stereo satellite images. Further, the line-of-sight (LOS) and azimuthal (AZI) directions far-field deformations by InSAR, the horizontal near-field deformation by optical pixel correlation based on pre- and post-earthquake GF-2/7 images, and the vertical deformation by differencing the pre- and post-earthquake DEMs, were incorporated to comprehensively solve the complete three-dimensional deformation fields of the 2022 Mw 6.6 Menyuan earthquake. Firstly, the vertical displacement field by differencing DEMs indicated there were significant vertical displacements about 1 m at the bend region, which was induced by the local compressive stress. In addition, at the epicenter on the middle segment of ruptured Lenglongling fault (LLLF), the maximum lifting occurred on the southern sides of the main and secondary faults exceeded 2 m. Secondly, the three-dimensional deformation fields solved by multiple deformation data demonstrated that the near-field deformation field calculated by optical pixel correlation method could quantify the displacements distributed over the rupture fault zone, which were not available on the InSAR deformation maps. Finally, the surface 2D strain derived from the displacement maps calculated by optical pixel correlation revealed high strain concentration on the rupture fault zone. Our study focused on the complete surface three-dimensional deformation estimation from multiple far- and near-field deformation data, and provided a new perspective for a deeper understanding of the characteristics of coseismic surface deformation and the rupture pattern of the fault.

With the continuous increase in ship size combined with the generally slower increase in the size of waterways, the need for the prediction of ship maneuvering in shallow waterways continues to attract attention from the international scientific community. Ship behavior in shallow water is relevant in seabed effects that result in changing the hydrodynamic forces acting on the ship. In this study, the maneuvering characteristics of container ship with four degrees of freedom in shallow water are analyzed. The Reynolds-averaged Navier Stokes approach in Ansys Fluent code is used to produce the maneuvering coefficients through the simulations of forward running, static drift, static heel, circular motion, the combined motions, and the pure roll motion of the Kriso container ship. The maneuvering characteristics of the ship are estimated for evaluating the ship behaviors in shallow water condition. The obtained results show that the roll has a significant decrease and the ship’s turning diameter has a significant increase when the ship operates in the shallow waterway. The prediction of maneuvering characteristics of the ship are in good agreement with those of free running model tests, indicating that the numerical simulation based on the Computational Fluid Dynamics method has compromising capability to predict the maneuvering derivatives and the 4-DOF ship maneuvering motion in shallow water as well.

The article presents a multiphase numerical computational fluid dynamics (CFD) for simulating hydrothermal liquefaction of sewage sludge in a continuous plug-flow reactor. The discrete particle method (DPM) was used to analyze the solid particle interaction in liquid-solid high shear flows to investigated coupling computational fluid dynamics (CFD) and discrete phase model (DPM) method. Increasing solid particles interaction were observed with increasing liquid velocity. The study examined the influence of parameters such as flow rate, temperature and residence time on the efficiency of bio-oil production. Furthermore, the trajectories of the solid particles are affected significantly by the ash flow and the solid particle size, thus the trajectories of particles resembles the streamlines. A parametric study is performed to investigate the effect of slurry flow rate, temperature and external heat transfer coefficient on the biocrude oil productions performance. This simulation data will be used in the future to design and scale up a large-scale HTL reactor.

Heat transfer enhancement and hydrodynamic characteristics of nanofluid use in a micro heat exchanger is investigated for cooling electronic components working in hot climatic conditions. The cooling fluid employed was water and TiO2 nanoparticles at concentrations of 1% and 5%, Reynolds numbers ranged from 400 to 2000, and inlet temperatures ranged between 35ºC and 65ºC. At a nanofluid inlet temperature of 55ºC and a nanoparticle concentration of 1%, the Nusselt number increased by 23% up to 54% as the Reynolds number varied between 400 and 2000. At a nanoparticle concentration of 5%, the percentages which correspondingly enhanced the Nusselt number were 32% and 63%. The temperature of the electronic heating component decreased by 4.6-5.2ºC when the nanofluid concentration was increased from 0 to 5% at a Reynolds number of 400 and a nanofluid inlet temperature of 35ºC. Small increments in pressure drop of about 6% and 13% were observed at nanofluid concentrations of 1% and 5%, respectively. With nanoparticle concentrations of 1% and 5%, a Reynolds number of 2000 and a nanofluid inlet temperature of 35ºC, performance evaluation criterion (PEC) values of 1.36 and 1.45 were obtained. When the nanofluid inlet temperature was increased to 65ºC, the PEC parameter decreased to 1.02-1.10 for both concentrations.

Timber-concrete composite (TCC) systems join the positive aspects of engineered wood products (good seismic behaviour, low thermal conductivity, environmental sustainability, good behaviour under fire if appropriately designed) with those of concrete (high thermal inertia, durability, excellent fire resistance). TCC facades are typically composed by an internal insulated timber-frame wall and an external concrete slab, separated by a ventilated air cavity. However, there is very limited knowledge concerning the performance of TCC facades, especially for what concerns their thermal behaviour. The present paper deals with the development and optimisation of a 2D CFD model for the analysis of TCC ventilated façades thermal behaviour. The model is calibrated and validated against the experimental data collected during the annual monitoring of a real TCC ventilated envelope in the north of Italy. Also, a new solver algorithm is developed to significantly speed up the simulation. The final model can be used for the time-efficient analysis and optimisation of the thermal performance of TCC ventilated facades, as well as other ventilated facades with external massive cladding, avoiding the expensive and time-consuming construction of mock-ups, or the use of comparably slow (conventional) CFD solvers that are less suitable for optimization studies.