37 void ignitionSite::findIgnitionCells(
const fvMesh&
mesh)
44 label ignCell = mesh.findCell(location_);
50 scalar radius = diameter_/2.0;
56 cellVolumes_[0] = vols[ignCell];
58 scalar minDist = GREAT;
59 label nearestCell = 0;
64 scalar
dist =
mag(centres[celli] - location_);
72 if (dist < radius && celli != ignCell)
75 cellVolumes_.
setSize(nIgnCells+1);
77 cells_[nIgnCells] = celli;
78 cellVolumes_[nIgnCells] = vols[celli];
86 Pout<<
"Found ignition cells:" <<
endl << cells_ <<
endl;
97 const_cast<ignitionSite&
>(*this).findIgnitionCells(mesh_);
107 scalar curTime = db_.
value();
112 (curTime - deltaT >= time_)
114 (curTime - deltaT < time_ +
max(duration_, deltaT) + SMALL)
121 scalar curTime = db_.
value();
124 return(curTime - deltaT >= time_);
132 location_ = is.location_;
133 diameter_ = is.diameter_;
135 duration_ = is.duration_;
136 strength_ = is.strength_;
138 cellVolumes_ = is.cellVolumes_;