Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh Weissenberg number simulations with incompressible Smoothed Particle Hydrodynamics and the log-conformation formulation
105
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Viscoelastic flows occur widely, and numerical simulations of them are important for a range of industrial applications. Simulations of viscoelastic flows are more challenging than their Newtonian counterparts due to the presence of exponential gradients in polymeric stress fields, which can lead to catastrophic instabilities if not carefully handled. A key development to overcome this issue is the log-conformation formulation, which has been applied to a range of numerical methods, but not previously applied to Smoothed Particle Hydrodynamics (SPH). Here we present a 2D incompressible SPH algorithm for viscoelastic flows which, for the first time, incorporates a log-conformation formulation with an elasto-viscous stress splitting (EVSS) technique. The resulting scheme enables simulations of flows at high Weissenberg numbers (accurate up to Wi=85 for Poiseuille flow). The method is robust, and able to handle both internal and free-surface flows, and a range of linear and non-linear constitutive models. Several test cases are considerd included flow past a periodic array of cylinders and jet buckling. This presents a significant step change in capabilties compared to previous SPH algorithms for viscoelastic flows, and has the potential to simulate a wide range of new and challenging applications.
rate research
Read More
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.
In this study we investigated the capabilities of the mesh-free, Lagrangian particle method (Smoothed Particle Hydrodynamics, SPH) to simulate the detailed hydrodynamic processes generated by both spilling and plunging breaking waves within the surf zone. The weakly-compressible SPH code DualSPHysics was applied to simulate wave breaking over two distinct bathymetric profiles (a plane beach and fringing reef) and compared to experimental flume measurements of waves, flows, and mean water levels. Despite the simulations spanning very different wave breaking conditions (including an extreme case with violently plunging waves on an effectively dry reef slope), the model was able to reproduce a wide range of relevant surf zone hydrodynamic processes using a fixed set of numerical parameters. This included accurate predictions of the nonlinear evolution of wave shapes (e.g., asymmetry and skewness properties), rates of wave dissipation within the surf zone, and wave setup distributions. By using this mesh-free approach, the model was able to resolve the critical crest region within the breaking waves, which provided robust predictions of the wave-induced mass fluxes within the surf zone responsible for the undertow. Within this breaking crest region, the model results capture how the potential energy of the organized wave motion is initially converted to kinetic energy and then dissipated, which reproduces the distribution of wave forces responsible for wave setup generation across the surf zone. Overall, the results reveal how the mesh-free SPH approach can accurately reproduce the detailed wave breaking processes with comparable skill to state-of-the-art mesh-based Computational Fluid Dynamics (CFD) models, and thus can be applied to provide valuable new physical insight into surf zone dynamics.
The standard formulation of the smoothed particle hydrodynamics (SPH) assumes that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is underestimated. As a result, the pressure of the low (high) density side is over (under) estimated. Thus, unphysical repulsive force appears at the contact discontinuity, resulting in the effective surface tension. This tension suppresses fluid instabilities. In this paper, we present a new formulation of SPH, which does not require the differentiability of density. Instead of the mass density, we adopt the internal energy density (pressure), and its arbitrary function, which are smoothed quantities at the contact discontinuity, as the volume element used for the kernel integration. We call this new formulation density independent SPH (DISPH). It handles the contact discontinuity without numerical problems. The results of standard tests such as the shock tube, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, point like explosion, and blob tests are all very favorable to DISPH. We conclude that DISPH solved most of known difficulties of the standard SPH, without introducing additional numerical diffusion or breaking the exact force symmetry or energy conservation. Our new SPH includes the formulation proposed by Ritchie & Thomas (2001) as a special case. Our formulation can be extended to handle a non-ideal gas easily.
At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.
In this paper, we present a new formulation of smoothed particle hydrodynamics (SPH), which, unlike the standard SPH (SSPH), is well-behaved at the contact discontinuity. The SSPH scheme cannot handle discontinuities in density (e.g. the contact discontinuity and the free surface), because it requires that the density of fluid is positive and continuous everywhere. Thus there is inconsistency in the formulation of the SSPH scheme at discontinuities of the fluid density. To solve this problem, we introduce a new quantity associated with particles and density of that quantity. This density evolves through the usual continuity equation with an additional artificial diffusion term, in order to guarantee the continuity of density. We use this density or pseudo density, instead of the mass density, to formulate our SPH scheme. We call our new method as SPH with smoothed pseudo-density (SPSPH). We show that our new scheme is physically consistent and can handle discontinuities quite well.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
