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] problem with restarting

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 4, 2022, 15:37
Default problem with restarting
  #1
Neb
Member
 
Join Date: Mar 2020
Posts: 59
Rep Power: 4
Neb is on a distinguished road
Hi,
I am using solids4foam on openfoam-7. I need to restart my simulation but I get the following error.
I also have another question: in restarting the tutorial I notice that the forces exchanged in the fluid interface are different from those of the solid interface. Is this normal to happen?
Thanks for help.
S.M.

Code:
Selecting physicsModel fluidSolidInteraction

Selecting fluidSolidInterface method fixedRelaxation

Selecting fluidModel pimpleFluid
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading g from constant directory
Using dynamicCode for patch Inlet on field U at line 28 in "/home/prova/processor0/2/fluid/U.boundaryField.Inlet"
Using dynamicCode for patch Outlet on field p at line 32 in "/home/prova/processor0/2/fluid/p.boundaryField.Outlet"
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon

PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Constructing face velocity Uf

Selecting solidModel linearGeometry
Selecting dynamicFvMesh staticFvMesh
Creating solidTraction boundary condition
[10] 
[10] 
[10] --> FOAM FATAL IO ERROR: 
[10] compound has already been transferred from token
    on line 51 the empty compound of type List<vector>
[10] 
[10] file: /home/prova/processor10/2/solid/DD.boundaryField.Inlet_solido.value at line 51.
[10] 
[10]     From function Foam::token::compound& Foam::token::transferCompoundToken(const Foam::Istream&)
[10]     in file lnInclude/token.C at line 99.
[10] 
FOAM parallel run aborting
[10] 
[10] #0  Foam::error::printStack(Foam::Ostream&)[11] 
[11] 
[11] --> FOAM FATAL IO ERROR: 
[11] compound has already been transferred from token
    on line 51 the empty compound of type List<vector>
[11] 
[11] file: /home/prova/processor11/2/solid/DD.boundaryField.Inlet_solido.value at line 51.
[11] 
[11]     From function Foam::token::compound& Foam::token::transferCompoundToken(const Foam::Istream&)
[11]     in file lnInclude/token.C at line 99.
[11] 
FOAM parallel run aborting
[11] 
[11] #0  Foam::error::printStack(Foam::Ostream&) addr2line failed
[10] #1  Foam::IOerror::abort() addr2line failed
[10] #2  Foam::token::transferCompoundToken(Foam::Istream const&) addr2line failed
[10] #3  Foam::Istream& Foam::operator>><Foam::Vector<double> >(Foam::Istream&, Foam::List<Foam::Vector<double> >&) addr2line failed
[10] #4  Foam::Field<Foam::Vector<double> >::Field(Foam::word const&, Foam::dictionary const&, int) addr2line failed
[10] #5  Foam::fixedDisplacementZeroShearFvPatchVectorField::fixedDisplacementZeroShearFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) addr2line failed
[10] #6  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::fixedDisplacementZeroShearFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) addr2line failed
[10] #7  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&,

Last edited by Neb; July 5, 2022 at 05:05.
Neb is offline   Reply With Quote

Old   July 7, 2022, 12:32
Default
  #2
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,026
Rep Power: 32
bigphil will become famous soon enoughbigphil will become famous soon enough
What is the boundary condition for Inlet_solido in 0/DD?
In addition, what does the boundary condition look like in processor10/2/solid/DD?
bigphil is offline   Reply With Quote

Old   July 7, 2022, 12:49
Default
  #3
Neb
Member
 
Join Date: Mar 2020
Posts: 59
Rep Power: 4
Neb is on a distinguished road
Hi,
thank you so much for your reply.
I've been working on this simulation for months and can't seem to solve.

Quote:
Originally Posted by bigphil View Post
What is the boundary condition for Inlet_solido in 0/DD?

Code:
dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    Outlet_solido
    {
        type            fixedDisplacementZeroShear;
     //   type            fixedDisplacement;   
        refValue        uniform (0 0 0);
        value           uniform (0 0 0);
    }
    Inlet_solido
    {
        type            fixedDisplacementZeroShear;
     //   type            fixedDisplacement;
        refValue        uniform (0 0 0);
        value           uniform (0 0 0);
    }
    wall_fluid
    {
        type            solidTraction;
        traction        uniform ( 0 0 0 );
        pressure        uniform 0;  //wss medio
        value           uniform (0 0 0);
    }
    wall_solid
    {
        type            solidTraction;
        traction        uniform ( 0 0 0 );
        pressure        uniform 0;   //pressione addominale
        value           uniform (0 0 0);
    }
}

Quote:
Originally Posted by bigphil View Post
In addition, what does the boundary condition look like in processor10/2/solid/DD?
I didn't understand if you ask me this:
Code:
dimensions      [0 1 0 0 0 0 0];

