De Vahl Davis natural convection test - Revision history

E6WikiAdmin at 19:54, 22 October 2022

2022-10-22T19:54:09Z

← 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]
	}}		}}

2022-10-22T19:41:50Z

2022-10-22T19:41:39Z

← Older revision		Revision as of 19:41, 22 October 2022
Line 2:		Line 2:

	{{Box-round\|title= Related papers \|		{{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]

			[https://e6.ijs.si/ParallelAndDistributedSystems/publications/1905659.pdf G. Kosec, M. Založnik, B. Šarler, H. Combeau; A meshless approach towards solution of macrosegregation phenomena, Computers, materials & continua : CMC, vol. 22, 2011 ]

			[http://comms.ijs.si/~gkosec/data/papers/27339815.pdf G. Kosec, M. Depolli, A. Rashkovska, R. Trobec; Super linear speedup in a local parallel meshless solution of thermo-fluid problem, Computers & Structures, vol. 133, 2014]

			[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]

	}}		}}

2022-10-22T19:40:39Z

← 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=

2020-10-07T11:16:25Z

Dead link fix

← Older revision		Revision as of 11:16, 7 October 2020
Line 25:		Line 25:

	=Code=		=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==		==Explicit ACM method with CBS looks==

2019-10-28T21:41:35Z

Explicit pressure correction

← Older revision		Revision as of 21:41, 28 October 2019
Line 70:		Line 70:
	<syntaxhighlight lang="c++" line>		<syntaxhighlight lang="c++" line>
	for (int step = 0; step <= O.t_steps; ++step) {		for (int step = 0; step <= O.t_steps; ++step) {

	// Explicit Navier-Stokes computed on whole domain, including boundaries		// Explicit Navier-Stokes computed on whole domain, including boundaries
	// without pressure		// without pressure
Line 83:		Line 82:
	rhs_pressure(N) = 0; // = 0 part of the regularization equation		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: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));		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);		VecXd solution = solver_p.solve(rhs_pressure);
Line 89:		Line 87:
	VecXd P_c = solution.head(N);		VecXd P_c = solution.head(N);
	for ( int i = interior) v_2[c] -= O.dt / O.rho * op.grad(P_c, c);		for ( int i = interior) v_2[c] -= O.dt / O.rho * op.grad(P_c, c);

	v_2[boundary] = 0; // force boundary conditions		v_2[boundary] = 0; // force boundary conditions

	//heat transport		//heat transport
	for (auto c : interior) {		for (auto c : interior) {

2019-10-28T21:41:12Z

Explicit pressure correction

← Older revision		Revision as of 21:41, 28 October 2019
Line 91:		Line 91:

	v_2[boundary] = 0; // force boundary conditions		v_2[boundary] = 0; // force boundary conditions

			//heat transport
			for (auto c : interior) {
			T2[c] = T1[c] + O.dt * O.lam / O.rho / O.c_p * op.lap(T1, c) -
			O.dt * v1[c].transpose() * op.grad(T1, c);
			}
			for (auto c : top) T2[c] = op.neumann(T2, c, vec_t{0, -1}, 0.0);
			for (auto c : bottom) T2[c] = op.neumann(T2, c, vec_t{0, 1}, 0.0);
	}		}
	</syntaxhighlight>		</syntaxhighlight>

2018-03-12T09:52:09Z

Implicit

← Older revision		Revision as of 09:52, 12 March 2018
Line 105:		Line 105:
	op.valuevec(M_velocity, i, 1 / O.dt_i);		op.valuevec(M_velocity, i, 1 / O.dt_i);
	op.lapvec(M_velocity, i, -O.mu / O.rho);		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		for (int i : boundary) op.valuevec(M_velocity, i, 1); //sets velocity to 0

2018-03-12T09:50:54Z

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

← 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]]~~

2018-03-12T09:47:18Z

Full implicit

← Older revision		Revision as of 09:47, 12 March 2018
Line 94:		Line 94:
	</syntaxhighlight>		</syntaxhighlight>

	== ~~Full implicit~~ ==		== Implicit ==
	<syntaxhighlight lang="c++" line>		<syntaxhighlight lang="c++" line>
	for (int step = 0; step <= O.t_steps_i; ++step) {		for (int step = 0; step <= O.t_steps_i; ++step) {