fppimenta / rheotool Goto Github PK
View Code? Open in Web Editor NEWToolbox to simulate GNF and viscoelastic fluid flows in OpenFOAM®
License: GNU General Public License v3.0
Toolbox to simulate GNF and viscoelastic fluid flows in OpenFOAM®
License: GNU General Public License v3.0
Dear Francisco Pimenta
I want to use both DynamicRefineFvMesh and dynamicMotionSolverFvMesh techniques together in rheoInterFoam solver. How can I write these functions together in DynamicMeshDict file? I tried many solutions, but all of the solutions used only one of the techniques, however I want use both of them together. How can I solve my problem?
Best regards
I tried to compile rheoTool with OpenFoam v4.1 on CentOS 7.x, but it failed in the linking stage for solvers/rheoFoam and solvers/rheoInterFoam:
-lm -o /mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/bin/rheoFoam
/usr/bin/ld: warning: libEDFModels.so, needed by /mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/lib/libBCRheoTool.so, not found (try using -rpath or -rpath-link)
/mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/lib/libpostProcessingRheoTool.so: undefined reference to Foam::EDFEquation::FK_' /mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/lib/libBCRheoTool.so: undefined reference to
Foam::EDFEquation::eK_'
/mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/lib/libBCRheoTool.so: undefined reference to `Foam::EDFEquation::kbK_'
collect2: error: ld returned 1 exit status
make: *** [/mnt/data1/users/software/OpenFoam/rheoTool/of40/src/build/bin/rheoFoam] Error 1
It appears in solvers/rheoFoam/Make/options and solvers/rheoInterFoam/Make/options EDFModels library link option and include path are not specified.
Adding the path and libray linkng options has fixed the problem:
cp solvers/rheoEFoam/Make/options solvers/rheoFoam/Make/options
sed -i.back -e 's!postProcessing/postProcUtils/lnInclude!postProcessing/postProcUtils/lnInclude -I../../libs/EDFModels/lnInclude!' -e 's!-lpostProcessingRheoTool!-lpostProcessingRheoTool -L$(FOAM_USER_LIBBIN) -lEDFModels!' solvers/rheoInterFoam/Make/options
Hello
First of all, thank you very much for developing and sharing this tool. The documentation is very well written and makes it easy to use.
I am working on the cylinder in the channel and when running the tutorial with of60, I found Cd=117.649 instead of 117.357 as stated in the doc (and in literature).
I tried to decrease the dt value and I found :
Delta Time | Drag Coefficient |
---|---|
0,01 | 117,64941909 |
0,001 | 117,665944642 |
0,0001 | 117,773294144 |
I also tried to refine the mesh (roughly 10% finer in each direction, for every region with dt=0.01) and found Cd=117.615 .
Do you know what may explain such difference ? May it be my rheotool installation ?
Thanks you for your help !
PS: To reproduce the results, I have set up a docker based on the of60 container containing rheotool that can be started by executing the the following script in the case directory:
openfoam6-rheotool.txt
Note that once in the container, environment variables must be sourced by executing source /home/rheotool/envRheotool
Hello guys. I have been reading the user guide but it was not clear to me: is there any solver for Dynamic Oscilatory Shear cases?
Thanks in advance.
Hello,
I have been trying to use the Hypre solver with some of the test examples but it's not working properly. OpenFOAM solvers work fine, while the Hypre solver or preconditioners from the Hypre package are not "registered" in the selection table. I use OF 6.0, I checked all Path settings and did some Petsc tests, also there were no compile issues. The only adjustment I made, is that I don't use libatlas as it makes no sense due to a BIOS/hardware shortfall of my PC.
E.g. case: slipSmoluchowski
Adjustments to fvSolution file:
/*
"(p|phiE|psi)"
{
solver GAMG;
tolerance 1e-8;
relTol 0.;
smoother DIC;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
*/
"(p|phiE|psi)"
{
solverType hypreSolver;
saveSystem true;
solver GMRES;
tolerance 1e-8;
relTol 0.01;
maxIter 1000;
updatePrecondFrequency 1;
updateMatrixCoeffs true;
preconditioner
{
preconditioner BoomerAMG;
maxIter 1;
}
}
Error: ...
Build : 6
Exec : rheoEFoam -parallel
Date : Dec 16 2019
Time : 16:30:56
Host : "sammy"
PID : 16406
I/O : uncollated
Case : /home/klaus/OpenFOAM/ThirdParty-6/rheoTool/of60/tutorials/rheoEFoam/EKmixer/slipSmoluchowski
nProcs : 2
Slaves : 1("sammy.16407")
Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
No MRF models present
No finite volume options present
Using hypreSolver solver for p
Using openFoamSolver solver for U
Selecting constitutive equation Newtonian
Selecting EDF model: slipSmoluchowski
Using dynamicCode for patch inletWest on field phiE at line 26 in "/home/klaus/OpenFOAM/ThirdParty-6/rheoTool/of60/tutorials/rheoEFoam/EKmixer/slipSmoluchowski/processor0/0/phiE.boundaryField.inletWest"
Using dynamicCode for patch inletNorth on field phiE at line 43 in "/home/klaus/OpenFOAM/ThirdParty-6/rheoTool/of60/tutorials/rheoEFoam/EKmixer/slipSmoluchowski/processor0/0/phiE.boundaryField.inletNorth"
SIMPLE: No convergence criteria found
Courant Number mean: 0 max: 0
Starting time loop
Courant Number mean: 0 max: 0
Time = 0.002
Inner iteration: 0
GAMG: Solving for Ux, Initial residual = 0.999999999375, Final residual = 1.35808664139e-10, No Iterations 3
GAMG: Solving for Uy, Initial residual = 0.999999999687, Final residual = 2.03786098887e-09, No Iterations 4
Hypre:GMRES:AMG: Solving for p, Initial residual = 0.999999999795 Final residual = 0.999999999795, No Iterations 0
Hypre:GMRES:AMG: Solving for p, Initial residual = 0.999999999795 Final residual = 0.999999999795, No Iterations 0
Hypre:GMRES:AMG: Solving for p, Initial residual = 0.999999999795 Final residual = 0.999999999795, No Iterations 0
time step continuity errors : sum local = 0.0144155398411, global = -0.0144155397542, cumulative = -0.0144155397542
[1]
[1]
[1] --> FOAM FATAL IO ERROR:
[1] Unknown symmetric matrix solver GMRES
Valid symmetric matrix solvers are :
4
(
GAMG
PBiCGStab
PCG
smoothSolver
)
[1]
[1]
[1] file: IOstream.solvers.(p|phiE|psi) from line 0 to line 0.
[1]
[1] [0]
[0]
[0] --> FOAM FATAL IO ERROR:
[0] Unknown symmetric matrix solver GMRES
Valid symmetric matrix solvers are :
4
(
GAMG
PBiCGStab
PCG
smoothSolver
)
[0]
[0]
[0] file: /home/klaus/OpenFOAM/ThirdParty-6/rheoTool/of60/tutorials/rheoEFoam/EKmixer/slipSmoluchowski/system/fvSolution.solvers.(p|phiE|psi) from line 35 to line 48.
[0]
[0] From function From function static Foam::autoPtr<Foam::lduMatrix::solver> Foam::lduMatrix::solver::New(const Foam::word&, const Foam::lduMatrix&, const Foam::FieldField<Foam::Field, double>&, const Foam::FieldField<Foam::Field, double>&, const lduInterfaceFieldPtrsList&, const Foam::dictionary&)
[1] in file matrices/lduMatrix/lduMatrix/lduMatrixSolver.C at line 74.
[1]
FOAM parallel run exiting
[1]
static Foam::autoPtr<Foam::lduMatrix::solver> Foam::lduMatrix::solver::New(const Foam::word&, const Foam::lduMatrix&, const Foam::FieldField<Foam::Field, double>&, const Foam::FieldField<Foam::Field, double>&, const lduInterfaceFieldPtrsList&, const Foam::dictionary&)
[0] in file matrices/lduMatrix/lduMatrix/lduMatrixSolver.C at line 74.
[0]
FOAM parallel run exiting
[0]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[sammy:16401] 1 more process has sent help message help-mpi-api.txt / mpi-abort
[sammy:16401] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
What could cause the problem?
Klaus
Dear Francesco Pimenta,
I tried to compile rheoTool with OpenFoam v7 on Ubuntu 19.10 and while compiling the following error appeared:
constitutiveEqs/otherModels/ML-IKH/ML_IKH.C:254:41: error: ‘const Mesh’ {aka ‘const class Foam::fvMesh’} has no member named ‘solver’; did you mean ‘Foam::dictionary Foam::solution::solvers_’? (not accessible from this context)
254 | lambdaIEqn.solve(phi().mesh().solver("lambdaI"));
| ^~~~~~
In file included from /opt/openfoam7/src/finiteVolume/lnInclude/fvSolution.H:37,
from /opt/openfoam7/src/finiteVolume/lnInclude/fvMesh.H:56,
from /opt/openfoam7/src/finiteVolume/lnInclude/volMesh.H:42,
from /opt/openfoam7/src/finiteVolume/lnInclude/volFields.H:38,
from lnInclude/constitutiveEq.H:40,
from constitutiveEqs/otherModels/ML-IKH/ML_IKH.H:36,
from constitutiveEqs/otherModels/ML-IKH/ML_IKH.C:27:
/opt/openfoam7/src/OpenFOAM/lnInclude/solution.H:74:20: note: declared private here
74 | dictionary solvers_;
| ^~~~~~~~
make: *** [/opt/openfoam7/wmake/rules/General/transform:26: Make/linux64GccDPInt32Opt/constitutiveEqs/otherModels/ML-IKH/ML_IKH.o] Error 1
Could you please help me with this problem?
Best regards
Dear Francisco Pimenta
I monitored wall shear stress for two phase flows. but I got zero values for WSS. I added the calcWSS.c function to fvsolution. For one -phase flow the results were correct. however, this was not true for two phase flows which are should be run by rheointerFoam solver.
How can i solve my problem?
Best regards
when i want install rheotool on mac os i find a problem
Petsc has been downloaded successfully to directory:
/home/openfoam/ThirdParty/petsc-3.10.2
./installPetsc: 121: ./installPetsc: sudo: not found
./installPetsc: 122: ./installPetsc: sudo: not found
./installPetsc: 123: ./installPetsc: sudo: not found
./installPetsc: 124: ./installPetsc: sudo: not found
/usr/bin/env: 'python': No such file or directory
how can i fi it?pleaswe help
Running Allwmake with Openfoam 6 fails. It appears as if a class referenced isn't available.
constitutiveEqs/otherModels/ML-IKH/ML_IKH.C:255:41: error: ‘const Mesh {aka const class Foam::fvMesh}’ has no member named ‘solver’ lambdaIEqn.solve(phi().mesh().solver("lambdaI")); ^~~~~ make: *** [Make/linux64GccDPInt32Opt/constitutiveEqs/otherModels/ML-IKH/ML_IKH.o] Error 1
Prior to that point, all appears to move forward as it should. Suggestions on proceeding or is this a bug?
Dear Sir..
I am trying to simulate the case rheoEFoam/channelEDF/PowerLaw/NernstPlanck, in which boltzmannEquilibrium bc is applied to the walls for ciAnion and ciCation. However, i wanted to simulate the same case using zeroIonicFlux bc as it is more accurate for but unfortunately it does not give any solution or diverges. I have tried a lot to know where the problem is but could not find it.
I have tried some other cases also with zeroIonicFlux bc but nothing works.
Please help...
Hi:
I have a question about CFL number.The case of Nerstblanck in ICEO has a very small CFL number and the result gets worse if I increase the number.Normally,we set the CFL to 2or1,but why it's so small in our case?Is there any theory about it?
When trying to compile with the arm compilers, after a successful OpenFOAM v7 build, I get errors of the following type when trying to compile LMatrixs.C
in /lib/fvmb/
(the LMatrixs.C src compiles fine with gcc) :
LMatrixs.C:38:1: error: cannot define or redeclare 'registerDebugSwitch' here because namespace 'fvmb' does not enclose namespace 'RegisterDebugSwitch<Foam::fvmb::LMatrix<double> >' defineTemplateTypeNameAndDebug(LMatrix<scalar>, 0); ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/className.H:131:5: note: expanded from macro 'defineTemplateTypeNameAndDebug' defineTemplateTypeNameAndDebugWithName(Type, #Type, DebugSwitch) ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/className.H:126:5: note: expanded from macro 'defineTemplateTypeNameAndDebugWithName' defineTemplateDebugSwitchWithName(Type, Name, DebugSwitch) ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/defineDebugSwitch.H:119:5: note: expanded from macro 'defineTemplateDebugSwitchWithName' registerTemplateDebugSwitchWithName(Type, Name) ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/defineDebugSwitch.H:78:42: note: expanded from macro 'registerTemplateDebugSwitchWithName' Foam::RegisterDebugSwitch<Type>::registerDebugSwitch(Name) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^ LMatrixs.C:39:1: error: cannot define or redeclare 'registerDebugSwitch' here because namespace 'fvmb' does not enclose namespace 'RegisterDebugSwitch<Foam::fvmb::LMatrix<Foam::Vector<double> > >' '
Arm compiler version 20.0 (build number 51) (based on LLVM 9.0.1)
Hi Francisco,
I just submitted this to the cfd-online forum, but I'll also leave the question here:
All the files I will be talking about can be found here:
https://github.com/ptorres-uom/rheoI...am-drippingBUG
And some snapshots here:
https://imgur.com/a/tw9x9VD
Essentially I have an axisymmetric simulation of a channel dispensing fluid in an air chamber where I define the boundary conditions for alpha.water and velocity in the air chamber to be essentially zeroGradient in the outlet "walls", and pressure=0 in the same boundaries. I have done this for jets with no issues, but at lower speeds we have dripping instead of jetting and the droplet does not leave the domain smoothly - instead having a small "bounce" and then exploding all over the domain and creating nonsensical results. I have tried setting the lower boundary to alpha.water=0 and this does not solve the problem.
What do you think I might be doing wrong? Or can this be an issue with the package?
Thank you for any help,
Pedro Torres
Hi,
I hope you are keeping well.
Can Boger, Newtonian and shear thinning elastic fluids be modelled using rheoFoam. If not can the constitutiveProperties be changed to transportProperties of various fluids such as Newtonian (from the OpenFOAM user guide)?
Also, how can the drag coefficient/Strouhal number be plotted in post-processing.
Your help is much appreciated, thanks in advance.
When I try to download eigen using the downloadEigen script in of60, it gives a 404 error and the bitbucket link for the eigen repository is either moved or not working. Can somebody help me install Eigen for Rheotool?
This issue is not directly related to one of your solvers, but it's related to the interface to PETsc
Hi,
First of all, thanks for this marvelous piece of software you're maintaining!
I'm attempting to solve a coupled system (for two scalars T1, T2, to which I already have a working segregated solver):
From what I gather from rheoFoam.C
, I should:
autoPtr<coupledSolver> cs (new coupledSolver(...));
cs->insertMesh(mesh);
cs->insertField(T1);
cs->insertField(T2);
// Start time loop
// S1 and S2 are matrices for T2
// Define equations
auto T1Eqn( fvm::ddt(C1, T1) + fvc::(phi1G) ); // gravity treated explicitely
auto T2InT1Eqn ( - fvm::laplacian(-M1, T2) + S1 );
cs->insertEquation(T1.name(), T1.name(), T1Eqn);
cs->insertEquation(T1.name(), T2.name(), T2InT1Eqn);
// Define T2Eqn and T1InT2Eqn similarly and insert them into the matrix
auto T2Eqn( fvm::laplacian(-M2, T2) + fvc::(phi2G) + S2);
auto T1InT2Eqn ( - fvm::ddt(C2, T1) );
cs->insertEquation(T2.name(), T2.name(), T1Eqn);
cs->insertEquation(T2.name(), T1.name(), T1InT2Eqn);
// Then
cs->solve();
// End time loop
But this doesn't really work; So I decided to check here if the interface is general enough to support such circumstances.
ddt
terms in this way (Should it "just work"?)I'm new to PETsc. Appologies if this seems trivial
Hello everybody
I solved the problem of flow in a curved duct and captured the secondary flow patterns and I could see the deviation in the velocity profile using simpleFoam as reported by Bara in his thesis (*). However, when I used rheoFoam the velocity profile became symmetric (non-physical results). I checked many things, and I found that is because of the div(phi,U) scheme that was set on :GaussDefCmpw none. I changed it to GaussDefCmpw cubista and got the correct results. This scheme (GaussDefCmpw none) is used in many of the viscoelastic cases. Why does this scheme have such a effect on the results? why doesn't affect the accuracy of the results of viscoelastic cases?
Thanks
(*) An experimental and numerical study of the Dean problem: flow development towards two-dimensional multiple solutions
Sir, i am trying to solve some basic cases where electroosmosis with temperature dependency and the effects of joule heating in electroosmosis. How can i modify the solver to get the desired result or just i should add the temperature as added to icoFoam solver.
Also i have some problem with test cases of rheoEFoam. Can i get your mail id to communicate further.
Hi Francisco,
I've been trying to install your tool for a couple of days now but I can't seem to compile, I'm wondering if you can help? Here is what I have tried:
On Ubuntu 20.04
Installed using OpenFoam 7 from https://openfoam.org/news/ubuntu-20-04-lts/ using packages.
kaleido-user@djs-of7-rheotool-20-04-intel-n2-standard-16:~/rheoTool/of70/src$ ./Allwmake
Allwmake /home/kaleido-user/rheoTool/of70/src
Allwmake /home/kaleido-user/rheoTool/of70/src/libs
On Ubuntu 16.04
Installed using OpenFoam 7 from https://openfoam.org/news/ubuntu-20-04-lts/ using packages.
Failed with error:
kaleido-user@djs-of7-rheotool-16-04-intel-n2-standard-16:~/rheoTool/of70/src$ ./Allwmake Allwmake /home/kaleido-user/rheoTool/of70/src
Allwmake /home/kaleido-user/rheoTool/of70/src/libs
On Ubuntu 20.04
OpenFoam 7 compiled and installed from source. Instructions here: https://openfoamwiki.net/index.php/Installation/Linux/OpenFOAM-7/Ubuntu/20.04
kaleido-user@djs-of7-rheotool-20-04-intel-n2-standard-1-src:~/rheoTool/of70/src$ ./Allwmake
Allwmake /home/kaleido-user/rheoTool/of70/src
Allwmake /home/kaleido-user/rheoTool/of70/src/libs
Other information i've seen you ask for on other questions on here:
gcc -v
Using built-in specs.
COLLECT_GCC=gcc
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/9/lto-wrapper
OFFLOAD_TARGET_NAMES=nvptx-none:hsa
OFFLOAD_TARGET_DEFAULT=1
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu 9.3.0-17ubuntu1~20.04' --with-bugurl=file:///usr/share/doc/gcc-9/README.Bugs --enable-languages=c,ada,c++,go,brig,d,fortran,objc,obj-c++,gm2 --prefix=/usr --with-gcc-major-version-only --program-suffix=-9 --program-prefix=x86_64-linux-gnu- --enable-shared --enable-linker-build-id --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --libdir=/usr/lib --enable-nls --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object --disable-vtable-verify --enable-plugin --enable-default-pie --with-system-zlib --with-target-system-zlib=auto --enable-objc-gc=auto --enable-multiarch --disable-werror --with-arch-32=i686 --with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --with-tune=generic --enable-offload-targets=nvptx-none=/build/gcc-9-HskZEa/gcc-9-9.3.0/debian/tmp-nvptx/usr,hsa --without-cuda-driver --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 9.3.0 (Ubuntu 9.3.0-17ubuntu1~20.04)
gfortran -v
Using built-in specs.
COLLECT_GCC=gfortran
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/9/lto-wrapper
OFFLOAD_TARGET_NAMES=nvptx-none:hsa
OFFLOAD_TARGET_DEFAULT=1
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu 9.3.0-17ubuntu1~20.04' --with-bugurl=file:///usr/share/doc/gcc-9/README.Bugs --enable-languages=c,ada,c++,go,brig,d,fortran,objc,obj-c++,gm2 --prefix=/usr --with-gcc-major-version-only --program-suffix=-9 --program-prefix=x86_64-linux-gnu- --enable-shared --enable-linker-build-id --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --libdir=/usr/lib --enable-nls --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object --disable-vtable-verify --enable-plugin --enable-default-pie --with-system-zlib --with-target-system-zlib=auto --enable-objc-gc=auto --enable-multiarch --disable-werror --with-arch-32=i686 --with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --with-tune=generic --enable-offload-targets=nvptx-none=/build/gcc-9-HskZEa/gcc-9-9.3.0/debian/tmp-nvptx/usr,hsa --without-cuda-driver --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 9.3.0 (Ubuntu 9.3.0-17ubuntu1~20.04)
mpifort -v
Using built-in specs.
COLLECT_GCC=/usr/bin/gfortran
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/9/lto-wrapper
OFFLOAD_TARGET_NAMES=nvptx-none:hsa
OFFLOAD_TARGET_DEFAULT=1
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu 9.3.0-17ubuntu1~20.04' --with-bugurl=file:///usr/share/doc/gcc-9/README.Bugs --enable-languages=c,ada,c++,go,brig,d,fortran,objc,obj-c++,gm2 --prefix=/usr --with-gcc-major-version-only --program-suffix=-9 --program-prefix=x86_64-linux-gnu- --enable-shared --enable-linker-build-id --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --libdir=/usr/lib --enable-nls --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object --disable-vtable-verify --enable-plugin --enable-default-pie --with-system-zlib --with-target-system-zlib=auto --enable-objc-gc=auto --enable-multiarch --disable-werror --with-arch-32=i686 --with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --with-tune=generic --enable-offload-targets=nvptx-none=/build/gcc-9-HskZEa/gcc-9-9.3.0/debian/tmp-nvptx/usr,hsa --without-cuda-driver --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 9.3.0 (Ubuntu 9.3.0-17ubuntu1~20.04)
mpicc -v
Using built-in specs.
COLLECT_GCC=/usr/bin/gcc
COLLECT_LTO_WRAPPER=/usr/lib/gcc/x86_64-linux-gnu/9/lto-wrapper
OFFLOAD_TARGET_NAMES=nvptx-none:hsa
OFFLOAD_TARGET_DEFAULT=1
Target: x86_64-linux-gnu
Configured with: ../src/configure -v --with-pkgversion='Ubuntu 9.3.0-17ubuntu1~20.04' --with-bugurl=file:///usr/share/doc/gcc-9/README.Bugs --enable-languages=c,ada,c++,go,brig,d,fortran,objc,obj-c++,gm2 --prefix=/usr --with-gcc-major-version-only --program-suffix=-9 --program-prefix=x86_64-linux-gnu- --enable-shared --enable-linker-build-id --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --libdir=/usr/lib --enable-nls --enable-clocale=gnu --enable-libstdcxx-debug --enable-libstdcxx-time=yes --with-default-libstdcxx-abi=new --enable-gnu-unique-object --disable-vtable-verify --enable-plugin --enable-default-pie --with-system-zlib --with-target-system-zlib=auto --enable-objc-gc=auto --enable-multiarch --disable-werror --with-arch-32=i686 --with-abi=m64 --with-multilib-list=m32,m64,mx32 --enable-multilib --with-tune=generic --enable-offload-targets=nvptx-none=/build/gcc-9-HskZEa/gcc-9-9.3.0/debian/tmp-nvptx/usr,hsa --without-cuda-driver --enable-checking=release --build=x86_64-linux-gnu --host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 9.3.0 (Ubuntu 9.3.0-17ubuntu1~20.04)
which gcc
/usr/bin/gcc
which gfortran
/usr/bin/gfortran
which mpifort
/usr/bin/mpifort
which mpicc
/home/kaleido-user/OpenFOAM/OpenFOAM-7/bin/mpicc
which mpicxx
/usr/bin/mpicxx
Hi!
I have been trying to compile Rheotool onto my OpenFOAM v.7 on ubuntu 18.04. I get an error when doing ./InstallPetsc,
and my openFoam version is v.7 and there are only one version of openFoam on my ubuntu 18.04
also my gcc -v is 64 bit and gfortran is also 64 bit .
here are the error messages:
TESTING: checkFortranNameMangling from config.compilers
Unknown Fortran name mangling: Are you sure the C and Fortran compilers are compatible?
Perhaps one is 64 bit and one is 32 bit?
The docs seems to contradict themselves: 2.4 talks about setting EIGEN_REHO
, but 2.5 says Eigen must be installed in a specified directory.
In any case, compiling with EIGEN_REHO
set as suggested didn't work for me
wmake libso constitutiveEquations
g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -isystem/opt/gridware/libs/eigen/3.2.1 -I/constitutiveEqs/utils -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/finiteVolume/lnInclude -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/meshTools/lnInclude -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/transportModels/twoPhaseMixture/lnInclude -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/transportModels/interfaceProperties/lnInclude -IlnInclude -I. -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/OpenFOAM/lnInclude -I/opt/gridware/apps/gcc/openfoam/6/OpenFOAM-6/src/OSspecific/POSIX/lnInclude -fPIC -c constitutiveEqs/constitutiveEq/constitutiveEq.C -o Make/linux64GccDPInt32Opt/constitutiveEqs/constitutiveEq/constitutiveEq.o
constitutiveEqs/constitutiveEq/constitutiveEq.C:28:50: fatal error: Eigen/Dense: No such file or directory
#include <Eigen/Dense> // For eigen decomposition
^
compilation terminated.
make: *** [Make/linux64GccDPInt32Opt/constitutiveEqs/constitutiveEq/constitutiveEq.o] Error 1
but instead I used
export EIGEN_RHEO=/opt/gridware/libs/eigen/3.2.1/include/eigen3
Perhaps this is Eigen version specific, but it would be worth including.
Hi,
I hope you are keeping well.
I am trying to run the flow around a cylinder tutorial, I am following the steps but I am getting an error (when I run rheoFoam). The rest of the tutorials I tried worked.
Please see attached error. Can you advise on what I can do, your help is much appreciated.
Best Regards
I am facing the following problem:
"rheoFoam: error while loading shared libraries: libHYPRE-2.18.2.so: cannot open shared object file: No such file or directory"
Anyone has a clue on how may I solve it?
Thanks in advance.
In the User guide its says, with reference to the installation folder:
A good location for rheoTool is, for example, directory $WM PROJECT USER DIR, which is defined by default when OpenFOAM⃝R or foam-extend is installed.
It isn't clear how to define the installation folder. At the moment it gets put in $FOAM_USER_APPBIN
.
How do I tell Allwmake where to put rheoTool?
Thanks
sir, i am having some convergence issue in rheoEFoam and tried to change the ddt schemes but it seems other schemes are not working in rheoEFoam and the error comes. What needs to be done in that case?
Hi,
I am wondered if there is a solver similar to multiphaseinterfoam of openfoam in rheotool. In fact, I am trying to simulate 3-layer viscoelastic+Newtonian flow.
Thanks
Hello. I have a question regarding GaussDefCmpw
I'm trying to implement CUBISTA within OpenFOAM's original LimiterFunc
. The limiter is defined by
max(min(min(1.5*r, (3.0 + r)/4.0), 1.5), 0)
according to
A convergent and universally bounded interpolation scheme for the treatment of advection, MA Alves, 2003
But the benchmark result is very terrible. I'm just wondering that two methods are theoretically equivalent. Why the original LimiterFunc
won't work and I have to use GaussDefCmpw
you provided? Thanks.
The installation instructions don't make it obvious that I need to clone the repo in order to get the of60/src
directory.
Rheo-interfoam can only account for multi phase flow. How to add the electrical body force to the solver and surface charge to the drop in a multi-phase platform.
Hello,
I have compiled and used RheoTool v3 and v5 on my workstation with OpenFOAM v6 and serial and parallel runs it works well.
I have asked our High-Performance Computing team at my institution to install the rheoTool v5 for use with OpenFOAM v7 on our HPC system. The compilation appears to run without any errors being generated. However, we have run into an issue when running cases in parallel. A serial case will run fine and the log file will show the residuals values and timesteps while it is being solved, and fully solved fields (U, p, tau, theta etc.). However, when running in parallel it appears that simulation will not solve and the residuals produce nan values and all the fields also produce nan values. This has been for the tutorial cases for rheoFoam (Channel/Oldroyd-BLog). Log files for the serial case and the parallel case are attached.
Do you know of any issues regarding rheoTool v5 and OpenFOAM v7 with running in parallel?
Thanks.
Line 44 in fade0eb
Hi,
I hope you are keeping well.
I have been trying to compile Rheotool onto my OpenFOAM v.6 on ubuntu 18.04. I get an error when doing ./InstallPetsc, I get the configure is confusing and doesn't know what compiler to select and use (Shown in the pic).
Your help is much appreciated. (I am new to coding)
Hi,
I cannot find the contribution of psi potential in species equation defined in NernstPlanckCoupled. Also, there is not a Poisson equation for psi as in the segregated solver (laplacian(psi)==-rhoE). Am I missing something? (using rheoTool for OpenFOAM 6).
Kind regards,
syavash
Dear Fppimenta
I am interested in calculating the lift coefficient of a cylinder in a channel with a viscoelastic flow, depending on its y-position. I made some computations and got some results and to explain the results obtained, I then checked the Newtonian case for Re=0 (ok, I should have done this in first place !).
Unexpectedly, the lift coefficient for the Newtonian case doesn't remain null when moving the cylinder in the y direction. Even if at the center it is null (
I have doubled checked that the convective term is disabled in fvSchemes, tried multiple convergence test (in time, mesh, etc) and I am not able to figure out what is wrong.
Have you an idea ?
If you want to have a look, I set up a minimum working example here based on your tutorial:
https://cloud.univ-grenoble-alpes.fr/index.php/s/drbj62AKBwZMNBx
Thanks a lot in advance for any suggestion you may have !
Guillaume
Hello, it would be greatly helpful if someone solve the issue.
[Ubuntu 20.04 & OpenFOAM-dev]
While following the installation guide of rheoTool, I have successfully downloaded and installed Eigen and Petsc. And they've also been sourced in the bashrc.
Following the fourth step which is to compile the rheotool, there has been a following error
Make/linux64GccDPInt32Opt/files:37: *** missing separator. Stop.
Make/linux64GccDPInt32Opt/files:37: *** missing separator. Stop.
I have added missing separator "SOURCE += " before coupled/coupledSolver.C (line 37) in
rheoTool-master/of70/src/libs/sparseMatrixSolvers/Make/linux64GccDPInt32Opt/files.
Sequentially the error was corrected and the following error appeared
segregated/hypreSolver/hypreSolver.C:159:31: error: ‘Foam::coupledPolyPatch::transformType’ has not been declared
159 | coupledPolyPatch::transformType::ROTATIONAL
| ^~~~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C:176:31: error: ‘Foam::coupledPolyPatch::transformType’ has not been declared
176 | coupledPolyPatch::transformType::ROTATIONAL
| ^~~~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C:192:31: error: ‘Foam::coupledPolyPatch::transformType’ has not been declared
192 | coupledPolyPatch::transformType::ROTATIONAL
| ^~~~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C: In member function ‘void Foam::hypreSolver::computeAllocationHypre(Foam::fvMatrix&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&)’:
segregated/hypreSolver/hypreSolver.C:734:53: error: ‘const class Foam::cyclicAMIPolyPatch’ has no member named ‘neighbPatch’
734 | const cyclicAMIPolyPatch& neicamipp = camipp.neighbPatch();
| ^~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C: In member function ‘void Foam::hypreSolver::assembleHypreAbx(hypre_IJMatrix_struct*&, hypre_IJVector_struct*&, hypre_IJVector_struct*&, Foam::fvMatrix&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&, int, int)’:
segregated/hypreSolver/hypreSolver.C:1039:28: error: ‘Foam::coupledPolyPatch::transformType’ has not been declared
1039 | coupledPolyPatch::transformType::ROTATIONAL
| ^~~~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C:1090:51: error: ‘const class Foam::cyclicAMIPolyPatch’ has no member named ‘neighbPatch’
1090 | const cyclicAMIPolyPatch& neicamipp = camipp.neighbPatch();
| ^~~~~~~~~~~
segregated/hypreSolver/hypreSolver.C: In member function ‘void Foam::hypreSolver::assembleHypreBx(hypre_IJVector_struct*&, hypre_IJVector_struct*&, Foam::fvMatrix&, const Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&, int)’:
segregated/hypreSolver/hypreSolver.C:1287:30: error: ‘Foam::coupledPolyPatch::transformType’ has not been declared
1287 | coupledPolyPatch::transformType::ROTATIONAL
| ^~~~~~~~~~~~~
make: *** [/opt/openfoam-dev/wmake/rules/General/transform:26: Make/linux64GccDPInt32Opt/segregated/hypreSolver/hypreSolvers.o] Error 1
Thank you!!
Hi fppimenta ,
May you kindly assist me, I am using rheoTool. Trying to run the rheoFoam solver using the one of the tutorials provided GiesekusLog. In the solver provided they are simulating the fluid flow in a channel. I want to simulate this using a wedge, I tried using one of the blockMeshDict files from the rheoFoam GNF - Newtonian tutorial .It gives me the below errors all and I did was to change the geomerty from a channel to a wedge. I further tried to create my own circular pipe blockMesh but it still gives me the below error
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
No MRF models present
No finite volume options present
Using openFoamSolver solver for p
Using openFoamSolver solver for U
Selecting constitutive equation GiesekusLog
Selecting thermoFunction: Constant
Selecting thermoFunction: Constant
SIMPLE: No convergence criteria found
SIMPLE: No convergence criteria found
SIMPLE: Operating solver in transient mode with 1 outer corrector
SIMPLE: Operating solver in PISO mode
Courant Number mean: 0.166666666667 max: 20.0000000003
Starting time loop
Courant Number mean: 0.166666666667 max: 20.0000000003
Time = 0.05
Inner iteration: 0
DICPCG: Solving for Ux, Initial residual = 1, Final residual = 8.85012301813e-11, No Iterations 9
DICPCG: Solving for Uy, Initial residual = 1.68230678437e-35, Final residual = 1.68230678437e-35, No Iterations 0
DICPCG: Solving for Uz, Initial residual = 2.61013189937e-34, Final residual = 2.61013189937e-34, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 1.69922271366e-11, No Iterations 128
time step continuity errors : sum local = 1.53609117535e-09, global = 1.54235196027e-11, cumulative = 1.54235196027e-11
DILUPBiCG: Solving for thetaxx, Initial residual = 1, Final residual = 4.73377676548e-14, No Iterations 5
DILUPBiCG: Solving for thetaxy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for thetaxz, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for thetayy, Initial residual = 1, Final residual = 4.91394049142e-14, No Iterations 5
DILUPBiCG: Solving for thetayz, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for thetazz, Initial residual = 1, Final residual = 1.9270971297e-14, No Iterations 5
ExecutionTime = 0.08 s ClockTime = 0 s
Courant Number mean: 41.5991225201 max: 287.839559909
Time = 0.1
Inner iteration: 0
DICPCG: Solving for Ux, Initial residual = 0.34212118528, Final residual = 3.77858587337e-11, No Iterations 7
DICPCG: Solving for Uy, Initial residual = 0.184547365086, Final residual = 1.44042045095e-11, No Iterations 8
DICPCG: Solving for Uz, Initial residual = 0.737359389751, Final residual = 2.44027650758e-11, No Iterations 7
DICPCG: Solving for p, Initial residual = 0.00408777212477, Final residual = 8.09619199808e-12, No Iterations 126
time step continuity errors : sum local = 3.02672895644e-09, global = 2.47038699989e-11, cumulative = 4.01273896016e-11
DILUPBiCG: Solving for thetaxx, Initial residual = 0.260162746334, Final residual = 4.58120438555e-11, No Iterations 4
DILUPBiCG: Solving for thetaxy, Initial residual = 1, Final residual = 6.3196806269e-15, No Iterations 5
DILUPBiCG: Solving for thetaxz, Initial residual = 0.999904252247, Final residual = 3.19930763771e-11, No Iterations 3
DILUPBiCG: Solving for thetayy, Initial residual = 0.264513260497, Final residual = 1.77389571288e-14, No Iterations 5
DILUPBiCG: Solving for thetayz, Initial residual = 0.997980409006, Final residual = 7.08204831597e-11, No Iterations 4
DILUPBiCG: Solving for thetazz, Initial residual = 0.301114558716, Final residual = 8.78553338396e-12, No Iterations 4
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libpthread.so.0"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 expf64 in "/lib/x86_64-linux-gnu/libm.so.6"
#5 Foam::constitutiveEq::calcEig(Foam::GeometricField<Foam::SymmTensor, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam::Tensor, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#6 Foam::constitutiveEqs::GiesekusLog::correct() at ??:?
#7 ? in "/home/machine/OpenFOAM/machine-7/platforms/linux64GccDPInt32Opt/bin/rheoFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? in "/home/machine/OpenFOAM/machine-7/platforms/linux64GccDPInt32Opt/bin/rheoFoam"
Floating point exception (core dumped)
Dear Francisco Pimenta,
I want to compile Rheotool master version 3.0 on Ubuntu 18.04. I installed Eigen library successfully ,however I will be faced with this error:
Allwmake /home/kiakaren/rheoTool-master/of60/src
Allwmake /home/kiakaren/rheoTool-master/of60/src/libs
Could you please help me with this problem?
Best regards
Sasan kheirandish
Hello, I have a problem with the ./Allwmake command. When I run it, it gives this error:
Allwmake /home/jcbottino/rheoTool-master/of70/src
Allwmake /home/jcbottino/rheoTool-master/of70/src/libs
I use Ubuntu 20.04 and g++ 6.4.0
I get the following error compiling with arm compilers (which I know you don't support).
Wondering if you have seen anything similar before, or know what the fix could be?
In file included from operators/blockOperators.C:26:
In file included from operators/blockOperators.H:39:
In file included from operators/div.H:39:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/finiteVolume/lnInclude/volFields.H:37:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/GeometricFields.H:34:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/GeometricScalarField.H:38:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/GeometricField.H:42:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/regIOobject.H:41:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/IOobject.H:79:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/NamedEnum.H:38:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/HashTable.H:566:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/HashTable.C:30:
In file included from /lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/List.H:316:
/lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/src/OpenFOAM/lnInclude/List.C:211:26: error: no matching function for call to 'distance'
List<T>(first, last, std::distance(first, last))
^~~~~~~~~~~~~
operators/div.C:147:13: note: in instantiation of function template specialization 'Foam::List<long>::List<int>' requested here
labelList cT(6,-1);
^
/opt/allinea/20.0.0.0/gcc-9.2.0_Generic-AArch64_SUSE-12_aarch64-linux/lib/gcc/aarch64-linux-gnu/9.2.0/../../../../include/c++/9.2.0/bits/stl_iterator_base_funcs.h:138:5: note: candidate template ignored: substitution failure [with _InputIterator = int]: no type named 'difference_type' in 'std::iterator_traits<int>'
distance(_InputIterator __first, _InputIterator __last)
^
1 error generated.
make: *** [/lustre/home/ba-corourke/bin/XCI/BUILD_FOAM_TEST/ARM/OpenFOAM-7/wmake/rules/General/transform:26: Make/linuxAArch64ClangDPInt64Opt/operators/blockOperators.o] Error 1
Hello, fppimenta
Thank you for developing such a useful tool,
when I try to compile rheoTool in centos-6, I've had some problems.
I follow this link to compile OpenFOAM-6 in centos-6, everything works fine
https://openfoamwiki.net/index.php/Installation/Linux/OpenFOAM-6/CentOS_SL_RHEL#CentOS_6.10_.28without_SCL.29
Configuring PETSc to compile on your system
--with-fc=/home/lzh/OpenFOAM/ThirdParty-6/platforms/linux64Gcc48/openmpi-2.1.1/bin/mpifort is specified along with --with-mpi-dir=/home/lzh/OpenFOAM/ThirdParty-6/platforms/linux64Gcc48/openmpi-2.1.1 which implies using /home/lzh/OpenFOAM/ThirdParty-6/platforms/linux64Gcc48/openmpi-2.1.1/bin/mpif90.
configure is confused and does not know which compiler to select and use! Please specify either [mpi] compilers or --with-mpi-dir - but not both!
In most cases, specifying --with-mpi-dir - and not explicitly listing compilers could be preferable.
to solve this problem, I delete "with-mpi-dir ", the configure is look like this:
./configure --with-debugging=0 --download-hypre --download-parmetis --download-metis --download-ptscotch --download-mumps --download-scalapack --with-precision=double --with-shared-libraries=1 --with-scalar-type=real --with-fc=$mpiFort --with-cc=$mpiCC --with-cxx=$mpiCXX CXX_LINKER_FLAGS=-Wl,--no-as-needed CFLAGS="-g -O2 -fPIC -fstack-protector-strong -Wformat -Werror=format-security" CXXFLAGS="-g -O2 -fPIC -fstack-protector-strong -Wformat -Werror=format-security" FCFLAGS="-g -O2 -fstack-protector-strong" FFLAGS="-g -O2 -fstack-protector-strong" CPPFLAGS="-Wdate-time -D_FORTIFY_SOURCE=2" LDFLAGS="-Wl,-Bsymbolic-functions -Wl,-z,relro" MAKEFLAGS=w
C compiler you provided with -with-cc=/home/lzh/OpenFOAM/ThirdParty-6/platforms/linux64Gcc48/openmpi-2.1.1/bin/mpicc cannot be found or does not work.
Cannot compile C with /home/lzh/OpenFOAM/ThirdParty-6/platforms/linux64Gcc48/openmpi-2.1.1/bin/mpicc.
Could you give me some hint?
Hi,
I am performing transient simulations with rheoFoam for a viscoelastic fluid, and I have realized that residual control is not working. This is my fvSolution.SIMPLE:
SIMPLE
{
nInIter 20;
nNonOrthogonalCorrectors 0;
pRefCell `0;`
pRefValue 0;
residualControl
{
U 1e-3;
p 1e-4;
}
}
and the solvers seems to detect it, as sees in the log:
...
Reading field p
Reading field U
Reading/calculating face flux field phi
No MRF models present
No finite volume options present
Selecting constitutive equation Oldroyd-BLog
SIMPLE: Convergence criteria found
U: tolerance 0.001
p: tolerance 0.0001
Constructing face velocity Uf
...
but, it does not finish the iterations when the tolerance is reached:
...
DICPCG: Solving for p, Initial residual = 6.66936302471e-06, Final residual = 9.05055469255e-11, No Iterations 280
time step continuity errors : sum local = 3.60673177163e-18, global = -1.64336289958e-21, cumulative = -3.48595785849e-18
DILUPBiCG: Solving for thetaxx, Initial residual = 2.67170260934e-09, Final residual = 4.80165417222e-15, No Iterations 1
DILUPBiCG: Solving for thetaxy, Initial residual = 4.31444085111e-09, Final residual = 8.98125828349e-15, No Iterations 1
DILUPBiCG: Solving for thetaxz, Initial residual = 6.60873120169e-09, Final residual = 1.7853899454e-14, No Iterations 1
DILUPBiCG: Solving for thetayy, Initial residual = 7.34315939268e-09, Final residual = 5.2628551726e-14, No Iterations 1
DILUPBiCG: Solving for thetayz, Initial residual = 3.82597468662e-09, Final residual = 1.38512305502e-14, No Iterations 1
DILUPBiCG: Solving for thetazz, Initial residual = 2.20213545615e-09, Final residual = 1.7186490182e-15, No Iterations 1
Inner iteration: 17
DICPCG: Solving for Ux, Initial residual = 3.68764188913e-07, Final residual = 9.23093950545e-12, No Iterations 3
DICPCG: Solving for Uy, Initial residual = 1.62164633149e-06, Final residual = 4.55061515852e-11, No Iterations 3
DICPCG: Solving for Uz, Initial residual = 2.25684877287e-07, Final residual = 4.9095227096e-12, No Iterations 3
DICPCG: Solving for p, Initial residual = 6.66429510978e-06, Final residual = 9.83202783838e-11, No Iterations 279
time step continuity errors : sum local = 3.91810636528e-18, global = -3.04615734829e-21, cumulative = -3.48900401584e-18
DILUPBiCG: Solving for thetaxx, Initial residual = 2.65943929899e-09, Final residual = 4.80115918637e-15, No Iterations 1
DILUPBiCG: Solving for thetaxy, Initial residual = 4.28414014436e-09, Final residual = 9.03823690697e-15, No Iterations 1
DILUPBiCG: Solving for thetaxz, Initial residual = 6.57404526402e-09, Final residual = 1.76928126639e-14, No Iterations 1
DILUPBiCG: Solving for thetayy, Initial residual = 7.27898434841e-09, Final residual = 4.7993380484e-14, No Iterations 1
DILUPBiCG: Solving for thetayz, Initial residual = 3.79017453322e-09, Final residual = 1.35110271671e-14, No Iterations 1
DILUPBiCG: Solving for thetazz, Initial residual = 2.19071150407e-09, Final residual = 1.78510477553e-15, No Iterations 1
Inner iteration: 18
DICPCG: Solving for Ux, Initial residual = 3.68214260462e-07, Final residual = 9.21792845251e-12, No Iterations 3
DICPCG: Solving for Uy, Initial residual = 1.61651614552e-06, Final residual = 4.55098190599e
-11, No Iterations 3
DICPCG: Solving for Uz, Initial residual = 2.25129870065e-07, Final residual = 4.90005586699e-12, No Iterations 3
DICPCG: Solving for p, Initial residual = 6.68393364572e-06, Final residual = 9.62462003748e-11, No Iterations 279
time step continuity errors : sum local = 3.83544984885e-18, global = -2.81013028623e-21, cumulative = -3.49181414613e-18
DILUPBiCG: Solving for thetaxx, Initial residual = 2.64809261994e-09, Final residual = 4.80074325712e-15, No Iterations 1
DILUPBiCG: Solving for thetaxy, Initial residual = 4.25564669653e-09, Final residual = 9.0941442038e-15, No Iterations 1
DILUPBiCG: Solving for thetaxz, Initial residual = 6.54172230786e-09, Final residual = 1.75434175854e-14, No Iterations 1
DILUPBiCG: Solving for thetayy, Initial residual = 7.21914661713e-09, Final residual = 4.43877275503e-14, No Iterations 1
DILUPBiCG: Solving for thetayz, Initial residual = 3.7566295381e-09, Final residual = 1.32052941499e-14, No Iterations 1
DILUPBiCG: Solving for thetazz, Initial residual = 2.180107147e-09, Final residual = 1.8406693216e-15, No Iterations 1
Inner iteration: 19
DICPCG: Solving for Ux, Initial residual = 3.67688609184e-07, Final residual = 9.20532055853e-12, No Iterations 3
DICPCG: Solving for Uy, Initial residual = 1.61160230323e-06, Final residual = 4.55109222577e-11, No Iterations 3
DICPCG: Solving for Uz, Initial residual = 2.24604859284e-07, Final residual = 4.89087271634e-12, No Iterations 3
DICPCG: Solving for p, Initial residual = 6.72449513199e-06, Final residual = 9.33686887486e-11, No Iterations 279
...
Is that something I am doing wrong or is it a bug?
Btw, congratulations for the solvers and library. I think they are a superb work.
Robert
Hello everyone
May you kindly assist me, I am using Ubuntu 18.0 . After installation of PETSC, it prompted the installation successful. Now when run the rheomultiRegion tutorials they run successfully but when I change input files in one of the tutorials according to the problem I am working on I get the below error :
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 50152059.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
I am not sure on how to resolve this issue, may you kindly assist with guidance on how to resolve this .
Thanks
Hi. I am getting an error while compiling rheotool with OpenFOAM6, on running Allwmake. Pasting the error message here, I tried googling, but to no avail. Can somebody help me with this problem? Thanks.
In file included from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/List.H:309:0,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/HashTable.C:30,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/HashTable.H:566,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/NamedEnum.H:38,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/IOobject.H:79,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/regIOobject.H:41,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricField.H:42,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricScalarField.H:38,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricFields.H:34,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/finiteVolume/lnInclude/volFields.H:37,
from operators/div.H:39,
from operators/blockOperators.H:39,
from operators/blockOperators.C:26:
/home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/List.C: In instantiation of ‘Foam::List<T>::List(InputIterator, InputIterator) [with InputIterator = int; T = long int]’:
operators/div.C:147:20: required from here
/home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/List.C:211:39: error: no matching function for call to ‘distance(int&, int&)’
List<T>(first, last, std::distance(first, last))
~~~~~~~~~~~~~^~~~~~~~~~~~~
In file included from /usr/include/c++/7/bits/stl_algobase.h:66:0,
from /usr/include/c++/7/bits/char_traits.h:39,
from /usr/include/c++/7/string:40,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/string.H:51,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/word.H:42,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/fileName.H:48,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/IOobject.H:75,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/regIOobject.H:41,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricField.H:42,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricScalarField.H:38,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/GeometricFields.H:34,
from /home/thebhatman99/OpenFOAM/OpenFOAM-6/src/finiteVolume/lnInclude/volFields.H:37,
from operators/div.H:39,
from operators/blockOperators.H:39,
from operators/blockOperators.C:26:
/usr/include/c++/7/bits/stl_iterator_base_funcs.h:138:5: note: candidate: template<class _InputIterator> typename std::iterator_traits<_Iterator>::difference_type std::distance(_InputIterator, _InputIterator)
distance(_InputIterator __first, _InputIterator __last)
^~~~~~~~
/usr/include/c++/7/bits/stl_iterator_base_funcs.h:138:5: note: template argument deduction/substitution failed:
/usr/include/c++/7/bits/stl_iterator_base_funcs.h: In substitution of ‘template<class _InputIterator> typename std::iterator_traits<_Iterator>::difference_type std::distance(_InputIterator, _InputIterator) [with _InputIterator = int]’:
/home/thebhatman99/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/List.C:211:39: required from ‘Foam::List<T>::List(InputIterator, InputIterator) [with InputIterator = int; T = long int]’
operators/div.C:147:20: required from here
/usr/include/c++/7/bits/stl_iterator_base_funcs.h:138:5: error: no type named ‘difference_type’ in ‘struct std::iterator_traits<int>’
/home/thebhatman99/OpenFOAM/OpenFOAM-6/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt64Opt/operators/blockOperators.o' failed
make: *** [Make/linux64GccDPInt64Opt/operators/blockOperators.o] Error 1
Dear professor:
what is the definition of psi in rheoTool ?In EDL,zeta potential and surface potential may have a slightly different in some cases.
Dear Francisco Pimenta
I am using rheoInterFoam solver for viscoelastic multi-phase fluids. In constitutive properties that is in constant folder, there is a parameter named uTauCoupling , which can be set as true or false.
what does utaucoupling do in viscoelastic fluids?
Best regards
Any idea why a floating point exception is thrown when trying to solve the 1D case for the ion-exchange membrane when setting ddt scheme to steady-state?
Output:
Running rheoEFoam in solution1D and ouputing to null.
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#7 Foam::NernstPlanck::correct() at ??:?
#8 ? at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? at ??:?
Floating point exception (core dumped)
Sir, I think there is some incorrect definitions and dimensions in rheoTool appendix A. Please look into the matter and clarify those.
Waiting for your reply..
I have some question about the calculation of psi and EDL.I sent a email to you because I failed to open github but now is ok.
My case is to describe the thermophoresis(The motion of a spherical particle in response to a temperature gradient).I use the model Nernstplanck to be my basic model.In this basic model I add the temperature field and the electric body force and the dielectric force.My case shows perfect result with 3D mesh.
Now,i want to use 2D mesh with my case,because the number of grids becomes very large when ka>10(a:radius of sphericity).So I use wedge type to describe a sphere in the middle of a cube with 2D mesh.And I successfully get the result when ka=0.1.The relative error is within 1%.But when ka=10,it falsed.I can not get the right distribution of psi and EDL with my own model,so I do some test using the basic model Nernstplanck,it failed too and the result is as same with my own.In order to describe it as clearly as possible, I wrote a file to the attachment.
Case test.pdf
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.