 |
 |
 |
|
|
|
[Sponsors] | ||||
UDF to change angular velocity after every timestep based on the previous value |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 9, 2022, 08:13
|
UDF to change angular velocity after every timestep based on the previous value
|
#1 |
|
New Member
Join Date: Jul 2019
Posts: 14
Rep Power: 7  |
Hello,
I am trying to run a simulation where the angular velocity of a dynamic mesh zone changes every time step based on this equation of motion where moment_z is the moment on the two aerofoils defined as one wall. omega_prev + = ((1/I) * (moment_z - (alpha * omega_prev))) When I set up the dynamic mesh with the udf, the initialisation returns an "All compute processes received SIGFPE." error. I am not sure what the cause of this error is. Could someone more knowledgeable please help me solve this? This is the UDF I am using: Code: Code:
#include "udf.h"
real current_dt;
real omega_old;
real omega_new;
real moment_z;
static I = 0.0275; /* Moment of inertia */
static alpha = 0.00626;
static real omega_prev = 0.0;
DEFINE_CG_MOTION (rotational_motion1, dt, vel, omega, time, dtime)
{
#if !RP_HOST
current_dt = RP_Get_Integer("time-step");
if (current_dt ==1)
omega_prev = -128;
else
omega_prev = omega[2];
Message("Omega_old1 is %g",omega_prev);
Domain *d =Get_Domain(1);
Thread *t_object = Lookup_Thread(d, 2);
real moment [ND_ND], cg[ND_ND], force [ND_ND];
Compute_Force_And_Moment(d, t_object, cg, force, moment, TRUE);
moment_z= moment [2];
/* Equation of motion */
omega_prev += ((1/I) * (moment_z - (alpha * omega_prev)));
Message("Omega_new1 is %g",omega_prev);
omega[2] = omega_prev;
#endif
}
|
|
|
|
|
|
November 10, 2022, 01:43
|
|
#2 |
|
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 |
show compilation log
__________________
best regards ****************************** press LIKE if this message was helpful |
|
|
|
|
|
|
November 11, 2022, 06:48
|
Compilation log
|
#3 |
|
New Member
Join Date: Jul 2019
Posts: 14
Rep Power: 7 |
Dear Alexander,
Thank you for your response. Here is the compilation log for the udf. Code:
Working...
for d in lnamd64/[23]*; do \
( \
cd $d; \
rm -rf *.{c,h,cpp,hpp}; \
for f in ../../src/*.{c,h,cpp,hpp} ../../src/makefile ../../src/user.udf; do \
if [ -f $f -a ! -f `basename $f` ]; then \
echo "# linking to" $f "in" $d; \
ln -s $f .; \
fi; \
done; \
echo ""; \
echo "# building library in" $d; \
if [ "" = "1" ]; then \
echo "# using gcc64"; \
make ARCHC=gcc64 -k>makelog 2>&1; \
else \
if [ "" = "1" ]; then \
echo "# using gcc"; \
make ARCHC=gcc -k>makelog 2>&1; \
else \
make -k>makelog 2>&1; \
fi; \
fi;\
cat makelog; \
) \
done
# linking to ../../src/rotCG_motion1.c in lnamd64/2ddp_host
# linking to ../../src/rotCG_motion2.c in lnamd64/2ddp_host
# building library in lnamd64/2ddp_host
make[1]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
# Generating udf_names.c because of makefile rotCG_motion1.c
make[2]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
make libudf.so "CFLAGS=-D_lnamd64 -D_GNU_SOURCE -fpic -shared -ansi -Wall -O -DPTR_RESTRICT= " "LDFLAGS=-shared -lm"
make[3]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
# Compiling rotCG_motion1.o because of rotCG_motion1.c
cc -D_lnamd64 -D_GNU_SOURCE -fpic -shared -ansi -Wall -O -DPTR_RESTRICT= -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/`expr "\`pwd\`" : '.*/\(.*\)/[23].*'`/`basename "\`pwd\`"` -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/main -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/addon-wrapper -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/io -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/species -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/pbns -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/numerics -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/sphysics -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/storage -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/mphase -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/bc -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/models -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/material -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/amg -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/util -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/mesh -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/udf -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/ht -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dx -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/turbulence -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/parallel -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/etc -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/ue -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dpm -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dbns -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/cortex/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/client/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/tgrid/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/mpi_wrapper/include -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/mpi_wrapper/src -I. -c rotCG_motion1.c
rotCG_motion1.c:7:13: warning: I defined but not used [-Wunused-variable]
static real I = 0.0275; /* Moment of inertia */
^
rotCG_motion1.c:8:13: warning: alpha defined but not used [-Wunused-variable]
static real alpha = 0.00626;
^
rotCG_motion1.c:9:13: warning: omega_prev defined but not used [-Wunused-variable]
static real omega_prev = 0.0;
^
# Linking libudf.so because of udf_names.c rotCG_motion1.o
cc -shared -lm udf_names.o rotCG_motion1.o -o libudf.so
make[3]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
make[2]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
make[1]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host'
# linking to ../../src/rotCG_motion1.c in lnamd64/2ddp_node
# linking to ../../src/rotCG_motion2.c in lnamd64/2ddp_node
# building library in lnamd64/2ddp_node
make[1]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
# Generating udf_names.c because of makefile rotCG_motion1.c
make[2]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
make libudf.so "CFLAGS=-D_lnamd64 -D_GNU_SOURCE -fpic -shared -ansi -Wall -O -DPTR_RESTRICT= " "LDFLAGS=-shared -lm"
make[3]: Entering directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
# Compiling rotCG_motion1.o because of rotCG_motion1.c
cc -D_lnamd64 -D_GNU_SOURCE -fpic -shared -ansi -Wall -O -DPTR_RESTRICT= -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/`expr "\`pwd\`" : '.*/\(.*\)/[23].*'`/`basename "\`pwd\`"` -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/main -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/addon-wrapper -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/io -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/species -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/pbns -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/numerics -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/sphysics -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/storage -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/mphase -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/bc -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/models -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/material -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/amg -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/util -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/mesh -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/udf -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/ht -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dx -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/turbulence -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/parallel -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/etc -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/ue -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dpm -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/src/dbns -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/cortex/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/client/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/tgrid/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/src -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/mpi_wrapper/include -I/usr045/ansys_inc/v182/fluent/fluent18.2.0/multiport/mpi_wrapper/src -I. -c rotCG_motion1.c
# Linking libudf.so because of udf_names.c rotCG_motion1.o
cc -shared -lm udf_names.o rotCG_motion1.o -o libudf.so
make[3]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
make[2]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
make[1]: Leaving directory `/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node'
You can also see the 'log'-file in
the working directory for this compilation history
Done.
Opening library "/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf"...
Library "/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_host/libudf.so" opened
Opening library "/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf"...
Library "/usr045/people/greeshma45/Desktop/DFBI_Trial/CG_Motion/libudf/lnamd64/2ddp_node/libudf.so" opened
rotational_motion1
Done.
Opening
|
|
|
|
|
|
|
November 16, 2022, 21:58
|
|
#4 |
|
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 |
this is compilation log I got, when I tried to compile your code
fix these errors, put here revised version of code and describe problems if any Code:
..\..\src\debug.c(18) : error C2275: 'Domain' : illegal use of this type as an expression ..\..\src\debug.c(18) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(19) : error C2275: 'Thread' : illegal use of this type as an expression ..\..\src\debug.c(19) : error C2065: 't_object' : undeclared identifier ..\..\src\debug.c(19) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(19) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(19) : warning C4024: 'Lookup_Thread' : different types for formal and actual parameter 1 ..\..\src\debug.c(20) : error C2275: 'real' : illegal use of this type as an expression ..\..\src\debug.c(20) : error C2146: syntax error : missing ';' before identifier 'moment' ..\..\src\debug.c(20) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(20) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(20) : error C2065: 'cg' : undeclared identifier ..\..\src\debug.c(20) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(20) : error C2065: 'force' : undeclared identifier ..\..\src\debug.c(20) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(21) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(21) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(21) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 1 ..\..\src\debug.c(21) : error C2065: 't_object' : undeclared identifier ..\..\src\debug.c(21) : warning C4047: 'function' : 'Thread *' differs in levels of indirection from 'int' ..\..\src\debug.c(21) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 2 ..\..\src\debug.c(21) : error C2065: 'cg' : undeclared identifier ..\..\src\debug.c(21) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(21) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 3 ..\..\src\debug.c(21) : error C2065: 'force' : undeclared identifier ..\..\src\debug.c(21) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(21) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 4 ..\..\src\debug.c(21) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(21) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(21) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 5 ..\..\src\debug.c(22) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(22) : error C2109: subscript requires array or pointer type
__________________
best regards ****************************** press LIKE if this message was helpful |
|
|
|
|
|
|
November 17, 2022, 10:00
|
Edited UDF
|
#5 |
|
New Member
Join Date: Jul 2019
Posts: 14
Rep Power: 7 |
Hi Alexander,
Thank you for your response. I made a few changes to the UDF and the method. Instead of using the dynamic mesh method, now I am trying to use the UDF to change the rotational velocity in the Cell Zone conditions tab. I changed the UDF as follows. The rotor shown in the image should depend on the moments on the aerofoils (id 13 and 14). The code compiles without any errors for me. But it seems like the UDF is not being read as the first message I try to print does not appear on the console. There are no options to 'hook' the UDF other than in the Cell Zone conditions tab. Could you please help me solve this issue? Code:
#include "udf.h"
real current_dt;
real omega_old;
real omega_new;
real moment_z;
real moment_z2;
real moment_R;
static real I = 0.0275; /* Moment of inertia */
static real alpha = 0.00626;
static real omega_prev = -128.0;
DEFINE_ZONE_MOTION (rotational_motion1,omega,axis,origin,velocity,time,dtime)
{
#if !RP_HOST
Message("Begin udf1");
omega_prev = omega[2];
Message("Omega_old1 is %g",omega_prev);
Domain *d =Get_Domain(1);
Thread *t_object = Lookup_Thread(d, 13);
Thread *t_object2 = Lookup_Thread(d, 14);
real moment[ND_ND], cg[ND_ND], force[ND_ND];
Compute_Force_And_Moment(d, t_object, cg, force, moment, TRUE);
moment_z= moment[2];
Compute_Force_And_Moment(d, t_object2, cg, force, moment, TRUE);
moment_z2= moment[2];
/* Equation of motion */
moment_R = moment_z + moment_z2;
omega_prev += ((1/I) * (moment_R- (alpha * omega_prev)));
Message("Omega_new1 is %g",omega_prev);
omega[2] = omega_prev;
origin[0]= 0.0;
origin[1]= 0.325;
origin[2]= 0.0;
N3V_D(axis,=,0.0,0.0,1.0);
#endif
}
|
|
|
|
|
|
|
November 18, 2022, 00:27
|
|
#6 |
|
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 |
this is compilation log of your code "without errors"
Code:
..\..\src\debug.c(18) : error C2275: 'Domain' : illegal use of this type as an expression ..\..\src\debug.c(18) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(19) : error C2275: 'Thread' : illegal use of this type as an expression ..\..\src\debug.c(19) : error C2065: 't_object' : undeclared identifier ..\..\src\debug.c(19) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(19) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(19) : warning C4024: 'Lookup_Thread' : different types for formal and actual parameter 1 ..\..\src\debug.c(20) : error C2275: 'Thread' : illegal use of this type as an expression ..\..\src\debug.c(20) : error C2065: 't_object2' : undeclared identifier ..\..\src\debug.c(20) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(20) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(20) : warning C4024: 'Lookup_Thread' : different types for formal and actual parameter 1 ..\..\src\debug.c(21) : error C2275: 'real' : illegal use of this type as an expression ..\..\src\debug.c(21) : error C2146: syntax error : missing ';' before identifier 'moment' ..\..\src\debug.c(21) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(21) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(21) : error C2065: 'cg' : undeclared identifier ..\..\src\debug.c(21) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(21) : error C2065: 'force' : undeclared identifier ..\..\src\debug.c(21) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(23) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(23) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(23) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 1 ..\..\src\debug.c(23) : error C2065: 't_object' : undeclared identifier ..\..\src\debug.c(23) : warning C4047: 'function' : 'Thread *' differs in levels of indirection from 'int' ..\..\src\debug.c(23) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 2 ..\..\src\debug.c(23) : error C2065: 'cg' : undeclared identifier ..\..\src\debug.c(23) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(23) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 3 ..\..\src\debug.c(23) : error C2065: 'force' : undeclared identifier ..\..\src\debug.c(23) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(23) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 4 ..\..\src\debug.c(23) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(23) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(23) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 5 ..\..\src\debug.c(24) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(24) : error C2109: subscript requires array or pointer type ..\..\src\debug.c(25) : error C2065: 'd' : undeclared identifier ..\..\src\debug.c(25) : warning C4047: 'function' : 'Domain *' differs in levels of indirection from 'int' ..\..\src\debug.c(25) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 1 ..\..\src\debug.c(25) : error C2065: 't_object2' : undeclared identifier ..\..\src\debug.c(25) : warning C4047: 'function' : 'Thread *' differs in levels of indirection from 'int' ..\..\src\debug.c(25) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 2 ..\..\src\debug.c(25) : error C2065: 'cg' : undeclared identifier ..\..\src\debug.c(25) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(25) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 3 ..\..\src\debug.c(25) : error C2065: 'force' : undeclared identifier ..\..\src\debug.c(25) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(25) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 4 ..\..\src\debug.c(25) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(25) : warning C4047: 'function' : 'real *' differs in levels of indirection from 'int' ..\..\src\debug.c(25) : warning C4024: 'Compute_Force_And_Moment' : different types for formal and actual parameter 5 ..\..\src\debug.c(26) : error C2065: 'moment' : undeclared identifier ..\..\src\debug.c(26) : error C2109: subscript requires array or pointer type
__________________
best regards ****************************** press LIKE if this message was helpful |
|
|
|
|
|
|
November 18, 2022, 00:34
|
|
#7 |
|
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 |
Code:
#include "udf.h"
real current_dt;
real omega_old;
real omega_new;
real moment_z;
real moment_z2;
real moment_R;
static real I = 0.0275; /* Moment of inertia */
static real alpha = 0.00626;
static real omega_prev = -128.0;
DEFINE_ZONE_MOTION(rotational_motion1,omega,axis,origin,velocity,time,dtime)
{
#if !RP_HOST
Domain *d;
Thread *t_object,*t_object2;
real moment[ND_ND], cg[ND_ND], force[ND_ND];
d =Get_Domain(1);
Message0("Begin udf1");
omega_prev = omega[2];
Message0("Omega_old1 is %g",omega_prev);
t_object = Lookup_Thread(d, 13);
t_object2 = Lookup_Thread(d, 14);
Compute_Force_And_Moment(d, t_object, cg, force, moment, TRUE);
moment_z= moment[2];
Compute_Force_And_Moment(d, t_object2, cg, force, moment, TRUE);
moment_z2= moment[2];
/* Equation of motion */
moment_R = moment_z + moment_z2;
omega_prev += ((1/I) * (moment_R- (alpha * omega_prev)));
Message0("Omega_new1 is %g",omega_prev);
omega[2] = omega_prev;
origin[0]= 0.0;
origin[1]= 0.325;
origin[2]= 0.0;
N3V_D(axis,=,0.0,0.0,1.0);
#endif
return;
}
__________________
best regards ****************************** press LIKE if this message was helpful |
|
|
|
|
|
|
| Tags |
| rotating domain, udf |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| udf for one dimensional linear motion based on force | maccheese | Fluent UDF and Scheme Programming | 2 | September 1, 2019 03:18 |
| Fluent do not use my velocity field(by UDF) to solve energy equation | tangleiplus | Fluent UDF and Scheme Programming | 6 | January 21, 2019 22:28 |
| UDF velocity inlet. Different values (compile - Fluent) | asking | FLUENT | 0 | July 13, 2018 18:02 |
| adapt viscous resistance based on velocity / UDF | Häwimeddel | FLUENT | 2 | April 20, 2009 11:08 |
| UDF For Angular velocity for gear pump | Mustafa Ayad | FLUENT | 0 | December 30, 2008 03:52 |
Linear Mode
