laplacian(tensor,tensor) seg faults
Hello All,
In a recent attempt to implement a tensor based viscosity field I ran into a problem with taking a Laplacian involving two volTensorFields. To perform some more simplified testing I altered laplacianFoam with the following: Code:
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) Code:
|
Maybe this bugreport can help you out, which is very similar:
http://www.openfoam.com/mantisbt/view.php?id=305 Maybe for the solution you have to look in the same direction. |
FYI: bugs should be reported on the bug tracker, so that developers can see them.
OpenFOAM: http://www.openfoam.com/mantisbt/main_page.php OpenFOAM-ext: http://sourceforge.net/apps/mantisbt.../main_page.php Best, |
Perfect!
Thanks a lot Bernhard and Alberto. I didn't post it to the bug tracker outright because I figured I was doing something wrong. It appears that Henry W. fixed it in 2.0.x but the fix hasn't been propagated to the extend branch. I'll post a bug report on it now. I fixed it locally in 1.6-ext in the same manner that Henry did with no issues. |
laplacian(tensor, tensor)
Hello all! Just to let everybody know (although the thread is quite old already and this problem does not arise in newer versions of OF):
The solution of henry posted here http://www.openfoam.com/mantisbt/view.php?id=305 worked for me too Problem: Wanted to use the LRR-RSTM with the originally commented diffusion term in OF-1.6 (src/turbulenceModels/compressible/RAS/LRR/LRR.C) Code:
- fvm::laplacian(Ceps_*rho_*(k_/epsilon_)*R_, epsilon_) Solution: add the line Code:
makeFvLaplacianTypeScheme(SS, symmTensor, symmTensor) \ src/finiteVolume/finiteVolume/laplacianSchemes/laplacianScheme/ recompile the finiteVolume library and that's it! The only strange thing which I could not resolve until now is, that it works fine in serial with Gauss linearUpwind scheme for divSchemes but not in parallel. In parallel I had to change the divSchemes to Gauss Gamma |
Thanks for the extra info on the div schemes! The fact that some schemes aren't functioning in parallel for this is suspicious.
cheers! |
Greetings to all!
@Johannes: Can you provide a simple test case? And which OpenFOAM version/variant are you using? Best regards, Bruno |
Hi! Hello Bruno!
I use OpenFOAM-1.6. Well, with respect to a simple test case, actually I did not try this on a simple case. My case is a stationary flame (so not that simple). I think one can reproduce this problem with any simple case using the LRR model, but with the diffusion term that originally is commented in OF-1.6 Code:
- fvm::laplacian(Ceps_*rho_*(k_/epsilon_)*R_, epsilon_) Code:
- fvm::laplacian(DepsilonEff(), epsilon_) Code:
tmp<surfaceInterpolationScheme<GType> > tinterpGammaScheme_; It seems that I just used the wrong schemes. But as I sad, in this topic I'm not an expert. |
All times are GMT -4. The time now is 07:55. |