| pcg |
August 20, 2020 05:12 |
In CHeatSolver.hpp you need to declare the method
Code:
void Source_Residual(CGeometry *geometry,
CSolver **solver_container,
CNumerics **numerics_container,
CConfig *config,
unsigned short iMesh) override;
because that solver does not have it yet. If a solver has a method with this signature it will be automatically called during the spatial integration phase.
The implementation should be placed in CHeatSolver.cpp, something like:
Code:
void CHeatSolver::Source_Residual(CGeometry *geometry,
CSolver **solver_container,
CNumerics **numerics_container,
CConfig *config,
unsigned short iMesh) {
/*--- Current time. ---*/
const auto dt = config->GetDelta_UnstTimeND();
const su2double t = config->GetTimeIter() * dt;
for (auto iPoint = 0ul; iPoint < nPointDomain; iPoint++) {
/*--- Some geometrical properties you may/will need. ---*/
const auto V = geometry->nodes->GetVolume(iPoint);
const auto X = geometry->nodes->GetCoord(iPoint);
//// x = X[0], y = X[1], z = X[2] (z does not exist in 2D) ////
/*--- Temperature at this point. ---*/
const auto T = nodes->GetSolution(iPoint,0);
/*--- Compute the source term. ---*/
const su2double source = ...
/*--- Add the source to the residual. ---*/
LinSysRes[iPoint] += source;
/*--- If your source depends on temperature you can also update the,
* Jacobian matrix. Take care not to break diagonal dominance. ---*/
Jacobian.AddVal2Diag(iPoint, dsource_dT);
}
}
That should get you going, you need to recompile the code after making changes ( https://su2code.github.io/docs_v7/Bu...2-Linux-MacOS/) after building it the first time it's enough to run "./ninja -C build install".
If you need more help, join the Slack channel (follow the link above and click the # on the top bar) and come to one of the developers meetings (Wednesdays at 4pm CET, everyone is welcome and usually there is no agenda, so plenty of time for random questions). |
