|June 15, 2012, 08:15||
turbulent inlet boundary condition
Join Date: Mar 2009
Location: Bruxelles, Belgium
Posts: 47Rep Power: 8
I try to modify the boundary condition which reads a outside (measurement data) timeVaryingMappedFixedValue in order to superimpose perturbations to it. I'd like to use this BC for LES. Up to now I first added random perturbations normalized by the Reynolds stresses. The results were rather discouraging, the velocity field improved, however the pressure field became completely incoherent. Well, this was a well known phenomenon in incompressible solvers by my professor, so he suggested to do the following trick, before going to more sophisticated coherent synthetic turbulent inlet:
1. compute rotPerturbation = rot(perturbation)
2. perturbation = eps*(rotPerturbation), where eps is a scaling parameter to keep the perturbations in the desired range.
What one gains is that in this way the imposed velocity field stays divergence free (which is the property of incompressible flows, coming from the continuity equations), so pressure won’t go crazy.
The pressure not, me yes. So I took the aforementioned BC and renamed, modified. It was not working, fine it never does for me at first. Now there is only 2 lines stayed uncommented:
1. The one where I access the velocity filed: const volVectorField UField = db().lookupObject<volVectorField>("U"); – no problem with this.
2. The one where I compute the vorticity field: const volVectorField vorticity(fvc::curl(UField)); – if this line is there, the pressure field go crazy. Note, I’m not using the vorticity for anything, the rest of the code stays unchanged. I plotted the inlet velocity field, stays constant, I plotted the overall velocity field, nothing remarkable happens after one time step. If I plot the pressure field, it changes a lot.
Can anyone tell me why?
|boundary condition, incompressible flows, turbulence|
|Thread||Thread Starter||Forum||Replies||Last Post|
|Wind turbine simulation||Saturn||CFX||34||October 16, 2014 05:27|
|External Radiation Boundary Condition (Two sided wall), Grid Interface||CFD XUE||FLUENT||0||July 8, 2010 06:49|
|New topic on same subject - Flow around race car||Tudor Miron||CFX||15||April 2, 2004 06:18|
|Please help with flow around car modelling!||Tudor Miron||CFX||17||March 19, 2004 20:23|