|
[Sponsors] |
November 17, 2011, 15:34 |
laplacian(tensor,tensor) seg faults
|
#1 |
Senior Member
Kyle Mooney
Join Date: Jul 2009
Location: San Francisco, CA USA
Posts: 323
Rep Power: 18 |
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++) { solve ( fvm::ddt(tau) - fvm::laplacian(dTau, tau) ); } Code:
#0 0x00007ffff76e8850 in Foam::HashTable<Foam::tmp<Foam::fv::laplacianScheme<Foam::Vector<double>, Foam::Tensor<double> > > (*)(Foam::fvMesh const&, Foam::Istream&), Foam::word, Foam::string::hash>::find (this=0x0, key=...) at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/HashTable.C:177 #1 0x00007ffff76de993 in Foam::fv::laplacianScheme<Foam::Vector<double>, Foam::Tensor<double> >::New (mesh=..., schemeData=...) at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/laplacianScheme.C:73 #2 0x00007ffff76d646b in Foam::fvm::laplacian<Foam::Vector<double>, Foam::Tensor<double> > (gamma=..., vf=..., name=...) at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/fvmLaplacian.C:220 #3 0x00007ffff76cf35d in Foam::fvm::laplacian<Foam::Vector<double>, Foam::Tensor<double> > (tgamma=..., vf=..., name=...) at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/fvmLaplacian.C:233 #4 0x00007ffff76c99be in Foam::LPTT::divTau (this=0x4c9c430, U=...) at viscoelasticLaws/LPTT/LPTT.C:146 #5 0x00007ffff7790ad6 in Foam::multiMode::divTau(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const () from /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/lib/linux64GccDPOpt/libviscoelasticTransportModels.so #6 0x00007ffff76c6e89 in Foam::viscoelasticModel::divTau(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const () from /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/lib/linux64GccDPOpt/libviscoelasticTransportModels.so #7 0x000000000041bbf1 in main () Last edited by kmooney; November 17, 2011 at 17:43. |
|
November 18, 2011, 02:05 |
|
#2 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22 |
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. |
|
November 18, 2011, 12:33 |
|
#3 |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
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,
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
|
November 18, 2011, 12:36 |
|
#4 |
Senior Member
Kyle Mooney
Join Date: Jul 2009
Location: San Francisco, CA USA
Posts: 323
Rep Power: 18 |
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. |
|
November 26, 2013, 06:51 |
laplacian(tensor, tensor)
|
#5 |
New Member
Johannes
Join Date: Mar 2011
Location: Austria
Posts: 12
Rep Power: 15 |
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 Last edited by jfw_cfd; November 26, 2013 at 12:05. |
|
November 26, 2013, 09:54 |
|
#6 |
Senior Member
Kyle Mooney
Join Date: Jul 2009
Location: San Francisco, CA USA
Posts: 323
Rep Power: 18 |
Thanks for the extra info on the div schemes! The fact that some schemes aren't functioning in parallel for this is suspicious.
cheers! |
|
November 26, 2013, 17:31 |
|
#7 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,980
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
@Johannes: Can you provide a simple test case? And which OpenFOAM version/variant are you using? Best regards, Bruno
__________________
|
|
November 27, 2013, 03:13 |
|
#8 |
New Member
Johannes
Join Date: Mar 2011
Location: Austria
Posts: 12
Rep Power: 15 |
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_; tmp<snGradScheme<Type> > tsnGradScheme_; It seems that I just used the wrong schemes. But as I sad, in this topic I'm not an expert. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Memory,CPU , page faults problems | pk | FLUENT | 1 | October 27, 2005 04:27 |