[Nek5000-users] array reshape

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Jul 12 09:57:06 CDT 2018


Hi Philipp,

If I have 

common /avg/    uvag(lx1,ly1,lz1,lelt)


in one subroutine, then can I call all data stored in that array as 


common /avg/ uvag(lx1*ly1*lz1,lelt) 


in another subroutine? Will this subroutine be able to read the data stored in this 

(common /avg/ uvag(lx1*ly1*lz1,lelt) ) array?



> On 12 Jul 2018, at 15:23, nek5000-users-request at lists.mcs.anl.gov wrote:
> 
> Send Nek5000-users mailing list submissions to
> 	nek5000-users at lists.mcs.anl.gov
> 
> To subscribe or unsubscribe via the World Wide Web, visit
> 	https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> or, via email, send a message with subject or body 'help' to
> 	nek5000-users-request at lists.mcs.anl.gov
> 
> You can reach the person managing the list at
> 	nek5000-users-owner at lists.mcs.anl.gov
> 
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Nek5000-users digest..."
> 
> 
> Today's Topics:
> 
>   1. array reshape (nek5000-users at lists.mcs.anl.gov)
>   2. Re: array reshape (nek5000-users at lists.mcs.anl.gov)
>   3. Re: Number of Elements per Characteristic Length for DNS
>      (nek5000-users at lists.mcs.anl.gov)
>   4. Re: Number of Elements per Characteristic Length for DNS
>      (nek5000-users at lists.mcs.anl.gov)
>   5. Re: Number of Elements per Characteristic Length for DNS
>      (nek5000-users at lists.mcs.anl.gov)
> 
> 
> ----------------------------------------------------------------------
> 
> Message: 1
> Date: Thu, 12 Jul 2018 12:07:42 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
> 	<nek5000-users at lists.mcs.anl.gov>
> Subject: [Nek5000-users] array reshape
> Message-ID:
> 	<mailman.6007.1531397285.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="iso-8859-1"
> 
> Hi all,
> 
> 
> How can I reshape this array uvag(lx1,ly1,lz1,lelt) into abc(lx1*ly1*lz1,lelt)?
> 
> 
> 
> The following array reshape which works for fortran didn't work:
> 
> 
> reshape(uavg, (  \  lx1*ly1*lz1 , lelt \ ))
> 
> 
> 
> Any suggestions?
> 
> 
> Thanks.
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/85abe829/attachment-0001.html>
> 
> ------------------------------
> 
> Message: 2
> Date: Thu, 12 Jul 2018 14:13:18 +0200
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] array reshape
> Message-ID:
> 	<mailman.6009.1531397608.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset=windows-1252; format=flowed
> 
> I dont think you need to do anything as the first index is the faster 
> varying in Fortran.
> Philipp
> 
> On 2018-07-12 14:07, nek5000-users at lists.mcs.anl.gov wrote:
>> Hi all,
>> 
>> 
>> How can I reshape this array?uvag(lx1,ly1,lz1,lelt) into 
>> abc(lx1*ly1*lz1,lelt)?
>> 
>> 
>> 
>> The following array reshape which works for fortran didn't work:
>> 
>> 
>> reshape(uavg, (? \ lx1*ly1*lz1?,?lelt?\?))
>> 
>> 
>> 
>> Any suggestions?
>> 
>> 
>> Thanks.
>> 
>> 
>> 
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
> 
> 
> ------------------------------
> 
> Message: 3
> Date: Thu, 12 Jul 2018 08:56:18 -0400
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
> 	Length for DNS
> Message-ID:
> 	<mailman.6015.1531400181.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
> 
> Hi Philipp,
> 
> Thank you for getting back to me. Do you have any references that you could
> provide so that I can look into this further?
> 
> Thanks,
> 
> Tommy
> 
> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
> 
>> Hi,
>> typically, in wall bounded turbulence you would design the mesh based on
>> inner scaled quantities (i.e. based on the viscous length scale l*). This
>> is what we have done in the paper by El Khoury, in accordance to studies in
>> channel flow. TOwards the centre, one can then also use Kolmogorov units
>> (not sure whether Taylor units are so useful). I am quite convinced that
>> the El Khoury simulations are well-resolved DNS, which can also be seen
>> when looking at spectra.
>> 
>> Philipp
>> 
>> 
>> On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov wrote:
>> 
>>> Hello All,
>>> 
>>> I am currently investigating the turbulence in a pipe with periodic
>>> boundary conditions. I am working on building my mesh and was planning on
>>> using Re^(3/4) as the number of elements per characteristic length, but I
>>> have found differing methods in literature detailing the number of elements
>>> used for DNS. One paper in particular (G.K. El Khoury et al, Direct
>>> Numerical Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>> Numbers,  Flow Turbulence Combust, 2013) uses Nek5000 to simulate turbulent
>>> flow in a pipe at different Reynolds numbers. The number of elements
>>> calculated in this paper seem to be in between the LES and DNS suggested
>>> number of elements (element size between the Taylor microscale and the
>>> Kolmogorov scale). Turbulent Flows by Pope however suggests a number of
>>> elements that is slightly larger than Re^(3/4), and is dependent on the
>>> Taylor microscale Reynolds number.
>>> 
>>> Any suggestions on the correct way to determine the number of spectral
>>> elements in a characteristic length would be appreciated.
>>> 
>>> Regards,
>>> 
>>> Tommy
>>> 
>>> 
>>> 
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>> 
>>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/3ca13d06/attachment-0001.html>
> 
> ------------------------------
> 
> Message: 4
> Date: Thu, 12 Jul 2018 15:07:49 +0200
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
> 	Length for DNS
> Message-ID:
> 	<mailman.6017.1531400891.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset=utf-8; format=flowed
> 
> I guess the El Khoury paper might be a good starting point; we really 
> tried to make sure that the resolution is sufficient for good DNS. 
> Otherwise, there is a resolution discussion for spectral methods in the 
> paper by Coleman, Johnstone and Spalart:
> 
> https://aip.scitation.org/doi/abs/10.1063/1.3247176
> 
> Philipp
> 
> On 2018-07-12 14:56, nek5000-users at lists.mcs.anl.gov wrote:
>> Hi Philipp,
>> 
>> Thank you for getting back to me. Do you have any references that you 
>> could provide so that I can look into this further?
>> 
>> Thanks,
>> 
>> Tommy
>> 
>> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov 
>> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
>> 
>>    Hi,
>>    typically, in wall bounded turbulence you would design the mesh
>>    based on inner scaled quantities (i.e. based on the viscous length
>>    scale l*). This is what we have done in the paper by El Khoury, in
>>    accordance to studies in channel flow. TOwards the centre, one can
>>    then also use Kolmogorov units (not sure whether Taylor units are so
>>    useful). I am quite convinced that the El Khoury simulations are
>>    well-resolved DNS, which can also be seen when looking at spectra.
>> 
>>    Philipp
>> 
>> 
>>    On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov
>>    <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>> 
>>        Hello All,
>> 
>>        I am currently investigating the turbulence in a pipe with
>>        periodic boundary conditions. I am working on building my mesh
>>        and was planning on using Re^(3/4) as the number of elements per
>>        characteristic length, but I have found differing methods in
>>        literature detailing the number of elements used for DNS. One
>>        paper in particular (G.K. El Khoury et al, Direct Numerical
>>        Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>        Numbers,? Flow Turbulence Combust, 2013) uses Nek5000 to
>>        simulate turbulent flow in a pipe at different Reynolds numbers.
>>        The number of elements calculated in this paper seem to be in
>>        between the LES and DNS suggested number of elements (element
>>        size between the Taylor microscale and the Kolmogorov scale).
>>        Turbulent Flows by Pope however suggests a number of elements
>>        that is slightly larger than Re^(3/4), and is dependent on the
>>        Taylor microscale Reynolds number.
>> 
>>        Any suggestions on the correct way to determine the number of
>>        spectral elements in a characteristic length would be appreciated.
>> 
>>        Regards,
>> 
>>        Tommy
>> 
>> 
>> 
>>        _______________________________________________
>>        Nek5000-users mailing list
>>        Nek5000-users at lists.mcs.anl.gov
>>        <mailto:Nek5000-users at lists.mcs.anl.gov>
>>        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>        <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>> 
>>    _______________________________________________
>>    Nek5000-users mailing list
>>    Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov>
>>    https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>    <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>> 
>> 
>> 
>> 
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
> 
> 
> ------------------------------
> 
> Message: 5
> Date: Thu, 12 Jul 2018 10:22:15 -0400
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] Number of Elements per Characteristic
> 	Length for DNS
> Message-ID:
> 	<mailman.6023.1531405402.86639.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
> 
> Hello Philipp,
> 
> Thanks again for getting back to me and for sending these references.
> 
> Regards,
> 
> Tommy
> 
> On Thu, Jul 12, 2018 at 9:07 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
> 
>> I guess the El Khoury paper might be a good starting point; we really
>> tried to make sure that the resolution is sufficient for good DNS.
>> Otherwise, there is a resolution discussion for spectral methods in the
>> paper by Coleman, Johnstone and Spalart:
>> 
>> https://aip.scitation.org/doi/abs/10.1063/1.3247176
>> 
>> Philipp
>> 
>> On 2018-07-12 14:56, nek5000-users at lists.mcs.anl.gov wrote:
>> 
>>> Hi Philipp,
>>> 
>>> Thank you for getting back to me. Do you have any references that you
>>> could provide so that I can look into this further?
>>> 
>>> Thanks,
>>> 
>>> Tommy
>>> 
>>> On Thu, Jul 12, 2018 at 4:25 AM, <nek5000-users at lists.mcs.anl.gov
>>> <mailto:nek5000-users at lists.mcs.anl.gov>> wrote:
>>> 
>>>    Hi,
>>>    typically, in wall bounded turbulence you would design the mesh
>>>    based on inner scaled quantities (i.e. based on the viscous length
>>>    scale l*). This is what we have done in the paper by El Khoury, in
>>>    accordance to studies in channel flow. TOwards the centre, one can
>>>    then also use Kolmogorov units (not sure whether Taylor units are so
>>>    useful). I am quite convinced that the El Khoury simulations are
>>>    well-resolved DNS, which can also be seen when looking at spectra.
>>> 
>>>    Philipp
>>> 
>>> 
>>>    On 2018-07-11 20:58, nek5000-users at lists.mcs.anl.gov
>>>    <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>>> 
>>>        Hello All,
>>> 
>>>        I am currently investigating the turbulence in a pipe with
>>>        periodic boundary conditions. I am working on building my mesh
>>>        and was planning on using Re^(3/4) as the number of elements per
>>>        characteristic length, but I have found differing methods in
>>>        literature detailing the number of elements used for DNS. One
>>>        paper in particular (G.K. El Khoury et al, Direct Numerical
>>>        Simulation of Turbulent Pipe Flow at Moderately High Reynolds
>>>        Numbers,  Flow Turbulence Combust, 2013) uses Nek5000 to
>>>        simulate turbulent flow in a pipe at different Reynolds numbers.
>>>        The number of elements calculated in this paper seem to be in
>>>        between the LES and DNS suggested number of elements (element
>>>        size between the Taylor microscale and the Kolmogorov scale).
>>>        Turbulent Flows by Pope however suggests a number of elements
>>>        that is slightly larger than Re^(3/4), and is dependent on the
>>>        Taylor microscale Reynolds number.
>>> 
>>>        Any suggestions on the correct way to determine the number of
>>>        spectral elements in a characteristic length would be appreciated.
>>> 
>>>        Regards,
>>> 
>>>        Tommy
>>> 
>>> 
>>> 
>>>        _______________________________________________
>>>        Nek5000-users mailing list
>>>        Nek5000-users at lists.mcs.anl.gov
>>>        <mailto:Nek5000-users at lists.mcs.anl.gov>
>>>        https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>        <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>> 
>>>    _______________________________________________
>>>    Nek5000-users mailing list
>>>    Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mc
>>> s.anl.gov>
>>>    https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>    <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>> 
>>> 
>>> 
>>> 
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>> 
>>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180712/b5066fc5/attachment.html>
> 
> ------------------------------
> 
> Subject: Digest Footer
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> 
> 
> ------------------------------
> 
> End of Nek5000-users Digest, Vol 113, Issue 11
> **********************************************



More information about the Nek5000-users mailing list