[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