CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent UDF and Scheme Programming

Moving average using UDF

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 12, 2021, 18:36
Default Moving average using UDF
  #1
New Member
 
jo samuel
Join Date: Apr 2020
Posts: 3
Rep Power: 3
josamuel is on a distinguished road
Hi all,

I'm trying to obtain the pressure values at a particular interface inside the domain in a transient 2D model using UDF. This calculation happens for every iteration. And I also want the moving average(aka rolling average) of these obtained pressure values over few timesteps(not iterations). After obtaining this averaged value I use them back as an input into the simulation.

I have tried a code that I have pasted below. I was able to successfully obtain the pressure values at the interface for every iteration. But, UDF calculates the moving average of the obtained pressure values for every iteration, instead, I need it to be calculated over a few timesteps. Please help me where I went wrong.

Thanks in advance!




#include "udf.h"

#define NumTempsToAvg 6 /* this number is the number of timesteps taken for average */

real Pr;

int temps[NumTempsToAvg];

real sum = 0.0;

DEFINE_ADJUST(Pint, domain)

{
int i;
i= N_TIME; /* obtains the integer value of the timestep*/

real coord[ND_ND];
face_t f;
cell_t c0;
int ID = 23;
Thread *thread, *t0;
thread = Lookup_Thread(domain, ID);
begin_f_loop(f, thread)
{
c0 = F_C0(f,thread);
t0 = THREAD_T0(thread);
F_CENTROID(coord,c0,t0);

printf("x = %f\t y = %f\n\n ",coord[0], coord[1]);
Pr += F_P(c0,t0);

printf("Pr = %f\n\n", Pr);
}

end_f_loop(f,thread)

Pr=Pr/14; /* this averages over the 14 faces */
printf("Pr = %f\n", Pr);

for ( i= 1; i < NumTempsToAvg; i++)
{
temps[i] = Pr; /* array name is temps */
sum += temps[i];
printf("the sum is %f\n\n", sum);
printf("avg: %f\n", sum / NumTempsToAvg); /* moving average */
sum -= temps[i];
temps[i] = Pr;
sum += temps[i];
}

}
josamuel is offline   Reply With Quote

Old   June 17, 2021, 00:09
Default
  #2
New Member
 
Genevieve
Join Date: Jun 2021
Posts: 1
Rep Power: 0
gensqy is on a distinguished road
You can try using DEFINE_EXECUTE_AT_END, I believe it runs at the end of every timestep in a transient run.
gensqy is offline   Reply With Quote

Old   June 25, 2021, 14:36
Default
  #3
New Member
 
jo samuel
Join Date: Apr 2020
Posts: 3
Rep Power: 3
josamuel is on a distinguished road
Thanks Gensqy,

I used DEFINE_EXECUTE_AT_END and obtained values at the end of the timestep. Thanks for the idea and reply. However, my final goal of this UDF is to get the moving average. Unfortunately, it still ends up showing some errors. Any leads would be helpful.

thanks in advance.
josamuel is offline   Reply With Quote

Old   June 27, 2021, 21:41
Default
  #4
Senior Member
 
Alexander
Join Date: Apr 2013
Posts: 1,760
Rep Power: 28
AlexanderZ will become famous soon enoughAlexanderZ will become famous soon enough
show code you are using
show compilation log file
explain in details problems/errors which you have
__________________
best regards


******************************
press LIKE if this message was helpful
AlexanderZ is offline   Reply With Quote

Reply

Tags
fluent, moving average, transient 2d, udf

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Moving mesh udf problem vitali31 Fluent UDF and Scheme Programming 20 August 29, 2018 11:42
UDF for moving wall (oscillating) Sambhu.jayachandran Fluent UDF and Scheme Programming 3 March 12, 2018 03:23
UDF for a moving sphere in 2D Sufyan Abushaala Fluent UDF and Scheme Programming 0 February 19, 2014 06:19
Moving Reference Frame (MRF) with UDF Tobard FLUENT 0 April 14, 2011 15:18
Moving Reference frame - UDF - Moving mesh modisa FLUENT 0 April 18, 2008 13:31


All times are GMT -4. The time now is 18:57.