Simulation Of Fluid Flow Research Articles

Numerical Analysis on the Aerodynamic Characteristics of SUV Car Model Install with Vortex Generator

The investigation of automotive aerodynamics involves analysing several forces operating on a car while driving on the road, such as drag and lift forces. The flow separation at the vehicle's rear end is one of the primary causes of aerodynamic drag for automobile vehicles. It is feasible to improve fuel efficiency by lowering the drag force. The study focuses on the influence of a vortex generator (VG) on the aerodynamics of a sport utility vehicle (SUV) car. The study aims to simulate fluid flow analysis for an SUV car that uses VG and without VG, as well as to assess the impact of a different configuration of VG and a varying number and fillet radius of VG. The number of VG are 3, 5, and 9. The different fillet radius of VG are 5, 10, and 15 mm. Using the Reynold-Averaged Navier Stokes Equation (RANS) in the numerical simulation, the Reynolds number at the computational domain is 1.1391 × 107 and 1.4808 × 107, which is determined by the height of the model and the freestream velocity. The results show that aerodynamic characteristics are significantly influenced by the number of VG and various size radius fillets of VG. From the result, 9 number of VG and 5 mm fillet radius obtained the lowest value of coefficient of drag, Cd compared with the others which is Cd = 0.3747 for 27.78 m/s and Cd = 0.5031 for 33.33 m/s, respectively. Furthermore, the analysis of flow structures suggested the locations of vortex formation and wake turbulence at the rear of the vehicle. In contrast, the 9 number of VG with a 5 mm radius fillet of VG emerged as the most suitable VG for this scenario, exhibiting a Cd value closest to the base model and the lowest Cd value.

Enhanced Thermoelectric Generator Power Output Based on Electrodeposited Bi2Te3 Nanocomposites with Pt NPs-CNTs Through Multiphysics Simulation

Thermoelectric generators (TEG) possess the potential to transform unused heat into electrical energy, making thermoelectricity a viable alternative for addressing the energy issue. This study provides insight into the enhanced power density of a TEG device when embedded with bismuth telluride (Bi2Te3) nanocomposite. The TEG power density has been analysed under the fluid flow simulation condition by utilizing several Bi2Te3 nanocomposite materials with Pt NPs-CNTs inclusion. The impact of static and steady air flow conditions has been studied using a detailed computational fluid dynamics (CFD) parameter, with flow velocities of 0.01 m/s and 1 m/s. As a main part of the work, the power density of nanocomposites embedded as N-type legs was compared to that of pure material. The Bi2Te3 with hybrid nanocomposite produces the highest power density, precisely combining Bi2Te3 with platinum nanoparticles (Pt Nps) and single-wall carbon nanotube (SWCNT) materials (Bi2Te3/Pt-SWCNTs). At velocity of 1 m/s, the power density was calculated to be 156.85μW/cm2, which increased to 158.86 μW/cm2 at a fluid velocity of 1 m/s, marking an 88% increment compared to the TEG model using pristine Bi2Te3 and higher than that of the previous work about 4 times of increment.

Fluid flow simulation with an [formula omitted]-accelerated Boundary-Domain Integral Method

The development of new numerical methods for fluid flow simulations is challenging but such tools may help to understand flow problems better. Here, the Boundary-Domain Integral Method is applied to simulate laminar fluid flow governed by a dimensionless velocity–vorticity formulation of the Navier–Stokes equation. The Reynolds number is chosen in all examples small enough to ensure laminar flow conditions. The false transient approach is utilised to improve stability.As all boundary element methods, the Boundary-Domain Integral Method has a quadratic complexity. Here, the H2-methodology is applied to obtain an almost linear complexity. This acceleration technique is not only applied to the boundary only part but more important to the domain related part of the formulation. The application of the H2-methodology does not allow to use the rigid body method to determine the singular integrals and the integral free term as done until now. It is shown how to apply the technique of Guigiani and Gigante to handle the strongly singular integrals in this application. Further, a parametric study shows the influence of the introduced approximation parameters. For this purpose the example of a lid driven cavity is utilised. The second example demonstrates the performance of the proposed method by simulating the Hagen–Poiseuille flow in a pipe. The third example considers the flow around a rigid cylinder to show the behaviour of the method for an unstructured grid. All examples show that the proposed method results in an almost linear complexity as the mathematical analysis promisses.

Matlab implementation of pore morphology method for modeling liquid residue in porous media with heterogeneous wettabilities

