Skip to content

Reply To: Using setBouzidiAdeDirichlet BC properly

#9776
shota
Participant

Hi talco,

in the following a short example on how to use the bouzidiAdeBoundary:

void prepareADLattice( SuperLattice<T,RADDESCRIPTOR>& adLattice,
                       AdeUnitConverter<T,RADDESCRIPTOR> const& adConverter,
                       STLreader<T>& stlReader,
                       SuperGeometry<T,3>& superGeometry)
{
  auto bulkIndicator = superGeometry.getMaterialIndicator({1});
  
  // Boundary conditions
  ...
  setBouzidiBoundary<T,RADDESCRIPTOR,BouzidiAdeDirichletPostProcessor>(adLattice, superGeometry, 5, stlReader);

  // Initial conditions
  AnalyticalConst3D<T,T> tempA( 10. );
  Vector<T,3> velocityV;
  AnalyticalConst3D<T,T> uF(velocityV);

  ...
  adLattice.defineRho( superGeometry, 5, tempA);
  adLattice.iniEquilibrium(superGeometry, 5, tempA, uF);
  setBouzidiAdeDirichlet(adLattice, superGeometry, 5, tempA);

  adLattice.template setParameter<descriptors::OMEGA>(adConverter.getLatticeRelaxationFrequencyFromDiffusivity<RADDESCRIPTOR>(D));
  adLattice.initialize();
}

I hope this helps.

BR
Shota

  • This reply was modified 3 weeks, 2 days ago by shota.