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

Interfoam... free surface simulation urgent

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 14, 2009, 14:05
Unhappy Interfoam... free surface simulation urgent
  #1
New Member
 
lostin
Join Date: Jul 2009
Location: India
Posts: 12
Rep Power: 7
lostin4ever is on a distinguished road
Hi all..

I am trying to simulate a free motion of water in a rectangular tank. I have taken water(second phase) in some initial distorted form and now I want to simulate its free motion in tank. I am using interfoam for this.

After running interFoam I am getting this error...
Build : 1.6-53b7f692aa41
Exec : interFoam
Date : Nov 15 2009
Time : 00:08:29
Host : lostin
PID : 19839
Case : /home/akash/OpenFOAM/akash-1.6/run/multiphase/interFoam/laminar/8cm_unforced
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading field p

Reading field alpha1

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 7.12628e-12, No Iterations 1
time step continuity errors : sum local = 3.34444e-09, global = -1.69132e-27, cumulative = -1.69132e-27
Courant Number mean: 2.51417e-09 max: 4.04544e-07

Starting time loop

Courant Number mean: 5.02834e-08 max: 8.09087e-06
deltaT = 0.1
Time = 0.1

MULES: Solving for alpha1
Liquid phase volume fraction = 0.2 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.2 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.2 Min(alpha1) = -1.16775e-40 Max(alpha1) = 1


mesh flux field does not exists, is the mesh actually moving?

From function fvMesh:hi()
in file fvMesh/fvMeshGeometry.C at line 336.

FOAM exiting


I am searched a lot but could not find solution of this error. I am using 1.6 version.


Thanks in advance..
lostin4ever is offline   Reply With Quote

Old   November 15, 2009, 06:46
Default
  #2
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 264
Rep Power: 8
phsieh2005 is on a distinguished road
Hi,

If you send me your case, I can take a look at it.

Pei
phsieh2005@yahoo.com
phsieh2005 is offline   Reply With Quote

Old   November 15, 2009, 08:50
Default
  #3
fab
New Member
 
Join Date: Jun 2009
Posts: 13
Rep Power: 7
fab is on a distinguished road
hi guys, i have a problem with the surface tension implementation!
like the tutorial:free surface flow over a bump.
i have a film(polymersolution) over a platt with 2 diferent temperatures,length scale 50 mm,laminar other settings and the boundarys condition are the same like in the tutorials.
i have to get the surfacetension-curve as funktion of the temperature!
what i should to do? thank a lot
fab is offline   Reply With Quote

Old   November 20, 2009, 07:33
Default
  #4
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 264
Rep Power: 8
phsieh2005 is on a distinguished road
Hi, lostin4ever,

Since you are using interFoam, so, I assume that you are not running moving mesh.

In your 0/U, change wall type to fixedValue and in constant/dynamicMeshDict change it to staticMeFvMesh

you should be able to run your case after the two changes.

Pei
phsieh2005 is offline   Reply With Quote

Old   October 12, 2010, 08:29
Default compiling error
  #5
Member
 
mohsen kh
Join Date: Nov 2009
Posts: 41
Rep Power: 5
mohsenkh599 is an unknown quantity at this point
Hi every body,
I have a problem when I want to compile a model.
the error is:

************************************************** *********************
************************************************** *********************

