CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions > OpenFOAM CC Toolkits for Fluid-Structure Interaction

[solids4Foam] elasticWallPressure no work in parallel. BC robin.

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 16, 2022, 11:11
Default elasticWallPressure no work in parallel. Don't work movingWallPressure and Robin BC.
  #1
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Hi everyone,
I have a big problem. I have to do a simulation with elasticWallPressure and work the parallel, but it doesn't work in Openfoam 7. I also tried the robin3dTube tutorial and it doesn't work in parallel.

Instead the boundary conditions movingWallPressure and Robin give the following error:

Quote:
FOAM FATAL ERROR:

request for volVectorField ddt(U) from objectRegistry fluid failed
available objects of type volVectorField are

5
(
grad(p)
cellMotionU
U_0
U
HbyA
)
I'm trying to add volVectorField ddt(U) in pisoFoam.C but fail.

Thanks for your help

Last edited by Neb; December 21, 2022 at 10:22.
Neb is offline   Reply With Quote

Old   December 21, 2022, 11:14
Default
  #2
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Neb,

Try the latest version of solids4foam from https://solids4foam.github.io.

You will see that the Robin method is used in the fourth tutorial: https://solids4foam.github.io/tutorials/tutorial4.html

Philip
bigphil is online now   Reply With Quote

Old   December 21, 2022, 11:34
Default
  #3
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Quote:
Originally Posted by bigphil View Post
Hi Neb,

Try the latest version of solids4foam from https://solids4foam.github.io.

You will see that the Robin method is used in the fourth tutorial: https://solids4foam.github.io/tutorials/tutorial4.html

Philip
Thanks, I'll try to use the new version!
Neb is offline   Reply With Quote

Old   December 22, 2022, 05:50
Default
  #4
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Quote:
Originally Posted by bigphil View Post
Hi Neb,

Try the latest version of solids4foam from https://solids4foam.github.io.

You will see that the Robin method is used in the fourth tutorial: https://solids4foam.github.io/tutorials/tutorial4.html

Philip


Hi, I tried the 3dTube tutorial but it explodes after few steps. I notice a rapid growth of the forces exchanged at the interface up to enormous values.
I have used pimpleFluid with both Dirichlet-Neumann and Robin.
What's up? When I tested the tutorials I found errors.
I used Openfoam-9.

The movingWallPressure boundary condition produces the same error of old version, which does not find the field ddt(U).

Thank you
Attached Images
File Type: jpg error.jpg (124.3 KB, 4 views)
Attached Files
File Type: txt testLoopReport.txt (9.7 KB, 0 views)

Last edited by Neb; December 22, 2022 at 08:31.
Neb is offline   Reply With Quote

Old   December 22, 2022, 09:07
Default
  #5
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Neb,

At https://solids4foam.github.io/tutorials/tutorial4.html you will see that the Robin condition is only currently compatible with foam-extend-4.1 (as it requires changes to the pimpleFluid solver).

Currently, the Robin conditions do not work with OpenFOAM-9 and OpenFOAM-v2012. I have just created an issue at https://github.com/solids4foam/solids4foam/issues/26 so this does not get forgotten. The change should not be too difficult to make, but pimpleFluid/pimpleFoam looks a bit different between the different versions.
bigphil is online now   Reply With Quote

Old   December 22, 2022, 10:42
Default
  #6
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Quote:
Originally Posted by bigphil View Post
Hi Neb,

At https://solids4foam.github.io/tutorials/tutorial4.html you will see that the Robin condition is only currently compatible with foam-extend-4.1 (as it requires changes to the pimpleFluid solver).

Currently, the Robin conditions do not work with OpenFOAM-9 and OpenFOAM-v2012. I have just created an issue at https://github.com/solids4foam/solids4foam/issues/26 so this does not get forgotten. The change should not be too difficult to make, but pimpleFluid/pimpleFoam looks a bit different between the different versions.
Yes, now I read about it into the tutorial. I needed openfoam-9 to take advantage of codedFixedValue which is not there in foam-extend. Is there a similar functionality to codedFixedValue in foam-extend? I absolutely need it.
Neb is offline   Reply With Quote

