Ordering of the sequence of loops over the edges and nodes

paulzhang · March 27, 2016, 21:16

Hello,

This is Paul Zhang. I got a question on page 5-6 when I read the paper "Towards High-Performance Optimizations of the Unstructured Open-Source SU2 Suite", where I am a little confused about when to apply boundary conditions on the boundary node.

As far as I know, should the boundary conditions be imposed in the very beginning instead of step 3, so the convective and diffusive fluxes can be evaluated using that when we approach boundaries?

Thanks,
Paul

fpalacios · March 31, 2016, 00:54

In SU2 we the weak formulation (boundary conditions are imposed via the residual). and we typically add the missing part flux form the exterior to the interior of the domain at the end (after computing internal fluxes).
This subroutine

void CIntegration::Space_Integration(CGeometry *geometry,
CSolver **solver_container,
CNumerics **numerics,
CConfig *config, unsigned short iMesh,
unsigned short iRKStep,
unsigned short RunTime_EqSystem) {

in integration_structure.cpp

could be useful for you.

Best,
Francisco

March 27, 2016, 21:16	Ordering of the sequence of loops over the edges and nodes	#1
paulzhang Member Paul Zhang Join Date: Feb 2011 Posts: 44 Rep Power: 15	Hello, This is Paul Zhang. I got a question on page 5-6 when I read the paper "Towards High-Performance Optimizations of the Unstructured Open-Source SU2 Suite", where I am a little confused about when to apply boundary conditions on the boundary node. As far as I know, should the boundary conditions be imposed in the very beginning instead of step 3, so the convective and diffusive fluxes can be evaluated using that when we approach boundaries? Thanks, Paul

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[ICEM] Assigning different number of nodes to parallel edges	Rohit	ANSYS Meshing & Geometry	3	September 22, 2017 09:16
[snappyHexMesh] No layers in a small gap	bobburnquist	OpenFOAM Meshing & Mesh Conversion	6	August 26, 2015 09:38
[ICEM] Nodes on Edges...	RodriguezFatz	ANSYS Meshing & Geometry	2	November 22, 2012 13:02
[ICEM] Nodes on sides of the edges don't match	Attesz	ANSYS Meshing & Geometry	1	February 10, 2011 20:20
Displaying only edges and not nodes of the grid in contour views	val3312	FLUENT	1	May 26, 2009 02:27

March 31, 2016, 00:54		#2
fpalacios Super Moderator Francisco Palacios Join Date: Jan 2013 Location: Long Beach, CA Posts: 404 Rep Power: 15	In SU2 we the weak formulation (boundary conditions are imposed via the residual). and we typically add the missing part flux form the exterior to the interior of the domain at the end (after computing internal fluxes). This subroutine void CIntegration::Space_Integration(CGeometry geometry, CSolver solver_container, CNumerics numerics, CConfig config, unsigned short iMesh, unsigned short iRKStep, unsigned short RunTime_EqSystem) { in integration_structure.cpp could be useful for you. Best, Francisco