[daniel_BA]

Hi guys, first time I submit a thread to the forum. Although I've used many tips provided here.
I've been trying for a while now to impose a velocity-field-dependent momentum source in a simpleFOAM simulation by configuring fvOptions.
I searched many times for similar questions in the forum. Although some have helped me to configure a basic fvOptions file, for uniform, constant source, for example, none lead me to the final task.


configuring my source field through something like

const vector fonte (fonteX,fonteY,fonteZ);
uSource[i] -= fonte;


where fonteX, fonteY and fonteZ are scalar functions of the velocity components.

The computation carry on if I set these scalars are set as constants or even as function of the cell volume accessed through


const scalarField& Vol = mesh().V();


If I set fonteX, fonteY and fonteZ as functions of the velocity components, the simulation is stopped before the first iteration.

I'm trying to access the velocity field through the following lines:


const vectorField& Vel = eqn.psi();
const direction Xcmpt = 0;
const direction Ycmpt = 1;
const direction Zcmpt = 2;
const scalarField& U = Vel.component(Xcmpt);
const scalarField& V = Vel.component(Ycmpt);
const scalarField& W = Vel.component(Zcmpt);


and the message I get is


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 __kernel_sigreturn
#3 Foam::fv::sourceUFvOptionvectorSource::addSup(Foam ::fvMatrix<Foam::Vector<double> >&, int) at system/fvOptions.codedSource:32
#4 ? in ~/OpenFOAM/OpenFOAM-v1806/platforms/linuxGccDPInt32Opt/bin/simpleFoam
#5 main in ~/OpenFOAM/OpenFOAM-v1806/platforms/linuxGccDPInt32Opt/bin/simpleFoam
#6 __libc_start_main in /lib/i386-linux-gnu/libc.so.6
#7 ? in ~/OpenFOAM/OpenFOAM-v1806/platforms/linuxGccDPInt32Opt/bin/simpleFoam

[1]+ Segmentation fault simpleFoam >> run.log


I'm not experienced in C++ programming, but even contributions that demand knowledge about C++ would be very appreciated.


Best regards,
Daniel