<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<pre wrap="">Hi Paul,
Thanks for your answer and sorry for having been late in answering in return.
In fact, I wanted to define a usrdiv based on the temperature field for the perturbation topic that I opened some time ago (see <a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2013-October/002334.html">https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2013-October/002334.html</a>) but I now understand that my ideas were wrong. It is not possible to do what I wanted with usrdiv because the problem is fully coupled. I rewrite the set of perturbations equations that I would like to solve below. "k" is the wavelength of the perturbation field (vxp,vyp,tp) in the z direction, and "tp" is used to model the third velocity component "vzp" which is not computed in 2D. Thus the perturbation velocity field has 3 components but (VX,VY) is the baseflow and is 2D.
<i>continuity : d(vpx)/dx + d(vpy)/dy + k * tp=0
</i>
<i>momentum : d(vpx)/dt + VX dvpx/dx + VY dup/dy = - d(prp)/dx - vpx dVX/dx - vxp dVX/dy + nu (</i><i>d^2vxp/dx^2 + d^2vxp/dy^2 - k^2 * vxp )
</i>
<i> d(vpy)/dt + VX dvpy/dx + VY dvpy/dy = - d(prp)/dy - vpx dVY/dx - vyp dVY/dy + nu (</i><i>d^2(vyp)/dx^2 + d^2(vyp)/dy^2 - k^2 * vpy )
</i>
<i> d(tp)/dt + VX d(tp)/dx + VY d(tp)/dy = -k prp + nu ( d^2(tp)/dx^2 + d^2(tp)/dy^2 - k^2 * tp )</i>
in which "d" means partial derivative.
In my opinion the only solution would be to include the third velocity component in the 2D calculations of the perturbation field in nek which would be useful for any situation in which the baseflow is 2D but the perturbation velocity field is 3D with a harmonic evolution along the normal to the plane direction.
Best regards,
Vincent.
Message: 1
Date: Thu, 21 Nov 2013 13:27:48 -0600 (CST)
From: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
To: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
Subject: Re: [Nek5000-users] PN field to PN-2 field
Message-ID:
<a class="moz-txt-link-rfc2396E" href="mailto:mailman.6318.1385062069.3278.nek5000-users@lists.mcs.anl.gov"><mailman.6318.1385062069.3278.nek5000-users@lists.mcs.anl.gov></a>
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
Hi Vincent,
It looks like the routine map12 will do this.
It's to be called on an element by element basis --- thus:
integer e
do e=1,nelt
call map12(b(1,1,1,e),a(1,1,1,e))
enddo
would yield the desired result.
May I inquire as to the application? (There are times when
it might be better to use the transpose of the interpolation
operator.)
Paul
On Thu, 21 Nov 2013, <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> wrote:
</pre>
<blockquote type="cite" style="color: #000000;">
<pre wrap="">Hi Nek's,
I would like to copy a PN field to a PN-2 field, how can I do that ?
let's say "a(lx1,ly1,lz1,nelt)" is given and I want to interpolate "a" onto
a field "b(lx2,ly2,lz2,nelt)" and of course I use a PN-PN-2 method, not
splitting.
Vincent.
_______________________________________________
Nek5000-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
</pre>
</blockquote>
<br>
<pre wrap="">ield to PN-2 field</pre>
<br>
</body>
</html>