Computational Fluid Dynamics Method Research Articles

3-D full-field reconstruction of chemically reacting flow towards high-dimension conditions through machine learning

Monitoring chemically reacting flow is crucial for optimizing and controlling the conversion processes in various chemical engineering scenarios. These processes are often influenced by multiple factors, creating a high-dimensional condition space that imposes intensive computational cost to computational fluid dynamics (CFD) methods. We present a machine learning (ML) framework to reconstruct the 3-D scalar field preserving the accuracy and resolution of CFD modeling. We demonstrate the framework’s effectiveness through a case study on biomass combustion in an industrial-scale grate furnace, which exemplifies the challenges of modeling chemically reacting flow in complex geometric domains influenced non-linearly by multiple parameters. The framework implements three core principles: (1) a Reactor-Structure-Resembled Network (RSRNet) reflecting the operational logics of the reactor, (2) an incremental learning strategy based on random Latin Hypercube Sampling, and (3) the incorporation of self-extracted high-level features from data as soft constraints in the loss function. Ground truth data were calculated through a CFD model, with all cases sampled from a 11-variable condition space. Using only 240 randomly sampled cases (around 1 % total cases in the discretized condition space), RSRNet precisely replicated 1-D and 2-D scalar distributions through incremental training. Further reconstruction of 3-D temperature field only used 40 3-D CFD cases, resulting in the training and testing errors decreased to 7.94 × 10-4 and 1.55 × 10-3, respectively, approaching the level of the 2-D reconstruction, with the average error in temperature field predictions not exceeding 50 K. The generalization ability was validated using a separate test dataset with continuously changing excess air ratio and capacity ratio in wide ranges, including some extreme values, and the model successfully captured complex phenomena under unseen conditions.

Influence of the stratification parameters on the wake field of an underwater vehicle in a two-layer stratified flow

Because the ocean exhibits different stratification states in different seasons and regions, the effect of different stratification parameters on a submarine wake must be studied. In particular, to investigate such effect in a two-layer stratified flow, this study established a computational fluid dynamics (CFD) model based on the URANS equations, shear-stress transport (SST) k-ω turbulence model, and Eulerian multiphase flow model. The volume of fluid (VOF) method was adopted to implement the free-surface capture method. The working conditions of the SUBOFF model under an infinite diving depth and near the free surface were compared with those of the experiments, and the convergence of the time step and grid number under stratified flow conditions was verified. These validation steps proved the feasibility of the CFD method. Finally, the effects of the stratification parameters on the free surface waves, internal waves, and turbulent kinetic energy behind the vehicle were observed by varying the stratification parameters to achieve different working conditions. The results show that when an underwater vehicle sails in a two-layer stratified flow, the interstratified density ratio and depth of the internal wave surface affect its wake field. When the underwater vehicle sails in the water layer, with an increase in the interstratified density ratio, the free-surface roughness increases and the turbulent kinetic energy decreases. When the underwater vehicle sails in the salt water layer, the result is just the opposite.

Savonius wind turbine blade selection based on computational fluid dynamics

The Savonius vertical axis wind turbine is widely used because of its simple design and efficiency at low wind speeds. The next development is the mini multi-blade Savonius helical wind turbine. At this miniature size, the performance of a mini multi-blade helical Savonius turbine due to variations in the number of blades and its impact on the energy produced has not been investigated in depth. Therefore, adjustments to the size and number of blades need to be made to optimize the performance of this turbine. The research aims to obtain the performance of the Savonius type S wind turbine with 3 blades and 4 blades using the CFD (Computational Fluid Dynamics) method. The main dimensions of the simulation model are the shaft diameter of 0.008 m, the outer diameter of 0.24 m, the height of 0.4 m, and the number of blades 3 and 4. The inner diameter of the 3-blade turbine is 0.21 m while the inner diameter of the 4-blade turbine is also 0.21 m. Both turbine models are subjected to wind speeds of 4 m/s, 5 m/s, and 6 m/s. Based on the simulation results, a turbine with 4 blades has more optimal performance compared to the 3-blade type. At a wind speed of 6 m/s, the 4-blade Savonius helical turbine produces a rotation of 498,215 rpm and a maximum power of 6,129 watts. Meanwhile, the Savonius turbine with 3 blades at a wind speed of 6 m/s produces a rotation of 545,655 rpm and a maximum power of 4,390 watts. This difference in performance shows the importance of adjusting the number of blades in wind turbine design to achieve optimal efficiency. This research provides useful insights for the further development of mini helical Savonius wind turbines in various wind conditions

