|
[Sponsors] |
March 22, 2021, 09:47 |
raddi translating deformable cylinder UDF
|
#1 |
New Member
moh
Join Date: Jan 2014
Posts: 20
Rep Power: 12 |
hi every one
i wan simulate a raddi deformable cylinder with plunging motion i use this UDF for radial deformation DEFINE_GRID_MOTION(RADII, domain, dt, time, dtime) { Thread *tf = DT_THREAD (dt); face_t f; Node *v; real x_coor; real x_centre=0; real y_coor; real rayon; real freq = 3.183; real alpha; real temps; real R =0.25; int n; // FILE *fp = fopen(UDF_FILENAME, "r"); /* set deforming flag on adjacent cell zone */ SET_DEFORMING_THREAD_FLAG (THREAD_T0 (tf)); //temps = RP_Get_Real("flow-time"); temps = time+dtime; begin_f_loop (f, tf) { f_node_loop (f, tf, n) { v = F_NODE (f, tf, n); if (NODE_POS_NEED_UPDATE (v)) { NODE_POS_UPDATED (v); x_coor= NODE_X (v); y_coor= NODE_Y (v); alpha = 0.125; NODE_X (v) = x_coor- x_coor*(alpha*cos(2*M_PI*4*freq*temps)); NODE_Y (v) = y_coor- y_coor*(alpha*cos(2*M_PI*4*freq*temps)); } } } end_f_loop (f, tf); } for plunging motion i use zone motion udf DEFINE_ZONE_MOTION(motion_TRANS,omega,axis,origin, velocity,time,dtime) { /************************************************** ***/ real corde = 0.1; // real freq = 0.224; real freq = 4; real ho = 1; real ao = 0.84811*(M_PI/2); /*76.33°*/ real phy = M_PI/2; /************************************************** ***/ //if (time > (1/(4*freq))) //{ velocity[1] = (ho*corde*cos(2*M_PI*freq*(time+dtime)) - ho*corde*cos(2*M_PI*freq*time))/dtime; // origin[0] =0.0; // origin[1] = ho*corde*cos(2*M_PI*freq*(time+dtime)); //} return; } after simulation i generate animation file in TECPLOT i found Unphysical results my question is : are the combination between ZONE MOTION and GRID MOTION udfs appropriate for combined deflection translation motions or not please any help will be very appreciated |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Udf for 2D rotating cylinder | Rick 037 | Fluent UDF and Scheme Programming | 2 | August 31, 2020 08:59 |
How to implement time-varying boundary conditions in an enclosed cylinder | Alejandro-FA | OpenFOAM Pre-Processing | 0 | May 14, 2020 07:29 |
Rotating Cylinder UDF | greg_2666 | Fluent UDF and Scheme Programming | 0 | February 11, 2016 10:33 |
UDF for a rotationally oscillating cylinder | aria.ghiasi | Fluent UDF and Scheme Programming | 0 | December 21, 2015 07:20 |
UDF of natural convection in 3D cylinder | a_Sarlak | Fluent UDF and Scheme Programming | 0 | November 7, 2015 09:19 |