Обновление

Уточнение ГУ

MemMAPR_winda.vcxproj

@@ -111,7 +111,6 @@
       <SDLCheck>true</SDLCheck>
       <PreprocessorDefinitions>_DEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
       <ConformanceMode>true</ConformanceMode>
-      <AdditionalIncludeDirectories>./include</AdditionalIncludeDirectories>
     </ClCompile>
     <Link>
       <SubSystem>Console</SubSystem>

es10.plt

@@ -1,10 +0,0 @@
-set cbrange [0:100]
-set size ratio 0.8
-unset key
-
-set palette defined (0 0 0 1, 0.25 0 1 1, 0.5 0 1 0, 0.75 1 1 0, 1 1 0 0)
-
-do for [i=0:99]{
-plot 'explicit10.txt' u 1:2:3 index i w points pt 5 palette
-pause 1e-09}
-pause mouse

es5.plt

@@ -1,10 +0,0 @@
-set cbrange [0:100]
-set size ratio 0.8
-unset key
-
-set palette defined (0 0 0 1, 0.25 0 1 1, 0.5 0 1 0, 0.75 1 1 0, 1 1 0 0)
-
-do for [i=0:99]{
-plot 'explicit5.txt' u 1:2:3 index i w points pt 5 palette
-pause 1e-09}
-pause mouse

include/System.hpp

@@ -6,16 +6,18 @@
 class System{
 	Object& _obj;
 	Mesh _mesh;
-	double _a;
+	double _a1;
+	double _a2;
 	double _step;
 public:
-System(Object& obj, double step = 10., double a1 = 1.): _obj(obj), _mesh(obj, step), _a(a1), _step(step) {}
+System(Object& obj, double step = 10., double a1 = 1., double a2 = 1.): _obj(obj), _mesh(obj, step), _a1(a1), _a2(a2), _step(step) {}
 	void DefineBounds(int, int, int, int);
 	std::vector<std::vector<Node*>>& Nodes();
 	std::vector<Node*>& LineX();
 	std::vector<Node*>& LineY();
 	double step() const { return _step; }
-	double a() const { return _a; };
+	double a1() const { return _a1; };
+	double a2() const { return _a2; };
 };
 
 #endif

is10.plt

@@ -1,10 +0,0 @@
-set cbrange [0:100]
-set size ratio 0.8
-unset key
-
-set palette defined (0 0 0 1, 0.25 0 1 1, 0.5 0 1 0, 0.75 1 1 0, 1 1 0 0)
-
-do for [i=0:99]{
-plot 'implicit10.txt' u 1:2:3 index i w points pt 5 palette
-pause 1e-09}
-pause mouse

is5.plt

@@ -1,10 +0,0 @@
-set cbrange [0:100]
-set size ratio 0.8
-unset key
-
-set palette defined (0 0 0 1, 0.25 0 1 1, 0.5 0 1 0, 0.75 1 1 0, 1 1 0 0)
-
-do for [i=0:99]{
-plot 'implicit5.txt' u 1:2:3 index i w points pt 5 palette
-pause 1e-09}
-pause mouse

src/Node.cpp

@@ -7,31 +7,6 @@ double Node::T() const {
 		return 100.;
 
 	if (_btype == 2) {
-		if (!_left)
-			if (_right)
-				return _right->T() / (1 + Dist(_right));
-		if (!_right)
-			if (_left)
-				return _left->T() / (1 + Dist(_left));
-		if (!_above)
-			if (_below)
-				return _below->T() / (1 + Dist(_below));
-		if (!_below)
-			if (_above)
-				return _above->T() / (1 + Dist(_above));
-		if (_right && _left) {
-			if (_right->IsBound())
-				return _left->T() / (1 + Dist(_left));
-			return _right->T() / (1 + Dist(_right));
-		}
-		if (_above && _below) {
-			if (_above->IsBound())
-				return _below->T() / (1 + Dist(_below));
-			return _above->T() / (1 + Dist(_above));
-		}
-	}
-
-	if (_btype == 3) {
 		if (!_left)
 			if (_right)
 				return _right->T();
@@ -56,6 +31,31 @@ double Node::T() const {
 		}
 	}
 
+	if (_btype == 3) {
+		if (!_left)
+			if (_right)
+				return _right->T() / (1 + Dist(_right));
+		if (!_right)
+			if (_left)
+				return _left->T() / (1 + Dist(_left));
+		if (!_above)
+			if (_below)
+				return _below->T() / (1 + Dist(_below));
+		if (!_below)
+			if (_above)
+				return _above->T() / (1 + Dist(_above));
+		if (_right && _left) {
+			if (_right->IsBound())
+				return _left->T() / (1 + Dist(_left));
+			return _right->T() / (1 + Dist(_right));
+		}
+		if (_above && _below) {
+			if (_above->IsBound())
+				return _below->T() / (1 + Dist(_below));
+			return _above->T() / (1 + Dist(_above));
+		}
+	}
+
 	return _t;
 }
 

src/Solver.cpp

@@ -1,6 +1,6 @@
 #include "Solver.hpp"