An alternative to the popular volume-of-fluid (VOF) and lattice Boltzmann (LB) methods for simulating fluid flow in porous media is the purely geometric pore morphology method (PMM). The PMM does not require solving partial differential equations and is therefore, orders of magnitude faster than the VOF and LB methods in terms of CPU time (e.g., about a week on 100 cores for the VOF vs. about 30 mins on one core for the PMM in our case). The PMM employs voxel-based morphological operations to provide a quasi-static approximation for capillary flow in porous structures. This paper presents step-by-step instructions for how to implement the PMM in the MATLAB environment, and further to benchmark the method using the simulation results obtained from the VOF and other continuum-based numerical methods. More specifically, this paper presents PMM simulations for the intrusion of a non-wetting fluid in two- and three-dimensional porous media with heterogeneous contact angles. Detailed discussions are also presented for the pros and cons of the PMM and how its non-continuum nature can be used as an advantage to conduct two-phase flow simulations in geometries where a symmetry boundary condition cannot be applied.

Simulation of Fluid Flow through Orifice Meter

The modeling of fluid flow using an orifice meter is an engineering endeavor aimed at thoroughly investigating and optimizing fluid flow measurement in various industrial applications. Orifice meters are highly regarded for their accurate measurement of fluid flow rates as they pass through a specifically designed orifice plate, utilizing the principles of fluid dynamics. The primary objective of this project is to develop a computational fluid dynamics (CFD) model capable of simulating the intricate fluid flow through an orifice meter. The simulation encompasses detailed geometric features of the orifice meter, including the precise dimensions of the orifice plate, pipe diameter, and any relevant taper angles. Additionally, the model takes into account various fluid parameters, boundary conditions, and factors such as velocity and pressure at different points within the system. The comprehensive study visualizes and quantifies key fluid properties, such as velocity profiles, pressure differentials, and flow velocities, across the entire length of the orifice meter using CFD analysis. This in-depth analysis provides crucial insights into the dynamic performance of the orifice meter. The obtained information holds significant potential to revolutionize flow rate measurement accuracy, efficiency, and cost-effectiveness in diverse industries, including water management, chemical processing, and oil and gas. The project's methodology for simulating fluid flow through orifice meters can bring about substantial improvements in the understanding and optimization of fluid flow in industrial processes.

Hybrid CFD PINN FSI Simulation in Coronary Artery Trees

This paper presents a novel hybrid approach that integrates computational fluid dynamics (CFD), physics-informed neural networks (PINN), and fluid–structure interaction (FSI) methods to simulate fluid flow in stenotic coronary artery trees and predict fractional flow reserve (FFR) in areas of stenosis. The primary objective is to utilize a 1D PINN model to accurately predict outlet flow conditions, effectively addressing the challenges of measuring or estimating these conditions within complex arterial networks. Validation against traditional CFD methods demonstrates strong accuracy while embedding physics-based training to ensure compliance with fundamental fluid dynamics principles. The findings indicate that the hybrid CFD PINN FSI method generates realistic outflow boundary conditions crucial for diagnosing stenosis, requiring minimal input data. By seamlessly integrating initial conditions established by the 1D PINN into FSI simulations, this approach enables precise assessments of blood flow dynamics and FFR values in stenotic regions. This innovative application of 1D PINN not only distinguishes this methodology from conventional data-driven models that rely heavily on extensive datasets but also highlights its potential to enhance our understanding of hemodynamics in pathological states. Ultimately, this research paves the way for significant advancements in non-invasive diagnostic techniques in cardiology, improving clinical decision making and patient outcomes.

The minimal flow unit and origin of two-dimensional elasto-inertial turbulence

The research on elasto-inertial turbulence (EIT), a new type of turbulent flow, has reached the stage of identifying the minimal flow unit (MFU). On this issue, direct numerical simulations of FENE-P fluid flow in two-dimensional channels with variable sizes are conducted in this study. We demonstrate with the increase of channel length that the simulated flow experiences several different flow patterns, and there exists an MFU for EIT to be self-sustained. At Weissenberg number ( $Wi$ ) higher than the one required to excite EIT, when the channel length is relatively small, a steady arrowhead regime (SAR) flow structure and a laminar-like friction coefficient is achieved. However, as the channel length increases, the flow can fully develop into EIT characterized with high flow drag. Close to the size of the MFU, the simulated flow behaves intermittently between the SAR state with low drag and EIT state with high drag. The flow falling back to ‘laminar flow’ is caused by the insufficient channel size below the MFU. Furthermore, we give the relationship between the value of the MFU and the effective $Wi$ , and explain its physical reasons. Moreover, the intermittent flow regime obtained based on the MFU gives us an opportunity to look into the origin and exciting process of EIT. Through capturing the onset process of EIT, we observed that EIT originates from the sheet-like extension structure located near the wall, which is maybe related to the wall mode rather than the centre mode. The fracture and regeneration of this sheet-like structure is the key mechanism for the self-sustaining of EIT.

