CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Low Mixing time Problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 15, 2013, 11:21
Default Low Mixing time Problem
  #1
New Member
 
cholaole@gmail.com
Join Date: Mar 2011
Posts: 4
Rep Power: 15
Mavier is on a distinguished road
Hello everybody.

I am trying to reproduce the conditions form this paper: Verification and application of a mathematical model for the assessment of the effect of guiding walls on the hydraulic efficiency of chlorination tanks.

First i ran a steady-state simulation to serve as the starting point for my transient analisys.
I´m using the standard k-epsilon model.
I created a NaCl material with its standard characteristics.
The top of the tank is a symmetry boundary condition.
I create two separete simulations for the transient analisys: firste is the mixture (water+NaCl) entering my storage tank for 3 minutes (i use a 1s timestep); the second part is just water entering the domain (mass fraction of NaCl = 0, water as constraint).
I have to outlet pipes that i have put opening boundary to conditions (to avoid the wall problem at the outlet), both with medium intensity (5%) and opening Pres. and Dirn = 0 Pa
However the mixture is exiting the tank faster than expected, like the following graph shows:

My results are the ones in green (with a time step equal to 1s in the second part of the simulation); and the purple one, time step equal to 4s at the second part of the simulation.

where:
C = mean NaCl concentration at the outlet
Co = average tracer concentration
t = current time
T = theoretical detention time

I really would like some help if i should change my turbulence model or if i should try perhaps a structured mesh (i generated mine with the ansys mesh generator using inflation at the walls) or a different time step.