Study on oil film spreading characteristics on the jet lubricating tooth surface of aviation herringbone gear

Purpose During the oil lubrication process of aviation herringbone gears, lubricating oil is injected into and spreads around the tooth surface. Its spreading characteristics directly affect the lubrication state and transmission efficiency of gears. This study aims to investigate the effect of changing the injection parameters and gear operating conditions on the oil film deposition during injection lubrication. Design/methodology/approach The computational fluid dynamics method was used to establish the computational domain for the simulation and optimization using an orthogonal test. The oil film spreading and thickness of the tooth surface were calculated with the change in gear speed, fuel injection distance and fuel injection speed. Optimized oil injection lubrication parameters were obtained, and design ideas were provided for the oil injection lubrication analysis of aviation herringbone gears. Findings Under different operating conditions, such as varying the injection speed, gear rotational speed and injection distance, there were differences in the distribution of the lubricating oil film formed on the tooth surface. The thickness of the oil film decreases as the distance from the injection port increases. The thickness of the oil film deposition can be increased by increasing the injection velocity and reducing the gear rotation speed. The injection distance had a relatively small effect on the spreading thickness of the deposited oil film, affecting only the dynamic spreading process of oil film deposition. Furthermore, the simulation results are compared with the analytical calculation results, and finally the experimental verification is carried out. Originality/value Oil spray lubrication characteristics are important for achieving precise gear lubrication designs in the aerospace field. Peer review The peer review history for this article is available at: https://publons.com/publon/10.1108/ILT-08-2023-0250/

Research and Analysis of the Small Disturbance Equation in Subsonic, Transonic, and Supersonic Regimes

This paper explores the application of the Small Disturbance Equation (SDE) across subsonic, transonic, and supersonic flow regimes. Derived from the Euler and Navier-Stokes equations, the SDE offers an efficient framework for analyzing aerodynamic behaviors, particularly through the utilization of discretization techniques and iterative solving methods executed in Python. The study assesses the accuracy and limitations of the SDE in detailing essential flow characteristics, revealing that while the equation performs effectively in subsonic and transonic flows, it encounters challenges in supersonic regimes. Nonlinear effects such as shock waves significantly hinder its performance at high speeds. Compared with conventional computational fluid dynamics (CFD) methods, the SDE stands out in scenarios where computational efficiency is paramount. However, its limitations in handling high-speed flows must be carefully considered, highlighting the need for further refinement in its application to supersonic dynamics. This analysis suggests that while the SDE is beneficial for certain aerodynamic studies, its scope and utility are constrained by the inherent complexities of high-speed fluid dynamics.

Numerical Simulation of Resistance and Flow Field for Submarines near Ice Surface

When a submarine operates in polar regions, the polar environment inevitably impacts its resistance and flow field characteristics, especially when the submarine navigates near the ice surface. This paper investigates the hydrodynamic characteristics of a submarine sailing near the free water surface and the ice surface using computational fluid dynamics (CFD) methods. In order to quantify the impact of ice on the resistance and flow field characteristics of the submarine, the resistance coefficients are calculated for both near ice surface and free surface. The resistance, velocity field, and pressure distribution around the submarine at different depths and speeds are analyzed. The results indicate that the total resistance of the submarine sailing near the ice surface is lower than the free water surface. When the submarine is sailing near the ice surface, its total resistance coefficient decreases with increased submergence depth at a constant Froude number. At a fixed depth, the resistance coefficient also decreases as the Froude number increases. Additionally, when the dimensionless depth relative to the maximum hull diameter (D) exceeds 3.5, it has little effect on the resistance coefficient.

Engineering and Modeling of Water Flow via Computational Fluid Dynamics (CFD) and Modern Hydraulic Analysis Methods

The manuscripts presented in this Special Issue will both present engineering analyses of problems of contemporary importance in fluid mechanics and hydrodynamics and describe historical water systems and the technologies used to design and operate them by ancient Old and New World water engineers [...]

