CFD Online Discussion Forums - problems on dsmc FreeStream Inflow Field Patch setup

Good morning, Im perfoming an validation experiment for a paper on dsmcFoam+
one of the simulation requirements is a non-uniform, vector list (U, T, numerical density)
reading the source code i came across this boundary condition:

Code:

/*---------------------------------------------------------------------------*\

  =========                 |

  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox

   \\    /   O peration     |

    \\  /    A nd           | Copyright (C) 2016-2021 hyStrath

     \\/     M anipulation  |

-------------------------------------------------------------------------------

License

    This file is part of hyStrath, a derivative work of OpenFOAM.



    OpenFOAM is free software: you can redistribute it and/or modify it

    under the terms of the GNU General Public License as published by

    the Free Software Foundation, either version 3 of the License, or

    (at your option) any later version.



    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT

    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

    for more details.



    You should have received a copy of the GNU General Public License

    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.



Class

    dsmcFreeStreamInflowFieldPatch



Description



    Fully deletion boundary patch: atomisticMolecules get deleted instantly upon

    collision with a face on the boundary. (aka vacuum boundary condition)



    Selective deletion of one specie is possible. Other species get

    specularly reflected.



SourceFiles

    dsmcFreeStreamInflowFieldPatch.C





\*---------------------------------------------------------------------------*/



#ifndef dsmcFreeStreamInflowFieldPatch_H

#define dsmcFreeStreamInflowFieldPatch_H



#include "dsmcGeneralBoundary.H"

#include "vector.H"

#include "polyMesh.H"

#include "volFields.H"

#include "dsmcCloud.H"

#include "distribution.H"



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



namespace Foam

{



/*---------------------------------------------------------------------------*\

                           Class dsmcFreeStreamInflowFieldPatch Declaration

\*---------------------------------------------------------------------------*/



class dsmcFreeStreamInflowFieldPatch

:

    public dsmcGeneralBoundary

{



private:



    // Private data

        dictionary propsDict_;



    labelList typeIds_;



    //vector so different translational, rotational

    //vibrational temperatures can be set

    vectorField inletTemperatures_;



    vectorField inletVelocities_;



    // list of (one value) mass densities corresponding to the type Ids

    List <scalarField> numberDensities_;



    // list of species

    List <scalarField> accumulatedParcelsToInsert_;



    //- boundary temperatures

    volVectorField boundaryT_;



    //- boundary velocity

    volVectorField boundaryU_;



    //- list of inlet densities (one entry for each species)

    List<autoPtr<volScalarField> > boundaryNumberDensity_;



    void setProperties();



public:



    //- Runtime type information

    TypeName("dsmcFreeStreamInflowFieldPatch");



    // Constructors



        //- Construct from components

        dsmcFreeStreamInflowFieldPatch

        (

            Time& t,

            const polyMesh& mesh,

            dsmcCloud& cloud,

            const dictionary& dict

        );



    // Destructor



        ~dsmcFreeStreamInflowFieldPatch();





    // Member Functions



        void initialConfiguration();



        //- calculate any required properties

        void calculateProperties();



        //- control the mols

        void controlParcelsBeforeMove();

        void controlParcelsBeforeCollisions();

        void controlParcelsAfterCollisions();



        //- output data

        void output

        (

            const fileName& fixedPathName,

            const fileName& timePath

        );



        //- update properties from a modified dictionary

        void updateProperties(const dictionary&);



        //reset fields if mesh has changed

        void setNewBoundaryFields();







        // Access



};





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



} // End namespace Foam



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



#endif



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

unfortunately i wasn't able to run this condition with my code:

Code:

 boundary

    {

        generalBoundaryProperties

        {

        patchName                                                outflow;

        }



        boundaryModel                    dsmcFreeStreamInflowFieldPatch;



        dsmcFreeStreamInflowFieldPatchProperties

        {

                  typeIds                   (N2cold);                       

                              boundaryU

                              N2cold

                            100

                            (

                                (719.1304875090954 25.172575474152957 0)

                                (719.3627785226914 26.328806602546166 0)

                                (719.5950695362875 27.29233254287383 0)

                                (719.8273605498836 28.54491626529981 0)

                                (720.0596515634796 29.604794799660247 0)

                                (720.2919425770757 30.857378522086222 0) ...

                               



                            );  

                            



                  boundaryT

                              N2cold

                              100

                              (

                                (8.988987042333086 8.988987042333086)

                                (8.98943049818169 8.98943049818169)

                                (8.9907608657275 8.9907608657275)

                                (8.992638162153256 8.992638162153256) ...

              



  

                                );

                            



                            boundaryNumberDensity

                            N2cold

                            100

                            (

                            (6.804743383294108)

                            (6.80455480107654) ...

                          



                            );

          }

    }

the main errors are the following:

Code:

[0] cannot find file

[0] 

[0] file: /home/nanor4k/openfoam/hyStrath/run/hyStrath/dsmcFoam+/supersonicFlatPlate2/processor0/0/boundaryT at line 0.

[0] 

[0]     From function regIOobject::readStream()

[0]     in file db/regIOobject/regIOobjectRead.C at line 237.

and

Code:

--> FOAM FATAL IO ERROR: [5] attempt to read beyond EOF

How can i properly setup this boundary conditions with my vector data ?

best regards nanoR4K.