Dear Guillermo,

I think the interpolation on the homogeneous direction is not yet implemented, for that reason Nektar give you the value of the nearest homogeneous plane changing z = 0.075 to z = 0.05.

In theory the implementation could be easy because you have the values of the Fourier coefficients, but I tried for a while without results (I had to focus on another aspects of the code). Probably somebody else have solved this issue.

Please let me know if you can solve it. Or we can sum efforts to solve it.

Cheers,

On Thu, 4 Jul 2019 at 04:43, Guillermo <guillermo180395@gmail.com> wrote:

Dear all,

I am conducting a quasi 3D simulation, being the z-axis the homogeneous
direction, with an homogeneous length equal to 0.1 (LZ=0.1) and up to 3
Fourier modes (HomModesZ = 4).
I am interested on placing some history points along the homogeneous
direction, however I am facing several issues which I would like to clarify.

- First of all, I have made several tests and it seems that it is not
possible to parametrize the points coordinates. I mean, the coordinates
do not accept parameters defined on the parameters XML element. The
following example does not provide any output result:
<PARAM NAME="Points">
         0.0 LZ 0.0
</PARAM>
Is there any way to parametrize the coordinates? since it will be very
beneficial.

- I would like to place probes at LZ/4, LZ/2 and 3*LZ/4, thus I defined
the HistoryPoints filter as follows:
<PARAM NAME="Points">
         0.0 0.0 0.025
         0.0 0.0 0.05
         0.0 0.0 0.075
</PARAM>

However, I obtain the following 'error' message: Reseting History point
from z = 0.075 to z = 0.05. As not Fourier interpolation is required,
why do I get this error?

- At the FilterHistoryPoints.cpp file I see that it is possible to
specify the plane on which to evaluate the history point by using the
"OutputPlane" tag. Since I could not find any example, I would like to
know how to use this tag, more in precise, how to define the
z-coordinate of the points along with the output plane. Is the following
example correct?

<FILTER TYPE="HistoryPoints">
<PARAM NAME="OutputFile">TimeValues</PARAM>
<PARAM NAME="OutputFrequency">10</PARAM>
<PARAM NAME="OutputPlane">0,1,2,3</PARAM>
<PARAM NAME="Points">
     1 0.5 0
     1 0.5 0
     1 0.5 0
     1 0.5 0
</PARAM>
</FILTER>

I expect to obtain the measurements at the coordinates
(1,0.5,0),(1,0.5,0.025),(1,0.5,0.5),(1,0.5,0.075).

- Also at the FilterHistoryPoints.cpp file I have seen that there is
another tag, called "WaveSpace". If I implement it, will I obtain the
Fourier coefficients on these points? It is somehow confusing to me, so
I please you to bring some light into it, and how to use this field.

I am using nektar 4.3.3 version and running in parallel with 8 cores.
Let me know if you require any other information.

With best regards,
Guillermo

_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users