Here are my both CCL for the transient analysis:
Part one (with NaCl):
&replace FLOW: Flow Analysis 1
ANALYSIS TYPE:
Option = Transient
EXTERNAL SOLVER COUPLING:
Option = None
END
INITIAL TIME:
Option = Value
Time = 0 [s]
END
TIME DURATION:
Option = Total Time
Total Time = 180 [s]
END
TIME STEPS:
Option = Timesteps
Timesteps = 1 [s]
END
END
DOMAIN: Default Domain
Coord Frame = Coord 0
Domain Type = Fluid
Location = B41
BOUNDARY: Default Domain Default
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = F294.41,F35.41,F36.41,F37.41,F38.41,F43.41,F44.41, F45.41,F46.41,F47.41,F48.41,F49.41
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: inlet
Boundary Type = INLET
Interface Boundary = Off
Location = F293.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0.00234
Option = Mass Fraction
END
FLOW DIRECTION:
Option = Normal to Boundary Condition
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Mass Flow Rate = 752.4 [kg s^-1]
Option = Mass Flow Rate
END
TURBULENCE:
Option = High Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: outlet_1
Boundary Type = OPENING
Interface Boundary = Off
Location = F219.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0.0
Option = Mass Fraction
END
FLOW DIRECTION:
Option = Normal to Boundary Condition
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Option = Opening Pressure and Direction
Relative Pressure = 0 [Pa]
END
TURBULENCE:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: outlet_2
Boundary Type = OPENING
Interface Boundary = Off
Location = F220.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0.0
Option = Mass Fraction
END
FLOW DIRECTION:
Option = Normal to Boundary Condition
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Option = Opening Pressure and Direction
Relative Pressure = 0 [Pa]
END
TURBULENCE:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: simetria
Boundary Type = SYMMETRY
Location = F42.41
END
DOMAIN MODELS:
BUOYANCY MODEL:
Option = Non Buoyant
END
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
REFERENCE PRESSURE:
Reference Pressure = 1 [atm]
END
END
FLUID DEFINITION: Dentro do Reservatorio
Material = Water NaCl
Option = Material Library
MORPHOLOGY:
Option = Continuous Fluid
END
END
FLUID MODELS:
COMBUSTION MODEL:
Option = None
END
COMPONENT: Sodium Chloride
Kinematic Diffusivity = 1.5e-9 [m^2 s^-1]
Option = Transport Equation
END
COMPONENT: Water
Option = Constraint
END
HEAT TRANSFER MODEL:
Fluid Temperature = 25 [C]
Option = Isothermal
END
THERMAL RADIATION MODEL:
Option = None
END
TURBULENCE MODEL:
Option = k epsilon
END
TURBULENT WALL FUNCTIONS:
Option = Scalable
END
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic
END
COMPONENT: Sodium Chloride
Mass Fraction = 0.0
Option = Automatic with Value
END
STATIC PRESSURE:
Option = Automatic
END
TURBULENCE INITIAL CONDITIONS:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
OUTPUT CONTROL:
MONITOR OBJECTS:
MONITOR BALANCES:
Option = Full
END
MONITOR FORCES:
Option = Full
END
MONITOR PARTICLES:
Option = Full
END
MONITOR POINT: Concentracao_Outlet1
Expression Value = ConcOutlet1
Option = Expression
END
MONITOR POINT: Concentracao_Outlet2
Expression Value = ConcOutlet2
Option = Expression
END
MONITOR POINT: Mass_fraction_outlet1
Expression Value = MF1
Option = Expression
END
MONITOR POINT: Mass_fraction_outlet2
Expression Value = MF2
Option = Expression
END
MONITOR RESIDUALS:
Option = Full
END
MONITOR TOTALS:
Option = Full
END
END
RESULTS:
File Compression Level = Default
Option = Standard
END
TRANSIENT RESULTS: Concentracoes
File Compression Level = Default
Include Mesh = No
Option = Selected Variables
Output Variables List = Velocity,Velocity u,Velocity v,Velocity w,Water.Mass Concentration,Water.Mass Fraction,Sodium Chloride.Mass Fraction,Sodium Chloride.Mass Concentration
OUTPUT FREQUENCY:
Option = Every Timestep
END
END
END
SOLUTION UNITS:
Angle Units = [rad]
Length Units = [m]
Mass Units = [kg]
Solid Angle Units = [sr]
Temperature Units = [K]
Time Units = [s]
END
SOLVER CONTROL:
Turbulence Numerics = High Resolution
ADVECTION SCHEME:
Option = High Resolution
END
CONVERGENCE CONTROL:
Maximum Number of Coefficient Loops = 10
Minimum Number of Coefficient Loops = 1
Timescale Control = Coefficient Loops
END
CONVERGENCE CRITERIA:
Residual Target = 0.00001
Residual Type = RMS
END
EQUATION CLASS: continuity
ADVECTION SCHEME:
Option = High Resolution
END
END
TRANSIENT SCHEME:
Option = Second Order Backward Euler
TIMESTEP INITIALISATION:
Option = Automatic
END
END
END
END


