Simulation Report
Simulation Report
Contents
1. Global ............................................................................................................................................... 3
2. Component 1 ................................................................................................................................... 4
2.1.
Definitions ................................................................................................................................ 4
2.2.
Geometry 1 .............................................................................................................................. 4
2.3.
Materials .................................................................................................................................. 5
2.4.
Laminar Flow .......................................................................................................................... 10
2.5.
Mesh 1.................................................................................................................................... 23
3. Study 1 ........................................................................................................................................... 25
3.1.
Stationary ............................................................................................................................... 25
3.2.
Solver Configurations ............................................................................................................. 25
4. Results ............................................................................................................................................ 29
4.1.
Data Sets ................................................................................................................................ 29
4.2.
Plot Groups ............................................................................................................................ 30
Global
Name
Untitled.mph
Program
Unit system SI
Used products
COMSOL Multiphysics
Component 1
Unit system
SI
2.1 Definitions
2.1.1
Coordinate Systems
Boundary System 1
Coordinate system type Boundary system
Tag
sys1
Settings
to
Settings
Name
Value
2.2 Geometry 1
Geometry 1
Units
Length unit
mm
Property
Value
Space dimension
Number of domains
Number of boundaries 25
Number of vertices
2.2.1
25
Import 1 (imp1)
Name
Value
Geometry import
DXF file
Filename
2.3 Materials
2.3.1
Air
Air
Selection
Domain 1
5
Material parameters
Name
Value
Unit
Pa*s
rho(pA[1/Pa],T[1/K])[kg/m^3] kg/m^3
Basic Settings
Description
Value
Relative permeability
Relative permittivity
Dynamic viscosity
eta(T[1/K])[Pa*s]
1.4
Electrical conductivity
Cp(T[1/K])[J/(kg*K)]
Density
rho(pA[1/Pa], T[1/K])[kg/m^3]
Thermal conductivity
Speed of sound
cs(T[1/K])[m/s]
Functions
Piecewise
Cp
Piecewise
rho
Analytic
Piecewise
cs
Analytic
eta
Cp
cs
Refractive index Settings
Description
Value
Refractive index
Refractive index, imaginary part {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
2.3.2
Aluminum
Aluminum
Selection
Boundaries 125
Basic Settings
Description
Value
Relative permeability
900[J/(kg*K)]
Thermal conductivity
Electrical conductivity
Relative permittivity
Density
2700[kg/m^3]
Description
Value
0.33
Murnaghan Settings
Description
Value
Description
Value
Lam parameter
5.1e10[Pa]
Laminar Flow
Selection
Domain 1
Equations
Settings
Description
Value
10
Description
Value
Discretization of fluids
P1 + P1
Off
Compressibility
Off
None
1[atm]
Streamline diffusion
On
Crosswind diffusion
On
Isotropic diffusion
Off
Off
Used products
COMSOL Multiphysics
Variables
Name
Expression
Unit Description
Selection
spf.dz
Thickness
Domain 1
spf.pref
1[atm]
Pa
Domain 1
spf.pA
p+spf.pref
Pa
Absolute pressure
Domain 1
spf.nx
dnx
Boundaries 125
spf.ny
dny
spf.nz
Boundaries 125
spf.nxmesh
root.dnxmesh
Boundaries 125
spf.nymesh root.dnymesh 1
spf.nzmesh
11
Boundaries 125
2.4.1
Fluid Properties 1
Fluid Properties 1
Selection
Domain 1
Equations
Settings
Description
Value
Density
From material
Property
Density
Air
Basic
Basic
Variables
Name
Expression
Unit
Description
Selection
spf.rho
material.rho
kg/m^3
Density
Domain 1
12
Name
Expression
Unit
Description
Selection
spf.mu
material.mu
Pa*s
Dynamic viscosity
Domain 1
spf.sr
sqrt(0.5*(4*ux^2+2*(uy+
vx)^2+4*vy^2)+eps)
1/s
Shear rate
Domain 1
spf.divu
ux+vy
1/s
Divergence of velocity
field
Domain 1
spf.Fx
N/m^3
Volume force, x
component
Domain 1
spf.Fy
N/m^3
Volume force, y
component
Domain 1
spf.Fz
N/m^3
Volume force, z
component
Domain 1
spf.U
sqrt(u^2+v^2)
m/s
Velocity magnitude
Domain 1
spf.vorticityx
1/s
Vorticity field, x
component
Domain 1
spf.vorticityy
1/s
Vorticity field, y
component
Domain 1
spf.vorticityz
vx-uy
1/s
Vorticity field, z
component
Domain 1
spf.vort_magn
sqrt(spf.vorticityx^2+spf.
vorticityy^2+spf.vorticity
z^2)
1/s
Vorticity magnitude
Domain 1
spf.cellRe
0.25*spf.rho*sqrt(emetri
c(u,v)/emetric2)/spf.mu
Domain 1
spf.nu
spf.mu/spf.rho
m^2/s
Kinematic viscosity
Domain 1
spf.betaT
d(spf.rho,p)/spf.rho
1/Pa
Isothermal
compressibility
coefficient
Domain 1
spf.mu_eff
spf.mu+spf.muT
Pa*s
Dynamic viscosity
Domain 1
spf.muT
Pa*s
Turbulent dynamic
viscosity
Domain 1
spf.T_stressx
spf.K_stressxp*spf.nxmesh
N/m^2
Total stress, x
component
Boundaries
125
spf.T_stressy
spf.K_stressyp*spf.nymesh
N/m^2
Total stress, y
component
Boundaries
125
spf.T_stressz
spf.K_stresszp*spf.nzmesh
N/m^2
Total stress, z
component
Boundaries
125
spf.K_stressx
spf.mu_eff*(2*ux*spf.nx
mesh+(uy+vx)*spf.nymes
N/m^2
Viscous stress, x
Boundaries
13
Name
Expression
Unit
h2*spf.divu*spf.nxmesh/3
)
Description
Selection
component
125
spf.K_stressy
spf.mu_eff*((vx+uy)*spf.
nxmesh+2*vy*spf.nymes
h2*spf.divu*spf.nymesh/3
)
N/m^2
Viscous stress, y
component
Boundaries
125
spf.K_stressz
2*spf.divu*spf.mu_eff*s
pf.nzmesh/3
N/m^2
Viscous stress, z
component
Boundaries
125
spf.K_stress_tens
orxx
spf.mu_eff*(2*ux2*spf.divu/3)
N/m^2
Domain 1
spf.K_stress_tens
oryx
spf.mu_eff*(vx+uy)
N/m^2
Domain 1
spf.K_stress_tens
orzx
N/m^2
Domain 1
spf.K_stress_tens
orxy
spf.mu_eff*(uy+vx)
N/m^2
Domain 1
spf.K_stress_tens
oryy
spf.mu_eff*(2*vy2*spf.divu/3)
N/m^2
Domain 1
spf.K_stress_tens
orzy
N/m^2
Domain 1
spf.K_stress_tens
orxz
N/m^2
Domain 1
spf.K_stress_tens
oryz
N/m^2
Domain 1
spf.K_stress_tens
orzz
-2*spf.mu_eff*spf.divu/3
N/m^2
Domain 1
spf.K_stress_tens
or_testxx
spf.mu_eff*(2*test(ux)2*test(ux+vy)/3)
N/m^2
Domain 1
spf.K_stress_tens
or_testyx
spf.mu_eff*(test(vx)+test
(uy))
N/m^2
Domain 1
spf.K_stress_tens
or_testzx
N/m^2
Domain 1
spf.K_stress_tens
or_testxy
spf.mu_eff*(test(uy)+tes
t(vx))
N/m^2
Domain 1
spf.K_stress_tens
or_testyy
spf.mu_eff*(2*test(vy)2*test(ux+vy)/3)
N/m^2
Domain 1
14
Name
Expression
Unit
Description
Selection
spf.K_stress_tens
or_testzy
N/m^2
Domain 1
spf.K_stress_tens
or_testxz
N/m^2
Domain 1
spf.K_stress_tens
or_testyz
N/m^2
Domain 1
spf.K_stress_tens
or_testzz
2*spf.mu_eff*test(ux+vy
)/3
N/m^2
Domain 1
spf.upwind_helpx
m/s
Upwind term, x
component
Domain 1
spf.upwind_helpy
m/s
Upwind term, y
component
Domain 1
spf.upwind_helpz
m/s
Upwind term, z
component
Domain 1
spf.tau_vdxx
spf.mu*(2*ux2*(ux+vy)/3)
Pa
Strain rate, xx
component
Domain 1
spf.tau_vdyx
spf.mu*(vx+uy)
Pa
Strain rate, yx
component
Domain 1
spf.tau_vdzx
Pa
Strain rate, zx
component
Domain 1
spf.tau_vdxy
spf.mu*(uy+vx)
Pa
Strain rate, xy
component
Domain 1
spf.tau_vdyy
spf.mu*(2*vy2*(ux+vy)/3)
Pa
Strain rate, yy
component
Domain 1
spf.tau_vdzy
Pa
Strain rate, zy
component
Domain 1
spf.tau_vdxz
Pa
Strain rate, xz
component
Domain 1
spf.tau_vdyz
Pa
Strain rate, yz
component
Domain 1
spf.tau_vdzz
-2*spf.mu*(ux+vy)/3
Pa
Strain rate, zz
component
Domain 1
spf.Qvd
spf.tau_vdxx*ux+spf.tau
_vdxy*uy+spf.tau_vdyx*
vx+spf.tau_vdyy*vy
W/m^3
Viscous dissipation
Domain 1
spf.res_u
px+spf.rho*u*ux+spf.rho
*v*uy-(d(2*ux2*spf.divu/3,x)+d(uy+vx,
N/m^3
Equation residual
Domain 1
15
Name
Expression
Unit
Description
Selection
y))*spf.mu-spf.Fx
spf.res_v
spf.rho*u*vx+py+spf.rho
*v*vy(d(vx+uy,x)+d(2*vy2*spf.divu/3,y))*spf.muspf.Fy
N/m^3
Equation residual
Domain 1
spf.res_p
spf.rho*spf.divu+u*d(spf
.rho,x)+v*d(spf.rho,y)
kg/(m^3*s)
Pressure equation
residual
Domain 1
Shape functions
Name
Shape function
Unit
Description
Shape frame
Selection
Lagrange (Linear)
m/s
Velocity field, x
component
Material
Domain 1
Lagrange (Linear)
m/s
Velocity field, y
component
Material
Domain 1
Lagrange (Linear)
Pa
Pressure
Material
Domain 1
Weak expressions
Weak expression
Integration frame
Selection
(p-spf.K_stress_tensorxx)*test(ux)spf.K_stress_tensorxy*test(uy)spf.K_stress_tensoryx*test(vx)+(pspf.K_stress_tensoryy)*test(vy)
Material
Domain 1
spf.Fx*test(u)+spf.Fy*test(v)
Material
Domain 1
spf.rho*(-(ux*u+uy*v)*test(u)-(vx*u+vy*v)*test(v))
Material
Domain 1
-spf.rho*spf.divu*test(p)
Material
Domain 1
-(u*d(spf.rho,x)+v*d(spf.rho,y))*test(p)
Material
Domain 1
spf.streamlinens
Material
Domain 1
spf.crosswindns
Material
Domain 1
16
2.4.2
Wall 1
Wall 1
Selection
Equations
Settings
Description
Value
Boundary condition
No slip
Off
Constraint method
Elemental
Variables
Name
Expression
Unit
Description
Selection
spf.ubndx
m/s
Velocity at boundary, x
component
spf.ubndy
m/s
Velocity at boundary, y
component
spf.ubndz
m/s
Velocity at boundary, z
component
17
Shape functions
Constraint
Constraint force
Shape function
Selection
-u+spf.ubndx
test(-u)
Lagrange (Linear)
-v+spf.ubndy
test(-v)
Lagrange (Linear)
spf.ubndz
2.4.3
Initial Values 1
Initial Values 1
Selection
Domain 1
Settings
Description
Value
18
2.4.4
Inlet 1
Inlet 1
Selection
Boundary 13
Equations
Settings
Description
Value
Boundary condition
Velocity
Off
Constraint method
Elemental
Variables
Name
Expression
Unit
Description
Selection
spf.ubndx
-nojac(spf.nxmesh)*spf.U0in
m/s
Velocity at
boundary, x
component
Boundary 13
spf.ubndy
-nojac(spf.nymesh)*spf.U0in
m/s
Velocity at
boundary, y
Boundary 13
19
Name
Expression
Unit
Description
Selection
component
spf.ubndz
-nojac(spf.nzmesh)*spf.U0in
m/s
Velocity at
boundary, z
component
Boundary 13
spf.U0in
m/s
Normal inflow
velocity
Boundary 13
Shape functions
Constraint
Constraint force
Shape function
Selection
test(-v+spf.ubndy)
spf.ubndz
test(spf.ubndz)
2.4.5
Outlet 1
Outlet 1
Selection
Boundary 5
Equations
20
Settings
Description
Value
Boundary condition
Pressure
Pressure
Normal flow
On
Suppress backflow
Off
Off
Elemental
Variables
Name
Expression
Unit
Description
Selection
spf.p0
Pa
Pressure
Boundary 5
spf.un_here
u*nojac(spf.nxmesh)+v*nojac
(spf.nymesh)
m/s
Intermediate
variable
Boundary 5
spf.ut_herex
uspf.un_here*nojac(spf.nxmes
h)
m/s
Intermediate
variable, x
component
Boundary 5
spf.ut_herey
vspf.un_here*nojac(spf.nymes
h)
m/s
Intermediate
variable, y
component
Boundary 5
spf.ut_herez
spf.un_here*nojac(spf.nzmes
h)
m/s
Intermediate
variable, z
component
Boundary 5
spf.un_there
m/s
Intermediate
variable
Boundary 5
spf.ut_therex
spf.un_there*nojac(spf.nxmes
h)
m/s
Intermediate
variable, x
component
Boundary 5
spf.ut_therey
spf.un_there*nojac(spf.nymes
h)
m/s
Intermediate
variable, y
component
Boundary 5
spf.ut_therez
spf.un_there*nojac(spf.nzmes
h)
m/s
Intermediate
variable, z
component
Boundary 5
spf.KStressn_avx
spf.K_stress_tensorxx*spf.nx
mesh+spf.K_stress_tensorxy*
spf.nymesh+spf.K_stress_tens
orxz*spf.nzmesh
N/m^2
Average viscous
stress, x
component
Boundary 5
21
Name
Expression
Unit
Description
Selection
spf.KStressn_avy
spf.K_stress_tensoryx*spf.nx
mesh+spf.K_stress_tensoryy*
spf.nymesh+spf.K_stress_tens
oryz*spf.nzmesh
N/m^2
Average viscous
stress, y
component
Boundary 5
spf.KStressn_avz
spf.K_stress_tensorzx*spf.nx
mesh+spf.K_stress_tensorzy*
spf.nymesh+spf.K_stress_tens
orzz*spf.nzmesh
N/m^2
Average viscous
stress, z
component
Boundary 5
spf.KStressTestn_avx
spf.K_stress_tensor_testxx*sp
f.nxmesh+spf.K_stress_tensor
_testxy*spf.nymesh+spf.K_str
ess_tensor_testxz*spf.nzmes
h
N/m^2
Average viscous
stress, x
component
Boundary 5
spf.KStressTestn_avy
spf.K_stress_tensor_testyx*sp
f.nxmesh+spf.K_stress_tensor
_testyy*spf.nymesh+spf.K_str
ess_tensor_testyz*spf.nzmes
h
N/m^2
Average viscous
stress, y
component
Boundary 5
spf.KStressTestn_avz
spf.K_stress_tensor_testzx*sp
f.nxmesh+spf.K_stress_tensor
_testzy*spf.nymesh+spf.K_str
ess_tensor_testzz*spf.nzmesh
N/m^2
Average viscous
stress, z
component
Boundary 5
spf.ujumpx
spf.ut_herex-spf.ut_therex
m/s
Velocity jump, x
component
Boundary 5
spf.ujumpy
spf.ut_herey-spf.ut_therey
m/s
Velocity jump, y
component
Boundary 5
spf.ujumpz
spf.ut_herez-spf.ut_therez
m/s
Velocity jump, z
component
Boundary 5
spf.ct_here
24*nojac(down(spf.mu)+spf.
muT)*spf.meshVol/spf.meshV
olInt
Pa*s
Intermediate
variable
Boundary 5
spf.meshVolInt
down(meshvol)
Volume of
interior mesh
element
Boundary 5
spf.meshVol
meshvol
spf.sigma_dg_ns
4*spf.ct_here
Pa*s
spf.f0
spf.p0
N/m^2
Normal stress
Boundary 5
spf.out1.Mflow
spf.out1.intFlow(spf.rho*(spf.
nx*u+spf.ny*v))
kg/s
Mass flow
Global
Boundary 5
22
Boundary 5
Weak expressions
Weak expression
Integration frame
Selection
spf.KStressn_avx*test(spf.ut_herex)+spf.KStressn_avy*test(sp
f.ut_herey)+spf.KStressn_avz*test(spf.ut_herez)+spf.KStressT
estn_avx*spf.ujumpx+spf.KStressTestn_avy*spf.ujumpy+spf.K
StressTestn_avz*spf.ujumpzspf.sigma_dg_ns*spf.ujumpx*test(spf.ut_herex)spf.sigma_dg_ns*spf.ujumpy*test(spf.ut_herey)spf.sigma_dg_ns*spf.ujumpz*test(spf.ut_herez)
Material
Boundary 5
-spf.f0*(test(u)*spf.nxmesh+test(v)*spf.nymesh)
Material
Boundary 5
2.5 Mesh 1
Mesh statistics
Property
Value
0.9475
Triangular elements
2221
Quadrilateral elements
Edge elements
217
Vertex elements
25
Mesh 1
2.5.1
Size (size)
Settings
23
Name
Value
21.4
0.0722
Curvature factor
0.25
Finer
Selection
Remaining
Settings
Name
Value
Number of iterations 20
24
Study 1
3.1 Stationary
Study settings
Property
Value
Physics interface
Discretization
Geometry
Mesh
Solution 1
Name
Value
Use study
Study 1
Name
Value
Name
Value
Solution Zero
Values of variables not solved for
Name
Value
Solution Zero
25
Name
Value
Name
Value
Name
Value
26
27
LinErr
LinRe
1e-011 3.4e2.4e-012 3.8e2.3e-012 4.2e1.7e-011 4.8e2.1e-011 6.7e1.9e-011 8.6e3.2e-011 1.1e3.7e-012 1.2e3.9e-011 1.7e3.7e-011 2.9e5.4e-011 8.2e2.4e-011 3.4e6.8e-011 2.7e6.5e-011 1.8e1.7e-011 1.6e1.1e-010 3.4e9.5e-011 1.9e3.7e-011 1.5e7.1e-011 1.4e3.6e-011
5e-
Name
Value
Name
Value
0.01
Name
Value
Solver PARDISO
28
Results
Study 1/Solution 1
Solution
Name
Value
Solution
Solution 1
29
Velocity (spf)
4.2.2
Pressure (spf)
30