Hemodynamics simulation and analysis of left coronary artery aneurysms with concomitant stenosis

The hemodynamic parameters in arteries are difficult to measure non-invasively, and the analysis and prediction of hemodynamic parameters based on computational fluid dynamics (CFD) has become one of the important research hotspots in biomechanics. This article establishes 15 idealized left coronary artery bifurcation models with concomitant stenosis and aneurysm lesions, and uses CFD method to numerically simulate them, exploring the effects of left anterior descending branch (LAD) stenosis rate and curvature radius on the hemodynamics inside the aneurysm. This study compared models with different stenosis rates and curvature radii and found that as the stenosis rate increased, the oscillatory shear index (OSI) and relative residence time (RRT) showed a trend of increase; In addition, the decrease in curvature radius led to an increase in the degree of vascular curvature and an increased risk of vascular aneurysm rupture. Among them, when the stenosis rate was less than 60%, the impact of stenosis rate on aneurysm rupture was greater, and when the stenosis rate was greater than 60%, the impact of curvature radius was more significant. Based on the research results of this article, it can be concluded that by comprehensively considering the effects of stenosis rate and curvature radius on hemodynamic parameters, the risk of aneurysm rupture can be analyzed and predicted. This article uses CFD methods to deeply explore the effects of stenosis rate and curvature radius on the hemodynamics of aneurysms, providing new theoretical basis and prediction methods for the assessment of aneurysm rupture risk, which has important academic value and practical guidance significance.

Research on flame characteristics and coating growth process of titanium alloy in HVAF spraying

AbstractIn this work, a flow field model of WC–12Co powder sprayed by high‐velocity air fuel (HVAF) was established based on the computational fluid dynamics (CFD) method. The macro–micro coupled thermodynamic model of coating multilayer growth process was established based on the birth and death element method. The temperature and velocity of the particles in the combustion flame were applied to the coating growth model, and the transient evolution of the coating temperature and thermal stress was revealed. The results show the temperature and thermal stress of coatings increase with the increasing the number of spraying layers, and their incremental gradient gradually decreases with the increasing the number of deposition layers. The peak temperature of the coating surface reaches 1494 K, and the peak thermal stress of the first coating reaches 2693 MPa. The maximum stress of spraying particles appears at the contact edge of particles. Further preparation of WC coating, through the hardness, friction, and wear, corrosion tests verify that WC–12Co coating can effectively provide the TC18 substrate performance. The microhardness of WC–12Co coating was approximately three times that of TC18 substrate. The average friction coefficient and corrosion rate of the coating are 0.15 and 0.06 mg/cm2/h lower than that of the substrate, respectively.

Numerical analysis of coupling effect between sloshing and motion responses of a KCS in uni- and bi-directional waves

The seakeeping behaviour can be severely affected by the coupling between sloshing and ship motions in waves. The sloshing behaviour of liquid-loaded ships in bi-directional waves differs significantly from that in uni-directional waves. Given that ships operate in short-crested or multi-directional waves during their whole service period, it is significant to gain a deeper understanding of the effect of wave directionality on the coupling between sloshing and ship motions. In this study, a KRISO container ship (KCS) model with type-B liquefied natural gas (LNG) cargo tanks is adopted for coupling analysis. The 2D long-crested regular wave is adopted as an uni-directional wave, and three dimensional (3-D) cross wave is adopted as a bi-directional wave. The computational fluid dynamics (CFD) method is applied for the sloshing and seakeeping simulation of the ship model in different wave fields. The differences in wave elevation, sloshing forces, and coupling effects between ship motions and sloshing are studied in these two wave fields by adjusting control parameters such as ship heading angle and filling ratio. The comparative results indicate that sloshing behaviour and coupled motion responses in cross waves should be paid enough attention during ship design and evaluation.

Damage analysis and assessment of concrete T-girder bridge based on fire scene numerical reconstruction

