A Fourier-based methodology without numerical diffusion for conducting dye simulations and particle residence time calculations

DOI10.1016/J.JCP.2023.112472MaRDI QIDQ6094734zbMATH OpenOpenAlexFDO

Authors Faisal Amlani, Heng Wei, Niema M. Pahlevan

Publication date 10 October 2023

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2112.05257

zbMATH Keywords

advection-diffusionFourier continuationhemodynamics dye simulationhigh-order methodsmathematical physiologyparticle residence time

Mathematics Subject Classification ID

Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems (65Mxx) Basic methods in fluid mechanics (76Mxx) Biological fluid mechanics (76Zxx)

Abstract: Dye experimentation is a widely used method in experimental fluid mechanics for flow analysis or for the study of the transport of particles within a fluid. This technique is particularly useful in biomedical diagnostic applications ranging from hemodynamic analysis of cardiovascular systems to ocular circulation. However, simulating dyes governed by convection-diffusion partial differential equations (PDEs) can also be a useful post-processing analysis approach for computational fluid dynamics (CFD) applications. Such simulations can be used to identify the relative significance of different spatial subregions in particular time intervals of interest in an unsteady flow field. Additionally, dye evolution is closely related to non-discrete particle residence time (PRT) calculations that are governed by similar PDEs. PRT is a widely used metric for various fluid dynamics applications (e.g., environmental fluids, biological flows) and is a well-accepted biomarker for cardiovascular diseases since it is linked to thrombus formation. This contribution introduces a pseudo-spectral method based on Fourier continuation (FC) for conducting dye simulations and non-discrete particle residence time calculations without numerical diffusion errors. Convergence and error analyses are performed with both manufactured and analytical solutions. The methodology is applied to three distinct physical/physiological cases: 1) flow over a two-dimensional (2D) cavity; 2) pulsatile flow in a simplified partially-grafted aortic dissection model; and 3) non-Newtonian blood flow in a Fontan graft. Although velocity data is provided in this work by numerical simulation, the proposed approach can also be applied to velocity data collected through experimental techniques such as from particle image velocimetry.

Recommendations

• scientific article; zbMATH DE number 1547793
• DPIV-driven flow simulation: a new computational paradigm
• Methods for evaluating fluid velocities in spectral simulations of turbulence
• A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow.
• PDE methods in flow simulation post processing

Cites work

• A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems
• A spectral FC solver for the compressible Navier-Stokes equations in general domains. I: Explicit time-stepping
• A stable high-order FC-based methodology for hemodynamic wave propagation
• Accurate, high-order representation of complex three-dimensional surfaces via Fourier continuation analysis
• An FC-based spectral solver for elastodynamic problems in general three-dimensional domains
• Application of the method of manufactured solutions to the verification of a pressure-based finite volume numerical scheme
• Computation of residence time in the simulation of pulsatile ventricular assist devices
• Dispersion analysis of finite element semidiscretizations of the two-dimensional wave equation
• Error dynamics: Beyond von Neumann analysis
• Exponentially-convergent strategies for defeating the Runge phenomenon for the approximation of non-periodic functions. I: Single-interval schemes
• Fourier continuation method for incompressible fluids with boundaries
• Global spectral analysis of multi-level time integration schemes: numerical properties for error analysis
• Global spectral analysis of three-time level integration schemes: focusing phenomenon
• High-order unconditionally stable FC-AD solvers for general smooth domains. I: Basic elements
• High-order unconditionally stable FC-AD solvers for general smooth domains. II: Elliptic, parabolic and hyperbolic PDEs; theoretical considerations
• Higher-order implicit-explicit multi-domain compressible Navier-Stokes solvers
• IMMERSED BOUNDARY METHODS
• Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers?
• LATTICE BOLTZMANN METHOD FOR FLUID FLOWS
• Lattice Boltzmann method. Fundamentals and engineering applications with computer codes
• Lattice Boltzmann simulation of advection-diffusion of chemicals and applications to blood flow
• Lattice-Boltzmann lattice-spring simulations of influence of deformable blockages on blood fluids in an elastic vessel
• Modeling a no-slip flow boundary with an external force field
• Multi-dimensional hybrid Fourier continuation-WENO solvers for conservation laws
• Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws
• On the significance of blood flow shear-rate-dependency in modeling of Fontan hemodynamics
• Simulation of flexible filaments in a uniform flow by the immersed boundary method
• Spatially dispersionless, unconditionally stable FC-AD solvers for variable-coefficient PDEs
• Spectral analysis of finite difference schemes for convection diffusion equation
• The immersed boundary method

This page was built for publication: A Fourier-based methodology without numerical diffusion for conducting dye simulations and particle residence time calculations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6094734)