CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   CFX (http://www.cfd-online.com/Forums/cfx/)
-   -   Error finding variable "THERMX" (http://www.cfd-online.com/Forums/cfx/115029-error-finding-variable-thermx.html)

sunilpatil March 22, 2013 04:25

Error finding variable "THERMX"
 
Hello Sir,
In my Transient Simulation of conjugate heat transfer problem i want to subtract heat transfer by convection (Q convection) from total radiation(Q) The equation is
Q on wall = alpa*Q incident-Q convection-Q night sky
The CEL Expression which i used is

LIBRARY:
CEL:
&replace EXPRESSIONS:
Dry Bulb = DRY BULB TEMP(t)
Tlower = areaAve(T)@WALL_SIDE1
Tupper = areaAve(T)@WALL_SIDE2
K = 0.0261 [W m^-1 K^-1]
alpa = 0.2
Kinematic Vis = 0.000016036 [m^2 s^-1]
L = 1.075 [m]
Pr = 0.707
q incident = GLOBAL RADIATION(t)
q night sky = night sky radiation(t)
Grashof = (g*beta*(Tupper-Dry Bulb)*L^3)/(Kinematic Vis^2)
Ray = Grashof*Pr
Nu L L = 0.27*(Ray^0.25)
Nu L U = 0.54*(Ray^0.25)
Nusselt = if(TSIDE1 > TSIDE2,Nu L L,Nu L U)
h bar = (Nusselt*K)/(L)
q con = h bar*(Tupper-Dry Bulb)
q Net = (alpa*q incident Cealing)-q con Cealing-q night sky
END
END
END



but before starting run it is giving me error


ERROR
| ERROR #001100279 has occurred in subroutine ErrAction. |
| Message: |
| NAME_MOD: Error finding variable "THERMX"


how to fix this problem sir
Thank you sir

ghorrocks March 22, 2013 05:10

CFX error messages can be cryptic sometimes. Have you specified a variabel called THERMX? You will have to post your whole CCL for us to help I suspect. Also a diagram of what you are modelling would help.

sunilpatil March 22, 2013 09:35

Error finding variable "THERMX"
 
Hello sir,
Thank you for your Reply,i have not used variable called THERMX in CFX pre variable in Present case i am just applying heat flux to top ceiling just to check weather the expression is working fine or not( my main aim is to subtract convection losses from global radiation with the help of expression. i have experimental values of global radiation)


CCL



&replace FLOW: Flow Analysis 1
ANALYSIS TYPE:
Option = Transient
EXTERNAL SOLVER COUPLING:
Option = None
END
INITIAL TIME:
Option = Automatic with Value
Time = 0 [s]
END
TIME DURATION:
Option = Total Time
Total Time = 86400 [s]
END
TIME STEPS:
Option = Timesteps
Timesteps = 120 [s]
END
END
DOMAIN INTERFACE: CEALING_ROOM
Boundary List1 = CEALING_ROOM Side 1
Boundary List2 = CEALING_ROOM Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = CEALING
Interface Region List1 = ROOM_CEILING
Interface Region List2 = CEILING_ROOM
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: CEALING_WALL
Boundary List1 = CEALING_WALL Side 1
Boundary List2 = CEALING_WALL Side 2
Filter Domain List1 = CEALING
Filter Domain List2 = WALL
Interface Region List1 = CEILING_WALL
Interface Region List2 = WALL_CEILING
Interface Type = Solid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: SLABE_ROOM
Boundary List1 = SLABE_ROOM Side 1
Boundary List2 = SLABE_ROOM Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = SLABE
Interface Region List1 = ROOM_BASE
Interface Region List2 = SLABE_ROOM
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: SLABE_WALL
Boundary List1 = SLABE_WALL Side 1
Boundary List2 = SLABE_WALL Side 2
Filter Domain List1 = SLABE
Filter Domain List2 = WALL
Interface Region List1 = SLABE_WALL
Interface Region List2 = WALL_SLABE
Interface Type = Solid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: WALL_ROOM_1
Boundary List1 = WALL_ROOM_1 Side 1
Boundary List2 = WALL_ROOM_1 Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = WALL
Interface Region List1 = ROOM_WALL_1,ROOM_WALL_2,ROOM_WALL_3,ROOM_WALL_4
Interface Region List2 = WALL_ROOM_1,WALL_ROOM_2,WALL_ROOM_3,WALL_ROOM_4
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN: CEALING
Coord Frame = Coord 0
Domain Type = Solid
Location = CEILING
BOUNDARY: Adiabatic
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = CEILING_CAVITY
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: CEALING_ROOM Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = CEILING_ROOM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: CEALING_WALL Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = CEILING_WALL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
END
END
BOUNDARY: TOP
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = CEILING_ATM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Heat Flux in = q con
Option = Heat Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1
Option = Opaque
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
SOLID DEFINITION: Solid 1
Material = CONCRETE
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Number of Histories = 10000
Option = Monte Carlo
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
END
END
DOMAIN: ROOM
Coord Frame = Coord 0
Domain Type = Fluid
Location = ROOM
BOUNDARY: CEALING_ROOM Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_CEILING
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: SLABE_ROOM Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_BASE
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: WALL_ROOM_1 Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_WALL_1,ROOM_WALL_2,ROOM_WALL_3,ROOM_WALL_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
DOMAIN MODELS:
BUOYANCY MODEL:
Buoyancy Reference Temperature = 25 [C]
Gravity X Component = 0 [m s^-2]
Gravity Y Component = -9.81 [m s^-2]
Gravity Z Component = 0 [m s^-2]
Option = Buoyant
BUOYANCY REFERENCE LOCATION:
Option = Automatic
END
END
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
REFERENCE PRESSURE:
Reference Pressure = 1 [atm]
END
END
FLUID DEFINITION: Fluid 1
Material = Air at 25 C
Option = Material Library
MORPHOLOGY:
Option = Continuous Fluid
END
END
FLUID MODELS:
COMBUSTION MODEL:
Option = None
END
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Option = Discrete Transfer
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
TURBULENCE MODEL:
Option = k epsilon
BUOYANCY TURBULENCE:
Option = None
END
END
TURBULENT WALL FUNCTIONS:
Option = Scalable
END
END
END
DOMAIN: SLABE
Coord Frame = Coord 0
Domain Type = Solid
Location = SLABE
BOUNDARY: ADIABATIC_WALL
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = SLABE_SIDE_WALLS,SLABE_SOIL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: SLABE_ROOM Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = SLABE_ROOM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: SLABE_WALL Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = SLABE_WALL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
SOLID DEFINITION: Solid 1
Material = CONCRETE
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Number of Histories = 10000
Option = Monte Carlo
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
END
END
DOMAIN: WALL
Coord Frame = Coord 0
Domain Type = Solid
Location = WALL
BOUNDARY: CEALING_WALL Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_CEILING
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: EAST_WALL
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_1
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: North Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_2
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: SLABE_WALL Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_SLABE
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: South Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: WALL_ROOM_1 Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_ROOM_1,WALL_ROOM_2,WALL_ROOM_3,WALL_ROOM_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: West_Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_3
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
TEMPERATURE:
Option = Automatic with Value
Temperature = 273 [K]
END
END
END
SOLID DEFINITION: Solid 1
Material = Brick Common
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Option = None
END
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = 0 [m s^-1]
END
RADIATION INTENSITY:
Option = Automatic
END
STATIC PRESSURE:
Option = Automatic with Value
Relative Pressure = 0 [Pa]
END
TEMPERATURE:
Option = Automatic with Value
Temperature = 298 [K]
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 RESIDUALS:
Option = Full
END
MONITOR TOTALS:
Option = Full
END
END
RESULTS:
File Compression Level = Default
Option = Standard
END
TRANSIENT RESULTS: Transient Results 1
File Compression Level = Default
Option = Standard
OUTPUT FREQUENCY:
Option = Time Interval
Time Interval = 1000 [s]
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 = 5
Minimum Number of Coefficient Loops = 1
Timescale Control = Coefficient Loops
END
CONVERGENCE CRITERIA:
Residual Target = 0.00001
Residual Type = RMS
END
TRANSIENT SCHEME:
Option = Second Order Backward Euler
TIMESTEP INITIALISATION:
Option = Automatic
END
END
END
END
#################################################





but if we apply Ra = constant in expression it will run fine but that is not a actual case

Thank you sir

sunilpatil March 22, 2013 09:37

Error finding variable "THERMX"
 
Hello sir,
Thank you for your Reply,i have not used variable called THERMX in CFX pre variable in Present case i am just applying heat flux to top ceiling just to check weather the expression is working fine or not( my main aim is to subtract convection losses from global radiation with the help of expression. i have experimental values of global radiation)
Image
http://postimg.org/image/wof4tkgoh/

CCL



&replace FLOW: Flow Analysis 1
ANALYSIS TYPE:
Option = Transient
EXTERNAL SOLVER COUPLING:
Option = None
END
INITIAL TIME:
Option = Automatic with Value
Time = 0 [s]
END
TIME DURATION:
Option = Total Time
Total Time = 86400 [s]
END
TIME STEPS:
Option = Timesteps
Timesteps = 120 [s]
END
END
DOMAIN INTERFACE: CEALING_ROOM
Boundary List1 = CEALING_ROOM Side 1
Boundary List2 = CEALING_ROOM Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = CEALING
Interface Region List1 = ROOM_CEILING
Interface Region List2 = CEILING_ROOM
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: CEALING_WALL
Boundary List1 = CEALING_WALL Side 1
Boundary List2 = CEALING_WALL Side 2
Filter Domain List1 = CEALING
Filter Domain List2 = WALL
Interface Region List1 = CEILING_WALL
Interface Region List2 = WALL_CEILING
Interface Type = Solid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: SLABE_ROOM
Boundary List1 = SLABE_ROOM Side 1
Boundary List2 = SLABE_ROOM Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = SLABE
Interface Region List1 = ROOM_BASE
Interface Region List2 = SLABE_ROOM
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: SLABE_WALL
Boundary List1 = SLABE_WALL Side 1
Boundary List2 = SLABE_WALL Side 2
Filter Domain List1 = SLABE
Filter Domain List2 = WALL
Interface Region List1 = SLABE_WALL
Interface Region List2 = WALL_SLABE
Interface Type = Solid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN INTERFACE: WALL_ROOM_1
Boundary List1 = WALL_ROOM_1 Side 1
Boundary List2 = WALL_ROOM_1 Side 2
Filter Domain List1 = ROOM
Filter Domain List2 = WALL
Interface Region List1 = ROOM_WALL_1,ROOM_WALL_2,ROOM_WALL_3,ROOM_WALL_4
Interface Region List2 = WALL_ROOM_1,WALL_ROOM_2,WALL_ROOM_3,WALL_ROOM_4
Interface Type = Fluid Solid
INTERFACE MODELS:
Option = General Connection
FRAME CHANGE:
Option = None
END
HEAT TRANSFER:
Option = Conservative Interface Flux
HEAT TRANSFER INTERFACE MODEL:
Option = None
END
END
PITCH CHANGE:
Option = None
END
END
MESH CONNECTION:
Option = Automatic
END
END
DOMAIN: CEALING
Coord Frame = Coord 0
Domain Type = Solid
Location = CEILING
BOUNDARY: Adiabatic
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = CEILING_CAVITY
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: CEALING_ROOM Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = CEILING_ROOM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: CEALING_WALL Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = CEILING_WALL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
END
END
BOUNDARY: TOP
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = CEILING_ATM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Heat Flux in = q con
Option = Heat Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1
Option = Opaque
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
SOLID DEFINITION: Solid 1
Material = CONCRETE
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Number of Histories = 10000
Option = Monte Carlo
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
END
END
DOMAIN: ROOM
Coord Frame = Coord 0
Domain Type = Fluid
Location = ROOM
BOUNDARY: CEALING_ROOM Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_CEILING
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: SLABE_ROOM Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_BASE
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
BOUNDARY: WALL_ROOM_1 Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = ROOM_WALL_1,ROOM_WALL_2,ROOM_WALL_3,ROOM_WALL_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
MASS AND MOMENTUM:
Option = No Slip Wall
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
END
DOMAIN MODELS:
BUOYANCY MODEL:
Buoyancy Reference Temperature = 25 [C]
Gravity X Component = 0 [m s^-2]
Gravity Y Component = -9.81 [m s^-2]
Gravity Z Component = 0 [m s^-2]
Option = Buoyant
BUOYANCY REFERENCE LOCATION:
Option = Automatic
END
END
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
REFERENCE PRESSURE:
Reference Pressure = 1 [atm]
END
END
FLUID DEFINITION: Fluid 1
Material = Air at 25 C
Option = Material Library
MORPHOLOGY:
Option = Continuous Fluid
END
END
FLUID MODELS:
COMBUSTION MODEL:
Option = None
END
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Option = Discrete Transfer
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
TURBULENCE MODEL:
Option = k epsilon
BUOYANCY TURBULENCE:
Option = None
END
END
TURBULENT WALL FUNCTIONS:
Option = Scalable
END
END
END
DOMAIN: SLABE
Coord Frame = Coord 0
Domain Type = Solid
Location = SLABE
BOUNDARY: ADIABATIC_WALL
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = SLABE_SIDE_WALLS,SLABE_SOIL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: SLABE_ROOM Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = SLABE_ROOM
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 0.8
Option = Opaque
END
END
END
BOUNDARY: SLABE_WALL Side 1
Boundary Type = INTERFACE
Interface Boundary = On
Location = SLABE_WALL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
THERMAL RADIATION:
Diffuse Fraction = 1.
Emissivity = 1.
Option = Opaque
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
SOLID DEFINITION: Solid 1
Material = CONCRETE
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Number of Histories = 10000
Option = Monte Carlo
Radiation Transfer Mode = Participating Media
SCATTERING MODEL:
Option = None
END
SPECTRAL MODEL:
Option = Gray
END
END
END
END
DOMAIN: WALL
Coord Frame = Coord 0
Domain Type = Solid
Location = WALL
BOUNDARY: CEALING_WALL Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_CEILING
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: EAST_WALL
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_1
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: North Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_2
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: SLABE_WALL Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_SLABE
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: South Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
BOUNDARY: WALL_ROOM_1 Side 2
Boundary Type = INTERFACE
Interface Boundary = On
Location = WALL_ROOM_1,WALL_ROOM_2,WALL_ROOM_3,WALL_ROOM_4
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Conservative Interface Flux
END
END
END
BOUNDARY: West_Wall
Boundary Type = WALL
Create Other Side = Off
Interface Boundary = Off
Location = WALL_CAVITY_3
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Option = Adiabatic
END
END
END
DOMAIN MODELS:
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
TEMPERATURE:
Option = Automatic with Value
Temperature = 273 [K]
END
END
END
SOLID DEFINITION: Solid 1
Material = Brick Common
Option = Material Library
MORPHOLOGY:
Option = Continuous Solid
END
END
SOLID MODELS:
HEAT TRANSFER MODEL:
Option = Thermal Energy
END
THERMAL RADIATION MODEL:
Option = None
END
END
END
INITIALISATION:
Option = Automatic
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = 0 [m s^-1]
END
RADIATION INTENSITY:
Option = Automatic
END
STATIC PRESSURE:
Option = Automatic with Value
Relative Pressure = 0 [Pa]
END
TEMPERATURE:
Option = Automatic with Value
Temperature = 298 [K]
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 RESIDUALS:
Option = Full
END
MONITOR TOTALS:
Option = Full
END
END
RESULTS:
File Compression Level = Default
Option = Standard
END
TRANSIENT RESULTS: Transient Results 1
File Compression Level = Default
Option = Standard
OUTPUT FREQUENCY:
Option = Time Interval
Time Interval = 1000 [s]
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 = 5
Minimum Number of Coefficient Loops = 1
Timescale Control = Coefficient Loops
END
CONVERGENCE CRITERIA:
Residual Target = 0.00001
Residual Type = RMS
END
TRANSIENT SCHEME:
Option = Second Order Backward Euler
TIMESTEP INITIALISATION:
Option = Automatic
END
END
END
END
#################################################




but if we apply Ra = constant in expression it will run fine but that is not a actual case

Thank you sir

sunilpatil March 22, 2013 09:41

Hello sir,
Here i attached CEL Expression And Variables
CEL EXPRESSION
LIBRARY:
CEL:
&replace EXPRESSIONS:
Avg Top Temp = areaAve(T)@CEILING_ATM
Differe = Avg Top Temp -Dry Bulb Temp
Dry Bulb Temp = DRY BULB TEMP(t)
Global Radiation = GLOBAL RADIATION ON CEALING(t)
Gr = (g*beta*(Avg Top Temp -Dry Bulb Temp )*L^3)/(visc^2)
Gr1 = if(Gr <= 0,-1*Gr,Gr)
K = 0.026 [W m^-1 K^-1]
L = 1.6 [m]
Nu = 0.5*(Ra1^0.25)
Pr = 0.71375
Ra = Gr1*Pr
Ra1 = if(Ra <= 0,-1*Ra,Ra)
h = (Nu*K)/L
q con = h*(Differe)
q inc = 0.2*Global Radiation
q net = abs(q inc - q con)
visc = 0.00001545 [m^2 s^-1]
END
END
END


USER FUNCTION
LIBRARY:
CEL:
&replace FUNCTION: DRY BULB TEMP
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [K]
INTERPOLATION DATA:
Data Pairs = 0,296.6,3600,296.6,7200,296.2,10800,295.4,14400,29 4.7,18000,293.9,21600,293.4,25200,294,28800,295.4, 32400,297.2,36000,299.3,39600,301.2,43200,303,4680 0,304.4,50400,305.6,54000,306.4,57600,306.8,61200, 306.2,64800,305,68400,303,72000,300.9,75600,299,79 200,297.9,82800,297.2,86400,297
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

LIBRARY:
CEL:
&replace FUNCTION: GLOBAL RADIATION ON CEALING
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [W m^-2]
INTERPOLATION DATA:
Data Pairs = 0,0,3600,0,7200,0,10800,0,14400,0,18000,0,21600,0, 25200,22.7,28800,142.6391808,32400,312.392491,3600 0,536.6750084,39600,777.0714132,43200,878.2510173, 46800,930.2914616,50400,944.7675523,54000,812.9697 022,57600,562.9522782,61200,303.4616159,64800,94.1 5321058,68400,0,72000,0,75600,0,79200,0,82800,0,86 400,0
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

LIBRARY:
CEL:
&replace FUNCTION: GLOBAL RADIATION ON EAST
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [W m^-2]
INTERPOLATION DATA:
Data Pairs = 0,0,3600,0,7200,0,10800,0,14400,0,18000,0,21600,0, 25200,26.25062504,28800,134.1803446,32400,613.6529 541,36000,709.9563769,39600,749.3350016,43200,585. 9693094,46800,382.3657474,50400,181.0886787,54000, 168.5290584,57600,136.8601059,61200,91.50844287,64 800,37.10802755,68400,0,72000,0,75600,0,79200,0,82 800,0,86400,0
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

LIBRARY:
CEL:
&replace FUNCTION: GLOBAL RADIATION ON NORTH
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [W m^-2]
INTERPOLATION DATA:
Data Pairs = 0,0,3600,0,7200,0,10800,0,14400,0,18000,0,21600,0, 25200,16.64501043,28800,127.7400598,32400,117.7366 103,36000,159.4982633,39600,193.2113825,43200,199. 1427246,46800,191.4867793,50400,181.0886787,54000, 168.5290584,57600,136.8601059,61200,91.50844287,64 800,37.10802755,68400,0,72000,0,75600,0,79200,0,82 800,0,86400,0
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

LIBRARY:
CEL:
&replace FUNCTION: GLOBAL RADIATION ON SOUTH
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [W m^-2]
INTERPOLATION DATA:
Data Pairs = 0,0,3600,0,7200,0,10800,0,14400,0,18000,0,21600,0, 25200,15.02644988,28800,19.11653453,32400,104.8584 968,36000,178.380534,39600,253.6895245,43200,287.9 649814,46800,299.1082256,50400,296.3887936,54000,2 62.4819492,57600,192.7043193,61200,110.7620703,648 00,38.61686584,68400,0,72000,0,75600,0,79200,0,828 00,0,86400,0
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

LIBRARY:
CEL:
&replace FUNCTION: GLOBAL RADIATION ON WEST
Argument Units = [s]
Option = Interpolation
Profile Function = Off
Result Units = [W m^-2]
INTERPOLATION DATA:
Data Pairs = 0,0,3600,0,7200,0,10800,0,14400,0,18000,0,21600,0, 25200,15.02644988,28800,19.11653453,32400,104.8584 968,36000,159.4982633,39600,193.2113825,43200,199. 1427246,46800,191.4867793,50400,206.4635465,54000, 382.6542143,57600,424.2491392,61200,306.7672572,64 800,112.6389917,68400,0,72000,0,75600,0,79200,0,82 800,0,86400,0
Extend Max = No
Extend Min = No
Option = One Dimensional
END
END
END
END

THANK YOU SIR

ghorrocks March 23, 2013 06:01

I see your functions like "GLOBAL RADIATION ON SOUTH" have spaces in the numbers - like "828 00". This will cause problems in CFX, but this maight be just an artefact of copying it into the forum.

If that is not the problem I would go through your CEL one function at a time. Replace a function with a constant value and see if it works. Then you should be able to find the function which is causing the problem.

sunilpatil March 23, 2013 06:49

Error finding variable "THERMX"
 
Hello Sir,
Thank you for your reply.i will check as you suggested by replacing function with constant
Thank you sir
Sunil Patil

liucheng860212 April 25, 2013 13:25

Quote:

Originally Posted by sunilpatil (Post 415841)
Hello Sir,
Thank you for your reply.i will check as you suggested by replacing function with constant
Thank you sir
Sunil Patil

Have you solved the problem yet? I got the same error saying " Error finding variable 'THERMX' " when I just do a thermo-hydrodynamic simulation with adiabatic walls.

PS: Well, I think I figured out. I got this error when I tried to monitor a force. I fixed it by cancling the user-defined monitor.

sunilpatil April 26, 2013 07:00

THERMX error
 
Hello liucheng860212,
ya i also did the same,i removed out put monitor of ceiling temperature it worked fine.
Thank you
Sunil patil


All times are GMT -4. The time now is 21:55.