[Pyrokrates]

Hello,


I would like to implement a udf for my uds scalar transport.
What I´m planning is to change the flux function as a function of a new calculated velocity (not the one of the fluid).


No I have two problems with following code:

#include "udf.h"

DEFINE_UDS_FLUX(nabla_v_C, f, t, i)
{
    real NV_VEC(psi), NV_VEC(A); // declaring vectors psi and A                  
    
    NV_D(psi, =, F_U(f,t), F_V(f,t), F_W(f,t)); // defining psi in terms of velocity field
    //NV_S(psi, *=, F_R(f,t)); // multiplying density to get psi vector

    //F_AREA(A,f,t); // face normal vector returned from F_AREA


    //return NV_DOT(psi,A); // dot product of the two returned


    return F_FLUX(f, t);

 }

1) If I use only the return of F_Flux(f, t), my program is working (code snippet from manual). But this is the calculation with fluid velocity... I want to change it like mentioned above.
2) When I tried the little more complex code from manual (all other lines in this function, I found out, that the NV_D line crashes my program instantly but I don´t know why...

The solution is initialized and I have a solution of 200 time steps using porous media.


The error in console and text file is following:

Error [node 999999] [time 3/29/22 17:10:12] Abnormal Exit!

Do you have any idea, why the code given by the manual does not work for me? And how can I replace the F_Flux function by my own to use my own velocity field?


Thanks in advance


Chris

[AlexanderZ]

code from manual should work well

how did you found out that the problem is in NV_D
code above is doing nothing but

return F_FLUX(f, t);

your log doesn;t make any sense as that error could mean whatever.
most likely the problem in your case not in udf