| 1. |
- Agrawal, Vishal, et al.
(författare)
-
An efficient isogeometric/finite-difference immersed boundary method for the fluid–structure interactions of slender flexible structures
- 2024
-
Ingår i: Computer Methods in Applied Mechanics and Engineering. - : Elsevier BV. - 0045-7825 .- 1879-2138. ; 418
-
Tidskriftsartikel (refereegranskat)abstract
- In this contribution, we present a robust and efficient computational framework capable of accurately capturing the dynamic motion and large deformation/deflection responses of highly-flexible rods interacting with an incompressible viscous flow. Within the partitioned approach, we adopt separate field solvers to compute the dynamics of the immersed structures and the evolution of the flow field over time, considering finite Reynolds numbers. We employ a geometrically exact, nonlinear Cosserat rod formulation in the context of the isogeometric analysis (IGA) technique to model the elastic responses of each rod in three dimensions (3D). The Navier–Stokes equations are resolved using a pressure projection method on a standard staggered Cartesian grid. The direct-forcing immersed boundary method is utilized for coupling the IGA-based structural solver with the finite-difference fluid solver. In order to fully exploit the accuracy of the IGA technique for FSI simulations, the proposed framework introduces a new procedure that decouples the resolution of the structural domain from the fluid grid. Uniformly distributed Lagrangian markers with density relative to the Eulerian grid are generated to communicate between Lagrangian and Eulerian grids consistently with IGA. We successfully validate the proposed computational framework against two- and three-dimensional FSI benchmarks involving flexible filaments undergoing large deflections/motions in an incompressible flow. We show that six times coarser structural mesh than the flow Eulerian grid delivers accurate results for classic benchmarks, leading to a major gain in computational efficiency. The simultaneous spatial and temporal convergence studies demonstrate the consistent performance of the proposed framework, showing that it conserves the order of the convergence, which is the same as that of the fluid solver.
|
|
| 2. |
- Crialesi-Esposito, Marco, et al.
(författare)
-
FluTAS : A GPU-accelerated finite difference code for multiphase flows
- 2023
-
Ingår i: Computer Physics Communications. - : Elsevier BV. - 0010-4655 .- 1879-2944. ; 284
-
Tidskriftsartikel (refereegranskat)abstract
- We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows. Program summary: Program Title: : Fluid Transport Accelerated Solver, FluTAS. CPC Library link to program files: https://doi.org/10.17632/tp6k8wky8m.1 Developer's repository link: https://github.com/Multiphysics-Flow-Solvers/FluTAS.git. Licensing provisions: MIT License. Programming language: Fortran 90, parallelized using MPI and slab/pencil decomposition, GPU accelerated using OpenACC directives. External libraries/routines: FFTW, cuFFT. Nature of problem: FluTAS is a GPU-accelerated numerical code tailored to perform interface resolved simulations of incompressible multiphase flows, optionally with heat transfer. The code combines a standard pressure correction algorithm with an algebraic volume of fluid method, MTHINC [1]. Solution method: the code employs a second-order-finite difference discretization and solves the two-fluid Navier-Stokes equation using a projection method. It can be run both on CPU-architectures and GPU-architectures.
|
|
| 3. |
- Crialesi-Esposito, Marco, et al.
(författare)
-
FluTAS: A GPU-accelerated finite difference code for multiphase flows
-
Ingår i: Computer Physics Communications. - 0010-4655 .- 1879-2944.
-
Tidskriftsartikel (refereegranskat)abstract
- We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows.
|
|
| 4. |
- Dalla Barba, Federico, et al.
(författare)
-
An interface capturing method for liquid-gas flows at low-Mach number
- 2021
-
Ingår i: Computers & Fluids. - : Elsevier Ltd. - 0045-7930 .- 1879-0747. ; 216
-
Tidskriftsartikel (refereegranskat)abstract
- Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the simulation of multiphase, viscous flows where a compressible and an incompressible phase interact in the low-Mach number regime. In this frame, acoustics are neglected but large density variations of the compressible phase can be accounted for as well as heat transfer, convection and diffusion processes. The problem is addressed in a fully Eulerian framework exploiting a low-Mach number asymptotic expansion of the Navier-Stokes equations. A Volume of Fluid approach (VOF) is used to capture the liquid-gas interface, built on top of a massive parallel solver, second order accurate both in time and space. The second-order-pressure term is treated implicitly and the resulting pressure equation is solved with the eigenexpansion method employing a robust and novel formulation. We provide a detailed and complete description of the theoretical approach together with information about the numerical technique and implementation details. Results of benchmarking tests are provided for five different test cases.
|
|
| 5. |
- Demou, Andreas, et al.
(författare)
-
A pressure-based diffuse interface method for low-Mach multiphase flows with mass transfer
- 2022
-
Ingår i: Journal of Computational Physics. - : Elsevier BV. - 0021-9991 .- 1090-2716. ; 448, s. 110730-
-
Tidskriftsartikel (refereegranskat)abstract
- This study presents a novel pressure-based methodology for the efficient numerical solution of a four-equation two-phase diffuse interface model. The proposed methodology has the potential to simulate low-Mach flows with mass transfer. In contrast to the classical conservative four-equation model formulation, the adopted set of equations features volume fraction, temperature, velocity and pressure as the primary variables. The model includes the effects of viscosity, surface tension, thermal conductivity and gravity, and has the ability to incorporate complex equations of state. Additionally, a Gibbs free energy relaxation procedure is used to model mass transfer. A key characteristic of the proposed methodology is the use of high performance and scalable solvers for the solution of the Helmholtz equation for the pressure, which drastically reduces the computational cost compared to analogous density-based approaches. We demonstrate the capabilities of the methodology to simulate flows with large density and viscosity ratios through extended verification against a range of different test cases. Finally, the potential of the methodology to tackle challenging phase change flows is demonstrated with the simulation of three-dimensional nucleate boiling. Superscript/Subscript Available
|
|
| 6. |
|
|
| 7. |
|
|
| 8. |
- Scapin, Nicolo, et al.
(författare)
-
A volume-of-fluid method for interface-resolved simulations of phase-changing two-fluid flows
- 2020
-
Ingår i: Journal of Computational Physics. - : ACADEMIC PRESS INC ELSEVIER SCIENCE. - 0021-9991 .- 1090-2716. ; 407
-
Tidskriftsartikel (refereegranskat)abstract
- We present a numerical method for interface-resolved simulations of evaporating two-fluid flows based on the volume-of-fluid (VoF) method. The method has been implemented in an efficient EFT-based two-fluid Navier-Stokes solver, using an algebraic VoF method for the interface representation, and extended with the transport equations of thermal energy and vaporized liquid mass for the single-component evaporating liquid in an inert gas. The conservation of vaporizing liquid and computation of the interfacial mass flux are performed with the aid of a reconstructed signed-distance field, which enables the use of well-established methods for phase change solvers based on level-set methods. The interface velocity is computed with a novel approach that ensures accurate mass conservation, by constructing a divergence-free extension of the liquid velocity field onto the entire domain. The resulting approach does not depend on the type of interface reconstruction (i.e. can be employed in both algebraic and geometrical VoF methods). We extensively verified and validated the overall method against several benchmark cases, and demonstrated its excellent mass conservation and good overall performance for simulating evaporating two-fluid flows in two and three dimensions.
|
|
| 9. |
|
|
| 10. |
- Scapin, Nicolo, et al.
(författare)
-
Evaporating Rayleigh-B\'enard convection : prediction of interface temperature and global heat transfer modulation
- 2022
-
Tidskriftsartikel (övrigt vetenskapligt/konstnärligt)abstract
- We propose an analytical model to estimate the interface temperature $\Theta_{\Gamma}$ and the Nusselt number $Nu$ for an evaporating two-layer Rayleigh-B\'enard configuration in statistically stationary conditions. The model is based on three assumptions: (i) the Grossmann-Lohse theory for thermal convection can be applied on the liquid and gas layers separately, (ii) the vapour content in the gas can be taken as the mean value at the gas-liquid interface and (iii) the bulk gas temperature can be determined neglecting the contributions of the thermal boundary layers. The resulting model can accommodate non-Oberbeck-Boussinesq effects in the liquid and the gas phases, as well as the variation of the liquid height due to evaporation. To obtain a simplified scaling between $Nu$ and the Rayleigh number $Ra$, we specify the model for the case of an Oberbeck-Boussinesq liquid and a gas phase with uniform properties except for the gas density and the vapour diffusion coefficient, which are functions of thermodynamic pressure, local temperature and vapour composition. We validate this simplified setting using direct numerical simulations for $Ra=10^6, 10^7$ and $10^8$ and for four values of the temperature differential $\varepsilon=0.05,0.10,0.15$ and $0.20$, which modulates the change of state variables in the gas layer. The proposed model agrees very well with the numerical simulations in the entire range of $Ra-\varepsilon$ investigated.
|
|
