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

mesh.surfaceInterpolation::deltaCoeffs()

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

Like Tree10Likes
  • 10 Post By santiagomarquezd

Reply
 
LinkBack Thread Tools Display Modes
Old   January 23, 2011, 09:49
Default mesh.surfaceInterpolation::deltaCoeffs()
  #1
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 10
ata is on a distinguished road
Hi FOAMers
can any one tell me what is the mesh.surfaceInterpolation::deltaCoeffs() in alphaCourantNo.H file of interFoam solver?
I mean what are mesh.surfaceInterpolation and deltaCoeffs()?
best regards
ata is offline   Reply With Quote

Old   January 24, 2011, 13:21
Default
  #2
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 16
santiagomarquezd will become famous soon enough
Hi, looking at surfaceInterpolation.C you have that deltaCoeffs are

Code:
00235     // Set local references to mesh data
00236     const volVectorField& C = mesh_.C();
00237     const unallocLabelList& owner = mesh_.owner();
00238     const unallocLabelList& neighbour = mesh_.neighbour();
00239     const surfaceVectorField& Sf = mesh_.Sf();
00240     const surfaceScalarField& magSf = mesh_.magSf();
00241 
00242     forAll(owner, facei)
00243     {
00244         vector delta = C[neighbour[facei]] - C[owner[facei]];
00245         vector unitArea = Sf[facei]/magSf[facei];
00246 
00247         // Standard cell-centre distance form
00248         //DeltaCoeffs[facei] = (unitArea & delta)/magSqr(delta);
00249 
00250         // Slightly under-relaxed form
00251         //DeltaCoeffs[facei] = 1.0/mag(delta);
00252 
00253         // More under-relaxed form
00254         //DeltaCoeffs[facei] = 1.0/(mag(unitArea & delta) + VSMALL);
00255 
00256         // Stabilised form for bad meshes
00257         DeltaCoeffs[facei] = 1.0/max(unitArea & delta, 0.05*mag(delta));
00258     }

delta=C_N-C_P
unitArea = Sf/|Sf|
deltaCoeffs = 1.0/max(unitArea & delta, 0.05*mag(delta))

where & is the dot product. Then they're the reciprocal of center to center distances projected over the face normal (with a stabilization factor) of cells' shared face.

Regards.
linch, kaifu, zhernadi and 7 others like this.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   January 24, 2011, 23:27
Default mesh.surfaceInterpolation::deltaCoeffs()
  #3
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 10
ata is on a distinguished road
Thank you very much again
ata 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



All times are GMT -4. The time now is 08:41.