De Vahl Davis natural convection test - Revision history

E6WikiAdmin at 19:54, 22 October 2022

2022-10-22T19:54:09Z

2022-10-22T19:41:50Z

2022-10-22T19:41:39Z

2022-10-22T19:40:39Z

2020-10-07T11:16:25Z

Dead link fix

2019-10-28T21:41:35Z

‎Explicit pressure correction

2019-10-28T21:41:12Z

‎Explicit pressure correction

2018-03-12T09:52:09Z

‎Implicit

2018-03-12T09:50:54Z

‎Quick comparison of explicit ACM, explicit pressure correction, and implicit methods

2018-03-12T09:47:18Z

‎Full implicit

← Older revision		Revision as of 19:54, 22 October 2022
Line 10:		Line 10:
	[https://e6.ijs.si/ParallelAndDistributedSystems/publications/3218939.pdf G. Kosec, B. Šarler; Simulation of macrosegregation with mesosegregates in binary metallic casts by a meshless method, Engineering analysis with boundary elements]		[https://e6.ijs.si/ParallelAndDistributedSystems/publications/3218939.pdf G. Kosec, B. Šarler; Simulation of macrosegregation with mesosegregates in binary metallic casts by a meshless method, Engineering analysis with boundary elements]

		+	[https://e6.ijs.si/ParallelAndDistributedSystems/publications/32388135.pdf M. Maksić, V. Djurica, A. Souvent, J. Slak, M. Depolli, G. Kosec; Cooling of overhead power lines due to the natural convection, International journal of electrical power & energy systems, 2019]
	}}		}}

@@ Line 2: / Line 2: @@
 {{Box-round|title= Related papers |
-[https://e6.ijs.si/ParallelAndDistributedSystems/publications/32230439.pdf J. Slak, G. Kosec; Adaptive radial basis function-generated finite differences method for contact problems, International journal for numerical methods in engineering, vol. 119, 2019 [DOI: 10.1002/nme.6067]]
+[https://e6.ijs.si/ParallelAndDistributedSystems/publications/992507.pdf G. Kosec, B. Šarler; Solution of thermo-fluid problems by collocation with local pressure correction, International journal of numerical methods for heat & fluid flow, vol.18, 2008]
 }}

← Older revision		Revision as of 19:40, 22 October 2022
Line 1:		Line 1:
	Click here to return back to [[Fluid Mechanics]]		Click here to return back to [[Fluid Mechanics]]
		+
		+	{{Box-round\|title= Related papers \|
		+	[https://e6.ijs.si/ParallelAndDistributedSystems/publications/32230439.pdf J. Slak, G. Kosec; Adaptive radial basis function-generated finite differences method for contact problems, International journal for numerical methods in engineering, vol. 119, 2019 [DOI: 10.1002/nme.6067]]
		+	}}
		+
		+

	=Intro=		=Intro=

@@ Line 25: / Line 25: @@
 =Code=
-Full examples can be found under the examples in the code repository [[Main Page]]).
+Full examples can be found under the examples in the [https://gitlab.com/e62Lab/medusa code repository].
 ==Explicit ACM method with CBS looks==

@@ Line 70: / Line 70: @@
 <syntaxhighlight lang="c++" line>
      for (int step = 0; step <= O.t_steps; ++step) {
          // Explicit Navier-Stokes computed on whole domain, including boundaries
          // without pressure
@@ Line 83: / Line 82: @@
          rhs_pressure(N) = 0; // = 0 part of the regularization equation
          for (int i:interior) rhs_pressure(c) = O.rho / O.dt * op.div(v_2, c);
          for (int i: boundary) rhs_pressure(c) = O.rho / O.dt * v_2[c].dot(domain.normal(c));
          VecXd solution = solver_p.solve(rhs_pressure);
@@ Line 89: / Line 87: @@
          VecXd P_c = solution.head(N);
          for ( int i = interior) v_2[c] -=  O.dt / O.rho * op.grad(P_c, c);
          v_2[boundary] = 0; // force boundary conditions
          //heat transport
          for (auto c : interior) {

@@ Line 91: / Line 91: @@
          v_2[boundary] = 0; // force boundary conditions
+}
 </syntaxhighlight>

@@ Line 105: / Line 105: @@
              op.valuevec(M_velocity, i, 1 / O.dt_i);
              op.lapvec(M_velocity, i, -O.mu / O.rho);
-             // op.gradvec(M_velocity, i, v_1[i]);
+             op.gradvec(M_velocity, i, v_1[i]);
+         }
          for (int i : boundary) op.valuevec(M_velocity, i, 1); //sets velocity to 0

← Older revision		Revision as of 09:50, 12 March 2018
Line 212:		Line 212:
	File:DVD_vMax_convRa8.png		File:DVD_vMax_convRa8.png
	</gallery>		</gallery>
−
−	~~== Quick comparison of explicit ACM, explicit pressure correction, and implicit methods ==~~
−
−	~~'''This section is not completed yet!'''~~
−
−	In below images you can see quick comparison of explicit stepping with ACM p-v coupling (EACM) [row 1 col 1], explicit stepping with pressure correction p-v coupling (EPC) [row 1 col 2], and fully implicit stepping with pressure correction p-v coupling (IPV) with 1000, 100, 10, 1 times bigger time step [rows 2-3] in comaprison to the EACM adn EPC. EACM and EPC agree well, also the implicit solution converges towards the EACM and EPC but only with the same time step as in explicit approaches.
−
−	~~[[File:ACM.png\|400px]] [[File:pp.png\|400px]]~~
−
−	~~[[File:implicit1.png\|400px]] [[File:implicit2.png\|400px]]~~
−
−	~~[[File:implicit3.png\|400px]] [[File:implicit5.png\|400px]]~~

@@ Line 94: / Line 94: @@
 </syntaxhighlight>
-== Full implicit ==
+== Implicit ==
 <syntaxhighlight lang="c++" line>
 for (int step = 0; step <= O.t_steps_i; ++step) {