|
[Sponsors] |
April 12, 2017, 08:59 |
Melting of sphere in box
|
#1 |
New Member
up
Join Date: Oct 2015
Posts: 10
Rep Power: 11 |
Hi,
I am simulating the melting of sphere in a box in 2D. The problem Specifications are: I.C. T_melt=1873 K T_sphere=300 K B.C. T_wall=1873 K vel_interface=(k_f*(dT/dr)-h*(T_melt-T_mpshell))/(rho*delta_H) The sphere_melt interface is moving with the 2nd boundary condition using dynamic mesh. The UDF for for dynamic mesh : #include "udf.h" #include "math.h" #include "dynamesh_tools.h" #define k_f 35 #define h 2000 #define T_bath 1873 #define T_mpshell 1423 #define rho 7000 #define delta_H 275753 DEFINE_CG_MOTION(vel_interface,dt,vel,omega,time,d time) { Thread *t; cell_t c; real x[ND_ND]; face_t f; t=DT_THREAD(dt); begin_f_loop(f,t) { F_CENTROID(x,f,t); vel[0]=(k_f*C_T_G(c,t)[0]-h*(T_bath-T_mpshell))/(rho*delta_H); vel[1]=(k_f*C_T_G(c,t)[1]-h*(T_bath-T_mpshell))/(rho*delta_H); } end_f_loop(f,t) } When I running the solution I get error : received a fatal signal (Segmentation fault) Please help me to get rid of this error. Thanks in advance link for the problem image: https://drive.google.com/file/d/0B3Z...ew?usp=sharing |
|
Tags |
dynamic meshing, solidification/melting, udf and programming |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] Sphere in a channel by snappyHexMesh | arsalan.dryi | OpenFOAM Meshing & Mesh Conversion | 0 | November 14, 2016 17:01 |
[Other] An sphere at the center of a box | babakflame | OpenFOAM Meshing & Mesh Conversion | 0 | June 29, 2016 16:26 |
Natural Convection around two spheres in a box - chtMultiRegioSimpleFoam | salvo-K61IC | OpenFOAM | 4 | January 16, 2015 14:27 |
Immersol Simulation of a Heated Box | Dong | Phoenics | 0 | March 2, 2006 22:20 |
meshing F1 front wing | Steve | FLUENT | 0 | April 17, 2003 13:37 |