Parallelization of an Unsteady ALE Solver with Deforming Mesh Using OpenACC

Algorithm 1

Two-dimensional parallel coefficient matrix assembly (2DPCMA) algorithm.

Input:  () Unstructured grid metrics (UGM) including the coordinates of the mesh nodes

     and the indices of the nodes that form the elements; adjacent data structure of grid

    () Row array, , which stores the starting positions of rows.

    () Inner node identifier and the GPU thread identifier .

Output:

     () Column array, , which stores the positions of the no-zero values in each row.

     () Value array, , which stores the no-zero values in each row.

     () Right hand vector of the linear system of equations, .

Procedure: Thread is responsible for the following tasks

     () initialize an integer counter: ;

     () initialize the elements of vector that is to calculate: ;

     () For each triangle that contains Do

     ()   Calculate torsional coefficients of ;

     ()   ;

     ()   ;

     () End For

     () For each neighbor of node Do

     ()   IF is a moving boundary node Then

     ()    For each triangle that contains edge Do

     ()     ;

     ()    End For

     ()   Else

     ()    ;

     ()    ;

     ()    ;

     ()    For each triangle that contains edge Do

     ()     ;

     ()     ;

     ()    End For

     ()   End IF

     () End For