Generating turbulence w/ forcing functions

JoshS · April 11, 2012, 15:05

Hey all, anyone have experience using forcing functions to generate turbulence? In my compressible LES code, I eventually get negative pressures when I use forcing. Are there any general rules people follow that I need to be aware of?

I've noticed some people use div-free forcing, but that's not a computational possibility with what I'm doing, as it's too expensive.

Fabian82 · April 13, 2012, 04:38

Usually, turbulence is forced by a statistical Ornstein-Uhlenbeck process. This methods allows to create forcing fields from the solenoidal (div free) to the fully compressible regime. As the forcing is usually only applied to low wave number modes in Fourier space, the computational overhead is very low in both cases.

The following references may be helpful to implement such a method:
http://dx.doi.org/10.1016/0045-7930(88)90013-8
http://nbn-resolving.de/urn/resolver...s2004012514889

April 11, 2012, 15:05	Generating turbulence w/ forcing functions	#1
JoshS New Member Josh Strodtbeck Join Date: Apr 2012 Posts: 6 Rep Power: 14	Hey all, anyone have experience using forcing functions to generate turbulence? In my compressible LES code, I eventually get negative pressures when I use forcing. Are there any general rules people follow that I need to be aware of? I've noticed some people use div-free forcing, but that's not a computational possibility with what I'm doing, as it's too expensive.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
High turbulence intensity	wenfengxie	FLUENT	1	April 16, 2010 12:00
question about turbulence model selection and sensitivity	karananand	Main CFD Forum	1	February 26, 2010 04:41
Centrifugal Pump and Turbulence Model	Michiel	CFX	12	January 25, 2010 03:20
turbulence modeling questions	llowen	Main CFD Forum	3	September 11, 1998 04:24

April 13, 2012, 04:38		#2
Fabian82 New Member Join Date: Mar 2012 Posts: 8 Rep Power: 14	Usually, turbulence is forced by a statistical Ornstein-Uhlenbeck process. This methods allows to create forcing fields from the solenoidal (div free) to the fully compressible regime. As the forcing is usually only applied to low wave number modes in Fourier space, the computational overhead is very low in both cases. The following references may be helpful to implement such a method: http://dx.doi.org/10.1016/0045-7930(88)90013-8 http://nbn-resolving.de/urn/resolver...s2004012514889