[openfoam_aero]

Dear Foamers,


Before this is deemed as an already solved issue, please let me explain the steps I have taken. I have a domain, that has the patches inlet/outlet, left/right and top/bottom walls with a cuboid placed inside such that the base of the cuboid is on the bottom wall. The remaining surfaces of the cuboid is not in contact with the domain. I need to apply cyclic patches to inlet and outlet, left and right sides and the top and bottom can be walls. I am using snappyHexMesh for meshing.


First, I tried this parallel error with cyclic BCs for pimpleDyMFoam and trouble in resuming running (check post #12) where i used cyclicAMI and transform translational along with separation vector. Works for a pair of cyclicAMI surfaces. But in my case, there are 2.



Secondly, I implemented the recommendation of creating topoSetDict, using blockMesh, snappyHexMesh, topoSet, createPatch where I set the entire domain to the patch 'wall', run blockMesh and snappyHexMesh and THEN create cyclicAMI patches. Seems to work fine, please find the checkMeshlog.dat. The order of executing the commands is

1. blockMesh
2. snappyHexMesh -overwrite
3. topoSet
4. createPatch -overwrite


similar to the post #5 in this link - Add / Remove faces from patch , create new


However, when I run icoFoam (or any other solver), I end up getting this:

--> FOAM FATAL ERROR:
Unable to find initial target face

From function bool Foam::AMIMethod::initialise(Foam::labelListList&, Foam::scalarListList&, Foam::labelListList&, Foam::scalarListList&, Foam::label&, Foam::label&)
in file AMIInterpolation/AMIInterpolation/AMIMethod/AMIMethod/AMIMethod.C at line 134.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::AMIMethod::initialise(Foam::List<Foam::List< int> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, int&, int&) at ??:?
#3 Foam::faceAreaWeightAMI::calculate(Foam::List<Foam ::List<int> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, int, int) at ??:?
#4 Foam::AMIInterpolation::update(Foam::PrimitivePatc h<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, bool) at ??:?
#5 Foam::AMIInterpolation::constructFromSurface(Foam: :PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::autoPtr<Foam::searchableSurface> const&, bool) at ??:?
#6 Foam::AMIInterpolation::AMIInterpolation(Foam::Pri mitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::autoPtr<Foam::searchableSurface> const&, Foam::faceAreaIntersect::triangulationMode const&, bool, Foam::AMIInterpolation::interpolationMethod const&, double, bool, bool) at ??:?
#7 Foam::cyclicAMIPolyPatch::resetAMI() const at ??:?
#8 Foam::cyclicAMIPolyPatch::AMIs() const at ??:?
#9 Foam::tmp<Foam::Field<double> > Foam::cyclicAMIPolyPatch::interpolate<double>(Foam ::Field<double> const&, Foam::UList<double> const&) const at ??:?
#10 Foam::cyclicAMIFvPatchField<double>:atchNeighbou rField() const at ??:?
#11 Foam::coupledFvPatchField<double>::evaluate(Foam:: UPstream::commsTypes) at ??:?
#12 Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPa tchField(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#13 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#14 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#15 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam:imensio nedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#16 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#17 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#18 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#19 ? in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
#20 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#21 ? in "/home/uttamcadambi/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/icoFoam"
Aborted (core dumped)

Please note that I have checked this link: CyclicAMI: Unable to find initial target face


I have changed 'noOrdering' to translational. I have used cyclicAMI everywhere. I am not able to debug this and it is really eating my head.


Please find attached the case files. I am also attaching the log for each step.



Any help on this will be largely appreciated.