[bss]

Hello,
I want to calculate 'curvature' in UDF. I have attached formulation of curvature.
'F' is volume fraction.
C_VOF_G(c,t) can calculate gradient of volume fraction
NV_MAG(C_VOF_G(c,t)) can calculate magnitude of volume fraction. I can calculate gradient of volume fraction but My problem is that I can't calculate divergence of volume fraction gradient in UDF.

Thank you in advance.

[anoop.k.unni]

Originally Posted by bss

Hello,
I want to calculate 'curvature' in UDF. I have attached formulation of curvature.
'F' is volume fraction.
C_VOF_G(c,t) can calculate gradient of volume fraction
NV_MAG(C_VOF_G(c,t)) can calculate magnitude of volume fraction. I can calculate gradient of volume fraction but My problem is that I can't calculate divergence of volume fraction gradient in UDF.

Thank you in advance.

Hi did you solve?

[pakk]

Going from the gradient to the divergence (in standard Cartesian coordinates) is super simple... Just sum them.

Maybe I don't understand the problem?

[anoop.k.unni]

Originally Posted by pakk

Going from the gradient to the divergence (in standard Cartesian coordinates) is super simple... Just sum them.

Maybe I don't understand the problem?

I have written a UDF for a multiphase (of) problem dealing with recoil pressure which has to be added as momentum source term. i have attached a figure which shows a definition of my problem and also udf. can you please check whether i have written in the right way because i am not able to run the case.
#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "mem.h"
#include "unsteady.h"


DEFINE_ADJUST(vol_gradient, domain)
{
Thread *t;
cell_t c;
real DVOFNX,DVOFNY,DVOFDNZ;
domain = Get_Domain(6);

/* Fill UDS with the variable. */

thread_loop_c (t,domain)
{
begin_c_loop (c,t)
{
C_UDSI(c,t,0) = C_VOF(c,t);

DVOFNX = C_UDSI_G(c,t,0)[0];
DVOFNY = C_UDSI_G(c,t,0)[1];
DVOFNZ = C_UDSI_G(c,t,0)[2];

C_UDMI(c,t,0) = DVOFNX;
C_UDMI(c,t,1) = DVOFNY;
C_UDMI(c,t,2) = DVOFNZ;

}
end_c_loop (c,t)
}

thread_loop_c (t,domain)
{
begin_c_loop (c,t)
{
C_UDMI(c,t,3) = NV_MAG(C_UDSI_G(c,t,0));
}
end_c_loop (c,t)
}

DEFINE_ON_DEMAND(store_gradient)
{
Domain *domain;
cell_t c;
Thread *t;
domain = Get_Domain(6);

/* Fill the UDM with magnitude of gradient. */

*thread_loop_c (t,domain)
{
begin_c_loop (c,t)
{
C_UDMI(c,t,4) = C_UDMI(c,t,0)/NV_MAG(C_UDSI_G(c,t,0));
C_UDMI(c,t,5) = C_UDMI(c,t,1)/NV_MAG(C_UDSI_G(c,t,0));
C_UDMI(c,t,6) = C_UDMI(c,t,2)/NV_MAG(C_UDSI_G(c,t,0));
}
end_c_loop (c,t)
}
}
}

DEFINE_SOURCE(x_recoil,c,t,dS,eqn)
{
real x[ND_ND];
real source;
real time;
real temp;
time=RP_Get_Real("flow-time"); //Acquire time from FLUENT solver
C_CENTROID(x,c,t);
temp=C_T(c,t); //Acquire the face centroid location
if (C_T(c,t) >= 3201.)
{
source=(54715.5*exp(13.14*(temp-(3200)/temp)))*(C_UDMI(c,t,4));
C_UDMI(c,t,7) = source;
dS[eqn] = 0;
}
else
{
source = dS[eqn] = 0.;
}
return source;
}
//same for y and z

[pakk]

Originally Posted by anoop.k.unni

can you please check whether i have written in the right way because i am not able to run the case.

"not able to run the case" can mean 20 things. Please be more detailed, that makes it much easier to help you.

[anoop.k.unni]

Originally Posted by pakk

Going from the gradient to the divergence (in standard Cartesian coordinates) is super simple... Just sum them.

Maybe I don't understand the problem?

Originally Posted by pakk

"not able to run the case" can mean 20 things. Please be more detailed, that makes it much easier to help you.

After initializing the case with adjust gradient hoking and it is showing' Received signal SIGSEGV. MPI Application rank 0 exited before MPI_Finalize()' So, I want to know whether is it because of any mistake in UDF or other (improper boundary conditions,patching issue)

[pakk]

Originally Posted by anoop.k.unni

After initializing the case with adjust gradient hoking and it is showing' Received signal SIGSEGV. MPI Application rank 0 exited before MPI_Finalize()' So, I want to know whether is it because of any mistake in UDF or other (improper boundary conditions,patching issue)

That is a good question. To find out if your error comes from the UDF or from somewhere else: run it without hooking the UDF! If you don't get an error, the problem was in your UDF.

If you do get a problem: debug! Simplify your UDF, and find out which line causes the error.