Hi Juan Diego, yes, first you create a random signal in Z with the appropriate spatial characteristics (in our case a low-pass filtered signal, but this could potentially be created in any other way). So the way you suggest should be fine, however, you just need to make sure that you signal does not depend on the mesh, i.e. give smaller scales close to element boundaries for instance. Regarding the second point. Essentially you just want to be able to restart your simulation. If you use the Lagrange-interpolants p as described in our paper, you would need to be able to re-create the same two h(z) you had when you stopped the simulation. That's what I meant. Best regards, Philipp On 2018-02-26 17:20, [email protected] wrote:

Dear Philipp,

Thank you for your response. So if I understand correctly, first you generate a random signal as a function of Z and then use a low-pass filter to obtain h^i(z). Thinking about how to implement this in Nek5000, I guess one could impose random numbers at each node (collocation point) and then apply some kind of filter. Is this correct?

Also, could you please clarify what you mean by: make sure that you have the same random numbers for the whole spanwise extent? Since h(z) is a function of z, I'm not sure I understand this statement.

Thank you very much for you help.

Juan Diego Colmenares

PhD Candidate

Department of Mechanical Engineering

University of New Mexico

_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users