[niko0807]

Dear Foamers,

I am currently working on a temperature dependent viscosity model for a modified solver based on simpleFoam and buoyantBuossinesqSimpleFoam.

I am planning on using the powerLaw viscosity as base for modification. Therefore, I have made a copy of it to:

OpenFOAM/user-2.3.x/src/transportModels/incompressible/viscosityModels/nuTemperature

I have changed the filenames to

nuTemperature.C

and

nuTemperature.H

Furthermore in Make/file I have

nuTemperature.C

LIB = $(FOAM_USER_LIBBIN)/libusernuTemperature

and in Make/options

EXE_INC = \
	-I$(LIB_SRC)/transportModels/incompressible/lnInclude/ \
	-I$(LIB_SRC)/finiteVolume/lnInclude 


LIB_LIBS = 
	-lfiniteVolume

Then I have run wmake and added

libs 			("libusernuTemperature.so");

to my test case controlDict, but when I run my solver I get the error

Build  : 2.3.x-6b18d98df2bc
Exec   : nuTemperatureSimpleFoam
Date   : Oct 14 2015
Time   : 18:05:36
Host   : "hulk2.iha.dk"
PID    : 34454
Case   : /home/dan/run/151009singleCylinder_tvsf_testCase
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Duplicate entry powerLaw in runtime selection table viscosityModel
#0	/usr/local/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam5error14safePrintStackERSo+0x27) [0x7f076e59a207]
#1	/usr/local/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libincompressibleTransportModels.so(_ZN4Foam14viscosityModel31adddictionaryConstructorToTableINS_15viscosityModels8powerLawEEC1ERKNS_4wordE+0xce) [0x7f076e91ce1e]
#2	/home/dan/OpenFOAM/dan-2.3.x/platforms/linux64GccDPOpt/lib/libusernuTemperature.so(+0x15ee9) [0x7f0769871ee9]
#3	/home/dan/OpenFOAM/dan-2.3.x/platforms/linux64GccDPOpt/lib/libusernuTemperature.so(+0x18d66) [0x7f0769874d66]
Create mesh for time = 0

Reading thermophysical properties

Reading field T

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model nuTemperature


--> FOAM FATAL ERROR: 
Unknown viscosityModel type nuTemperature

Valid viscosityModels are : 

5
(
BirdCarreau
CrossPowerLaw
HerschelBulkley
Newtonian
powerLaw
)


    From function viscosityModel::New(const volVectorField&, const surfaceScalarField&)
    in file viscosityModels/viscosityModel/viscosityModelNew.C at line 53.

FOAM exiting

^Z
[10]+  Stopped                 nuTemperatureSimpleFoam

Can someone please see where the problem is?

[alexeym]

Hi,

It is not enough just to rename file, you have to correct this piece of code:

...
public:

    //- Runtime type information
    TypeName("powerLaw");
...

in .H file (nuTemperature.H in your case). Put nuTemperature instead of powerLaw there.

UPD. And surely you have to rename class.

[TemC]

Hi foamers,

I'am working with simpleFoam, and I have implemented a new viscosity model based on the Herschel-Bulkley model.

Through the forum I found several disussions (and documents) about the implementation of a new viscosity model in OpenFOAM.

I followed all the tips and steps that are suggested there, but I still don't manage to run a case with my new viscosity model. I always got the error message posted by "niko08/07", suggesting that my new viscosity model is not valid.

Can somebody give a step by step workflow to handle this issue? Any help would be well appreciated.

Thanks in advance, and have a nice week.

Regards.