1 Numeca Russia User Meeting 2011 FINE TM /Open OpenLabs in FINE TM /Open: Create and customize your physical models Kilian Claramunt, Yingchen Li, Jan E. Anker, Nijso Beishuizen, Dirk Wunsch, Thomas Deconinck, Selvan Gnanakumaran
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 1/30
1
Numeca Russia
User Meeting 2011
FINETM /Open
OpenLabs in FINETM /Open:Create and customize your physical models
Kilian Claramunt, Yingchen Li, Jan E. Anker, Nijso Beishuizen,
Dirk Wunsch, Thomas Deconinck, Selvan Gnanakumaran
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 2/30
2 Numeca Russia User Meeting 2011
Motivation
Industrial applications usually involve multiple physical
phenomena
Users should be allowed to adapt and/or add modeling capabilities
in a flexible, efficient and user-friendly way ... … without programming Fortran or C++ code
Combustion, radiation,
pollutant formation, CHT,…
Multiphase flows, e.g. cavitation
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 3/30
3 Numeca Russia User Meeting 2011
OpenLabs
OpenLabs allows users to customize or add physical models in a
flexible and a user-friendly GUI
Users don’t need to care about programming details and code
structure Users benefit from NUMECA’s CFD industrial environment and
features (HPC, parallelization, meshing capabilities, advanced
numerical methods)
OpenLabs can be used in a wide variety of industrial and academic
applications Compared to source-coded models, CFD solutions are obtained
with identical computing and memory costs
Free access for all FINETM /Open community
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 4/30
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 5/30
5 Numeca Russia User Meeting 2011
OpenLabs’ functionalities in FINETM /Open 2.11 allow
addition/customization of models for
Initial field customization, e.g. initial free surface position for a VOF simulation
Boundary conditions, e.g. unsteady inlet boundary conditions
Turbulence modeling, e.g. realizable k-ε model, round/jet anomaly correction
Mass diffusion to track pollutant concentration
Steady or unsteady source terms, e.g. time-dependent heat source
… and many more!
OpenLabs
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 6/30
6 Numeca Russia User Meeting 2011
OpenLabs in FINETM /Open 2.12 will include
GUI to introduce the physical models
Library built by a simple click in the GUI and automatically loaded by theflow solver
The process is automatized and fully integrated in the FINETM/Open
environment
Customization of thermal and transport properties
Access to the combustion and radiation modules Multi-block with similar physical models; access to solid block customization
OpenLabs under Windows platform
OpenLabs
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 7/307 Numeca Russia User Meeting 2011
OpenLabs’ functionalities in FINETM /Open 2.12 will allow
addition/customization of models for
Combustion and radiation, e.g. radiation optical properties, soot formation
Multi-phase flows, e.g. droplet condensation model (currently being
implemented), cavitation models,…
Heat source term added to a solid block
… and many more!!!
OpenLabs
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 8/308 Numeca Russia User Meeting 2011
How?
Set up your test case withFINETM /Open
Add/customize
your physicalmodel withOpenLabs’
GUI
Create thelibrary with asimple click in
OpenLabs
Launch theflow solver of FINETM /Open
Analyze and
post-processthe results
withCFViewTM
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 9/309 Numeca Russia User Meeting 2011
OpenLabs - Examples of applications.
Managing initial conditions
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 10/3010 Numeca Russia User Meeting 2011
=>INITIAL_PROFILES
@ INITIAL_FIELD: Initial_VelocityX
->EXPRESSION:
->ExistingField: Vx
@ INITIAL_FIELD: Initial_VelocityY
->EXPRESSION:
->ExistingField: Vy
=>AUXTERMS
Imposing initial field for Vx and Vy to follow geometry initialization
Customize initial conditions.
Imposing initial velocity field
r
=>INITIAL_PROFILES
@ INITIAL_FIELD: Initial_VelocityX
->EXPRESSION: IF(xCoord<0.1) 1.0 \
ELSEIF(yCoord>2.65) 0.0 \
ELSE cos(asin((xCoord-0.1)/r))
->ExistingField: Vx
@ INITIAL_FIELD: Initial_VelocityY
->EXPRESSION:
->ExistingField: Vy
=>AUXTERMS
=>INITIAL_PROFILES
@ INITIAL_FIELD: Initial_VelocityX
->EXPRESSION: IF(xCoord<0.1) 1.0 \
ELSEIF(yCoord>2.65) 0.0 \
ELSE cos(asin((xCoord-0.1)/r))
->ExistingField: Vx
@ INITIAL_FIELD: Initial_VelocityY
->EXPRESSION: IF(xCoord<0.1) 0.0 \
ELSEIF(yCoord>2.65) 1.0 \
ELSE sin(asin((xCoord-0.1)/r))
->ExistingField: Vy
=>AUXTERMS
=>INITIAL_PROFILES
@ INITIAL_FIELD: Initial_VelocityX
->EXPRESSION: IF(xCoord<0.1) 1.0 \
ELSEIF(yCoord>2.65) 0.0 \
ELSE cos(asin((xCoord-0.1)/r))
->ExistingField: Vx
@ INITIAL_FIELD: Initial_VelocityY
->EXPRESSION: IF(xCoord<0.1) 0.0 \
ELSEIF(yCoord>2.65) 1.0 \
ELSE sin(asin((xCoord-0.1)/r))
->ExistingField: Vy
=>AUXTERMS
@ r=sqrt((xCoord-0.1)*(xCoord-0.1)+ (yCoord-2.65)*(yCoord-2.65))
Velocity field at
zero iterations
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 11/3011 Numeca Russia User Meeting 2011
OpenLabs - Examples of applications.
Managing boundary conditions
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 12/3012 Numeca Russia User Meeting 2011
To simulate the effect of upstream wake of the vane, the relative total
pressure is fitted and placed as boundary condition at the inlet of the rotor
The formula is dependent on the physical time t
Customize Boundary Conditions.
Unsteady inlet boundary to the GE-E3 blade
Stator
Po-bg total pressure at inlet, 223332.0 (Pa)
n number of blades which is 76 Ω rotation speed, 8283rpm
τ time period for one wake passage,
Rotor
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 13/3013 Numeca Russia User Meeting 2011
Customize Boundary Conditions.
Unsteady inlet boundary to the GE-E3 blade
=>AUXTERMS
@ n = 76
@ omega = 8283
@ P0_bg = 223331.0
@ PI = 3.1415926535
@ Period = 60.0/(omega*n)
Po-bg total pressure at inlet, 223332.0 (Pa)
n number of blades which is 76
Ω rotation speed, 8283rpm
τ time period for one wake passage,
=>AUXTERMS
@ n = 76
@ omega = 8283
@ P0_bg = 223331.0
@ PI = 3.1415926535
@ Period = 60.0/(omega*n)
=>CUSTOM_BOUNDARY_CONDITIONS
@ CUSTOMIZED_BOUNDARY_CONDITION: PtInlet
->EXPRESSION: P0_bg * (1.0 - 0.15 * pow((sin(n*tCoord/2+PI*Time/Period)),10))
->ExistingBC: "Absolute Total Pressure" , rotor_inlet
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 14/3014 Numeca Russia User Meeting 2011
Customize Boundary Conditions.
Unsteady inlet boundary to the GE-E3 blade
Simulation over one period
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 15/3015 Numeca Russia User Meeting 2011
OpenLabs - Examples of applications.
Managing additional transport equations
and source terms
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 16/3016 Numeca Russia User Meeting 2011
Diffusion of a pollutant in a room
Transport equation and source terms.
Mass diffusion to track pollutant concentration
Window
outlet
AC inletPOLLUTANT
Cold Temperature
Door inletAIR
Warm Temperature
The pollutant isheavier than air
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 17/3017 Numeca Russia User Meeting 2011
Modeling approach in an incompressible flow
Transport equation for the pollutant mass fraction
Density is considered constant except for the buoyancy term in themomentum equation (Boussinesq approx.)
The solutal buoyancy source term has to be added to the momentum equation
Transport equation and source terms.
Mass diffusion to track pollutant concentration
Thermal buoyancy
already in
FINETM /Open
Solutal buoyancy to
be added with
OpenLabs
βT is the thermal
expansion coeff.
βS is the solutal
expansion coeff.
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 18/3018 Numeca Russia User Meeting 2011
Transport equation and source terms.
Mass diffusion to track pollutant concentration
=>CONSTANTS
@ solutalExpansion = 0.14372@ Le = 1.0
@ SigmaY = 1.0
@ YRef = 0.0
@ gravityX = 0.0
@ gravityY = 0.0
@ gravityZ = -9.8
βS solutal expansion coeff.
Lewis number assumed unity
Reference mass fraction
Turbulent Prandtl number
Gravity
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 19/30
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 20/3020 Numeca Russia User Meeting 2011
Transport equation and source terms.
Mass diffusion to track pollutant concentration
=>SOURCETERMS
@ SOURCE: sourceToMomX
->EXPRESSION: Density * ( solutalExpansion*(Y-YRef)) * gravityX
->AddToExistingPde: MomentumXEquation
@ SOURCE: sourceToMomY
->EXPRESSION: Density * ( solutalExpansion*(Y-YRef)) * gravityY
->AddToExistingPde: MomentumYEquation
@ SOURCE: sourceToMomZ
->EXPRESSION: Density * ( solutalExpansion*(Y-YRef)) * gravityZ
->AddToExistingPde: MomentumZEquation
Source term to the
momentum equationThermal buoyancy
already in
FINETM /Open
Solutal buoyancy to
be added with
OpenLabs
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 21/30
21 Numeca Russia User Meeting 2011
OpenLabs.
Demonstration
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 22/30
22 Numeca Russia User Meeting 2011
OpenLabs.
Demonstration
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 23/30
23 Numeca Russia User Meeting 2011
OpenLabs.
Demonstration
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 24/30
24 Numeca Russia User Meeting 2011
Transport equation and source terms.
Mass diffusion to track pollutant concentration
Y
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 25/30
25 Numeca Russia User Meeting 2011
OpenLabs - Other applications
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 26/30
26 Numeca Russia User Meeting 2011
Initial free surface position for an unsteady VOF simulation
Broken Dam: water column; dam removed at start of the analysis
Evolution of the wave shape with time
The initial position of the free surface is prescribed with OpenLabs
OpenLabs.
Other applications
Thickness of the interface
Distance function
Distance function
C = 0
C = 1
(xLim,yLim)
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 27/30
27 Numeca Russia User Meeting 2011
Turbulence modeling
The Low-Re Yang-Shih k-ε turbulence model is available in FINETM/Open
As a validation exercise, the same model has been fully added with OpenLabs
Identical results are obtained using FINETM/Open and OpenLabs models
Results are shown for the turbulent flat plate: profiles of Vx and k at the outlet
OpenLabs.
Other applications
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 28/30
28 Numeca Russia User Meeting 2011
Soot formation model
Combustion and radiation models are already available in FINETM/Open
Soot strongly effects the radiation absorption coefficient and therefore the
radiative heat transfer
Soot formation models have been added with OpenLabs
The transport equation for the soot mass fraction with
the nucleation and oxidation source terms is added
The absorption coefficient of the radiation model is is
modified to take into account soot
OpenLabs.
Other applications (cont.)
7/28/2019 Create and Customize Your Physical Models
http://slidepdf.com/reader/full/create-and-customize-your-physical-models 29/30
29 Numeca Russia User Meeting 2011
OpenLabs libraries will be progressively shared with all FINETM /Open
community with a variety of functionalities and applications:
You are very welcome to share your experience and extend the
OpenLabs library !!!!
OpenLabs.
Conclusions
Type of functionality
How to customize initial
solutionsHow to customize boundary
conditions
How to add/modify source
terms in a transport equation
How to add transport equations
…
Application / Physical
modeling
Turbulence modeling
Combustion modeling
Radiation modeling
Multiphase modeling
BC customization
…