Fire is a sporadic disaster of concrete bridges, with diverse fire environments and complex damage mechanisms. Accurately evaluating the damage situation of concrete bridges after a fire is exceedingly challenging. This study formulates a damage analysis and assessment method based on the step-by-step and progressively deepening working principle. The method relies on fire scene numerical reconstruction and encompasses key technical aspects, including bridge detection and analysis during the fire incident, fire scene numerical reconstruction, and subsequent bridge damage assessment. Building upon these principles, the study utilizes results from the detection and analysis of the concrete T-girder bridge during a fire incident to establish Computational Fluid Dynamics (CFD) and Finite Element Method (FEM) models for the numerical reconstruction of the fire scene. These models enable the calculation of varying temperature distributions and the evolution of the bridge under fire. Compared with the parameters obtained through the ISO834 method, the numerical reconstruction approaches not only enhances the accuracy of replicating the bridge combustion process but also enables the extraction of temperature field distribution patterns within the bridge fire space and its concrete components.

CFD Analysis of Takeoff from a Water Surface for an Insect-Scale Aerial/Aquatic Robot

To develop an insect-scale aerial/aquatic robot, we analyzed takeoff mechanisms to counteract surface tension, such as paddling, slapping, and clap-and-fling. Because a diving beetle, Eretes griseus, takes off directly from the water surface, a flapping-wing robot is promising as an alternative to a drone with multiple rotary wings. In this study, we first investigated diving beetle flight with a three-dimensional high-speed camera system and analyzed the motion characteristics. Subsequently, we developed a computational fluid dynamics method that tracked the water surface using a volume of fluid method, reproduced the motion with a multibody model, treated the deformation of the elastic membrane wing with the phase delay of the joint angle functions, and simulated takeoff, that is, the transition from water to air, and hovering near the water surface. The simulation result showed that during the transition, the slapping motion exerted the maximum and average lift per unit of body weight of 18 and 9.2, respectively, while those of paddling produced 0.46 and 0.23, respectively. The water surface effect improved the lift by 25% at the normalized height of less than 0.44 and disappeared at a height greater than 0.7. During hovering, while the clap-and-fling motion improved lift by 2.6% and the water surface effect was 9.8%, the synergy effect was 22%. In addition, the former enhanced it significantly after the fling, while the latter was remarkable during the wing acceleration phase. In contrast to ground effects, flapping reduced the water level and caused the ripples, dynamically changing the water surface effect.

Effect of different PV modules surface energies and surface types on the particle deposition characteristics and PV efficiency

The particle deposition behavior of solar photovoltaic (PV) modules and its effect on PV efficiency were numerically investigated using computational fluid dynamics (CFD) methods. The surface flow field of PV modules using the shear stress transfer (SST) k-ω model with a user-defined function (UDF) for the development of entrance boundaries is predicted. A UDF coupled discrete phase model (DPM) that takes into account wall bounce and hydrophobic properties is used to predict particulate deposition in the flow field. Mesh independence and average pressure coefficients are verified. The effects of particle diameter, boundary type, surface energy and surface type on particle deposition of PV panel and PV efficiency are investigated. The results show that the smaller the surface energy γ, the smaller the particle deposition rate. At dp = 150 μm, γ = 0.001J/m2, it showed that the most significant reduction of 427.3 % in particle deposition rate compared to the trap boundary type, and γ = 0.001J/m2 showed a 210 % reduction in particle deposition rate compared to γ = 0.1J/m2. However, at dp ≤ 50 μm, the change in surface energy has little effect on the deposition rate. The effect of convex PV modules on particle deposition is more obvious at low surface energy (γ ≤ 0.001J/m2) compared to planar. The reduction in PV efficiency caused by particle deposition at different surface energies is predicted using an empirical model of PV output reductions developed by the researcher.

Numerical Study of the Combustion Process in the Vertical Heating Flue of Air Staging Coke Oven

To investigate the combustion process and reduce Nitric Oxide (NO) emissions in the vertical heating flue of air-staged coke ovens, a three-dimensional computational fluid dynamics method was applied to simulate the combustion process. The model integrates the k-ε turbulence model with a multi-component transport combustion model. The impact of air staging on the flow field and NO emissions in the vertical fire chamber was assessed through comparative validation with experimental data. The impact of air staging on the flow field and NO emissions in the vertical fire chamber was assessed through comparative validation with experimental data. Based on this research, the effects of the excess air coefficient and air inlet distribution ratio on NO emission levels at the flue gas outlet were further investigated. Analysis of the flow field structure, temperature at the center cross-section, component concentration, and NO emission levels indicates that as the excess air coefficient increases, the NO emission levels at the flue gas outlet initially decrease and then increase, accompanied by corresponding changes in outlet temperature. At an air excess factor of 1.3 and an air inlet distribution ratio of 7:3, NO emission levels are at their lowest—53% lower than those in a conventional coke oven—and the temperature distribution in the riser channel is more uniform. These results provide a theoretical foundation for designing the air-staged coke oven standing fire channel structure.

