CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Isotropic Turbulence IC / Rogallo-procedure

Register Blogs Community New Posts Updated Threads Search

Like Tree9Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 20, 2013, 12:55
Default
  #41
New Member
 
Join Date: Jun 2012
Posts: 14
Rep Power: 14
bigeddy is on a distinguished road
Hi LNA,

Just dropping a quick question, when you had "small strange oscillations at the large wave numbers" i.e. that little weird peak, do you have data readily available to plot with that version of your code, I was hoping you could obtain this function just so I can confirm something im not seeing in my DNS, the function is;

E(kappa)*kappa^(5/3)/epsilon^(2/3)=C_kol

i.e. a plot of C_kol vs kappa*eta with logarithmic x and y axis, Its a compensated kolmogorov energy spectrum function

If you could upload a plot that would be awesome!!
Thanks in advance!
bigeddy is offline   Reply With Quote

Old   October 16, 2013, 08:49
Default
  #42
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Hello

Thanks for telling us about the correction. I wanted to know how could one use the same code (for 3D) to generate a 2D velocity field.

Thank you.

Regards
anzillo is offline   Reply With Quote

Old   October 27, 2013, 17:43
Default
  #43
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 557
Rep Power: 20
cfdnewbie is on a distinguished road
I'm not sure I understand your question...are you asking about the exact solution to the 2D problem?
cfdnewbie is offline   Reply With Quote

Old   October 28, 2013, 05:02
Default
  #44
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Hello !

Thank you for your response. What I meant was if Rogallo's procedure could be modified to generate a 2D velocity field that is homogeneous and isotropic, for a given energy spectrum.

I managed to do that with a little modification. However, the velocity field is not divergence free so to speak, the divergence is about 10^(-3) which is okay with regard to the solution's accuracy. But I would be more content if I could get the divergence to say 10^(-5).

Regards
Dhruv
anzillo is offline   Reply With Quote

Old   December 14, 2013, 22:05
Default
  #45
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Hello Everyone

Does anyone know why does the velocity field lose its solenoidal character after being converted to physical space from wavenumber space by using an IFFT?

I used Rogallo's procedure to generate the field (took care of the 'minus' sign) and made sure that the field matrix is conjugate symmetric before using the IFFT (all done in MATLAB).

I would be grateful to receive some help.

Thank you !

Regards
Dhruv
anzillo is offline   Reply With Quote

Old   December 15, 2013, 08:08
Default
  #46
Senior Member
 
cfdnewbie
Join Date: Mar 2010
Posts: 557
Rep Power: 20
cfdnewbie is on a distinguished road
It all depends on how you compute the divergence in physical space. If you are not using a spectral method to do that (I am assuming that you are not), then you won't get a discrete divergence free condition. The way you are computing derivatives / divergence in both spaces (wave and physical) is not compatible, so you cannot expect symmetries etc. to hold.
FMDenaro likes this.
cfdnewbie is offline   Reply With Quote

Old   December 15, 2013, 09:22
Default
  #47
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Thank you for your prompt reply.

So, one can never convert a the spectral field generated with Rogallo's code into a divergence free velocity field in physical space?

That is why most authors simulate turbulence with spectral approach and mention that the field is divergence free (in spectral space)?
anzillo is offline   Reply With Quote

Old   December 15, 2013, 09:45
Default
  #48
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by anzillo View Post
Thank you for your prompt reply.

So, one can never convert a the spectral field generated with Rogallo's code into a divergence free velocity field in physical space?

That is why most authors simulate turbulence with spectral approach and mention that the field is divergence free (in spectral space)?

as previously said, you can compute in physical space the continuos velocity field using the same spectral reconstrution that derive from your computation
FMDenaro is offline   Reply With Quote

Old   December 15, 2013, 09:54
Default
  #49
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
as previously said, you can compute in physical space the continuos velocity field using the same spectral reconstrution that derive from your computation
Hello Filippo

Thank you for your reply and apologies for repeating my question again. However, I did not understand the meaning of"

"compute in physical space the continuos velocity field using the same spectral reconstrution that derive from your computation"

I would glad if you could elaborate the meaning of the above sentence. I am not familiar with spectral related things because my code is meant for physical space.

Thank you !
anzillo is offline   Reply With Quote

Old   December 15, 2013, 10:18
Default
  #50
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by anzillo View Post
Hello Filippo

Thank you for your reply and apologies for repeating my question again. However, I did not understand the meaning of"

"compute in physical space the continuos velocity field using the same spectral reconstrution that derive from your computation"

I would glad if you could elaborate the meaning of the above sentence. I am not familiar with spectral related things because my code is meant for physical space.

Thank you !

when you use a spectral code, you solve for the time evolution of the Fourier coefficients but that implies that you have the possibility to compute the velocity field by using the same coefficients in a discrete sum of a finite number of wavenumber components Ui(k)*epx(i*k*x). This will be your "continuous" velocity field in physical space.
FMDenaro is offline   Reply With Quote

Old   December 15, 2013, 11:39
Default
  #51
New Member
 
Dhruv Mehta
Join Date: Jun 2012
Posts: 22
Rep Power: 14
anzillo is on a distinguished road
Okay I got it now. Thank you
anzillo is offline   Reply With Quote

Old   November 24, 2017, 02:18
Default Hello,sorry to bother you
  #52
New Member
 
Join Date: Nov 2017
Posts: 1
Rep Power: 0
H.Zheng is on a distinguished road
Quote:
Originally Posted by cfdnewbie View Post
Hello bigeddy,
thanks for your offer. I did get the code working, but I'm still a little bit confused about one thing:

If you follow Rogallo's procedure and create the fourier coefficients in wave space, they will not be fully isotropic. In the x and y plane they will be, but not in the z direction. If I remember correctly, thats by design (by enforcing div u =0) implicitely, but it really confused me as to why people call it "ICs for isotropic flow" when it isn't in the z-direction.

Anyway, I assumed at the time that was ok (people also call the Taylor Green vortex isotropic, but it clearly isn't in the z-direction).....

But I'm currently looking into HIT again, but now for compressible HIT... do you know what's the standard procedure to generate compressible ICs? something like "compressible Rogallo"? Any references / paper would be greatly appreciated!


Thanks again!
could you give me the Rogallo code that initialize isotropic turbulence, I have not found it on other place, thank you!
H.Zheng is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
help on isotropic turbulence initial condition Hall Main CFD Forum 8 April 28, 2011 17:55
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 10:52
Help needed on homogeneous isotropic turbulence an Guoping Xia Main CFD Forum 0 March 12, 2006 22:54
Isotropic Turbulence gorka Main CFD Forum 2 June 30, 2003 16:41


All times are GMT -4. The time now is 16:12.