CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Tank sloshing in fluent diverging...need help

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

LinkBack Thread Tools Display Modes
Old   November 11, 2009, 05:26
Unhappy Tank sloshing in fluent diverging...need help
New Member
Join Date: Jul 2009
Location: India
Posts: 12
Rep Power: 8
lostin4ever is on a distinguished road
I am trying to simulate a sloshing tank problem using vof and 3ddp. To mimic the lateral motion of tank I have defined a udf of source term in motion direction (n my case it is Y dirn) . I have got results for a time period 0.44 and it seems to be okay. When I am trying to run he simulation for higher time period ( 1.33 ,1.25) my continuity residual is diverging after (1.8-2 sec ) simulation. here are my solver and other settings.

For different time I need to change my Udf time everytime that I am doing.I have attached my udf.

thanks for help in advance.

Reading "D:\tank\tanknearfirstmodet=1.25\tanknearfirstmode t=1.25.cas"...
61250 hexahedral cells, zone 2, binary.
9450 quadrilateral wall faces, zone 3, binary.
179025 quadrilateral interior faces, zone 5, binary.
66096 nodes, binary.
66096 node flags, binary.

Opening library "libudf"...
Library "libudf\ntx86\3ddp\libudf.dll" opened

fluid (phase-1)
tankwall (phase-1)
default-interior (phase-1)
fluid (phase-2)
tankwall (phase-2)
default-interior (phase-2)
shell conduction zones,

Version: 3d, dp, pbns, vof, lam, unsteady (3d, double precision, pressure-based, VOF, laminar, unsteady)
Release: 6.3.26


Model Settings
Space 3D
Time Unsteady, 1st-Order Implicit
Viscous Laminar
Heat Transfer Disabled
Solidification and Melting Disabled
Species Transport Disabled
Coupled Dispersed Phase Disabled
Pollutants Disabled
Pollutants Disabled
Soot Disabled

Boundary Conditions


name id type
fluid 2 fluid
tankwall 3 wall
default-interior 5 interior

Boundary Conditions


Condition Value
Material Name water-liquid
Specify source terms? yes
Source Terms ((x-momentum) (y-momentum ((profile udf cell_y_source::libudf) (inactive . #f) (constant . 0))) (z-momentum))
Specify fixed values? no
Local Coordinate System for Fixed Velocities no
Fixed Values ((x-velocity (inactive . #f) (constant . 0) (profile )) (y-velocity (inactive . #f) (constant . 0) (profile )) (z-velocity (inactive . #f) (constant . 0) (profile )))
Motion Type 0
X-Velocity Of Zone (m/s) 0
Y-Velocity Of Zone (m/s) 0
Z-Velocity Of Zone (m/s) 0
Rotation speed (rad/s) 0
X-Origin of Rotation-Axis (m) 0
Y-Origin of Rotation-Axis (m) 0
Z-Origin of Rotation-Axis (m) 0
X-Component of Rotation-Axis 0
Y-Component of Rotation-Axis 0
Z-Component of Rotation-Axis 1
Deactivated Thread no
Porous zone? no
Porosity 1


Condition Value
Enable shell conduction? no
Wall Motion 0
Shear Boundary Condition 0
Define wall motion relative to adjacent cell zone? yes
Apply a rotational velocity to this wall? no
Velocity Magnitude (m/s) 0
X-Component of Wall Translation 1
Y-Component of Wall Translation 0
Z-Component of Wall Translation 0
Define wall velocity components? no
X-Component of Wall Translation (m/s) 0
Y-Component of Wall Translation (m/s) 0
Z-Component of Wall Translation (m/s) 0
Rotation Speed (rad/s) 0
X-Position of Rotation-Axis Origin (m) 0
Y-Position of Rotation-Axis Origin (m) 0
Z-Position of Rotation-Axis Origin (m) 0
X-Component of Rotation-Axis Direction 0
Y-Component of Rotation-Axis Direction 0
Z-Component of Rotation-Axis Direction 1
X-component of shear stress (pascal) 0
Y-component of shear stress (pascal) 0
Z-component of shear stress (pascal) 0
Specularity Coefficient 0


Condition Value

Solver Controls


Equation Solved
Flow yes
Volume Fraction yes


Numeric Enabled
Absolute Velocity Formulation yes

Unsteady Calculation Parameters

Time Step (s) 0.0099999998
Max. Iterations Per Time Step 20

Non-Iterative Solver Factors

Variable Max. Corrections Correction Tolerance Residual Tolerance Relaxation Factor
Pressure 10 0.25 9.9999997e-05 1
Momentum 5 0.050000001 9.9999997e-05 1
Volume Fraction 5 0.050000001 9.9999997e-05 1

Linear Solver

Solver Termination Residual Reduction
Variable Type Criterion Tolerance
Pressure F-Cycle 0.1
X-Momentum Flexible 0.1 0.69999999
Y-Momentum Flexible 0.1 0.69999999
Z-Momentum Flexible 0.1 0.69999999
Volume Fraction Flexible 0.1 0.69999999

Pressure-Velocity Coupling

Parameter Value
Neighbour Correction 1

Discretization Scheme

Variable Scheme
Pressure PRESTO!
Momentum Second Order Upwind
Volume Fraction Second Order Upwind

Solution Limits

Quantity Limit
Minimum Absolute Pressure 1
Maximum Absolute Pressure 5e+10
Minimum Temperature 1
Maximum Temperature 5000
Attached Files
File Type: c tank_t=1.25.c (341 Bytes, 6 views)
lostin4ever is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Stopping a Fluent batch job AND saving the data! Possible? Volker Pawlik FLUENT 7 May 22, 2014 22:30
solving ocean wave with Fluent or CFX? gholamghar Main CFD Forum 1 March 21, 2009 13:49
Sloshing Boundary Conditions in Fluent Prabodh FLUENT 0 April 16, 2006 05:52
Integration of a Custom C++ Model into FLUENT Syed Haider FLUENT 1 February 20, 2006 05:49
Advanced Turbulence Modeling in Fluent, Realizable k-epsilon Model Jonas Larsson FLUENT 5 March 13, 2000 04:27

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