Design and Performance of High-Capacity Magnesium–Air Battery for Power Generator System

Efforts to achieve carbon neutrality, which aims to reduce the net carbon emissions to zero by decreasing carbon emissions from human activities and increasing carbon absorption, are actively underway. Additionally, the search for clean energy alternatives to fossil fuels has become a global research trend. This paper presents research on metal–air batteries, focusing on the development of energy supply technologies that do not generate carbon emissions during power generation and require less space for power generation compared to existing renewable energy sources. The proposed Mg–air battery (MAB) in this study uses magnesium as the metal anode and theoretically offers a maximum open-circuit voltage of 3.1 V and a high energy density of 6.8 kWh/kg. While previous research has primarily focused on designing small-capacity cells and maximizing the performance of metal anodes, this study differentiates itself by designing a large-capacity MAB cell and optimizing its electrical performance. For the large-capacity cell design, the weight, shape, and size of the anode were designed based on MAB performance factors, and research was conducted on manufacturing methods to optimize the performance of the air cathode. Furthermore, to enhance usability and extend the lifespan of the MAB cell, it was designed to allow electrolyte circulation, and the electrolyte circulation performance was verified through simulations of fluid flow within the cell. Based on the study of the power performance of the newly designed large-capacity MAB cell, the feasibility of constructing a kW-class system using multiple Mg–air battery cell stacks was confirmed.

Optimizing production well geometry in the Utah FORGE geothermal project using machine learning and fluid flow modeling

This study addresses the critical challenge of optimizing well placement in Enhanced Geothermal Systems (EGS), specifically within the framework of the Utah FORGE geothermal project, as part of the 2023 Society of Petroleum Engineers (SPE) Geothermal Datathon. Effective well placement is essential for enhancing geothermal production efficiency and maximizing resource utilization. We employed a discrete fracture network (DFN) modeling approach, utilizing the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm from the Scikit-Learn library to analyze microseismic event location data. Through rigorous simulations conducted in the open-source GeoDT fluid flow simulator, we identified an optimal production well configuration characterized by a spacing of 400 m, an injection rate of 0.03 m³/s, and alignment parameters that significantly improve thermal recovery. The results indicate a projected net present value (NPV) of $75 million over a 20-year operational horizon, underscoring the economic potential of optimized well placement strategies. This study offers valuable insights for the operation of the FORGE geothermal site. More importantly, it exclusively utilizes open-source tools, enhancing accessibility and adaptability for the broader geothermal community.

A conceptual and numerical model of fluid flow and heat transport in the Topusko hydrothermal system

Comprehensive understanding of hydrothermal systems is often obtained through the integration of conceptual and numerical modelling. This integrated approach provides a structured framework for the reconstruction and quantification of fluid dynamics in the reservoir, thereby facilitating informed decision-making for sustainable utilisation and environmental protection of the hydrothermal system. In this study, an updated conceptual model of the Topusko hydrothermal system (THS; central Croatia) is proposed based on structural, geochemical, and hydrogeological analyses. The stratigraphic sequence and the structural framework of the THS were defined based on geological maps and field investigations. As depicted by hydrochemical and isotope analyses, thermal waters in Topusko (temperatures of up to 65 °C) are of meteoric origin and circulate in a carbonate aquifer. The THS receives diffuse recharge approximately 13 km S from Topusko, where Triassic carbonates crop out. Gravity-driven regional groundwater circulation is favoured by regional thrusts that tectonically uplifted Palaeozoic low permeable rocks. These structures confine the fluid flow in the permeable, fractured and karstified, Triassic carbonates favouring the northward circulation of the water. A regional anticline lifts the aquifer closer to the surface in Topusko. Open fractures in the anticline hinge zone increase the fracturing and permeability field of the aquifer promoting the quick upwelling of the thermal water resulting in the Topusko thermal springs. Numerical simulations of fluid flow and heat transport corroborate the proposed conceptual model. In particular, a thermal anomaly was modelled in the Topusko subsurface with temperature values of 31.3 °C and 59.5 °C at the surface and at the base of the thermal aquifer, respectively, approaching the field observations. These findings showed that the circulation of Topusko thermal water is influenced by regional and local geological structures suggesting that the enhanced permeability field in the discharge area enables the formation of the natural thermal springs.

