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

turbulent schmidt number with scalarTransportFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 25, 2013, 02:38
Default turbulent schmidt number with scalarTransportFoam
  #1
New Member
 
dongjin kim
Join Date: Jul 2013
Posts: 4
Rep Power: 3
dgenekim is on a distinguished road
Hi Foamers,

I need a solver which can deal with turbulent diffusivity instead of constant DT.
I have tried to modify the original scalarTransportFoam code as:

In createField.H, I added

===================
Info<< "Reading field nut\n" << endl;

volScalarField nut
(
IOobject
(
"nut",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
=========================

and

=========================
Info<< "Reading diffusivity ScT\n" << endl;

dimensionedScalar ScT
(
transportProperties.lookup("ScT")
);
=========================

and modified the scalarTransportFoam.C as:

=========================
solve
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(nut/ScT, T)
==
fvOptions(T)
========================

The compiling was done without error message. However, when I ran a test simulation I got the error message as:

========================
--> FOAM FATAL IO ERROR:
"ill defined primitiveEntry starting at keyword 'laplacian(nut' on line 37 and ending at line 60"

file: /home/dongjin/OpenFOAM/dongjin-2.2.1/practice/crossflowmixingchannel/scalarTransportFoam/system/fvSchemes at line 60.

From function primitiveEntry::readEntry(const dictionary&, Istream&)
in file lnInclude/IOerror.C at line 132.

FOAM exiting
=======================

Before running a simulation, I copied the 'nut' file from the previous simpleFoam simulation under the '0' folder, and modified the transportProperties and fvScheme files as:

===================
laplacianSchemes
{
default none;
laplacian(nut/ScT,T) Gauss linear corrected;
}
====================
and
====================
ScT ScT [ 0 0 0 0 0 0 0 ] 0.15;
=====================

Please reply me why I got the above mentioned error?
dgenekim is offline   Reply With Quote

Old   July 25, 2013, 02:54
Default
  #2
New Member
 
dongjin kim
Join Date: Jul 2013
Posts: 4
Rep Power: 3
dgenekim is on a distinguished road
I could make the run further by changing fvScheme as:

from
laplacian(nut/ScT,T) Gauss linear corrected;
to
laplacian(nut|ScT,T) Gauss linear corrected;

but now I got:

=================================
file: /home/dongjin/OpenFOAM/dongjin-2.2.1/practice/crossflowmixingchannel/scalarTransportFoam/0/nut.boundaryField.upperWall from line 2906 to line 2910.

From function fvPatchField<Type>::New(const fvPatch&, const DimensionedField<Type, volMesh>&, const dictionary&)
in file /opt/openfoam221/src/finiteVolume/lnInclude/fvPatchFieldNew.C at line 136.

FOAM exiting
===============================

Any idea?
dgenekim is offline   Reply With Quote

Old   July 25, 2013, 03:14
Default
  #3
New Member
 
dongjin kim
Join Date: Jul 2013
Posts: 4
Rep Power: 3
dgenekim is on a distinguished road
In nut file, line number from 2906 to 2910 contains:

2903 }
2904 upperWall
2905 {
2906 type nutkWallFunction;
2907 Cmu 0.09;
2908 kappa 0.41;
2909 E 9.8;
2910 value nonuniform List<scalar>
dgenekim is offline   Reply With Quote

Old   July 25, 2013, 03:49
Default
  #4
New Member
 
dongjin kim
Join Date: Jul 2013
Posts: 4
Rep Power: 3
dgenekim is on a distinguished road
I just make the run successful by simply modifying nut file as:

from
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value nonuniform List<scalar>

to
type calculated;
value nonuniform List<scalar>

and by fixing fvScheme file as:

from
laplacian(nut|ScT,T) Gauss linear corrected;
to
laplacian((nut|ScT),T) Gauss linear corrected;


I will leave this posting for other users information ^^
dgenekim 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
Turbulent Schmidt Number John FLUENT 8 June 1, 2015 05:04
decomposePar pointfield flying OpenFOAM Running, Solving & CFD 28 December 30, 2013 16:05
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 2 November 25, 2012 09:54
DecomposePar unequal number of shared faces maka OpenFOAM Pre-Processing 6 August 12, 2010 09:01
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36


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