|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
yahya rezvani
Join Date: Aug 2019
Posts: 18
Rep Power: 7 ![]() |
Hi everyone!
I'm trying to simulate incompressible material flow in friction welding. im use udf for viscosity and heat transfer in pin surface. Boundary conditions: velocity inlet, pressure outlet, slip condition walls. heat input by udf on copled via system method in surface of pin. energy equation is on. pressure base solver used in steady state. After 3 iterations, Divergence detected in AMG solver: pressure coupled Stabilizing temperature to enhance linear solver robustness. temperature limited to 1.000000e+00 in 1695 cells on zone 4 in domain 1 temperature limited to 5.000000e+03 in 3961 cells on zone 4 in domain 1 Divergence detected in AMG solver: pressure coupled Error at host: floati Code:
iter continuity x-velocity y-velocity z-velocity energy time/iter temperature limited to 2.500000e+02 in 2034 cells on zone 4 in domain 1 temperature limited to 8.000000e+03 in 2 cells on zone 4 in domain 1 1 1.0000e+00 2.5067e-03 3.3226e-03 6.4718e-03 3.3263e-03 0:04:18 19 Stabilizing pressure coupled to enhance linear solver robustness. temperature limited to 2.500000e+02 in 111 cells on zone 4 in domain 1 temperature limited to 8.000000e+03 in 81 cells on zone 4 in domain 1 2 1.0000e+00 1.1354e-03 1.5119e-03 2.2842e-03 5.9263e-03 0:06:37 18 Stabilizing pressure coupled to enhance linear solver robustness. Stabilizing pressure coupled using GMRES to enhance linear solver robustness. Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Stabilizing pressure coupled to enhance linear solver robustness. temperature limited to 2.500000e+02 in 3164 cells on zone 4 in domain 1 temperature limited to 8.000000e+03 in 583 cells on zone 4 in domain 1 3 1.0000e+00 8.4188e-03 3.7937e-03 1.3550e-02 5.5499e+00 0:14:58 17 E Error at Node 1: floating point excreptionror at Node 0: floating point exception Stabilizing pressure coupled to enhance linear solver robustness. Stabilizing pressure coupled using GMRES to enhance linear solver robustness. Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Experiencing convergence difficulties - temporarily relaxing and trying again... Stabilizing pressure coupled to enhance linear solver robustness. Stabilizing pressure coupled using GMRES to enhance linear solver robustness. Experiencing convergence difficulties - temporarily relaxing and trying again... Divergence detected in AMG solver: pressure coupled Stabilizing temperature to enhance linear solver robustness. temperature limited to 2.500000e+02 in 900 cells on zone 4 in domain 1 temperature limited to 8.000000e+03 in 2187 cells on zone 4 in domain 1 Divergence detected in AMG solver: pressure coupled Error at host: floating point exception Best Regards. Last edited by rezvani; January 29, 2021 at 07:35. Reason: attach udf. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 ![]() ![]() |
there is no way to help you based on information you've provided
show your UDF if you think problem comes from it change variables inside UDF to constant values and try to run case again. if it doesn't work, it means, problem is in your settings not in UDF
__________________
best regards ****************************** press LIKE if this message was helpful |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Join Date: Nov 2013
Posts: 1,965
Rep Power: 26 ![]() |
Initialize, don't run, plot the viscosity. Is it what you expected?
|
|
![]() |
![]() |
![]() |
![]() |
#4 | |
New Member
yahya rezvani
Join Date: Aug 2019
Posts: 18
Rep Power: 7 ![]() |
Quote:
divergence may be produced by heat generation udf. temperature (heat flux) dont travel into flow Code:
#include "udf.h" real x[ND_ND], A[ND_ND], c[ND_ND], rcal[ND_ND], Atotal[ND_ND]; /* this will hold the position vector */ face_t f; real r, q, qt, y; DEFINE_PROFILE(heatfx,t,i) { int abst; Atotal[0]=0; Atotal[1]= 0; Atotal[2]=0; c[0]=0; c[1]= 0; c[2]=0.006; qt=0; begin_f_loop(f,t) { F_CENTROID(x,f,t); NV_VV(rcal,=,x,-,c); r=NV_MAG(rcal); if(r<0.0059&&r>0) { q=(3*2272*r)/(2*3.14*(pow(0.0059,3))); /*//w/m2*/ F_PROFILE(f,t,i)=q; F_AREA(A,f,t); qt=qt+q*NV_MAG(A); /*printf("qAreas=%f%f%f%f\n",q,A[0],A[1],A[2]); */ NV_VV(Atotal,=,Atotal,+,A); } else { F_PROFILE(f,t,i)=-30*(F_T(f,t)-298);/*//w/m2*/ y=F_T(f,t); } } end_f_loop(f,t) /*printf("\nQtdiff=%f\n",qt); */ } Last edited by rezvani; January 27, 2021 at 12:03. |
||
![]() |
![]() |
![]() |
![]() |
#5 | |
Senior Member
Join Date: Nov 2013
Posts: 1,965
Rep Power: 26 ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
yahya rezvani
Join Date: Aug 2019
Posts: 18
Rep Power: 7 ![]() |
this is viscosity of solid metal in 350 k . i try chose higer temperature like 800k for smaller viscosity.
Last edited by rezvani; January 29, 2021 at 07:33. |
|
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
Alexander
Join Date: Apr 2013
Posts: 2,363
Rep Power: 34 ![]() ![]() |
your code seems to be correct.
you may try to check values, which you apply. you may try to simplify your case for testing, until you will get what you expect you may try to use UDMs to store and plot your source and material properties
__________________
best regards ****************************** press LIKE if this message was helpful |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Generate convective pressure fluctuation | Bananenflanke | CFX | 10 | May 12, 2021 18:33 |
PEMFC model with FLUENT | brahimchoice | FLUENT | 22 | April 19, 2020 15:44 |
sonicFoam - pressure driven pipe: flow continuity violation and waveTransmissive BC | Endel | OpenFOAM Running, Solving & CFD | 3 | September 11, 2014 16:29 |
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues | michele | OpenFOAM Meshing & Mesh Conversion | 2 | July 15, 2005 04:15 |
Hydrostatic pressure in 2-phase flow modeling (long) | DS & HB | Main CFD Forum | 0 | January 8, 2000 15:00 |