|
[Sponsors] |
July 10, 2012, 12:43 |
UDS Flux Function Problem
|
#1 |
New Member
Jun Li
Join Date: Jul 2012
Posts: 9
Rep Power: 14 |
Hey all,
Recently I want to use Fluent's user define scalars(uds) to solve the nucleation of equiax grain in a ingot cast problem with Fluent Eulerian Mutiphase Model. The control equation for the nucleation of equiaxed grain as follow: N/t+(uN)=source Therefor, for both the Unsteady Function and Flux Function should be defined by udf. My code is shown: DEFINE_UDS_UNSTEADY(unsteady_terms,cell_index,cell _thread,i,apu,su) { real cell_volume, physical_dt, RHO_cell, vof, vof_new, phi_old; cell_volume = C_VOLUME(cell_index,cell_thread); physical_dt = RP_Get_Real("physical-time-step"); RHO_cell = C_R(cell_index,cell_thread); vof = C_VOF_M1(cell_index,cell_thread); vof_new = C_VOF(cell_index,cell_thread); /* IMPLICIT PART OF THE TIME DERIVATIVE DISCRETIZATION */ *apu = -(cell_volume)/physical_dt; /* EXPLICIT PART OF THE TIME DERIVATIVE DISCRETIZATION */ phi_old = C_STORAGE_R(cell_index,cell_thread,SV_UDSI_M1(i)); *su = (cell_volume*phi_old)/physical_dt; } /*-------------------------*/ /* CONVECTIVE (FLUX) TERMS */ /*-------------------------*/ DEFINE_UDS_FLUX(convective_flux_terms,face_index,f ace_thread,i) { Thread *t0, *t1 = NULL; cell_t c0, c1 = -1; real NV_VEC(psi_vec), NV_VEC(A), vof, flux = 0.0, rho = 2; c0 = F_C0(face_index,face_thread); t0 = F_C0_THREAD(face_index,face_thread); vof = C_VOF(c0, t0); c0 = F_C0(face_index,face_thread); t0 = F_C0_THREAD(face_index,face_thread); /* if(PRINCIPAL_FACE_P(face_index,face_thread)) */ { F_AREA(A,face_index,face_thread); /* NORMAL VECTOR OF THE FACE */ if(BOUNDARY_FACE_THREAD_P(face_thread)) { NV_DS(psi_vec,=,F_U(face_index,face_thread),F_V(fa ce_index,face_thread),F_W(face_index,face_thread), *,1.0); flux = NV_DOT(psi_vec,A); /* FLUX THROUGH THE FACE */ } else { c1 = F_C1(face_index,face_thread); t1 = F_C1_THREAD(face_index,face_thread); NV_DS(psi_vec, =,C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1.0); NV_DS(psi_vec,+=,C_U(c1,t1),C_V(c1,t1),C_W(c1,t1), *,1.0); flux = NV_DOT(psi_vec,A)/2.0; /* AVERAGE FLUX THROUGH THE FACE */ } return (flux); } } Only momentum equation is solved. For these code I done two different cases: Case I, only the melt (liquid) region is considered. For this case the simulation result is reasonable; Case II, both the melt (liquid) and model (solid) are considered which means this is a liquid solid couple problem. For this case, the simulation cannot running. The error message emerged even at the first iteration: Error: FLUENT received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart to continue. 4. Report error to your distributor. Error Object: () I'd appreciate your help. Thanks in advance. |
|
July 11, 2012, 03:54 |
|
#2 |
New Member
Jun Li
Join Date: Jul 2012
Posts: 9
Rep Power: 14 |
|
|
July 11, 2012, 14:14 |
|
#3 |
New Member
Jun Li
Join Date: Jul 2012
Posts: 9
Rep Power: 14 |
|
|
July 11, 2012, 14:15 |
|
#4 |
New Member
Jun Li
Join Date: Jul 2012
Posts: 9
Rep Power: 14 |
Please help me!
|
|
July 13, 2012, 05:29 |
|
#5 |
New Member
Jun Li
Join Date: Jul 2012
Posts: 9
Rep Power: 14 |
|
|
September 19, 2013, 06:13 |
|
#6 |
Member
Georgy
Join Date: Apr 2011
Location: Russia
Posts: 32
Rep Power: 15 |
sunbird04,
Did you use Windows? I often see such message on the Windows OS. It looks like on the Linux OS, this error message appears less frequently. |
|
April 28, 2014, 03:58 |
|
#7 |
New Member
Tsaba Nagy
Join Date: Nov 2011
Location: Miskolc, Hungary
Posts: 10
Rep Power: 15 |
This message is usually for a UDM number problem (for example). If you have 5 UDM-s in your UDF, but in the Fluent only 4 is set, you will have this error message. The UDF and UDS numbers are started with zero, so if your last UDM is "C_UDMI(c,t,4)", you have 5 UDMs. Check all the setups. Usually you can see the name of the desired UDF, which is providing the error.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDS Flux Function Choice | Yang Chung | FLUENT | 1 | November 17, 2015 02:14 |
Problem with function Compute_Force_And_Moment | zedas | FLUENT | 1 | October 15, 2012 01:36 |
OpenFOAM static build on Cray XT5 | asaijo | OpenFOAM Installation | 9 | April 6, 2011 13:21 |
Help with UDS FLUX | Andrew Garrard | FLUENT | 2 | February 14, 2005 06:47 |
ODD: no UDS flux across interior boundary | Bob | FLUENT | 1 | February 25, 2002 20:15 |