Part 2:
&replace FLOW: Flow Analysis 1
ANALYSIS TYPE:
Option = Transient
EXTERNAL SOLVER COUPLING:
Option = None
END
INITIAL TIME:
Option = Value
Time = 181 [s]
END
TIME DURATION:
Option = Total Time
Total Time = 5000 [s]
END
TIME STEPS:
Option = Timesteps
Timesteps = 1 [s]
END
END
DOMAIN: Default Domain
Coord Frame = Coord 0
Domain Type = Fluid
Location = B41
BOUNDARY: Default Domain Default
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = F294.41,F35.41,F36.41,F37.41,F38.41,F43.41,F44.41, F45.41,F46.41,F47.41,F48.41,F49.41
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: inlet
Boundary Type = INLET
Interface Boundary = Off
Location = F293.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0
Option = Mass Fraction
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Normal Speed = 2.66 [m s^-1]
Option = Normal Speed
END
TURBULENCE:
Option = High Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: outlet_1
Boundary Type = OPENING
Interface Boundary = Off
Location = F219.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0.00026494
Option = Mass Fraction
END
FLOW DIRECTION:
Option = Normal to Boundary Condition
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Option = Opening Pressure and Direction
Relative Pressure = 0 [Pa]
END
TURBULENCE:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: outlet_2
Boundary Type = OPENING
Interface Boundary = Off
Location = F220.41
BOUNDARY CONDITIONS:
COMPONENT: Sodium Chloride
Mass Fraction = 0.00028418
Option = Mass Fraction
END
FLOW DIRECTION:
Option = Normal to Boundary Condition
END
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Option = Opening Pressure and Direction
Relative Pressure = 0 [Pa]
END
TURBULENCE:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
BOUNDARY: simetria
Boundary Type = SYMMETRY
Location = F42.41
END
DOMAIN MODELS:
BUOYANCY MODEL:
Option = Non Buoyant
END
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
REFERENCE PRESSURE:
Reference Pressure = 1 [atm]
END
END
FLUID DEFINITION: Dentro do Reservatorio
Material = Water NaCl
Option = Material Library
MORPHOLOGY:
Option = Continuous Fluid
END
END
FLUID MODELS:
COMBUSTION MODEL:
Option = None
END
COMPONENT: Sodium Chloride
Kinematic Diffusivity = 1.5e-9 [m^2 s^-1]
Option = Transport Equation
END
COMPONENT: Water
Option = Constraint
END
HEAT TRANSFER MODEL:
Fluid Temperature = 25 [C]
Option = Isothermal
END
THERMAL RADIATION MODEL:
Option = None
END
TURBULENCE MODEL:
Option = k epsilon
END
TURBULENT WALL FUNCTIONS:
Option = Scalable
END
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic
END
COMPONENT: Sodium Chloride
Mass Fraction = 0.0
Option = Automatic with Value
END
STATIC PRESSURE:
Option = Automatic
END
TURBULENCE INITIAL CONDITIONS:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
END
OUTPUT CONTROL:
MONITOR OBJECTS:
MONITOR BALANCES:
Option = Full
END
MONITOR FORCES:
Option = Full
END
MONITOR PARTICLES:
Option = Full
END
MONITOR POINT: Concentracao_Outlet1
Expression Value = ConcOutlet1
Option = Expression
END
MONITOR POINT: Concentracao_Outlet2
Expression Value = ConcOutlet2
Option = Expression
END
MONITOR POINT: Mass_fraction_outlet1
Expression Value = MF1
Option = Expression
END
MONITOR POINT: Mass_fraction_outlet2
Expression Value = MF2
Option = Expression
END
MONITOR RESIDUALS:
Option = Full
END
MONITOR TOTALS:
Option = Full
END
END
RESULTS:
File Compression Level = Default
Option = Standard
END
TRANSIENT RESULTS: Concentracoes
File Compression Level = Default
Include Mesh = No
Option = Selected Variables
Output Variables List = Velocity,Velocity u,Velocity v,Velocity w,Water.Mass Concentration,Water.Mass Fraction,Sodium Chloride.Mass Fraction,Sodium Chloride.Mass Concentration
OUTPUT FREQUENCY:
Option = Every Timestep
END
END
END
SOLUTION UNITS:
Angle Units = [rad]
Length Units = [m]
Mass Units = [kg]
Solid Angle Units = [sr]
Temperature Units = [K]
Time Units = [s]
END
SOLVER CONTROL:
Turbulence Numerics = High Resolution
ADVECTION SCHEME:
Option = High Resolution
END
CONVERGENCE CONTROL:
Maximum Number of Coefficient Loops = 10
Minimum Number of Coefficient Loops = 1
Timescale Control = Coefficient Loops
END
CONVERGENCE CRITERIA:
Residual Target = 0.00001
Residual Type = RMS
END
EQUATION CLASS: continuity
ADVECTION SCHEME:
Option = High Resolution
END
END
TRANSIENT SCHEME:
Option = Second Order Backward Euler
TIMESTEP INITIALISATION:
Option = Automatic
END
END
END
END




Thank you all in advance. Best regards

