jakob.schratter/SciFEM-Project_CoffeeMugSimulation_Celebic-Schratter
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

add another 67 degree break condition

Browse source
This commit is contained in:
dino.celebic 2026-01-26 17:24:33 +01:00
commit c39a976ebd
5 changed files with 17 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

5
mgrid_2/main.cpp
View file

@ -83,9 +83,11 @@ int main(int argc, char **argv )
auto t3 = system_clock::now(); // start timer
double average_cup_temperature = u0_mug;
double percentage_temp_reached = 0.0;
double time_count = 0;
while (average_cup_temperature < 67.0)
// while (percentage_temp_reached < 60.0)
//for (int step = 0; step < steps; ++step)
{
vector<double> G(Mdt.Nrows(), 0.0);
@ -100,8 +102,9 @@ int main(int argc, char **argv )
JacobiSolve(SK, H, uv); // solve: (M/dt + K + C) * u_{n+1} = F + M/dt * u_{n}
// ----- SK ----- ------ H -------
average_cup_temperature = mesh_c.AverageVectorFunction_perSubdomain(uv, 0);
tie(average_cup_temperature, percentage_temp_reached) = mesh_c.AverageVectorFunction_perSubdomain(uv, 0);
cout << "Average cup temperature: " << average_cup_temperature << " after " << time_count << " seconds. " << endl;
cout << "% of mug elements reached temperature 67º: " << percentage_temp_reached << endl;
time_count += dt;
}
auto t4 = system_clock::now(); // stop timer