The computation of compressible and incompressible recirculating flows by a non-iterative implicit scheme
Abstract
The PISO algorithm, which is presented in a companion paper, is a non-iterative method for solving the implicity discretised, time-dependent, fluid flow equations. The algorithm is here applied in conjunction with a finite-volume technique employing a backward temporal difference scheme to the computation of compressible and incompressible flow cases. The results of calculations are compared with similar ones obtained with an existing iterative method. It is shown that for time-evolving flows the splitting error of PISO is negligibly small at the level of time-step required to eliminate the temporal truncation error, while the avoidance of iteration results in a substantial reduction in computing effort over that required by iterative methods. It is also demonstrated that PISO is stable for fairly large time steps, which renders it useful for steady-state calculations as well.
References (15)
- R.I. Issa
J. Comput. Phys.
(1985) - F.H. Harlow et al.
Phys. Fluids
(1965) - F.H. Harlow et al.
J. Comput. Phys.
(1971) - A.D. Gosman et al.
- L.S. Caretto et al.
- S.V. Patankar
Numerical Heat Transfer and Fluid Flow
(1980) - W.R. Briley et al.
J. Comput. Phys.
(1977)
Cited by (459)
Numerical simulation and spray model development of liquid ammonia injection under diesel-engine conditions
2024, EnergyAmmonia is regarded as a promising alternative fuel in internal combustion engine for its potential to reduce carbon emissions. Injection process of liquid ammonia will affect air-fuel mixture distribution and therefore combustion efficiency and NOx emissions. However, the experimental studies of liquid ammonia injection under engine conditions are limited, and the key impact factors of spray evolution are still unclear. This study aims to evaluate the applicability of numerical models for liquid ammonia spray, and numerically investigate spray characteristics under diesel-engine conditions (ambient temperature > 800K, ambient pressure > 2 MPa). Large eddy simulation coupled with Lagrangian particle tracking method were conducted, and the results show that the non equilibrium - Frössling phase change model coupled with Schiller-Naumann drag model can provide a good prediction of spray penetration. Liquid penetration increases with ambient pressure decreasing; vapor penetration increases with ambient density decreasing, but insensitive to ambient temperature. Temperature - mixture fraction distribution is sensitive to ambient temperature and superheat degree. Besides, a three-stage 0D model of spray tip penetration under diesel-engine conditions is developed firstly, which can reflect the key impact factors of ammonia spray evolution and provide theoretical support to optimize fuel injection strategy of liquid ammonia-fueled engines.
Accelerating flow simulations in the built environment by using the fast fluid dynamics initializer
2024, Building and EnvironmentThis study introduces the Fast fluid dynamics (FFD) based initializer as a novel approach to accelerate steady-state flow simulations. Three modes - Coarse-Mesh Mode (CMM), Large Timestep Mode (LTM), and Hybrid Mode (HM) - are developed on the OpenFOAM platform for efficiently simulating airflow in the built environment. All three modes exhibit faster computations than standard initialization methods by enabling quicker attainment of desired flow states and reducing the required matrix iterations. These initializers demonstrate better performance in flow problems characterized by higher Reynolds numbers. The CMM initializer achieves a nearly six-fold increase in speed by mapping initial solutions from a coarse mesh to a fine mesh. The LTM initializer allows for larger timesteps during initialization, resulting in a two-fold speed improvement in fluid simulation. The HM initializer combines the advantages of both CMM and LTM. However, it does not surpass the speed-up performance of the previous two modes. Apart from the benchmark test case, the proposed initializers can be applied to other indoor and outdoor flow scenarios. Furthermore, the simulation process is further streamlined by automating the execution of the initializer using scripts, eliminating the need for manual step-by-step inputs.
Acoustic streaming and cavitation effect on microstructure evolution of Nb-Si-Ti alloy treated by ultrasound
2024, International Journal of Heat and Mass TransferHeat and mass transfer play a critical role in grain refinement by ultrasound treatment. In order to reveal the mechanism of ultrasonic in Nb-Si molten metal, a 3-D model was simulated on the acoustic, temperature and flow fields with different acoustic frequencies and displacement amplitudes by Comsol Multiphysics. Simulation results show that a new harmonic wave appears with frequency increasing, and the acoustic streaming effect is significant at 20 kHz. As the displacement amplitude increases from 10 μm to 20 μm, the maximum velocity varies from 0.01 m/s to 0.03 m/s. The distributions of three fields are dependent on frequency, and the displacement amplitude determines acoustic pressure value and fluid velocity. The displacement amplitude has a positive effect on the acoustic streaming effect improvement in the melt. In the UST ingots with different displacement amplitudes, the size of the broken dendrites decreases from 57.5 μm to 10 μm with the displacement amplitude lower than 16 μm, by 82.6%. γ-Nb5Si3 phase increases and the primary Nbss phases are elongated, coarsened and spheroidized with the increase of the displacement amplitudes. The simulation results of the flow directions can be verified by the elongated Nbss phase. The acoustic streaming effect promotes heat and mass transfer and the sluggish eutectoid reaction in Nb-Si-Ti alloy during the ultrasound treatment.
Pore-scale numerical analysis of fluid flows in compressed polyurethane foams with a workflow of open-cell foams modeling
2024, Journal of Fluids and StructureseX-Poro-HydroDynamic (XPHD) lubrication presents a different scientific approach to dealing with tribological problems. It is an innovative inter- and multidisciplinary research topic which offers a promising sliding solution for various applications, such as thrust bearings, various guide components and in terms of load capacity and damping. XPHD lubrication is a lubrication mechanism of biomimetic inspiration which features an additional parameter to the system “the porous media”. It consists of self-sustained fluid films generated within compressible porous layers imbibed with liquids in replacement for using the fluid film only as in the classic lubrication system. Soft and porous structures imbibed with liquids generate a high load support under compression. The load support is generated through the resistance to flow inside the porous material. During compression, the resistance to flow and load support increases the greater the compression rate. The main objective of this work is then to understand the behavior of the fluid flow inside the porous structures when subjected to axial compression stress. In the scientific literature, the works studying the flow in compressible materials are essentially experimental because of their very complex geometrical shape, the CFD (Computational Fluid Dynamics) simulations offer an economical solution to study the performance of this new concept of lubrication. To create the geometry, the morphological structures of foam samples are reconstructed at different levels of compression rates from 3D (Three Dimensional) X-ray microtomography. This is achieved by using the commercial software Avizo that allows to process 3D images and create 3D meshes suitable for numerical simulations. The numerical simulations of flows will be performed with the solver IcoFoam of the toolbox OpenFOAM for incompressible laminar flows, making it possible to study the pressure drop in these porous structures. The performed simulations were made with a polyurethane foam of 96% porosity using five compression rates for creating the different structures. The analysis of the numerical simulations shows the impact of the polyurethane foam compression on different key parameters such as the decrease in the permeability as function of the compression rate, the anisotropy of the flow within the compressible structure and the actual increase in the tortuosity generated by the compression of the foam and the variation of the porosity.
Oscillation of an ultrasonically driven gas bubble in an asymmetric confined domain
2024, International Journal of Mechanical SciencesThis study numerically investigates the behaviour of an ultrasonically driven gas bubble between two parallel rigid circular walls with a cylindrical micro-indentation in one wall using the OpenFOAM software. The primary objective is to determine the conditions that facilitate the removal of particulate contamination from the indentation using the bubble jet. To accomplish this objective, various values of the indentation diameter were considered, ranging from 0 to 23. It was found that, when , during the collapse phase, the bubble splits into two smaller sub-bubbles, within each of which a liquid jet emerges pointing towards the adjacent wall. Once the lower sub-bubble transforms into a toroidal shape, the high-speed liquid jet penetrates the indentation. This allows for the effective removal of particulate contamination present within the indentation specifically for cases where . On the other hand, the lower sub-bubble completely vanishes when , rendering the removal of contamination impossible. It is important to note that for large values of , there may exist situations in which numerical instability arises, resulting in the creation of artificial gas-filled regions within the liquid. We effectively resolve this issue through the proper treatment of the computational grid.
An implicit splitting scheme with characteristic boundary conditions for compressible reactive flows on unstructured grids
2024, Journal of Computational and Applied MathematicsThe computation of combustion instabilities requires efficient numerical solvers. For low Mach number compressible flows splitting schemes provide an adequate family of solution strategies. Here, the set of governing equations is separated into several subsystems, e.g. into an advective and acoustic subsystem. Advancements for such an algorithm are presented improving accuracy and stability. Moreover a new, gradient based approach is given for coupling the splitting scheme with characteristic boundary conditions on collocated unstructured grids. A large number of academic test cases is chosen to demonstrate the accuracy and performance of the proposed improved numerical scheme. In addition, the improved solver is tested in a complex, non-reactive LES of a swirl combustor on an unstructured grid.