From: <jcch@ma...>  20060123 15:50:59

Hi, I am developing a RungeKutta discontinuous Galerkin (RKDG) solver to study compressible flow in turbomachines. By using this great library I was able to code the serial version of the 2D and the 3D solvers in less than an month. To compute steadystate in a turbomachines (3D) with a single row of blades we only need a domain with a single blade with appropriated periodic conditions. These periodic condition are imposed by applying a rotation of ( 360 deg / num_blades ) to the conservative variables and conservative fluxes. These periodic boundaries are geometrically equal but have different surface meshes. The periodic boundary conditions are handled quite easily in the serial version of the code because the RKDG only needs the numerical solution from the periodic neighbors at the surface integration points. I coded this by finding which element of the other periodic surface contains the desired surface integration point. The problem is when I try to perform the parallelization of the code. The parallel code that finds the periodic neighbors was almost equal to the serial version because the entire mesh is passed to all the processors. My BIG PROBLEM is how to get the solution at elements that contains the integration points in the periodic surfaces that belong to other processors but are not included the send_list of the processor? Best regards, João Henriques 