Old   December 22, 2022, 10:50
Default
  #7
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
foam-extend-4.1 (or foam-extend-5.0) does not seem to have a codeFixedValue or similar, although it does have codedFunctionObject.

I see the following options:
  • Create a custom boundary condition for foam-extend, which hard-codes the behaviour you want
  • Use groovyBC, although I have not checked if it is compatible.
  • In soldis4foam, port the Robin-related changes in pimpleFluid.foamextend.C to pimpleFluid.foundation.C
bigphil is online now   Reply With Quote

Old   January 10, 2023, 06:29
Default
  #8
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Quote:
Originally Posted by bigphil View Post
foam-extend-4.1 (or foam-extend-5.0) does not seem to have a codeFixedValue or similar, although it does have codedFunctionObject.

I see the following options:
  • Create a custom boundary condition for foam-extend, which hard-codes the behaviour you want
  • Use groovyBC, although I have not checked if it is compatible.
  • In soldis4foam, port the Robin-related changes in pimpleFluid.foamextend.C to pimpleFluid.foundation.C
Hi, I'm trying to transfer the fixes from the two solvers like you suggested. I've fixed tons of bugs, but there are a few left that I can't fix. I'm not very experienced with this language. Could I have help?


Update:
I fixed other errors and I have only one left on the cumulativeContErr_(0) function. But I realized that this error also appears in the original version of solids4Foam. I attach the log file log.Allwmake_solids4foam.txt.
Attached Files
File Type: txt log.Allwmake.txt (34.8 KB, 0 views)
File Type: c pimpleFluid.foundation.C (18.3 KB, 0 views)
File Type: h pimpleFluid.foundation.H (4.6 KB, 0 views)
File Type: txt log.Allwmake_solids4foam.txt (41.3 KB, 0 views)

Last edited by Neb; January 12, 2023 at 11:24.
Neb is offline   Reply With Quote

Old   January 31, 2023, 04:35
Default
  #9
Neb
Member
 
Join Date: Mar 2020
Posts: 66
Rep Power: 6
Neb is on a distinguished road
Quote:
Originally Posted by bigphil View Post
Hi Neb,

At https://solids4foam.github.io/tutorials/tutorial4.html you will see that the Robin condition is only currently compatible with foam-extend-4.1 (as it requires changes to the pimpleFluid solver).

Currently, the Robin conditions do not work with OpenFOAM-9 and OpenFOAM-v2012. I have just created an issue at https://github.com/solids4foam/solids4foam/issues/26 so this does not get forgotten. The change should not be too difficult to make, but pimpleFluid/pimpleFoam looks a bit different between the different versions.
Hello,
I'm trying my simulation on foam-extend. It gives me this error:

Selecting physicsModel fluidSolidInterface
Selecting fluidSolidInterface fixedRelaxation

Selecting fluidModel pimpleFluid
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading g from constant directory
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model LaunderSharmaKE

PIMPLE: Operating solver in PISO mode

Selecting solidModel linearGeometry
Selecting dynamicFvMesh staticFvMesh


--> FOAM FATAL IO ERROR:
unexpected class name faceCompactList expected faceList
while reading object faces

file: /mnt/c/Sim/constant/solid/polyMesh/faces at line 15.

From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 136.

FOAM exiting

The model is produced with ANSA v 21. With openFoam 7, 9 there are no problems with this model.

Thank you

Last edited by Neb; January 31, 2023 at 09:24.
Neb is offline   Reply With Quote

Old   February 14, 2023, 05:04
Default
  #10
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Does this mesh work with other solves in foam-extend? For example, does it work if you use laplacianFoam?
bigphil is online now   Reply With Quote

Old   February 14, 2023, 05:10
Default
  #11
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Neb,

Can you please post the top 20-30 lines of the solid/polyMesh/faces file?

Philip
bigphil is online now   Reply With Quote

