Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh-order simulations of isothermal flows using the local anisotropic basis function method (LABFM)
91
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Mesh-free methods have significant potential for simulations of flows in complex geometries, with the difficulties of domain discretisation greatly reduced. However, many mesh-free methods are limited to low order accuracy. In order to compete with conventional mesh-based methods, high order accuracy is essential. The Local Anisotropic Basis Function Method (LABFM) is a mesh-free method introduced in King et al., J. Comput. Phys. 415:109549 (2020), which enables the construction of highly accurate difference operators on disordered node discretisations. Here, we introduce a number of developments to LABFM, in the areas of basis function construction, stencil optimisation, stabilisation, variable resolution, and high order boundary conditions. With these developments, direct numerical simulations of the Navier Stokes equations are possible at extremely high order (up to 10th order in characteristic node spacing internally). We numerically solve the isothermal compressible Navier Stokes equations for a range of geometries: periodic and channel flows, flows past a cylinder, and porous media. Excellent agreement is seen with analytical solutions, published numerical results (using a spectral element method), and experiments. The potential of the method for direct numerical simulations in complex geometries is demonstrated with simulations of subsonic and transonic flows through an inhomogeneous porous media at pore Reynolds numbers up to Re=968.
rate research
Read More
Mesh-free methods have significant potential for simulations in complex geometries, as the time consuming process of mesh-generation is avoided. Smoothed Particle Hydrodynamics (SPH) is the most widely used mesh-free method, but suffers from a lack of consistency. High order, consistent, and local (using compact computational stencils) mesh-free methods are particularly desirable. Here we present a novel framework for generating local high order difference operators for arbitrary node distributions, referred to as the Local Anisotropic Basis Function Method (LABFM). Weights are constructed from linear sums of anisotropic basis functions (ABFs), chosen to ensure exact reproduction of polynomial fields up to a given order. The ABFs are based on a fundamental Radial Basis Function (RBF), and the choice of fundamental RBF has small effect on accuracy, but influences stability. LABFM is able to generate high order difference operators with compact computational stencils (4th order with 25 nodes, 8th order with 60 nodes in two dimensions). At domain boundaries (with incomplete support) LABFM automatically provides one-sided differences of the same order as the internal scheme, up to 4th order. We use the method to solve elliptic, parabolic and mixed hyperbolic-parabolic PDEs, showing up to 8th order convergence. The inclusion of hyperviscosity is straightforward, and can effectively provide stability when solving hyperbolic problems. LABFM is a promising new mesh-free method for the numerical solution of PDEs in complex geometries. The method is highly scalable, and for Eulerian schemes, the computational efficiency is competitive with RBF-FD for a given accuracy. A particularly attractive feature is that in the low order limit, LABFM collapses to SPH, and there is potential for Arbitrary Lagrangian-Eulerian schemes with natural adaptivity of resolution and accuracy.
In this numerical study, an original approach to simulate non-isothermal viscoelastic fluid flows at high Weissenberg numbers is presented. Stable computations over a wide range of Weissenberg numbers are assured by using the root conformation approach in a finite volume framework on general unstructured meshes. The numerical stabilization framework is extended to consider thermo-rheological properties in Oldroyd-B type viscoelastic fluids. The temperature dependence of the viscoelastic fluid is modeled with the time-temperature superposition principle. Both Arrhenius and WLF shift factors can be chosen, depending on the flow characteristics. The internal energy balance takes into account both energy and entropy elasticity. Partitioning is achieved by a constant split factor. An analytical solution of the balance equations in planar channel flow is derived to verify the results of the main field variables and to estimate the numerical error. The more complex entry flow of a polyisobutylene-based polymer solution in an axisymmetric 4:1 contraction is studied and compared to experimental data from the literature. We demonstrate the stability of the method in the experimentally relevant range of high Weissenberg numbers. The results at different imposed wall temperatures, as well as Weissenberg numbers, are found to be in good agreement with experimental data. Furthermore, the division between energy and entropy elasticity is investigated in detail with regard to the experimental setup.
We present numerical simulations of three-dimensional thermal convective flows in a cubic cell at high Rayleigh number using thermal lattice Boltzmann (LB) method. The thermal LB model is based on double distribution function approach, which consists of a D3Q19 model for the Navier-Stokes equations to simulate fluid flows and a D3Q7 model for the convection-diffusion equation to simulate heat transfer. Relaxation parameters are adjusted to achieve the isotropy of the fourth-order error term in the thermal LB model. Two types of thermal convective flows are considered: one is laminar thermal convection in side-heated convection cell, which is heated from one vertical side and cooled from the other vertical side; while the other is turbulent thermal convection in Rayleigh-Benard convection cell, which is heated from the bottom and cooled from the top. In side-heated convection cell, steady results of hydrodynamic quantities and Nusselt numbers are presented at Rayleigh numbers of $10^6$ and $10^7$, and Prandtl number of 0.71, where the mesh sizes are up to $257^3$; in Rayleigh-Benard convection cell, statistical averaged results of Reynolds and Nusselt numbers, as well as kinetic and thermal energy dissipation rates are presented at Rayleigh numbers of $10^6$, $3times 10^6$, and $10^7$, and Prandtl numbers of 0.7 and 7, where the nodes within thermal boundary layer are around 8. Compared with existing benchmark data obtained by other methods, the present LB model can give consistent results.
This work presents a new multiphase SPH model that includes the shifting algorithm and a variable smoothing length formalism to simulate multi-phase flows with accuracy and proper interphase management. The implementation was performed in the DualSPHysics code and validated for different canonical experiments, such as the single-phase and multiphase Poiseuille and Couette test cases. The method is accurate even for the multiphase case for which two phases are simulated. The shifting algorithm and the variable smoothing length formalism has been applied in the multiphase SPH model to improve the numerical results at the interphase even when it is highly deformed and non-linear effects become important. The obtained accuracy in the validation tests and the good interphase definition in the instability cases indicate an important improvement in the numerical results compared with single-phase and multiphase models where the shifting algorithm and the variable smoothing length formalism are not applied.
The nonlinear and nonlocal coupling of vorticity and strain-rate constitutes a major hindrance in understanding the self-amplification of velocity gradients in turbulent fluid flows. Utilizing highly-resolved direct numerical simulations of isotropic turbulence in periodic domains of up to $12288^3$ grid points, and Taylor-scale Reynolds number $R_lambda$ in the range $140-1300$, we investigate this non-locality by decomposing the strain-rate tensor into local and non-local contributions obtained through Biot-Savart integration of vorticity in a sphere of radius $R$. We find that vorticity is predominantly amplified by the non-local strain coming beyond a characteristic scale size, which varies as a simple power-law of vorticity magnitude. The underlying dynamics preferentially align vorticity with the most extensive eigenvector of non-local strain. The remaining local strain aligns vorticity with the intermediate eigenvector and does not contribute significantly to amplification; instead it surprisingly attenuates intense vorticity, leading to breakdown of the observed power-law and ultimately also the scale-invariance of vorticity amplification, with important implications for prevailing intermittency theories.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