Sensitivity analysis of a dual-continuum model system for integrated CO2 sequestration and geothermal extraction in a fractured reservoir

Depleted hydrocarbon reservoirs are considered as the most feasible option for CO2 geological sequestration and utilization. Most of the hydrocarbon reservoirs are naturally fractured. Simulation of fluid flow and heat transfer in these fractured formations remains a significant challenge in reservoir engineering. In this study, a dual-continuum model is developed to simulate integrated CO2 sequestration and CO2-circulated geothermal extraction in a fractured reservoir block in North Oman. The high-dimensional sensitivity of key parameters controlling CO2-brine flow and heat transfer in this matrix-fracture system is quantitatively evaluated by an efficient surrogate modeling approach. The surrogate models are constructed and validated based on a suite of physics-based model simulations. It is found that fracture permeability dominates the CO2 injectivity, storage, circulation and associated geothermal extraction. Response surface analysis shows that the flow area density between matrix-fracture and matrix block length controls the flux interaction between matrix and fracture formations. In contrast, the fracture aperture shows negligible influence in the dual-continuum modeling system. Particularly, sensitivity varying with locations on the response surface is analyzed for defined performance indicators.

Modelling of Fluid Permeability at the Interface of the Metal-to-Metal Sealing Surface.

This paper presents a method for modelling the permeability of fluid at the interface formed between flat parallel plates and the sharp-edged ridges of a metal gasket. This work was divided into three stages. In the first stage, numerical calculations simulating the deformation (compression of the gasket) were performed. The calculations were carried out using thermomechanical static analysis with commercial software. The purpose of these calculations was to determine the contact area of the gasket ridges with the plates, the deformation of the gasket ridges, and the reaction force resulting from the degree of compression of the gasket. In the second part of this work, analytical calculations were performed to estimate the tightness level. The analytical model proposed in this paper was based on Darcy's equation, simulating fluid flow through a ring-shaped porous layer. The analytical model also took into account the shape of the roughness profile of the sealed surfaces. A mathematical Ausloos-Berman function based on fractal theory was used to represent it. In the last part of this work, experimental tests were carried out to determine the actual fluid permeability and thus verify the numerical and analytical calculations.

Differential quadrature lattice Boltzmann method for simulating fluid flows

In this study, a novel method for simulating incompressible flows using the lattice Boltzmann method is presented, aiming to improve robustness. This method is enhanced by the Qadyan numerical method. In the Qadyan method, Partial differential equations are solved using semi-discrete schemes combined with the differential quadrature method. To discretize the spatial derivatives of the lattice Boltzmann equation, the upwind differential quadrature method is employed, while the first-order forward difference scheme and/or fourth-order Runge-Kutta method handle the temporal term. The decision to exclude more complex methods stems from their requirement for refined computational grids. This work focuses on achievable results with similar types of uniform grids. The proposed scheme introduces and evaluates a novel method for solving both steady and unsteady problems. The present numerical method is validated by solving five benchmark problems, including heat diffusion on a slab, Stokes’ first problem, Taylor-Green vortex flow, flow around a flat plate, and flow in a lid-driven cavity. Reasonable agreements are obtained between the solutions obtained using this method and those from analytical/other numerical approaches. The Qadyan formulation delivers more accurate results compared to the finite-difference lattice Boltzmann method, (FDLBM), without requiring an increase in the number of grids. Notably, the novel Qadyan method achieves this increased accuracy without introducing additional complexity to the standard lattice Boltzmann method or resorting to non-uniform grids. This is possible due to the nature of the differential quadrature method, which utilizes more combinations of candidate stencils. The Qadyan method has a higher computational cost in comparison with the FDLBM.

On the influence of matrix flow in the hydraulic permeability of rough-walled fractures