Old   June 3, 2023, 19:51
Default
  #12
New Member
 
Join Date: Jan 2018
Posts: 18
Rep Power: 8
Zhi Cheng is on a distinguished road
Hi, many thanks for your contribution!
I tried to run the microBeamInCrossFlow tutorial case in parallel. This case works well in serial, however, we are faced with following problems while in parallel running:


Interpolating point values using AMI
[1] [9]

[1] [9]
[1]
[9] --> FOAM FATAL ERROR: --> FOAM FATAL ERROR: (openfoam-(openfoam-20122012))
[1]
[9] bad size -1241566988bad size -990378326
[1]

[9] [1]
From Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][9] From
Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][1]
in file [9] /cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C at line in file 100/cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C.
at line [1] 100.
FOAM parallel run aborting


I noticed some others are also being faced with the same issues. Do you have any suggestions on this?

Best,
Zhi
Zhi Cheng is offline   Reply With Quote

Old   June 3, 2023, 19:52
Default
  #13
New Member
 
Join Date: Jan 2018
Posts: 18
Rep Power: 8
Zhi Cheng is on a distinguished road
Hi, many thanks for your contribution!
I tried to run the microBeamInCrossFlow tutorial case in parallel. This case works well in serial, however, we are faced with following problems while in parallel running:


Interpolating point values using AMI
[1] [9]

[1] [9]
[1]
[9] --> FOAM FATAL ERROR: --> FOAM FATAL ERROR: (openfoam-(openfoam-20122012))
[1]
[9] bad size -1241566988bad size -990378326
[1]

[9] [1]
From Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][9] From
Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][1]
in file [9] /cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C at line in file 100/cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C.
at line [1] 100.
FOAM parallel run aborting


I noticed some others are also being faced with the same issues. Do you have any suggestions on this?

Best,
Zhi

Quote:
Originally Posted by bigphil View Post
Hi Neb,

Can you please post the top 20-30 lines of the solid/polyMesh/faces file?

Philip
Zhi Cheng is offline   Reply With Quote

Old   June 6, 2023, 08:49
Default
  #14
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,089
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by Zhi Cheng View Post
Hi, many thanks for your contribution!
I tried to run the microBeamInCrossFlow tutorial case in parallel. This case works well in serial, however, we are faced with following problems while in parallel running:


Interpolating point values using AMI
[1] [9]

[1] [9]
[1]
[9] --> FOAM FATAL ERROR: --> FOAM FATAL ERROR: (openfoam-(openfoam-20122012))
[1]
[9] bad size -1241566988bad size -990378326
[1]

[9] [1]
From Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][9] From
Foam::List<T>::List(Foam::label) [with T = Foam::Vector<double>; Foam::label = int][1]
in file [9] /cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C at line in file 100/cvmfs/soft.computecanada.ca/easybuild/software/2020/avx2/MPI/gcc9/openmpi4/openfoam/v2012/OpenFOAM-v2012/src/OpenFOAM/lnInclude/List.C.
at line [1] 100.
FOAM parallel run aborting


I noticed some others are also being faced with the same issues. Do you have any suggestions on this?

Best,
Zhi
Feel free to add an issue at https://github.com/solids4foam/solids4foam; please also include details of which version of OpenFOAM you are using.

Separately, you may like to check if the other tutorials work in parallel for you.
bigphil is online now   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Make my solver to work in parallel gionni OpenFOAM Programming & Development 1 August 14, 2021 15:14
myBC does not work in parallel uosilos OpenFOAM Running, Solving & CFD 0 July 20, 2021 11:21
chtMultiRegionSimpleFoam: crash on parallel run student666 OpenFOAM Running, Solving & CFD 3 April 20, 2017 11:05
simpleFoam parallel AndrewMortimer OpenFOAM Running, Solving & CFD 12 August 7, 2015 18:45
Cases with small length scale work fine on a single processor but fail in parallel adona058 OpenFOAM Bugs 5 April 17, 2009 04:41


All times are GMT -4. The time now is 11:34.