CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

perturbU in twoPhaseEulerFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   April 28, 2011, 11:18
Default perturbU in twoPhaseEulerFoam
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 78
Rep Power: 6
cheng1988sjtu is on a distinguished road
Hi Foamers,

I want to check that the twoPhaseEulerFoam is reliable, and one of the test case is that I set the domain as one fluid (i.e. the beta=1 everywhere).

I modified the code to include pressure gradient which can give a designed Re_tau, and If I'm not wrong, when I turn on the k-epsilon model, the wall-function is automatically on. I set the y_plus=50 at the first grid.

OK, Maybe too much information. Now the problem is that when I want to use perturbU to introduce some perturbation, I get stuck. The following is the Perturb.C file. perturbU.C

I think the problem is Ub.write(), since if I let it print out Ub[celli], it goes well, however, I just can't have the new Ub file in 0/, so what should I do to let the Ub.write() work?

Any hint will be highly appreciated!
cheng1988sjtu is offline   Reply With Quote


Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Something wrong in UEqns.H within twoPhaseEulerFoam cheng1988sjtu OpenFOAM 2 June 24, 2011 10:48
Stratified flow using twoPhaseEulerFoam karthik1414 OpenFOAM 4 May 4, 2011 05:12
stratified horizontal two phase flow usinfg twoPhaseEulerFoam karthik1414 OpenFOAM 0 April 12, 2011 09:57
problems in Two Phase flow using twoPhaseEulerFoam with OpenFoam 1.6 raagh77 OpenFOAM Running, Solving & CFD 0 March 6, 2010 06:11
TwoPhaseEulerFoam bed tutorial case stable in 1.5, crashes in 1.6 hemph OpenFOAM 3 December 5, 2009 05:19

All times are GMT -4. The time now is 21:26.