Fractured porous media models may account for the influence of matrix flow in fracture permeability estimates through the slip boundary condition effect, i.e. permeable walls allow fluid to flow along them without viscous resistance. However, matrix flow may also create fluid pathways that unlock further fracture flow capacity, thus increasing its overall apparent permeability. We investigate this effect through numerical experiments by simulating fluid flow in a single fracture-matrix system with varying matrix permeability and calculating the apparent fracture permeability. We observe a consistent increase in fracture permeability with matrix permeability, which is yet underestimated by analytical models that factor in slip boundary conditions. Our study highlights the critical role of the surrounding porous matrix in accurately determining fracture permeability to prevent underestimations of the flow capacity of fractured porous media.

A hybrid multiscale model for fluid flow in fractured rocks using homogenization method with discrete fracture networks

Fluid flow in subsurface tight reservoirs containing pores, microcracks and macrocracks is notably influenced by the characteristics of macro/micro-cracks. A novel hybrid multiscale model is proposed to address the response of macrocracks and pores/microcracks in different spatial scales. Specifically, an equivalent macroscopic model (EMM) deduced from locally periodic representative element volume (REV) is developed using the asymptotic homogenization method to represent the poroelastic behavior of porous medium with microcracks. Simultaneously, the macrocracks are modeled explicitly using the discrete fracture model (DFM), where the hydraulic properties of cracks influenced by fluid pressure gradient is represented by the nonlinear opening/closure behavior. The obtained hybrid model takes into account the heterogeneous nature of fractured rock masses containing pores, micro/macro-cracks, which is fundamental to describe fluid flow behavior in fracture-matrix system. Specialized finite elements, regular meshing technique and adaptive time stepping algorithm are adopted to improve the computational efficiency. The hybrid multiscale model is firstly validated step by step to demonstrate the accuracy and then used to simulate fluid flow in fractured rock reservoir, shedding light on the underlying mechanisms of the enhanced flow capacity resulting from microcrack distribution, connectivity, and macrocrack stimulation.

Coupled CFD-DEM modeling of surface erosion in granular soils: Simulation of erosion function apparatus experiments

This study introduces a numerical modeling approach that couples the computational fluid dynamics (CFD) with the discrete element method (DEM) to simulate grain-scale soil erosion processes induced by water flows. In this modeling framework, CFD simulates fluid flows by solving the volume-averaged Navier-Stokes equations, and uses the k-ω turbulent model for turbulent flows. Simultaneously, DEM computes the displacement of solid particles by incorporating the fluid-particle interactions driven by fluid flows while adhering to Newton's laws of motion. These interactions encompass drag force, buoyancy force, pressure-gradient force, and viscous force exerted by fluid flows and acting on the particles. The coupled CFD-DEM modeling adeptly replicates soil erosion processes, demonstrating good alignment with results obtained from laboratory erosion function apparatus (EFA) tests. In particular, the DEM facilitates the estimation of shear stress acting on the soil surface based on fluid-particle interaction forces, which has been roughly approximated by empirical or semi-empirical models. This study underscores the capability of coupled CFD-DEM in providing valuable insights into the grain-scale behavior of soil particles subjected to fluid flows, with the potential for extension to address soil erosion and fines migration.

Research Progress of Three-dimensional Geological Modeling and Multi-field Coupling of Fractured Oil and Gas Reservoirs

Fractured oil and gas reservoirs play an important role in oil and gas exploration and development due to their complex fracture network structure. As the main fluid channel in the reservoir, fractures have a significant impact on the migration, accumulation and exploitation efficiency of oil and gas. This paper summarizes the key geological characteristics of fractured reservoirs, including fracture types, distribution rules and main controlling factors of fracture development, and discusses the influence of fractures on reservoir physical properties. Furthermore, the deterministic and stochastic modeling methods of three-dimensional geological modeling of fractured reservoirs and the latest progress of discrete fracture network modeling technology are introduced in detail. In addition, this paper also discusses in detail the application of multi-field coupling model in the numerical simulation of fluid flow in fractured reservoirs, including continuous medium model and discrete medium model, and their advantages and disadvantages in simulating fluid flow in fractured reservoirs. Finally, this paper summarizes the research progress of multi-physics coupling model of fractured reservoirs, and points out the direction and challenges of future research. Through these studies, it can provide theoretical support and technical guidance for oil and gas exploration and development, in order to achieve more effective development and management of fractured reservoirs.

Optimized design of Helical-Finned Double Pipe heat exchangers via numerical simulation and Artificial Intelligence

