|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
anand sudhi
Join Date: Sep 2015
Posts: 16
Rep Power: 11 ![]() |
Hi,
the improved version of limited interpolation Schemes for vectors like limitedLinearV, linearUpwindV are different from the the original versions (linearUpwind and limitedLinear). If anyone can tell me the extra modification done to linearUpwindV to account for direction of the field. I understand there is a new variable 'maxCorr' in the code. But i am not able to understand its implementation. Also does this modification make linearUpwindV bounded with out the need of specifying limited gradients? Any comments will be very much appreciated. Thank you Anand |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
anand sudhi
Join Date: Sep 2015
Posts: 16
Rep Power: 11 ![]() |
forAll(faceFlux, facei)
{ vector maxCorr; if (faceFlux[facei] > 0.0) { maxCorr = (1.0 - w[facei])*(vf[nei[facei]] - vf[own[facei]]); sfCorr[facei] = (Cf[facei] - C[own[facei]]) & gradVf[own[facei]]; } else { maxCorr = w[facei]*(vf[own[facei]] - vf[nei[facei]]); sfCorr[facei] = (Cf[facei] - C[nei[facei]]) & gradVf[nei[facei]]; } scalar sfCorrs = magSqr(sfCorr[facei]); scalar maxCorrs = sfCorr[facei] & maxCorr; if (sfCorrs > 0) { if (maxCorrs < 0) { sfCorr[facei] = vector::zero; } else if (sfCorrs > maxCorrs) { sfCorr[facei] *= maxCorrs/(sfCorrs + VSMALL); } } else if (sfCorrs < 0) { if (maxCorrs > 0) { sfCorr[facei] = vector::zero; } else if (sfCorrs < maxCorrs) { sfCorr[facei] *= maxCorrs/(sfCorrs - VSMALL); } } } |
|
![]() |
![]() |
![]() |
Tags |
limited, linearupwind, linearupwindv |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
A few simple questions about linearUpwind and limitedLinear | chegdan | OpenFOAM | 28 | March 30, 2024 02:44 |
strange processor boundary behavior with linearUpwindV | akimbrell | OpenFOAM Bugs | 19 | April 2, 2018 07:42 |
error on linearUpwind | immortality | OpenFOAM Running, Solving & CFD | 3 | May 22, 2013 18:20 |
rhoCentralFoam linearUpwind | Andy_bm | OpenFOAM Running, Solving & CFD | 8 | November 14, 2011 11:34 |
linearUpwind scheme in OpenFOAM 2.0.1 ??? | cabul | OpenFOAM | 8 | November 9, 2011 06:57 |