Complex Geometry Generation

[FernandoIGuerra]

Hi all,
Currently I am implementing a channel with rectangular obstacle in 2D,. Accordingly, I am employing cylinder to achieve so.

void
EBDDT_Geometry::build(
  const amrex::Geometry& geom, const int max_coarsening_level)
{

//channel with obstacles
// first obstacle is a d/2 distance of the right end wall
//height of the obstacle is a d/4
//width of each obstacle is d/16
//distance between obstacle is d
amrex::Real diameter = 2.0;
amrex::Real heigth = 1.2;
// the cylinder input is the radio. Hence the width required is divided into 2 to get
// the wodth
amrex::Real width = diameter/(4*64) ;

  amrex::EB2::CylinderIF cylinder1(width, heigth, 1, {AMREX_D_DECL(1.0,0.0,0.0)}, false);
   amrex::EB2::CylinderIF cylinder2(width, heigth, 1, {AMREX_D_DECL(5.0,0.0,0.0)}, false);
  amrex::EB2::CylinderIF cylinder3(width, heigth, 1, {AMREX_D_DECL(9,0.0,0.0)}, false);
  amrex::EB2::CylinderIF cylinder4(width, heigth, 1, {AMREX_D_DECL(13,0.0,0.0)}, false); 
  amrex::EB2::CylinderIF cylinder5(width, heigth, 1, {AMREX_D_DECL(17,0.0,0.0)}, false);

        
  auto polys1 = amrex::EB2::makeUnion(cylinder1, cylinder2);
  auto polys2 = amrex::EB2::makeUnion(polys1, cylinder3);
  auto polys3 = amrex::EB2::makeUnion(polys2, cylinder4);
  auto polys4 = amrex::EB2::makeUnion(polys3, cylinder5);

  auto gshop = amrex::EB2::makeShop(polys4);
  amrex::EB2::Build(gshop, geom, max_coarsening_level, max_coarsening_level);
}

but by running the following error arises:

Starting to call amrex_probinit ... 
Successfully run amrex_probinit
amrex::Abort::53::pc_fill_sv_ebg: zero apnorm !!!
amrex::Abort::62::pc_fill_sv_ebg: zero apnorm !!!
amrex::Abort::63::pc_fill_sv_ebg: zero apnorm !!!

which indicates that the normal surface is 0. How can I fix this issue. Thanks in adavance

Fernando

[marchdf]

Hi, thanks for reaching out! The easiest way to debug this is to make the geometry as simple as possible, make sure you have the geometry you expect and then add on bit by bit.

I am a little confused in that you talked about having 1 cylinder in your message but I see 5 defined and used.