-void Solver::SolveImplicit(System& program, double tstop) const {
-	std::ofstream output(filename_impl);
+void Solver::SolveExplicit(System& program, double tstop) const {
+	std::ofstream output(filename_expl);
 	//output << "t x y T" << std::endl;
 
 	for (double t = 0.0; t < tstop; t += delta) {
@@ -50,8 +50,8 @@ void Solver::SolveImplicit(System& program, double tstop) const {
 	}
 }
 
-void Solver::SolveExplicit(System& sys, double tstop) const {
-	std::ofstream output(filename_expl);
+void Solver::SolveImplicit(System& sys, double tstop) const {
+	std::ofstream output(filename_impl);
 	//output << "t x y T" << std::endl;
 
 	for (double t = 0.; t < tstop; t += delta) {
@@ -66,7 +66,7 @@ void Solver::SolveExplicit(System& sys, double tstop) const {
 
 					double tx = (node->r()->T() - 2 * node->T() + node->l()->T()) / pow(sys.step(), 2);
 					double ty = (node->u()->T() - 2 * node->T() + node->d()->T()) / pow(sys.step(), 2);
-					double t = delta * (tx + ty) * sys.a() + node->T();
+					double t = delta * (tx + ty) + node->T();
 					node->SetT(t);
 				}
 			}
@@ -85,17 +85,17 @@ void Solver::SolveLine(System& sys, std::vector<Node*>& n) const {
 	/* Защита от нуля */
 	if (mu2 == 0.)
 		mu2 = .1;
-	double val2 = -(2 * sys.a()) / (pow(sys.step(), 2)) - 1 / delta;
-	double val1 = sys.a() / (pow(sys.step(), 2));
+	double val2 = -(2 * sys.a1()) / (pow(sys.step(), 2)) - 1 / delta;
+	double val1 = sys.a1() / (pow(sys.step(), 2));
 	std::vector<std::vector<double>> _Temperature(size);
 	std::vector<double> right(size);
 	for (int i = 0; i < _Temperature.size(); i++)
 		_Temperature[i].resize(3, 0.0);
 
-	_Temperature[0][0] = -(2 * sys.a()) / (mu1 * pow(sys.step(), 2)) - 1 / delta; /* Первый узел по X */
-	_Temperature[0][1] = (2 * sys.a()) / ((mu1 + 1) * pow(sys.step(), 2)); /* Первый узел по Y */
-	_Temperature.back()[1] = (2 * sys.a()) / ((mu2 + 1) * pow(sys.step(), 2));
-	_Temperature.back()[2] = -(2 * sys.a()) / (mu2 * pow(sys.step(), 2)) - 1 / delta;
+	_Temperature[0][0] = -(2 * sys.a1()) / (mu1 * pow(sys.step(), 2)) - 1 / delta; /* Первый узел по X */
+	_Temperature[0][1] = (2 * sys.a1()) / ((mu1 + 1) * pow(sys.step(), 2)); /* Первый узел по Y */
+	_Temperature.back()[1] = (2 * sys.a1()) / ((mu2 + 1) * pow(sys.step(), 2));
+	_Temperature.back()[2] = -(2 * sys.a1()) / (mu2 * pow(sys.step(), 2)) - 1 / delta;
 
 	for (int i = 1; i < size - 1; i++) {
 		_Temperature[i][0] = val1;
@@ -105,8 +105,8 @@ void Solver::SolveLine(System& sys, std::vector<Node*>& n) const {
 
 	for (int i = 0; i < right.size(); i++)
 		right[i] = -n[i + 1]->T() / delta;
-	right.front() += -(2 * sys.a() * n.front()->T()) / (mu1 * (mu1 + 1) * pow(sys.step(), 2));
-	right.back() += -(2 * sys.a() * n.back()->T()) / (mu2 * (mu2 + 1) * pow(sys.step(), 2));
+	right.front() += -(2 * sys.a1() * n.front()->T()) / (mu1 * (mu1 + 1) * pow(sys.step(), 2));
+	right.back() += -(2 * sys.a1() * n.back()->T()) / (mu2 * (mu2 + 1) * pow(sys.step(), 2));
 	std::vector<double> tmps = ThomasMethod(_Temperature, right);
 	for (int i = 0; i < tmps.size(); i++)
 		n[i + 1]->SetT(tmps[i]);

src/main.cpp

@@ -17,7 +17,7 @@
 #define HOLE_Y 255.
 #define HOLE_RADIUS 50.
 
-#define CONDUCTIVITY 90. // Теплопроводность материала
+#define CONDUCTIVITY 50. // Теплопроводность материала
 
 void visualize(std::ofstream& file, std::string filename, int time_end) {
 	file << "set cbrange [" << 0 << ":" << 100 << "]" << std::endl;
@@ -36,12 +36,12 @@ int main()
 	3 - конвекция
 	4 - отсутствует
 	*/
-	int left = 1;
-	int top = 3;
-	int right = 1;
-	int bottom = 3;
+	int left = 3;
+	int top = 1;
+	int right = 3;
+	int bottom = 1;
 	int arc_bound = 3;
-	int hole_bound = 2;
+	int hole_bound = 1;
 
 	double step_5 = 5;
 	double step_10 = 10;
@@ -83,19 +83,19 @@ int main()
 	slv10.SolveExplicit(explicit10, time_end);
 	slv10.SolveImplicit(implicit10, time_end);
 
-	std::ofstream script("es10.plt");
+	std::ofstream script("ixplicit10.plt");
 	visualize(script, "explicit10.txt", time_end);
 	script.close();
 
-	script.open("is10.plt");
+	script.open("implicit10.plt");
 	visualize(script, "implicit10.txt", time_end);
 	script.close();
 
-	script.open("es5.plt");
+	script.open("explicit5.plt");
 	visualize(script, "explicit5.txt", time_end);
 	script.close();
 
-	script.open("is5.plt");
+	script.open("explicit5.plt");
 	visualize(script, "implicit5.txt", time_end);
 	script.close();
 	return 0;