Reply To: Using setBouzidiAdeDirichlet BC properly
OpenLB – Open Source Lattice Boltzmann Code › Forums › on Lattice Boltzmann Methods › General Topics › Using setBouzidiAdeDirichlet BC properly › Reply To: Using setBouzidiAdeDirichlet BC properly
January 21, 2025 at 4:36 pm
#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.