[banji]

Hi Foamers,

Please, is it possible to use the groovyBC to specify an open boundary condition. I am trying to simulate a die swell problem and need to specify a traction free surface on the lateral boundaries.

Thanks in anticipation.

[gschaider]

Originally Posted by banji

Hi Foamers,

Please, is it possible to use the groovyBC to specify an open boundary condition. I am trying to simulate a die swell problem and need to specify a traction free surface on the lateral boundaries.

Thanks in anticipation.

You should be a BIT more specific about the mathematical formulation of "traction free", "lateral boundaries" etc. Never expect people to put more time into an answer than you obviously took to formulate the question

[banji]

Originally Posted by gschaider

You should be a BIT more specific about the mathematical formulation of "traction free", "lateral boundaries" etc. Never expect people to put more time into an answer than you obviously took to formulate the question

Thanks Bernhard for the observation and reply.

Attached is a rough picture of the computational domain and some comments. In particular, what i would like to know is if i can specify the following bc in the attachment with groovybc.

'n' and 'm' denote the unit normal and tangent vectors to the surface and \tau is the total stress tensor.

Thanks.

[gschaider]

Originally Posted by banji

Thanks Bernhard for the observation and reply.

Attached is a rough picture of the computational domain and some comments. In particular, what i would like to know is if i can specify the following bc in the attachment with groovybc.

'n' and 'm' denote the unit normal and tangent vectors to the surface and \tau is the total stress tensor.

Thanks.

OK. We're getting nearer. It's a multiphase problem. That wasn't really clear from your initial question. Which solver are you going to use?

Still the sketch confuses me: are you talking about the boundary condition where you've written "free surface" or of the BC on the free surface between liquid and air (the bulb in your sketch)?

Anyway: BC for the stress tensor can usually not be directly enforced but by enforcing a BC on the other fields that ensures that the stress is 0 (something like zeroGradient for U or so)

[banji]

Originally Posted by gschaider

OK. We're getting nearer. It's a multiphase problem. That wasn't really clear from your initial question. Which solver are you going to use?

Still the sketch confuses me: are you talking about the boundary condition where you've written "free surface" or of the BC on the free surface between liquid and air (the bulb in your sketch)?

Anyway: BC for the stress tensor can usually not be directly enforced but by enforcing a BC on the other fields that ensures that the stress is 0 (something like zeroGradient for U or so)

I am using interFoam and neglecting surface tension effects. It's the free surface between liquid and air.

Is it possible to set a pde as a BC (in cases when it is difficult to obtain analytical solution of the variable on the remote field) and then groovyBC solves for the variable to be used as a BC?

[gschaider]

Originally Posted by banji

I am using interFoam and neglecting surface tension effects. It's the free surface between liquid and air.

Is it possible to set a pde as a BC (in cases when it is difficult to obtain analytical solution of the variable on the remote field) and then groovyBC solves for the variable to be used as a BC?

For the free interface between liquid and air you can't set a BC with a VOF-solver. They can only be set on the boundaries of the grid (thus the name): where in your sketch you wrote "free surface". Your nomenclature confuses me more and more

I'm not sure what you mean with "set a PDE as a BC". But generally I'd say "no". Usually the BC can be expressed in terms of values and gradients

[banji]

Originally Posted by gschaider

For the free interface between liquid and air you can't set a BC with a VOF-solver. They can only be set on the boundaries of the grid (thus the name): where in your sketch you wrote "free surface". Your nomenclature confuses me more and more

I'm not sure what you mean with "set a PDE as a BC". But generally I'd say "no". Usually the BC can be expressed in terms of values and gradients

You are correct, Bernhard. Maybe the snapshots attached might make my problem clearer.

The flow domain -

Inlet and Outlet bc -
I had no problem with this one.

free surface stress condition -

simplified form of the stress condition -

Here is a link to the paper .

Thanks for the support for the fine utility.

[gschaider]

Originally Posted by banji

You are correct, Bernhard. Maybe the snapshots attached might make my problem clearer.

The flow domain -

Inlet and Outlet bc -
I had no problem with this one.

free surface stress condition -

simplified form of the stress condition -

Here is a link to the paper .

Thanks for the support for the fine utility.

3.3 is in interFoam not a boundary condition but modelled on the interior. You can't influence that in it.

[banji]

Thank you!

[barati]

Hi fomers
i want to define U Field and T field in two rooms with BuoyantBoussinesqSimpleFoam solver when immplementation the BC and run model the error into log file has plesae Help me
error is :
Reading field T

--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for T on ceiling therefore using 160{0}
Reading field p_rgh

--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for p_rgh on ceiling therefore using 160{0}
--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for p_rgh on TwallR therefore using 96{0}
Reading field U

--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for U on ceiling therefore using 160{(0 0 0)}
--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for U on TwallR therefore using 96{(0 0 0)}
Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian

[barati]

hi fomers
i want to attach file which immplemantaion into groovyBC
why error and not to do this model please help me

[nimasam]

you should add keyword value to your syntax

type groovyBC;
valueExpression "vector(0,pos().x > 8 && pos().x < 10.2 && pos().y ==3.27 ? -5 : 0, 0)";
gradientExpression "vector(0, 0, 0)";// consider the rest of celing insulated
fractionExpression "1";
value uniform (0 -0.05 0);//should be considered

[barati]

thank a lot for note
please let me to check it and tell you if error to be clean

[S.M.H]

Originally Posted by nimasam

you should add keyword value to your syntax

hi i have same problem

what i must do to fix this problem?

thanks