Last edited by Mavier; April 15, 2013 at 12:09.
Mavier is offline   Reply With Quote

Old   April 15, 2013, 20:05
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Can you confirm no multi-phase stuff happens - so no droplets, bubbles, interfaces or anything like that? Both the wafer and NaCl dissolve in each other.

Your results are close so you do not seem to be too far off. The FAQ discusses general accuracy issues: http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
ghorrocks is offline   Reply With Quote

Old   April 16, 2013, 17:04
Default
  #3
New Member
 
cholaole@gmail.com
Join Date: Mar 2011
Posts: 4
Rep Power: 15
Mavier is on a distinguished road
Thank you for the reply Glenn.

My mixing conditions are the following:

I created a PURE SUBSTANCE for the NaCl (that i called Sodium Chloride) with the following values:
Molar Mass: 58.44 g/mol
Density: 2.165 g/cm³

Specific Heat Capacity: 629.53 J/(kg.K)

Then i created a VARIABLE COMPOSITION MIXTURE that consists of:
Water + Sodium Chloride as na Ideal Mixture that i called "Water NaCl".

Then at the domain i used this mixture of "Water NaCl"
In the componet i specified the following:
Sodium Chlorides as a Transport Equation with kinematic diffusity of: 1.5e-9 m²/s (the NaCl diffusity in water at 25°C).
Water as constraint.

At the inlet i´m using a Mass Fraction for the Sodium Chloride of: 0.00234

Turbulence is High intensity (10%)


There are a couple of things that the autor specifies in the paper that i don´t know where to put it, he mentions the following: The turbulent energy (k) and its dissipation rate (e) are assumed to be uniform, with values corresponding to an eddy viscosity at the inlet approximately 100 times the molecular viscosity of water. The tracer concentration, C, is assumed to be uniformly distributed at the inlet. At the outlet pipes the pressure is specified and the vertical gradients of k, e and C are set equal to zero.

Is thera a palce in CFX that i can change the values for the turbulent energy, and dissipation rate.

Anyway, thanks for the response again Glenn, best regards.
Mavier is offline   Reply With Quote

Old   April 16, 2013, 20:27
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The assumption that the k and e are uniform seems like a brave one to me. The whole point of CFD is to model these terms so you do not need to make assumptions like this. But if you want to reproduce their results then you can probably use a laminar flow model and increase the viscosity to be the sum of the molecular and turbulent viscosities. You would also need to increase the diffusivities of the heat and mass fraction equations as well.
ghorrocks is offline   Reply With Quote

Old   April 29, 2013, 00:42
Default
  #5
New Member
 
cholaole@gmail.com
Join Date: Mar 2011
Posts: 4
Rep Power: 15
Mavier is on a distinguished road
Thanlks for the answer Glenn.

I tried some alternatives here but the concentration is still leaving the tank early, do you think this could be mesh realted Glenn ?
I´m using a non-structured mesh, perhaps i should try using a structured one from icem.

Also, i tried to simulate it using the SST turbulence model, but i couldn´t get it to converge during the steady state simulation and i think this this could also be due to the unstructured mesh that i´m using.

Best Regards.
Mavier is offline   Reply With Quote

Old   April 29, 2013, 01:00
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,655
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
FAQ1: http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
FAQ2: http://www.cfd-online.com/Wiki/Ansys...gence_criteria
ghorrocks is offline   Reply With Quote

Reply

Tags
cfx, mixing, mixing time, water storage tank

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
TimeVaryingMappedFixedValue irishdave OpenFOAM Running, Solving & CFD 32 June 16, 2021 07:55
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Adaptive time stepping problem for periodic problem daysley FLUENT 7 April 14, 2013 13:27
How to find a good time scale strategy? StefanG CFX 19 June 8, 2012 09:41
convergence problem when use pisoFoam, LES for wind tunnel case Forrest_Lei OpenFOAM 3 July 19, 2011 07:00


All times are GMT -4. The time now is 06:13.