Transport Velocity Formulation

Some main principles and consequences of the transport velocity.

Full theoretical description of the transport velocity is beyond the scope of this manual, so some main principles and consequences of the transport velocity will be discussed. For full theoretical derivation and analysis, refer to the work of Adami. 1

Accuracy of the SPH method heavily relies on the ability of the code to accurately reconstruct the Shepard coefficient and provide full support to the particles. In reality the value of the Shepard coefficient will be ≈ 1, but rarely exactly 1. If you analyze the SPH method, you can easily understand that the accuracy of the method is actually directly related to the particle distribution. If the particles are uniformly ordered, the reconstruction of the variable fields will be accurate. If there are excesses, such as overly-packed or overly-sparse particle distributions, this will negatively reflect on the accuracy of the solution. It would; therefore, be ideal if you can keep the particles ordered as uniformly as possible without sacrificing computational time.

The transport velocity does precisely this. The numerical formulation of the transport velocity introduces a correction method to otherwise normal SPH velocity computation in a case-independent manner, while preserving the physicality of the solution. The dual correction is applied exclusively through the momentum equation, as opposed to the traditional background pressure approach which is explicitly appearing in the quasi-incompressible equation of state. The correction comes directly from the computation of the transport velocity u ˜ , or more precisely the time advancement of the transport velocity:

u i ˜ ( t + Δ t ) = u i ( t ) + Δ t ( d ˜ u i d t 1 ρ i p c )

Where, p c is the corrective pressure field, usually set to be equivalent to the initial pressure of the simulation p 0 .

These corrections actively maintain particle order which has a number of beneficial influences on the numerical behavior of the code.

The magnitude of the p c pressure directly influences the strength of the correction. The higher the p c value, the more vigorous will be the correction attempt. You should keep this in mind, as specifying the p c value too high, for example, p c = 10 p 0 , can lead to excessive correction force and in these cases the time step must be appropriately reduced.

The command is actually a coefficient with which you multiply the p 0 value and therefore determine the p c correction pressure (by default it is set to 1.0).

p 0 = p c _ f a c t o r * p 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGWbWaaSbaaSqaaiaaicdaaeqaaOGaeyypa0JaamiCamaaBaaa leaacaWGJbaabeaakiaac+facaWGMbGaamyyaiaadogacaWG0bGaam 4BaiaadkhacaGGQaGaamiCamaaBaaaleaacaaIWaaabeaaaaa@46BE@

In more graphical terms, the transport velocity formulation automatically detects particle vacuum and attempts to populate it with particles. As mentioned, this has profoundly beneficial influence in multiphase simulations, but in single phase simulation where you have intentionally left a large portion of the domain empty (particle vacuum), the use of the transport velocity could actually have detrimental effects. The reason is precisely because transport velocity is seeking for particle vacuum and tries to fill it, which would in single phase cases result in a pop-corn like behavior of the free surface. This is something to be avoided and therefore in single phase cases it is strongly recommended to turn the transport velocity off.

Artificial Particle Displacement (APD)

nanoFluidX has two schemes to address this. The transport velocity is used in combination with the Weighted interaction scheme, whereas with the Riemann interaction scheme, you rely on the Artificial Particle Displacement (APD) method.2

The APD method relies on adding a small correction vector to the position of the particle:
δ r i = β j = 1 N r i j r i j 3 r i , o 2 u max Δ t MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiTdqMabC OCayaalaWaaSbaaSqaaiaadMgaaeqaaOGaeyypa0JaeqOSdi2aaabC aeaadaWcaaqaaiqahkhagaWcamaaBaaaleaacaWGPbGaamOAaaqaba aakeaacaWGYbWaa0baaSqaaiaadMgacaWGQbaabaGaaG4maaaaaaaa baGaamOAaiabg2da9iaaigdaaeaacaWGobaaniabggHiLdGccaWGYb Waa0baaSqaaiaadMgacaGGSaGaam4BaaqaaiaaikdaaaGccaaMc8Ua amyDamaaBaaaleaaciGGTbGaaiyyaiaacIhaaeqaaOGaeuiLdqKaam iDaaaa@558C@
Where,
β
Numerical coefficient, typically of the value 0.01 to 0.05 for multiphase flows with bounded domain and recommended as 0 for single phase (free surface; open domain) flows
r i , o MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOCamaaBa aaleaacaWGPbGaaiilaiaad+gaaeqaaaaa@39A8@
Average particle distance between a particle and its neighbors
u max MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyDamaaBa aaleaaciGGTbGaaiyyaiaacIhaaeqaaaaa@39ED@
Maximum expected velocity in the domain (reference velocity)
j MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOAaaaa@36E2@ index
Refers to the standard SPH vector subtraction between the owner particle and its neighbors
Resulting particle fields display better order and provide smoother solutions without noticeable diffusion or detriment to the results. However, it is worth noting that pushing the APD coefficient β beyond the recommended values is likely to introduce negative consequences into the results – most notably numerical diffusion.
Important: Be cautious when changing the values of the APD coefficient ( β ).
1 S. Adami, Modeling and Simulation of Multiphase Phenomena with Smoothed Particle Hydrodynamics, Chair of Aerodynamics and Fluid Mechanics, Technical University of Munich, 2014
2 SPH Modelling of Long-term Sway-Sloshing Motion in a Rectangular Tank, Oezbulut, M.; Tofighi, N.; Yildiz, M.; Goeren Oe.; The Twenty-fifth International Ocean and Polar Engineering Conference, Kona, Hawaii, USA, 2015