• OpenLB - Open source lattice Boltzmann code

June 18, 2022 at 6:08 pm #6658

Keymaster

It is hard to tell what the problem is in this case without a full listing or at least a backtrace with debug symbols enabled. So the following is just guesswork:

As the code worked in 1.4 and fails in 1.5 I suspect that the issue is in the refactoring of the old dynamics definition to the new style. e.g. the replacement of your commented out


ForcedBGKdynamics<T, NSDESCRIPTOR> NSbulkDynamics(
  converter.getLatticeRelaxationFrequency(),
  instances::getBulkMomenta<T,NSDESCRIPTOR>());

AdvectionDiffusionBGKdynamics<T, TDESCRIPTOR> TbulkDynamics (
  converter.getLatticeThermalRelaxationFrequency(),
  instances::getAdvectionDiffusionBulkMomenta<T,TDESCRIPTOR>());

// (NS,AD)lattice.defineDynamics(...) in prepareLattice

in the prepareLattice function. Can you provide the relevant section / a full listing of prepareLattice?

It should look similar to:


NSlattice.defineDynamics<ForcedBGKdynamics>(...);
NSlattice.setParameter<descriptors::OMEGA>(converter.getLatticeRelaxationFrequency());

respectively


ADlattice.defineDynamics<AdvectionDiffusionBGKdynamics>(...);
ADlattice.setParameter<descriptors::OMEGA>(converter.getLatticeThermalRelaxationFrequency());

I hope this helps.

Reply To: Segmentation fault for OpenLB 1.5