Heat exchangers are essential in industrial applications, where optimizing their design offers significant environmental and economic benefits. However, traditional methods of enhancing heat transfer, such as adding fins, often lead to increased friction and pressure drop. This study aims to optimize the design of helical-finned double-pipe heat exchangers to improve heat transfer efficiency while minimizing frictional losses. Numerical simulations, combined with Artificial Neural Networks (ANNs), were used to model the impact of design parameters (Reynolds number, number of fins, fin height, and twist pitch) on performance. The k-ω shear-stress transport turbulence model and finite volume method were used for fluid flow and heat transfer simulations, while Genetic Algorithms (GAs) were employed for optimization. ANN models demonstrated high predictive accuracy (R2≈1), and the optimized designs achieved a relative Nusselt number of 2.01, a relative friction factor of 1.01, and a Performance Evaluation Criterion (PEC) of 1.43. Notably, fins did not always improve thermal efficiency, illustrating the complexity of optimizing heat exchanger designs. This study presents an approach by integrating ANN and GA, providing an effective strategy for improving heat exchanger performance.

Estimation of added effects and their frequency dependence in various fluid–structure interaction problems

This paper focuses on a fluid–structure interaction topic—the determination of added effects caused by fluid forces acting on a body, considering the standard linear equation of motion. We present various problems that assume small-displacement oscillations of single and multiple bodies in inviscid irrotational (potential) flow or viscous incompressible flow in both closed domain and external flow. For inviscid flow, effects of geometric parameters on the added effects were studied. The presented results extend results known from the literature. For viscous flow, frequency dependence of the added effects was studied for a wide range of frequency. The added effects were computed from data from numerical simulations of fluid flow, where the body oscillations were modeled using the dynamic mesh approach. Effects of the phase shift caused by the dynamic mesh were addressed. The added mass was compared with the corresponding value determined for inviscid flow where applicable. The results show strong dependence of the added effects on many parameters, making their proper computation challenging even for simplified cases.

DualSPHysics+: An enhanced DualSPHysics with improvements in accuracy, energy conservation and resolution of the continuity equation

This paper presents an enhanced version of the well-known SPH (Smoothed Particle Hydrodynamics) open-source code DualSPHysics for the simulation of free-surface fluid flows, leading to the DualSPHysics+ code. The enhancements are made through incorporation of several schemes with respect to stability, accuracy and energy/volume conservation issues in simulating incompressible free-surface fluid flows within the weakly compressible SPH formalism. The Optimized Particle Shifting (OPS) scheme is implemented to improve the accuracy of particle shifting vectors in the free-surface region. To mitigate energy dissipation and maintain consistency, the artificial viscosity in δ-SPH is substituted with a Riemann stabilization term, leading to the δR-SPH. The Velocity divergence Error Mitigating (VEM) and Volume Conservation Shifting (VCS) schemes are adopted in DualSPHysics+ to mitigate the velocity divergence error and improve the volume conservation, and hence to enhance the resolution of the continuity equation. To further reduce both the instantaneous and accumulated errors in velocity divergence, a Hyperbolic/Parabolic Divergence Cleaning (HPDC) scheme is incorporated in addition to the VEM scheme. The implementations of the introduced schemes on both CPU and GPU-based versions of the DualSPHysics+ code along with details on the compilation, running and computational performance are presented. Validations in terms of accuracy, energy conservation and convergence of DualSPHysics+ are shown via several relevant benchmarks. It is demonstrated that a better velocity divergence error cleaning in both instantaneous and accumulated errors can be achieved by the combination of VEM and HPDC. Meanwhile, the excessive energy dissipation by the artificial viscosity is shown to be suppressed by adopting the Riemann stabilization term. Enhanced resolution of the continuity equation along with improved energy conservation of DualSPHysics+ advance the SPH-based simulation of incompressible free-surface fluid flows. Program SummaryProgram title: DualSPHysics+.CPC Library link to program fileshttps://doi.org/10.17632/xnrfv9pgb5.1.Licensing provisions: GNU Lesser General Public License (LGPL).Programming language: C++, CUDA.External dependencies: DualSPHysics (https://dual.sphysics.org).Nature of problem: Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) method and the open-source code DualSPHysics have been widely applied to simulate free-surface fluid flows. Both the general WCSPH method and the more specific DualSPHysics need further improvements in several aspects, including spurious pressure fluctuations, non-conservation of volume and excessive energy dissipation, to enhance the accuracy and stability of simulations.Solution method: DualSPHysics+ implements a set of numerical schemes to enhance the overall accuracy, divergence-free velocity field, invariant density field and energy conservation in the simulation of free-surface fluid flows based on DualSPHysics.