[shaoyx]

Hi Everyone,

I am trying to use codestream to initialize internal velocity and here is my 0/U file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

//internalField	uniform (2 0 0);

internalField   #codeStream
{
    codeInclude
    #{
	#include "fvCFD.H"
    #};

    codeOptions
    #{
	-I$(LIB_SRC)/finiteVolume/lnInclude \
	-I$(LIB_SRC)/meshTools/lnInclude
    #};
	
    codeLibs
    #{
	-lmeshTools \
	-lfiniteVolume
    #};
	
    code

    #{
	const IOdictionary& d = static_cast<const IOdictionary&>(dict);
	const fvMesh& mesh = refCast<const fvMesh>(d.db());

	vectorField T(mesh.nCells());

	const vectorField& CC = mesh.C(); //cell center 


	forAll(CC,cellI)
	{
	    	scalar x = CC[cellI].x();
	    	scalar y = CC[cellI].y();
		
	    	T[cellI] = vector(1-pow(y,2),-sin(x),0);
	}
	T.writeEntry("", os);
    #};
};

boundaryField
{
    frontAndBack
    {
        type            	empty;
    }
    walls
    {
        type            	noSlip;
    }
    outlet
    {
	type			inletOutlet;
	value			$internalField;
	inletValue		uniform (0 0 0);
    }
    inlet
    {
	type			groovyBC;
	variables		"yp=pts().y;para=1-pow(yp,2);";

//	variables		"yp=pts().y;minY=min(yp);maxY=max(yp);para=-(maxY-pos().y)*(pos().y-minY)/(pow(maxY-minY,2))*normal();";

	valueExpression		"para";
	value			uniform (0 0 0);
    }
}

// ************************************************************************* //

But I get errors when running simpleFoam

--> FOAM FATAL IO ERROR: 
compound has already been transfered from token
    on line 0 the empty compound of type List<vector>

file: /home/yanxuan/OpenFOAM/yanxuan-5.0/run/2d_poi/0/U.boundaryField.outlet.value at line 0.

    From function Foam::token::compound& Foam::token::transferCompoundToken(const Foam::Istream&)
    in file lnInclude/token.C at line 99.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::IOerror::abort() at ??:?
#2  Foam::token::transferCompoundToken(Foam::Istream const&) at ??:?
#3  Foam::Istream& Foam::operator>><Foam::Vector<double> >(Foam::Istream&, Foam::List<Foam::Vector<double> >&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#4  Foam::Field<Foam::Vector<double> >::Field(Foam::word const&, Foam::dictionary const&, int) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#5  Foam::inletOutletFvPatchField<Foam::Vector<double> >::inletOutletFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#6  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::inletOutletFvPatchField<Foam::Vector<double> > >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#7  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#8  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#9  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#10  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#11  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#12  ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
Aborted (core dumped)

What is the problem here and what I should do? Thank you!

Charlotte

[simrego]

Hi!


The outlet value is your problem as you can see from the error. Instead of "$internalField;" set to "uniform (0 0 0);" for example. The code will correct this based on the inletOutlet BC.

[shaoyx]

Thanks for your reply! It is working well now!