viscoelasticLaws/mohsen/mohsen.C:159: error: no match for ‘operator==’ in ‘Foam:perator+(const Foam::tmp<Foam::fvMatrix<Type> >&, const Foam::tmp<Foam::fvMatrix<Type> >&) [with Type = Foam::SymmTensor<double>](((const Foam::tmp<Foam::fvMatrix<Foam::SymmTensor<double> > >&)((const Foam::tmp<Foam::fvMatrix<Foam::SymmTensor<double> > >*)(& Foam::fvm::div(const Foam::surfaceScalarField&, Foam::GeometricField<Type, Foam::fvPatchField, Foam::volMesh>&) [with Type = Foam::SymmTensor<double>](((Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>&)(&((Foam::mohsen*)this)->Foam::mohsen::a_))))))) == Foam:perator+(const Foam::tmp<Foam::GeometricField<TypeR, PatchField, GeoMesh> >&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Type1 = Foam::Tensor<double>, Type2 = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator*(const Foam::dimensioned<double>&, const Foam::tmp<Foam::GeometricField<Type, PatchField, GeoMesh> >&) [with Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator&(const Foam::tmp<Foam::GeometricField<TypeR, PatchField, GeoMesh> >&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Type1 = Foam::SymmTensor<double>, Type2 = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator-(const Foam::dimensioned<Type>&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Form = Foam::SymmTensor<double>, Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >&)((const Foam::tmp<Foam::GeometricField<Foam::SymmTensor<do uble>, Foam::fvPatchField, Foam::volMesh> >*)(& Foam:perator*(const Foam::scalar&, const Foam::GeometricField<Type, PatchField, GeoMesh>&) [with Type = Foam::SymmTensor<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh](((const Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>&)((const Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>*)(&((Foam::mohsen*)this)->Foam::mohsen::a_))))))))))))))))))))’

************************************************** ************
************************************************** ************

and the model is :

// Velocity gradient tensor
volTensorField L = fvc::grad( U() );

// Convected derivate term
// volTensorField C = tau_ & L;

// vorticity term
volTensorField Vor = L - L.T();

// Twice the rate of deformation tensor
volSymmTensorField twoD = twoSymm( L );


// fc
scalar fc = 0.5; //1 - 27 * det( a_ );

//double dot term

// volSymmTensorField mk = -( (2 / 35) * ( 1 - fc ) * twoD - (2 / 7) * ( a_ & twoD + twoD & a_ + (1 / 2) * tr( twoD & a_ ) * I_ ) - fc * tr( twoD & a_ ) * a_ );

// Stress transport equation
tmp<fvSymmTensorMatrix> aEqn
(
fvm::ddt(a_)
+ fvm::div(phi(), a_)
==
( 1 / 2 ) * ( Vor & a_ - a_ & Vor ) + ( 1 / 2 ) * keisi_ * ( twoD & a_ + a_ & twoD + 2 * ( ( (2 / 35) * ( 1 - fc ) * twoD - (2 / 7) * ( 1 - fc ) * ( a_ & twoD + twoD & a_ + (1 / 2) * tr( twoD & a_ ) * I_ ) - fc * tr( twoD & a_ ) * a_ ) ) + 2 * ci_ * twoD & ( I_ - 3 * a_ )
//+ twoSymm( C )
//- zeta_ / 2 * ( (tau_ & twoD) + (twoD & tau_) )
//- fvm::Sp( epsilon_ / etaP_ * tr(tau_) + 1/lambda_, tau_ )
);

aEqn().relax();
solve(aEqn);
// Viscoelastic stress
tau_ = (etaP_ / keisi_) * (a_ - I_);
}

************************************************** ****************
************************************************** ****************

where "a" , "tau" and "I" are symmetric tensors and "Vor" is asymmetric tensor. and the others are constant parameters.


Best
mohsenkh599 is offline   Reply With Quote

Reply

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
Thermocapillary free surface flow zakifoam OpenFOAM Running, Solving & CFD 6 August 24, 2014 03:24
Tensions on free surface interFoam sega OpenFOAM Running, Solving & CFD 0 April 29, 2008 14:26
VOF simulation free surface caused by gas inject rensb FLUENT 4 July 20, 2006 04:45
Variable Density - Free Surface with FIDAP Vitaliy Pavlyk FLUENT 7 May 2, 2000 15:56
Modeling Free Surface Flows Elliot Schwartz Main CFD Forum 5 August 25, 1998 21:03


All times are GMT -4. The time now is 15:42.