Hydrodynamic Performance Improvement of a Tirhandil Yacht by Stern Form Modifications

This paper presents a comprehensive investigation to improve the hydrodynamic performance of a Tirhandil hull form by modification efforts on the stern region. The form improvement approach combines computational fluid dynamics (CFD) methods with computer-aided design (CAD) systems. The design process for the reference and modified models was carried out by using CAD systems. The hydrodynamic characteristics of the reference hull form were evaluated by employing CFD methods and it was determined that form improvements should be concentrated on the stern region. The modification process was conducted by considering constraints on the design variables in the stern region and the main dimensions of the reference model. A grid independence study was performed to evaluate various grid structures to determine the optimal mesh configuration for the numerical analyses. The SST k-Omega turbulence model was used for the numerical analyses to simulate turbulence structure around the hull form. Achieving around a 13.4% reduction in the total resistance coefficient, the modified model also exhibited decreased wave amplitudes, smoother wave transitions, and a significant reduction or cancellation of shoulder and stern waves.

Numerical study of a chemical looping combustion system: Effects of fuel type and particle size distributions on the operating performance

Chemical looping combustion (CLC) using biomass offers a promising pathway for achieving negative carbon emissions. However, the different properties of coal and biomass affect reaction performance, which has not been fully explored in the CLC system. Additionally, the effects of particle size distributions (PSDs) of solid fuels and oxygen carriers (OCs) are inadequately addressed. In this study, a 3D reactive model using the computational particle fluid dynamics (CPFD) method is developed to investigate the effects of PSDs of both solid fuels and oxygen carriers on operating performance with coal and biomass as fuels. Results indicate that the average residence time in FR of coal and biomass particles is 0.86 s and 0.68 s, respectively, but biomass outperforms coal in reaction efficiency. For coal-fueled CLC, the average residence time is 0.46 s for coal particles with diameters of 100–200 μm, whereas it reaches 0.86 s for coal diameters of 300–700 μm. However, the reaction performance is better with smaller coal sizes due to char gasification being the rate-limiting step. In comparison, for biomass-fueled CLC, the reaction efficiency remains consistent across various particle sizes under the specified conditions due to the rapid release of volatile components and the low content of char. Additionally, it is found that the OC conversion ratio is mainly determined by their residence time in the zones with relatively higher concentrations of the combustible gas in FR. Employing OCs with overly wide PSDs such as 300–800 μm can compromise the material seal above the AR, increase the risk of gas leakage and decrease the separation efficiency of the inertial separator.

CFD-DEM simulation and experimental investigations on continuous hydraulic sampling of deep-sea polymetallic nodules

The physical sampling to grasp the distribution of deep-sea polymetallic nodules is essential for efficient mining. A hydraulic sampling scheme is proposed that balances low environmental disturbance, high sampling precision, broad coverage, and statistical data viability to address the limitations of existing methods. This paper utilizes a coupled approach of Computational Fluid Dynamics (CFD) and Discrete Element Method (DEM) to explore innovative solutions for precise and continuous hydraulic sampling of deep-sea polymetallic nodules. Focusing on the optimization of the double-row nozzle sampling mechanism design, simulations are conducted to analyze the dynamic characteristics of solid-liquid two-phase flow during the hydraulic sampling process under various nozzle configurations and operational water pressures of the pump. Laboratory experiments with the hydraulic sampling mechanism validate the accuracy of the simulation and the effectiveness of the optimization suggestions. The configuration is finalized with a front nozzle diameter of 10 × φ9 mm, a rear nozzle diameter of 10 × φ7 mm, and an initial working water pressure for the pump of 80 kPa. Subsequently, integrated driving-sampling experiments are conducted using a self-propelled tracked vehicle, validating the feasibility of the proposed configuration.

Advancing packaging technology: computational fluid dynamics modeling for capillary underfill encapsulant in multi-chip heterogenous packages

