34 void Foam::ignitionSite::findIgnitionCells(
const fvMesh&
mesh)
47 scalar radius = diameter_/2.0;
53 cellVolumes_[0] = vols[ignCell];
55 scalar minDist = GREAT;
60 scalar dist =
mag(centres[celli] - location_);
67 if (dist < radius && celli != ignCell)
70 cellVolumes_.
setSize(nIgnCells+1);
72 cells_[nIgnCells] = celli;
73 cellVolumes_[nIgnCells] = vols[celli];
81 Pout<<
"Found ignition cells:" <<
endl << cells_ <<
endl;
90 if (mesh_.changing() && timeIndex_ != db_.timeIndex())
92 const_cast<ignitionSite&
>(*this).findIgnitionCells(mesh_);
94 timeIndex_ = db_.timeIndex();
102 scalar curTime = db_.value();
103 scalar deltaT = db_.deltaTValue();
107 (curTime - deltaT >= time_)
109 (curTime - deltaT < time_ +
max(duration_, deltaT) + SMALL)
116 scalar curTime = db_.value();
117 scalar deltaT = db_.deltaTValue();
119 return(curTime - deltaT >= time_);
127 location_ = is.location_;
128 diameter_ = is.diameter_;
130 duration_ = is.duration_;
131 strength_ = is.strength_;
133 cellVolumes_ = is.cellVolumes_;