internalField   nonuniform List<vector> 
5337
In the meantime, I'm doing other tests. May I ask if there is a way to get a .dat file of AMI weights? I am trying to write a new function object in controlDict with coded, but I cannot access the weightsSum values ​​of AMIInterpolation.
Neb is offline   Reply With Quote

Old   July 8, 2022, 07:08
Default
  #4
Neb
Member
 
Join Date: Mar 2020
Posts: 59
Rep Power: 4
Neb is on a distinguished road
I'm trying to simulate flow in an artery. The mechanical characteristics of the wall are:
-poisson 0.45;
-young 1.5 * 10 ^ 6;
-density 2000;
After a lot of testing I noticed the following:
- if I run simulations based on the 3dTube tutorial, entering my boundary conditions, everything works. Also the restart.
-if I do the same simulation on the straight pipe I created, everything works. Also the restart.
- If I work with a pipe with a bend, the simulation is very slow, with 100000 with 100000 initial external cycles and the restart does not work. I can reduce the tolerance on the external FSI loop and finish the simulation faster. I had thought of using this simulation as the initial condition of a new one and bring the FSI tolerance to 10 ^ -2, but the restart doesn't work and I'm stuck.
I would have liked to see the AMI weights in a file, because my suspicion is that the wall is too deformable and for some reason, the solver struggles with coupling. But these are all hypotheses that I can't verify and after thousands (yes, I think I've done them even more) simulations, I can't go on.

If anyone can help me understand what happens, I would be very happy.
Neb is offline   Reply With Quote

Old   July 13, 2022, 04:57
Default
  #5
Neb
Member
 
Join Date: Mar 2020
Posts: 59
Rep Power: 4
Neb is on a distinguished road
Quote:
Originally Posted by Neb View Post
I'm trying to simulate flow in an artery. The mechanical characteristics of the wall are:
-poisson 0.45;
-young 1.5 * 10 ^ 6;
-density 2000;
After a lot of testing I noticed the following:
- if I run simulations based on the 3dTube tutorial, entering my boundary conditions, everything works. Also the restart.
-if I do the same simulation on the straight pipe I created, everything works. Also the restart.
- If I work with a pipe with a bend, the simulation is very slow, with 100000 with 100000 initial external cycles and the restart does not work. I can reduce the tolerance on the external FSI loop and finish the simulation faster. I had thought of using this simulation as the initial condition of a new one and bring the FSI tolerance to 10 ^ -2, but the restart doesn't work and I'm stuck.
I would have liked to see the AMI weights in a file, because my suspicion is that the wall is too deformable and for some reason, the solver struggles with coupling. But these are all hypotheses that I can't verify and after thousands (yes, I think I've done them even more) simulations, I can't go on.

If anyone can help me understand what happens, I would be very happy.
update:
I restart a simulation with icoFluid, RAS laminar, changing the solid boundary condition from fixedDispalcementZeroShear to fixedDispalcement.
If I try with icoFluid, RAS laminar and fixedDispalcementZeroShear, it always gives the above mentioned error.
If instead I try the parallel restart I get the following error:

Code:
Selecting mechanical law linearElastic
additionalMeshCorrection: false
Selecting interfaceToInterfaceMapping AMI
Creating pointHistory function object.
[16] History point ID: 2415
[16] History point coordinates: (0.0981215 0.0781573 0.07)
[16] Reference point coordinates: (0.09 0.01519 0.07)
Courant Number mean: 0.0164553 max: 0.159312 velocity magnitude: 0.113218
deltaT = 0.00119048
Time = 1

Setting traction on solid interfaces
Interpolating face values using AMI
Create AMI zone-to-zone interpolator
AMI: Creating addressing and weights between 712000 source faces and 712000 target faces
AMI: using globalPolyPatch
[2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] Unable to set source and target faces
[2] 
[2]     From function void Foam::newFaceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const boolList&, Foam::labelList&, const Foam::DynamicList<int>&, bool) const [with SourcePatch = Foam::PrimitivePatch<Foam::List<Foam::face>, Foam::Field<Foam::Vector<double> > >; TargetPatch = Foam::PrimitivePatch<Foam::List<Foam::face>, Foam::Field<Foam::Vector<double> > >; Foam::label = int; Foam::boolList = Foam::List<bool>; Foam::labelList = Foam::List<int>]
[2]     in file lnInclude/newFaceAreaWeightAMI.C at line 294.
[2] 
FOAM parallel run aborting
Could there be a bug?
Neb is offline   Reply With Quote

Reply

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
BuoyantBoussinesqSimpleFoam_Facing problem Mondal131211 OpenFOAM Running, Solving & CFD 1 April 10, 2019 20:41
Is COMSOL Multi Physics is suitable to solve complex flow problem? steve lee Main CFD Forum 7 October 28, 2018 22:02
Mesh& steptime independant: conduction-convection problem Fati1 Main CFD Forum 1 October 28, 2018 14:52
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13


All times are GMT -4. The time now is 12:35.