[Neb]

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.

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&,

[bigphil]

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?

[Neb]

Hi,
thank you so much for your reply.
I've been working on this simulation for months and can't seem to solve.

Originally Posted by bigphil

What is the boundary condition for Inlet_solido in 0/DD?

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);
    }
}

Originally Posted by bigphil

In addition, what does the boundary condition look like in processor10/2/solid/DD?

I didn't understand if you ask me this:

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]

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]

Originally Posted by Neb

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:

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?