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] integration of particle Solver CFDEMcoupling

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 18, 2019, 07:57
Default integration of particle Solver CFDEMcoupling
  #1
Member
 
Patrick Höhn
Join Date: Mar 2018
Location: Celle, Germany
Posts: 53
Blog Entries: 2
Rep Power: 8
hoehnp is on a distinguished road
Hello together,


while integrating CFDEMcoupling as an additional fluidModel in solids4Foam, I face the bug as described in https://www.cfdem.com/forums/bug-fin...using-parallel and polyMesh::findCell slow in parallel, but meshSearch::findCell sometimes fails


I suggested I replaced meshSearch::findCell by polyMesh::findCell. Sofar the code seems to compile and run. Is there anything I should still consider in this replacement?


Thanks a lot everyone.


Patrick
hoehnp is offline   Reply With Quote

Old   October 4, 2019, 04:47
Default reconstructPar does not reconstruct lagrangian fields in regions
  #2
Member
 
Patrick Höhn
Join Date: Mar 2018
Location: Celle, Germany
Posts: 53
Blog Entries: 2
Rep Power: 8
hoehnp is on a distinguished road
Hello everyone,
during the development I am facing a second issue now. solids4foam uses unlike the original FSI toolkit regions for the fluid and the solid phase. When trying to reconstruct a decomposed problem, I found that lagrangian fields are no longer recognized. This seems to be caused by the fact, that the IOobjects in the lagrangian parts of reconstructPar include the region "fluid" in the path, while the actual folder structure would require a region "".

One of course could start completely rewriting this part of reconstructPar. Before doing this I wonder, if anybody has faced and solved similar problems before?

Thanks a lot everyone
Patrick
hoehnp is offline   Reply With Quote

Old   October 12, 2019, 11:19
Default
  #3
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 Patrick,

Yep, it seems that you need to edit reconstructPar such that it treats lagrangian fields the same whether they are in a region or not.

Philip
bigphil is offline   Reply With Quote

Old   October 15, 2019, 06:43
Default Implementation of particle forces and mesh rotation
  #4
Member
 
Patrick Höhn
Join Date: Mar 2018
Location: Celle, Germany
Posts: 53
Blog Entries: 2
Rep Power: 8
hoehnp is on a distinguished road
Hi Philip,


thank you for your answer. I am now facing a few more questions during my further implementations:
1. If I want to introduce the particle forces in the fluidSolidInterface, from my understanding I have to do this in Foam::fluidSolidInterface::transferFacesZoneToZone . Is this assumption correct?
2. How could I handle a rotation of the solid domain? Rotate the solid mesh and also put a rotation boundary on the inner patch of the fluid mesh? Only rotate the fluid mesh?
3. Could I also get write access to your bitbucket repository to share my results so far?


Best Regards,
Patrick
hoehnp is offline   Reply With Quote

Old   October 15, 2019, 11:46
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 Patrick,

See my replies below:

Quote:
1. If I want to introduce the particle forces in the fluidSolidInterface, from my understanding I have to do this in Foam::fluidSolidInterface::transferFacesZoneToZone . Is this assumption correct?
Sort of. The transferFacesZoneToZone function maps a field on the fluid interface to the solid interface. You can see that this function is called once in the fluidSolidInterface class in the void Foam::fluidSolidInterface::updateForce() function at line 1694. You will see that it is mapping the fluidZoneTotalTraction field to the solidZoneTotalTraction field.

For you, I would suggest that you add your particle forces to this fluidZoneTotalTraction field before it is mapped to solidZoneTotalTraction; that way, you do not need to modify the transferFacesZoneToZone function.

One thing to be careful of is that fluidZoneTotalTraction and solidZoneTotalTraction are "zone" fields, meaning that this field is the same on every processor in parallel. The globalPolyPatch class is used to transfer between "patch" fields (local to each processor) and global "zone" field.

Quote:
2. How could I handle a rotation of the solid domain? Rotate the solid mesh and also put a rotation boundary on the inner patch of the fluid mesh? Only rotate the fluid mesh?
The solid domain uses a Lagrangian approach, and as such, its motion comes from solving its governing equation (remember the solid mesh is not independent of the solid motion/displacement/velocity field). Consequently, to impose a rotation on the solid, it should be done through appropriate boundary conditions on the displacement field e.g. a rotation boundary condition at the axle patch of your solid, like the fixedRotation boundary condition.

Then for the fluid, depending on the magnitude of the rotation, you may need to define a rotating region, where I guess you know the approximate rotational velocity in advance, and the fluid mesh motion solver would need to account for the interface motion relative to this rotation.

Quote:
3. Could I also get write access to your bitbucket repository to share my results so far?
Sure. If you have some working tutorials showing new features, you can send me an email describing them and I can grant you permission.

Best,
Philip
bigphil is offline   Reply With Quote

Reply

Tags
cfdem, cfdemcoupling, liggghts, solids4foam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Modify Particle Velocity in Lagrangian Solver AustinG OpenFOAM 0 August 23, 2018 17:38
Does OpenFoam have an appropriate solver for particle collision into a water droplet? Hossein1 OpenFOAM Running, Solving & CFD 3 October 10, 2015 08:28
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
reactingFoam: ODE Solver KRR4 "Too many integration steps" pajofego OpenFOAM 3 May 18, 2013 15:15
Two independent time integration loops in one solver ziemowitzima OpenFOAM Running, Solving & CFD 0 January 22, 2013 14:08


All times are GMT -4. The time now is 01:36.