Purpose The purpose of this study is to investigate the dynamics of capillary underfill flow (CUF) in flip-chip packaging, particularly in a multi-chip configuration. The study aims to understand how various parameters, such as chip-to-chip spacing (S12), chip thickness (tc) and others, affect the underfill flow process. By using computational fluid dynamics (CFD) simulations and experimental studies, the goal is to provide insights into understanding the dynamics of CUF in heterogeneous electronic packaging. Design/methodology/approach The paper introduces a CFD analysis and experimental study on CUF in a multi-chip configuration, aiming to understand underfill flow dynamics. A 3D geometry models of multi-chip arrangement are created using computer-aided design (CAD) software. After that, the CAD models are meshed and simulated in Ansys Fluent using incompressible and non-Newtonian fluid properties. The study maintains S12 of 2.86 and tc of 22.29 between experimental and simulation data for results validation. Next, a various of S12 values (1.14, 2.86, 5.71, 8.57, 14.29 and 20) which focus on tc of 22.29 have been investigated. Further studies have been conduct using S12 of 5.71 and tc of 8.00, 14.29 and 22.29. Findings Results show a strong correlation between simulation and experiment which validate the correctness and robustness of simulation. Further parameter’s studies using simulation for various of S12 indicated that higher S12 values lead to faster flow. This effect is due to large underfill weight from reservoir able to flow into S12 region which contributed to higher mass momentum movement. Furthermore, the effect of various of tc shows that the thicker the chip the faster the underfill to flow in S12 region. Research limitations/implications The intentional exclusion of solder bump pattern arrangements from the experiment and simulation may limit the study's ability to fully understand the impact of solder bump patterns on underfill flow. Therefore, more parameters can be investigated such as solder bump pattern, underfill weight and dispense pattern in the future using CFD. Practical implications The manuscript provides a comprehensive examination of the contributions of CFD to the advancement of knowledge regarding CUF phenomena in heterogeneous electronic packaging assemblies. Moreover, it delineates the utilization of CFD methodologies to assess the influence of chip-to-chip spacing (S12) and the thickness of the chip (tc) on the underfill flow characteristics. Originality/value This paper fulfills an identified need of computational fluid dynamics method to study capillary underfill flow dynamics in heterogenous electronic packaging.

Numerical Analysis of Hydrogen Behavior Inside Hydrogen Storage Cylinders under Rapid Refueling Conditions Based on Different Shapes of Hydrogen Inlet Ports

In order to investigate the effects of different shapes of hydrogen inlet ports on the behavioral characteristics of hydrogen in Type IV hydrogen storage cylinders under rapid refueling conditions, a mathematical model of hydrogen temperature rise and a three-dimensional numerical analysis model were developed. The rectangular, hexagonal, triangular, Reuleaux triangular, circular, elliptical and conical inlet ports were researched by using computational fluid dynamics methods. The results showed that, for the same refueling flow rate and cross-sectional area, the hydrogen temperature inside a cylinder with a rectangular inlet port is higher and the jet tilt angle is larger than for a hexagonal port, while the thermal stratification phenomenon is not obvious. The hydrogen temperature inside a cylinder with a triangular inlet port is lower than that with a Reuleaux triangle port and the jet tilt angle is larger, and neither has significant thermal stratification. The hydrogen temperature inside a cylinder with a circular inlet port is higher than that with an ellipse port, the jets are not tilted on either one, and the phenomenon of thermal stratification is prominent. Further analysis indicated that enlarging the cross-sectional area and increasing the refueling flow rate results in a higher hydrogen temperature and intensified thermal stratification and an upward-angled jet can effectively reduce or eliminate thermal stratification.

Comparative Analysis of NURBS and Finite Element Method in Computational Fluid Dynamics Applications: Case Study on NACA 2412 Airfoil Aerodynamics

In this research, an attempt was made to employ the Non-Uniform Rational B-Splines (NURBS) method for a challenging computational fluid dynamics (CFD) problem of aerodynamics around NACA 2412 airfoils. The comparison was carried out thoroughly by using the same boundary conditions and geometry, comparing NURBS to standard FEM implementations. Our study was interested in demonstrating the foreseeable functionalities of NURBS for solving complex CFD problems and conducting a comparative effectiveness performance evaluation between them with traditional FEM methodologies.