• CASD: Close MONITOR_POINT window 5. Choose SINGLE_FIELD_SCALAR from SCENARIO menu • CASD: Define which variables to report at monitors 6. Click on , drag mouse pushing MOUSE+LEFT across all monitors, • CASD: Log pressure at all three transducers 7. Repeat for and • CASD: Log pressure impulse and dynamic pressure, too 8. Click and choose SINGLE_FIELD_3D from SCENARIO menu • CASD: Define variables for contour plots 9. Click on , CTRL-, CTRL-, • CASD: Pressure, flame and velocity vectors. CTRL needed to select more than one (NB! deselect when using the scroll bar) 10. Choose SIMULATION in SCENARIO menu OR ∗ scen sim • CASD: Choose output and simulation parameters 11. Click on , enter 50 , • CASD: Increase number of contour plots, return to main menu 12. Click on GAS_COMP... in SCENARIO menu OR ∗ scen gas_c • CASD: Define gas cloud loc., size, comp. and concentration 13. Click on , 0 0 0 • CASD: Position of bounding box describing gas cloud FLACS v9.0 User’s Manual
2.6 Introductory example
25
14. Click on , 25.6 8 8 • CASD: Dimension of gas cloud equals module dimensions 15. Click on , 91.7 7 1.3 • CASD: Gas composition is defined 16. Click on 1.05 0 • CASD: Slightly rich gas mixture is chosen ER=1.05 17. Click on IGNITION in SCENARIO menu 12.5 4.1 4.25 OR ∗ scen ign pos 12.5 4.1 4.25 OK • CASD: Define location of ignition (12.5, 4.1, 4.25) 18. Choose SAVE from the FILE menu • CASD: Save all files, ready to run flacs 19. Minimize CASD • CASD: Leave CASD for now, can be activated easily
Figure 2.15: Adding monitoring points
FLACS v9.0 User’s Manual
26
Getting started
Figure 2.16: Choosing variables for 3D output
Figure 2.17: Adding a gas cloud and choosing the gas composition
FLACS v9.0 User’s Manual
2.6 Introductory example
2.6.4
27
Start FLACS simulation
Select the job in Run Manager and click simulate (if job not visible, use add directory or if directory is already added, right click and rescan), check how the simulation starts up (click log file)
Figure 2.18: Running a simulation in the FLACS Runmanager
2.6.5
Study results in post prosessor Flowvis
Use Run Manager → Tools → Flowvis (or click the FLACS post-processor icon) 1. choose ADD from Page menu (or CTRL+a) • FLOWVIS: Prepare first page 2. click MOUSE+RIGHT, choose PLOT_TYPE and SCALAR_TIME plot • FLOWVIS: Plotting of time histories of variables 3. choose 100001 and P with MOUSE+LEFT, select all 3 monitors (drag mouse) • FLOWVIS: Plot pressure time history at all monitors 4. • FLOWVIS: if simulation is running rescan will update plot 5. Choose MODIFY in the Page menu (or CTRL+m), enter 1 2 • FLOWVIS: divide page into 2 plots 6. Click at lower frame, then MOUSE+RIGHT, PLOT_TYPE, ANNOTATION_ST (or CTRL+0) • FLOWVIS: show numerical values from pressure plots FLACS v9.0 User’s Manual
28
Getting started 7. ADD page and do the same for the DRAG and PIMP variables 8. Choose ADD in Page menu (or CTRL+a), click MOUSE+RIGHT, PLOT_TYPE, 2D... (or CTRL+2) • FLOWVIS: prepare 2D contour plot 9. Choose 100001, P, click • FLOWVIS: contour plot of pressure
10. click MOUSE+RIGHT, choose PLOT_DOMAIN, change k-index to 5 • FLOWVIS: choose XY-cut plane through ignition 11. Click MOUSE+RIGHT, choose VARIABLE_APPEARANCE change Value Range Setting to Fixed • FLOWVIS: choose a user-defined fixed scale for all time steps 12. Choose Min. Value as 0.05 and Max. Value as 2.0 • FLOWVIS: define the scale
Figure 2.19: Showing pressure-time curves with annotation in Flowvis
FLACS v9.0 User’s Manual
2.6 Introductory example
29
Figure 2.20: 2D cutplane plot showing over-pressures
Figure 2.21: Setting plot domain for a volume plot
Time steps can now be changed moving the bottom scroll bar to the right, page can be varied using the right scroll bar. 1. Repeat this method for PROD and VVEC variables (these can be plotted on the same plot) FLACS v9.0 User’s Manual
30
Getting started • FLOWVIS: visualize flame and velocity vectors
Try to show PRESSURE and PROD on the same page using PAGE MODIFY (use a fixed scale for PROD from 0.15 to 0.2 and change Min. Color Index to 9 and Max to 10) Now that you are familiar with Flowvis, try the volume plot menu to study the development of flame (PROD) and pressure Use PLOT DOMAIN to narrow the view window and see below the ceiling
2.6.6
Study the effect of ignition location
Enter CASD, open the 100001.caj job-file, save this as a new job number e.g. 100002.caj Change ignition location in order to study how pressures may vary with different ignition locations End ignition (0.5, 4.1, 4.25), (job number 100002) Your own assumed worst-case location (job number 100003) Report highest pressure achieved on monitor point Make animation of either 2D or volume plots using the export menu (with all timesteps)
FLACS v9.0 User’s Manual
Chapter 3
CASD
32
CASD
The preprocessor CASD for the CFD simulator FLACS is used to prepare the input data, or job data , that defines a FLACS simulation: geometry model, computational grid, porosites, and scenario description. CASD is an acronym for Computer Aided Scenario Design. CASD 4 released in 1994, use X11 graphics, but a new version is available based on QT CASD 5 released in 2001, use Open Inventor graphics CASD 6 released in 2008, use QT and Coin 3D graphics This manual describes CASD 6, but the general functionality of CASD 6 is in principle the same for CASD 4 and CASD 5. CASD 6 is fully backward compatible with CASD 4 and CASD 5.
3.1
Overview
This section provides a general overview of the functionality in CASD.
3.1.1
Starting CASD
Users start CASD by clicking the CASD icon in the run manager window:
Figure 3.1: The CASD desktop icon
or alternatively by executing the command: > run9 casd6
on the command line in Linux.
3.1.2
CASD command line options
The following options can be given when starting CASD on the command line: Option -macro macro file name -numMat maximum number of materials -numObj maximum number of objects -numAsis maxmimum number of assemblies/instances -stackAsis maxmimu number of nested assembly levels -noLock
FLACS v9.0 User’s Manual
Description Read input from specified macro file Default is 50 Default is 10000 Default is 3500 Default is 8 Turns of locking on the database files. Must not be used if more than one user accesses the database simultaneously. This option speeds up the database operations significantly.
3.1 Overview
-display and others
33
Linux: options accepted by X Table 3.1: CASD command line options
Example: Linux: run9 casd -numObj 20000 -numAsis 20000 -noLock
Windows: casd -numObj 20000 -numAsis 20000 -noLock
Alternatively the options can be set permanently in the FLACS Runmanager, Options→Preferences. This will only apply if CASD is started from the Runmanager.
3.1.3
The main window in CASD
Starting CASD 6 opens the main window.
Figure 3.2: The main window in CASD
The main window is divided into the following parts: • The menu bar FLACS v9.0 User’s Manual
34
CASD • The icon bar • The command input field • The geometry window(s) • The status field
These parts are described in the following subsections.
3.1.4
The menu bar
The menu bar contains the following menus: • File • Geometry • Grid • Porosities • Scenario • Block • View • Options • Macro • Help The options on the various menus are described in separate sections in this chapter.
3.1.5
The icon bar
The icon bar contains the following toolbars: • Main toolbar, provides shortcuts to several of the commands on the meny bar: – New, Open, Save, Save as, Import, and Result on the File menu. – Database icon on the Geometry menu. – Calculate and Verify porosities on the Porosity menu. • Graphics toolbar, controls various features of the geometry window(s). – – – –
View splitting. Rectangle zoom. Spinning (toggle on/off). Highlighting option, from filled only (0) to various degrees of contour highlighting (1-5).
• Drawing toolbar, opens the plan drawing dialog box: – Specifying file names for texture (e.g. drawings). – File formats: PNG, JPEG, GIF, TIFF FLACS v9.0 User’s Manual
3.1 Overview
3.1.6
35
The command input field
The command input field represents an alternative interface between the user and CASD, in addition to the regular menus on the menu bar. The control input field contains a scrollable command history list, and a current command context indicator (left side). The user controls the command history list from the keyboard: • UP: retrieves the previous line from the command history list • DOWN: retrieves the next line from the command history list • RETURN: processes the content of the command input field Hence, the user can choose whether to use a menu options on the menu bar, e.g: File→Exit→Yes (to exit and save) or to execute, after typing or retrieving, the following command in the command input field:
∗ file exit yes yes Command line input will in many situations be the most efficient way to work with CASD, and other sections in this chapter present additional examples on how to use this feature. Examples: Using the command input field in CASD • Select a box primitive in an object. The following command moves the box to (2, 2, 2), and would cause the properties dialog to be shown – ∗ edit properties 2 2 – This is because the position is not completely specified. The user does not have to specify all parameters, but must include all values for the parameter specified. • If the user wants to edit one of the last parameters in the dialog, it is not necessary to specify all the parameters in front. The parameter name can be used to indicate which parameter to edit – ∗ edit properties size 2 2 2 vol_por 0.5 • The user can also supply the answer to a question in the input field. To delete an assembly/instance, CASD will ask to confirm the operation. To avoid the question dialog, type the following command – ∗ geometry delete yes – or shorter: ∗ ge de y • To direct the output from a list to a file, append the file name after the list command. For instance, to list geometries in the database, enter the following command, which will create the text file outfile.txt – ∗ geometry list outfile.txt
3.1.7
The graphical area
The graphical area in the main window displays the geometry and the computational grid. In addition to the options on the View menu, there are several ways of manipulating the view: • • • •
Rotation: MOUSE+LEFT Panning: CTRL+MOUSE+LEFT Zoom: MOUSE+SCROLL Rectangle zoom: MOUSE+RIGHT+SELECT FLACS v9.0 User’s Manual
36
CASD • Splitting and closing views: MOUSE+RIGHT+SELECT
The use of these features are quite intuitive, and they will not be described in more detail in this manual.
3.1.8
The message area
The message area in the main window contains information concerning the active database, project, geometry, grid, and units.
3.1.9
Files in CASD
CASD stores job data on a set of files. For the arbitrary job number 010100, the most important files are: • Header file, 010100.caj: ASCII file created by CASD; defines the co, cg, and cm files used by CASD. • Geometry file, co010100.dat3: binary file created by CASD; contains a list of primitives from a CASD database that define the geometry; used by Porcalc and Flowvis. • Grid file, cg010100.dat3: binary file created by CASD; defines the computational mesh; used by CASD, Flacs, and Flowvis. • Porosity file, cp010100.dat3: binary file created by Porcalc (typically from the Grid menu in CASD); defines the porosities for each grid cell; used by Flacs and Flowvis. • Polygon file, cm010100.dat3: binary file created by CASD; defines the polygon model; used by Flowvis (if the file exists). • Scenario file, cs010100.dat3: ASCII file created by CASD; defines the general scenario (monitor points, output variables, fuel region, pressure relief panels, ignition position, etc.); used by CASD, Flacs, and Flowvis. The grid-file is also called the obstruction file, or co-file, and is not a direct input to the simulation; it is however used by Porcalc when generating the porosity file. The File menu in the main window contains commands for creating, opening, and saving the various job files. See section Files in FLACS for further information.
3.1.10
Working with geometries in CASD
To implement the geometry model in CASD can often be the most time consuming part of a project. For modern process facilities it may be possible to import a geometry from an existing CAD model, but for many installations the geometry must be constructed manually from drawings, photographs, etc. A large projects, such as a full probabililistic analysis, can involve hundreds of CFD simulations, and each simulation will typically produce 10-15 different files. Hence, it is very important to organize the files in a well-structured manner. The building blocks in a CASD geometry are instances of objects. The structure within an object is a so-called Constructive Solid Geometry (CSG) model, where simple solid primitives (boxes and cylinders) are combined by Boolean operators (unions and left differences). Objects in CASD can be either global or local. Several geometries can contain instances of the same global object, whereas a local object can only be included in the geometry where it was created. It is generally recommended to use global objects, and avoid the use of local objects. FLACS v9.0 User’s Manual
3.1 Overview
37
The list of information required to implement a typical process facility, such as an offshore oil platform or an onshore process plant, is quite extensive: • • • • • • • •
Plot plan Sectional drawings Piping plan HVAC layout Cable trays layout Framing plans Cladding Deck plan
Most FLACS users find it convenient to define standardized axis directions, and the following convention is used by GexCon for typical process facilities: • East-West along the x-axis, with positive x towards the east. • North-South along the y-axis, with positive y towards the north. • Up-Down along the z-axis, with positive z pointing upwards. This results in a conventional right handed coordinate system, where the lower south-western corner of the facility coincides with the origin (0,0,0). Each object in a CASD database is assigned a material property, and each ’material’ is assigned a colour hue from the 0-360° colour circle. Many FLACS users find it convenient to assign certain hues to various structural elements, and the following convention is used by GexCon for typical process facilities. Hue 0 30 60 120 180 200 220 250 300
Colour Red Orange
Description solid walls and decks pressure relief and and louvred panels Yellow grated decks Green anticipated congestion Cyan equipment Light blue structure Medium Blue secondary structure Dark Blue piping Pink equipment Table 3.2: Colour convention used by GexCon
A standardized colour scheme makes it more straightforward to review geometries from old projects.
3.1.11
About congestion, confinement, and vents
In order to have a good representation of the effect of obstacles it is important that they are well represented geometrically by the chosen grid. In most practical situations it will not be possible to represent the smaller obstacles on the grid, these should still be included since they may be treated by proper sub-grid models. Larger obstacles like the floor (or the ground), the ceiling, the walls and larger equipment will be resolved on-grid. This means that they will be adjusted to match the grid lines. FLACS v9.0 User’s Manual
38
CASD
The most challenging geometry to represent properly is repeated obstacles of the same size and spacing as the chosen grid resolution, in such cases the user should consider to change the grid to achieve a better representation. If this type of geometry is dominant it is of vital importance for the accuracy of the result that the representation is good enough. In cases where such a geometry is not dominant one may pay less attention to how it is represented. For normal offshore modules there will be a range of subgrid sized obstacles which are more or less randomly distributed in space. In many experimental setups one will find repeated obstacles of the same size. The basic research on gas explosions past many years now has focused on the effect of obstacle arrays, perhaps to a greater extent than on the effect of more realistic geometries. Both categories are important in order to be able to validate tools like FLACS. It is important to represent the vent openings of a semi-confined geometry properly. If obstacles close to the outer boundaries are adjusted to match the grid, the effective vent area may be affected. In order to verify that the representation of the vent openings is as good as possible the user should check the porosity fields (using CASD or Flowvis).
3.2 3.2.1
File menu New
Shortcut CTRL+N Starts a new simulation job. The New command in the File menu creates a new empty job. If there were unsaved changes to the current job, a dialog box is displayed, asking about saving the changes.
3.2.2
Open
Shortcut: CTRL+O This command opens an existing set of simulation files. The default selection is defined in a ∗.caj file. The Open command in the File menu opens an existing job. If you enter the file name in the command input field, the path must be encapsulated in apostrophes, for instance:
∗ open "../../Test/000000.caj" If you select the command from the menu bar, or if no name is specified in the command input field, the Open dialog box is displayed, allowing you to specify a path and file name to open. By default, the file filter is initiated for selecting CASD job header files (type ∗.caj). But you may also select a geometry file (type co∗.dat3). CASD will then open all files with the same job number. If a geometry is open (in the database), the filter string will be constructed from the project and geometry numbers. It is not possible to open a job that is not compatible with the open project and geometry numbers. If there were unsaved changes to the current job, a dialog box is displayed, asking about saving the changes. The geometry file is not read when a geometry is open in the database. If no geometry is open FLACS v9.0 User’s Manual
3.3 Geometry menu
39
in the database, CASD will display the contents of the geometry file in the graphic area after successful open. The contents of the geometry file can be edited using the Edit File command in the Geometry menu, see section Geometry menu.
3.2.3
Save
Shortcut: CTRL+S Saves the current simulation job (i.e. the various files that define the job). The Save command in the File menu saves the current job.
3.2.4
Save as
Shortcut: CTRL+SHIFT+S The Save As command saves the current job under a new (user-defined) name (job number).
3.2.5
Import
Imports certain specifications from another simulation job (e.g. grid file, scenario file, etc.).
3.2.6
Exit
Shortcut: CTRL+Q Exits the CASD software.
3.3
Geometry menu
CASD stores the geometry in a database, and on the geometry file (co-file). The commands in the Geometry menu in the main window, except the Edit File command, are available when connected to a database. The Save and Save As commands in the File menu writes the geometry to the geometry file. The building blocks in a CASD geometry are instances of objects. Objects can be global or local. Several geometries can contain instances of the same global object, while a local object only can be included in the geometry where it was created. Instances can be grouped under assemblies. Several levels of assemblies can be created. Each instance and assembly has a transformation matrix. The position, scale, and orientation of an instance is the result of the matrices on all levels above the instance, in addition to the matrix for the instance itself. Each geometry is a member of a project. The project is the top level in the CASD data structures. A project can own a number of geometries. Instances and assemblies can be made invisible and visible using the following commands: CTRL+I Make the selected assembly/instance invisible CTRL+SHIFT+I Make the selected assembly/instance visible. FLACS v9.0 User’s Manual
40
CASD
Use the Position command in the Geometry menu to change the position of the selected assembly or instance.
3.3.1
Geometry Database
The first option on the Geometry menu in CASD opens the Database dialog box.
Figure 3.3: The geometry database window in CASD
In the Database dialog box the user can: • • • • •
Create a new database, project, geometry, or object. Connect to or save an existing database. Open or save existing, projects, geometries, or objects. Insert instances in a geometry. Define new materials or edit existing materials.
3.3.1.1
Geometry tab
On the Geometry tab the user can create, open and manipulate projects and geometries. Projects can be renamed and deleted, geometries can be renamed, copied and deleted. 3.3.1.2
Objects tab
The New Object button in the Database dialog box opens the Object window. 3.3.1.3
Materials tab
Each object in a CASD database is assigned a material property, and each ’material’ is assigned a colour hue from the 0-360° colour circle. To define a new material click the New Material button. The new material is defined by a name and a hue, a value between 0 and 360. FLACS v9.0 User’s Manual
3.3 Geometry menu
3.3.2
41
Creating a CASD database
To create a database choose Geometry→Database or type ∗ geometry database. The Geometry Database window is shown. Click the Connect button. A file selection dialog box is displayed. Move to the directory where the database should be created, and write the name of the database, e.g. my_database.db. Alternatively the database can be created using the command input: ∗ database create my_database.db, which will create a database in the current directory. If the Geometry Database window is not open, choose Geometry→Database. Use the New Project button to create a new project, or the Open Project button to open an existing project. When a project is opened, a new geometry can be created clicking the New Geometry button, or open an existing geometry clicking the Open Geometry button. When an existing geometry is opened, the assembly/instance structure and all objects and materials used are loaded into the CASD program. If the geometry contains many assemblies/instances, you may get an error message indicating that there were not room enough in the CASD data structures. See section CASD command line options for information on how you can use command line options to allocate more memory for these structures.
3.3.3
Connecting to a database
To create a new database, see section Creating a CASD database. To connect to an existing database choose Geometry→Database or type ∗ geometry database. The Geometry Database window is shown. Click the Connect button. A file selection dialog box is displayed. Select the CASD_DB file on the database directory you want to connect to. If you enter the file name in the command input field, the path must be encapsulated in apostrophes, for instance:
∗ database connect "MyCasdDB/CASD_DB"
3.3.4
Creating a new or opening an existing object
You can create a new object clicking the New Object button on the Objects tab in the Geometry Database window, or open an existing object using the Open button. When you have completed the New or Open Object command, the object window is displayed.
3.3.5
Selecting a node and a subtree
At any time, a part of the binary tree is selected. It may be a single node, or a subtree containing several nodes. If a subtree is selected, the top node is referred to as the selected node. In the postfix string, the top node is the rightmost node in the subtree. The selected subtree is highlighted in the graphic window, and underlined in the message area. There are two different methods for selecting a subtree. 1. Click MOUSE+LEFT while pointing at a primitive. If several primitives are hit, they are placed on a stack (list). Only one primitive is selected at a time. Press CTRL+TAB command to parse this stack. 2. Use the following commands: FLACS v9.0 User’s Manual
42
CASD CTRL+L Select the previous instance CTRL+R Select the next instance
3.3.6
Maintaining a CASD database
The dbfutil program is available for creating and maintaining CASD file databases. Linux: run9 dbfutil database command [option]
Windows: dbfutil database command [option]
The usage of this program is described in table Using the the dbfutil program. Make sure that no other users are connected to the database when you execute this program. Command create destroy force dellock
restoredep
restorehead
FLACS v9.0 User’s Manual
Description Create database Destroy database Destroy database, override any errors Delete all locks. Use this command if files in the database are still locked after a crash in CASD Restore dependencies. For each object in the database, there is a file containing a list of all geometries that contain instances of the object. (Executing the Information command in the File menu in the Object dialog lists the contents of this file.) This file is used for determining if the object can be deleted when you execute the Delete Object command in the Database menu. CASD updates these files when required. But if a problem should occur for some reason, the restoredep command might help. It updates the file mentioned above for all objects in the database. Restore header files. This command resets the process log file for the database. This file contains a list of (CASD) processes currently connected to the database.
3.3 Geometry menu list
43 List the content of all table files, e.g. list O lists all objects: P List the content of all project table files. O List the content of all object table files. M List the content of all material table files. G List the content of all geometry table files. L List the content of all local object table files. U List the content of all objects-used table files. A List the content of all asis table files. Table 3.3: Using the dbfutil program
We strongly recommend that you make backups of your databases on a regular basis.
3.3.7
Local objects
Local objects consist simply of one box or one cylinder. Use local objects to define entities like walls, floors etc. Define global objects for more complicated things. The name of a local object must start with an underscore character (_). The Local Object command in the Geometry menu creates a local object, and one instance of it. You can of course create several instances of the local object using the Instance command. The Local Object command has two sub choices, Box and Cylinder. Select the appropriate primitive type. CASD will first ask for the material name. Enter the name of an existing material. The material decides the colour of the object. If you haven’t defined any materials, use the New Material command in the Geometry Database window to create one. CASD will then ask for the sizes and porosities for the primitive. CASD creates an instance of the object in (0, 0, 0). Use the Position or Translate command to move it to the correct position. You can use the Properties command to edit material, sizes and porosities for a local object. The Rename command changes the name of the object.
3.3.8
Global objects
A global object is edited in a separate object window. All the commands described in this chapter refers to the menus in the object window. Global objects can have instances in several geometries. The structure within a global object is a constructive solid geometry (CSG) model where simple solid primitives are combined by means of Boolean set operations. The primitives and operations are nodes in a binary tree where the leaves are primitives and the internal nodes are operations. Boxes, cylinders, ellipsoids, general truncated cones (GTC) and complex polyhedrons (CP8) are FLACS v9.0 User’s Manual
44
CASD
the primitive types supported. The box primitive includes planes as a special case. Available operation types are union and difference.
Warning: Only boxes and cylinders should be used in by default, but ellipsoids, general truncated cones and complex polyhedrons can be used in special cases. These latter primitive types have the following important limitations: • No subgrid models, thus not contribution to turbulence and drag force • Porosity calculation takes a long time for these primitive types. There should be no more than 100-200 of these primitives in any given geometry
Figure 3.4: Supported primitive types
A root is a subtree that is not part of another subtree. The object typically contains several roots during editing. But it must contain only one root when it is saved. The postfix string represents a way of visualising the binary tree defining the object. The postfix string for the open object is displayed in the message area in the object window. The selected subtree is highlighted. A material is assigned to each object. The material decides the colour of the object. FLACS v9.0 User’s Manual
3.3 Geometry menu
45
Figure 3.5: The binary tree for an objects, and the corresponding postfix string
3.3.9
Assembly
Opens a dialog box where the user can specify an assembly of several instances. 3.3.9.1
Adding an assembly
Assemblies represents a way to group the instances in complicated geometries. The Assembly command in the Geometry menu adds an assembly to the geometry. CASD will ask for the assembly name. You must enter a name that doesn’t exist on the same level, see below. The assembly is placed in (0, 0, 0). You can transform an assembly in the same way as an instance. All geometries contains at least one assembly, called the top assembly. That assembly can not be deleted. When you create an assembly, it is placed in the geometry structure depending on what was selected on forehand. If an instance was selected, the new assembly is placed after that instance under the same assembly. If an assembly was selected, the new assembly is placed under that assembly. You can later rename the assembly using the Rename command. 3.3.9.2
Selecting an assembly or instance
The selected instance, or all the instances in the selected assembly, are highlighted in the graphic window. The name of the selection is written in the message area. The name is concatenated from the geometry name, the names of all assemblies above the selected assembly/instance, and the name of the selected assembly/instance. Each level is separated by a period (.). An example is shown below. Current Geometry Selection:
M24.A1.COOLER-2 FLACS v9.0 User’s Manual
46
CASD
Here, M24 is the geometry name, A1 is an assembly. The last part of the string is the lowest level. In this example, it is an instance, identified by the object name, COOLER, and the instance number. There are three different methods for selecting an assembly or instance. The first method is to select from the graphic window. To do this, click MOUSE+LEFT while pointing at the instance. If several instances are hit, they are placed on a stack. Only one instance is selected at a time. In CASD4 use the CTRL+TAB command to parse this stack. The second method is to use the following commands: • Select the parent assembly: Press CTRL+U. • Select the child assembly/instance: Press CTRL+D. • Select the assembly/instance name: Press CTRL+F. You are asked to enter the concatenated name to select. • Select the previous assembly/instance on the same level: Press CTRL+L. • Select the next assembly/instance on the same level: Press CTRL+R. The third method is to use the List command in the Geometry menu to pop up a list of the contents of the open geometry. You can use the mouse to select from the list.
3.3.10
Instance
Creates an instance in the current geometry and/or assembly. 3.3.10.1
Adding an instance
To add an instance of an object, use the Instance command in the Geometry menu. CASD will ask for the object name. You must enter the name of an existing object. The instance is placed in (0, 0, 0). Use the Position or Translate command to move it to the correct position. Alternatively the Instance button on the Objects tab in Geometry Datbase dialog can be used. When a new instance is created, it is placed in the geometry structure depending on what was selected on forehand. If an instance was selected, the new instance is placed after that instance under the same assembly. If an assembly was selected, the new instance is placed under that assembly.
3.3.11
Local object
Creates a local object in the current geometry.
3.3.12
Delete
Deletes either the currently selected instance, local object, or the current assembly (must be empty).
3.3.13
List
Lists all assemblies and instances in the current geometry, including modified positions. FLACS v9.0 User’s Manual
3.3 Geometry menu
3.3.14
47
Duplicate
Duplicates the selected instances in the current geometry.
3.3.15
Position
Defines the position of an instance.
3.3.16
Translate
Translates the current instance.
3.3.17
Rotate
The Rotate command rotates the selected assembly or instance. Note that CASD only accepts axis parallel geometry. That means that the rotation angle must be a multiples of 90 degrees.
3.3.18
Scale
Scales the current instance by a certain factor in each spatial direction
3.3.19
Matrix
Specifies the transformation matrix of the current instance. This command is normally not used directly, but is available for macro reading and writing.
3.3.20
Making an assembly or instance visible or invisible
Shortcut: CTRL+I CTRL+SHIFT+I This command lets the user make the current instance invisible/visible.
3.3.21
Select
Selects an instance in the current geometry through the following short cut options. See section Selecting an assembly or instance.
3.3.22
Substitute
Substitutes all instances of one object in the current geometry with instances of another object. The user specifies the name of the existing object and new objects.
3.3.23
Properties
Opens a dialog box where the user can observe and edit the properties of a local object. FLACS v9.0 User’s Manual
48
3.3.24
CASD
Rename
Opens a dialog box where the user can rename assemblies or local objects.
3.3.25
Object
Opens the currently selected object.
3.3.26
Edit file
The Edit File command in the Geometry menu makes it possible to edit the geometry file (co file) for the open job. This command is only available when no geometry is open in the database. The geometry is saved on the geometry file as one single object, when selecting Save in CASD. Upon the Edit File command, an object window is therefore shown for editing this object, if the geometry database is not available, or the user wants to make small modifications to the geometry outside of the database. Since the object structure lacks the assembly/instance mechanism, editing the geometry file directly without using the database is recommended only for geometries with a relatively small number of primitives. For geometries with many primitives, the postfix string is long and difficult to manage. Editing the geometry file for FLACS simulations may be advantageous when the user want to test the impact of small changes in the geometry on the simulation results. Note that there is no way to update the database from the geometry file.
3.4
Object window in CASD
The object window opens from the ’New Object’ button in the database dialog box. The object window opens from the database window. FLACS v9.0 User’s Manual
3.4 Object window in CASD
49
Figure 3.6: The object window in CASD
The message area in the object window shows the postfix string.
3.4.1
File menu in the Object window
The options on the file menu in the object window are explained below.
3.4.2
Save
If the user is editing an object in the database, the Save command in the File menu saves the object on the database. If the user is editing the geometry file, the changes are stored internally in the geometry database, and will be written to the file upon the Save and Save As commands in the File menu in the main window. Exiting from the object window without saving, the changes are lost. The object is stored only if it is consistent, that is if it has only one root. If the object is not consistent, an error message is displayed, and a Union or Left Difference should be added.
3.4.3
Information
The Information command in the File menu displays a list of all geometries containing instances of the open object. FLACS v9.0 User’s Manual
50
3.4.4
CASD
Exit
Upon the Exit command, CASD asks about saving the object, and then whether to exit from the object. If the answer is yes to the last question, the object window is closed.
3.4.5
Edit menu
The options on the edit menu in the object window are explained below. 3.4.5.1
Operations
The Operation command in the Edit menu changes the operation type if the selected node is an operation. 3.4.5.2
Properties
The Properties command in the Edit menu changes the primitive properties if the selected node is a primitive. If you have selected a subtree containing only one type of primitives, the Properties command can be used for changing one or more parameters for all these primitives. 3.4.5.3
Translate
Use the Translate command to translate the selected assembly or instance a specified distance in each axis direction. Use the Translate command in the Edit menu to translate the selected subtree a specified distance in each axis direction. 3.4.5.4
From To
Use the From To command to translate the subtree so that one specified position, the base point, is moved to another, the target point. A dialog box for specifying the two positions is displayed. A circle is displayed in the graphic window, indication the position being edited. CASD keeps a list of positions used in the object. By pressing CTRL+L or CTRL+R, you can parse this list. The coordinates in the dialog box is updated. 3.4.5.5
Rotate
The Rotate command rotates the selected subtree. You must specify a base point for the rotation, and the rotation angle. As for the From To command, you can parse the position list using the CTRL+L or CTRL+R commands. Note that CASD only accepts axis parallel geometry. That means that the rotation angle must be a multiple of 90 degrees. 3.4.5.6
Scale
The Scale command is only legal when an instance of a local object consisting of a box is selected. FLACS v9.0 User’s Manual
3.4 Object window in CASD
51
The Scale command scales the selected subtree. You must specify a base point for the scaling, and the scaling factor. You can parse the position list using the CTRL+L or CTRL+R commands
3.4.5.7
Delete
The Delete command in the Edit menu deletes the last node in the postfix string, the selected subtree or the current root. Note that if the postfix string for the object is consistent, it consists of only one root. Therefore deleting the current root deletes the entire object.
3.4.5.8
Mark
The Mark command is used in connection with the Substitute command. Select Mark command to mark the subtree to be substitued with the subtree selected when the Substitute command is selected.
3.4.5.9
Substitute
The Substitute command in the Geometry menu substitutes all instances of one object with instances of another object. You are asked to specify the two object names. The Substitute command in the Edit menu substitutes the selected subtree with another subtree. Use the Mark command to select the first subtree. The substitute command implies the following steps. (Let subtree 1 denote the first subtree and subtree 2 the second subtree.) 1. Make a copy of subtree 2 and give it a new identity, say subtree 3. 2. Delete subtree 1 from the postfix string. 3. Insert subtree 3 in the postfix string in the position where subtree 1 was situated.
3.4.5.10
Duplicate
The Duplicate command in the Geometry menu duplicates the selected instance. You are asked to enter the number of copies, and the distance between each copy in the three axis directions. Click on Ok, and a dialog box pops up for each copy, allowing you to edit the position. The Duplicate command in the Edit menu duplicates the selected sub tree. You are asked to enter the number of copies, and the distance between each copy in the three axis directions. Union operations are added automatically, so that the resulting sub tree includes the original one.
Creating pipe bundles Start with creating one cylinder with the appropriate diameter, length and direction. Use the Duplicate command in the Edit menu to duplicate the cylinder in one direction. Use the same command once more to duplicate the resulting row of cylinders in the other direction. If you need to change some parameters for all the cylinders, select the entire pipe bundle sub tree and use the Properties command. If you want to change the distances between the cylinders, this can be done by scaling the entire sub tree. Afterwards you can use the Properties command to reset the cylinder diameters and lengths. FLACS v9.0 User’s Manual
52
CASD
3.4.5.11
Material
The Material command in the Edit menu edits the material name for the object. You must enter an existing material name.
3.4.5.12
Matrix
The Matrix command was introduced to make it simple to create and run macros for creating geometries.
Warning: This command should normally not be used in interactive mode.
3.4.6
Add menu in the Object window
The options on the add menu in the object window are explained below.
3.4.6.1
Box
The Box command in the Add menu adds a box at the end of the postfix string. A dialog box for defining the box parameters is displayed.
3.4.6.2
Cylinder
The Cylinder command adds a cylinder at the end of the postfix string. A dialog box for defining the box parameters is displayed.
3.4.6.3
Ellipsoid
The Ellipsoid command adds an ellipsoid at the end of the postfix string. A dialog box for defining the ellipsoid parameters is displayed. Note warning about the use of ellipsoid.
3.4.6.4
CP8
The CP8 command adds a complex polyhedron at the end of the postfix string. A dialog box for defining the complex polyhedron parameters is displayed. Note warning about the use of complex polyhedron. FLACS v9.0 User’s Manual
3.4 Object window in CASD
53
Figure 3.7: Definition of a complex polyhedron
3.4.6.5
GTC
The GTC command adds a general truncated cone at the end of the postfix string. A dialog box for defining the general truncated cone parameters is displayed. Note warning about the use of general truncated cone.
Figure 3.8: Definition of a general truncated cone
3.4.6.6
Union
The Union command adds an union operation at the end of the postfix string. This command is only legal if the object contains at least two roots which can be connected by the operation.
3.4.6.7
Left Difference
The Left Difference command adds a difference operation at the end of the postfix string. This command is only legal if the object contains at least two roots which can be connected by the operation. If using CASD4, use the Shade command in the View menu, to see the result of the operation. Note that the right hand side operator of a difference operation must be a primitive.
3.4.6.8
Copy
The Copy command adds a copy of the selected sub tree at the end of the postfix string. FLACS v9.0 User’s Manual
54
CASD
3.4.6.9
Object
The Object command adds a copy of a specified object at the end of the postfix string.
3.4.7
Select menu in the Object window
The options on select file menu in the object window are explained below. 3.4.7.1
Previous
Shortcut: CTRL+L Selects the previous primitive or subtree 3.4.7.2
Next
Shortcut: CTRL+R Selects the next primitive or subtree 3.4.7.3
Stack
Shortcut: CTRL+TAB This command will parse (cycle through) the list of selected primitives or subtrees if more than one is selected.
3.4.8
View menu in the Object window
The options on the view menu in the object window are explained below. 3.4.8.1
Print
The Print menu allows exporting a screenshot of the CASD window into different formats: • Postscript • RGB • IV 3.4.8.2
Examiner Viewer and Fly viewer
The default and most widely used viewer is the Examiner viewer. The Fly viewer can be used to fly through the geometry. 3.4.8.3
The XY, XZ and YZ views
The option XY View and XZ View display a projection of the geometry in the XY and XZ planes respectively. The options YZ East View and YZ West View display a projection of the geometry in the YZ plane along the positive and negative Y-axis respectively. FLACS v9.0 User’s Manual
3.4 Object window in CASD
3.4.9
55
3D View
The 3D View option displays a default 3D view of the geometry. 3.4.9.1
Axis
The Axis option turns axis display on and off. 3.4.9.2
Maximize
The option Maximize maximizes the visible window to display the entire geometry and grid. 3.4.9.3
Grid Display
Three different options are available in the Grid Display menu: • Off: The grid is not displayed. Only the geometry would be displayed. • Working Direction: The grid would be displayed in the working direction only. • All Directions: The grid would be displayed in the three directions. 3.4.9.4
Annotation
The options in this menu are currently not used. 3.4.9.5
Draw Style
Different options are available in this menu: • • • •
Off: The geometry will not be displayed. Wireframe: Only the edges of the objects that compose the geometry would be displayed. Filled: Surfaces of the objects that compose the geometry would be displayed. Scenario Wireframe: Only the edges of scenario objects (for example, a fuel region) would be displayed. • Scenario Filled: Surfaces of scenario objects would be displayed. 3.4.9.6
LOD and Properties
The LOD (Level Of Details) and properties menus control the details of the geometry displayed.
3.4.10
Macro menu in the Object window
The options on the macro menu in the object window are explained below. for more infromation about CASD macros see section Macro menu. 3.4.10.1
Run
Read a macro file defining a single object. FLACS v9.0 User’s Manual
56
CASD
3.4.10.2
Record
Writes all commands given to a user specified macro file.
3.4.10.3
Write Object
Writes a macro file containing the current object.
3.5
Grid menu
The simulation volume is divided into a set of control volumes by three sets of grid planes, one in each axis direction. There is always a current grid working direction, and a selected region of grid planes in this direction. The current working direction is shown in the message area. The lines indicating the selected region is highlighted.
3.5.1
Simulation volume
The Simulation Volume command lets you change the simulation volume extent in all three directions. If you increase the volume, the original grid planes are kept, but one additional plane is added in each direction. If you decrease the volume, planes outside the new volume are deleted, and new planes are created on the volume borders.
3.5.2
Direction
The Direction command changes the working direction. Legal input is x, y or z. The Grid menu commands Region, Add, Position, Move, Delete, Smooth, Stretch and List affects the grid planes in the working direction.
3.5.3
Region
The Region command substitutes the selected grid planes by a new set of grid planes. CASD asks you to enter the new number of control volumes in the region.
3.5.4
Add
The Add command adds a new grid plane in the working direction. You are asked to enter the coordinate value for the new plane.
3.5.5
Position
The Position command lets you edit the position for the selected grid plane. FLACS v9.0 User’s Manual
3.5 Grid menu
3.5.6
57
Move
The Move command moves the selected grid planes a specified distance.
3.5.7
Delete
The Delete command deletes the selected grid planes.
3.5.8
Smooth
The Smooth command substitutes the selected grid planes by a new set of grid planes. For the Smooth command, the sizes of the control volumes at each end of the region is kept unchanged. The sizes of the control volumes between them varies gradually. This function is typically used when refining the grid around a leak.
3.5.9
Stretch
The Stretch commands substitutes the selected grid planes by a new set of grid planes. This is particularly useful when stretching the grid towards the outer boundaries. The Stretch command has two sub-choices: • Negative Direction (typically used at the boundaries at the negative end of the axis) • Positive Direction (typically used at the boundaries at the positive end of the axis) You must enter the size of the control volume at one end of the region, default is the current size. Then you must enter a factor by which the sizes of the control volumes in the specified direction increases/decreases. Attention: Note that stretching of the grid should be avoided in areas of the simulation domain where the main combustion is happening. The flame model in FLACS has been validated for cubical control volumes, thus the user should not stretch the grid in areas where accurate results are required. It is however good practice to stretch the grid towards the boundaries, to concerve simulation time and computer memory.
3.5.10
Information
The Information command displays status information about the defined grid, while the List command lists the grid coordinates in the working direction.
3.5.11
List
The Information command displays status information about the defined grid, while the List command lists the grid coordinates in the working direction. FLACS v9.0 User’s Manual
58
CASD
3.5.12
Display
The Display command turns grid display off, displays the grid in the working direction only, or displays the grid in all three directions.
3.5.13
Select
The selected region of grid planes is limited by two planes, the lower and upper limit. If only one plane is selected, the upper and lower limit is the same grid plane. Grid planes are selected using the following commands: • Lower boundary – Select the next grid plane: CTRL+RIGHT – Select the previous grid plane: CTRL+LEFT • Upper boundary – Select the next grid plane: CTRL+UP – Select the previous grid plane: CTRL+DOWN
3.5.14
Grid-related operations
3.5.15
Importing the grid from another job
Use the Import command in the File menu to import the grid from another job. If you enter the grid file name in the command input field, the path must be encapsulated in apostrophes, as described in section . If you select the command from the menu bar, or if no name is specified in the command input field, the Import dialog box is displayed, allowing you to specify the path and file name for a grid file. You will be asked to verify that the current grid is overwritten by the grid from the specified file.
3.5.16
Saving the grid
The Save and Save As commands in the File menu saves the grid, together with the rest of the job data. If the grid is changed, you will need to recalculate the porosities.
3.5.17
Grid-related utilities
FLACS is deleivered with a command line tool for creating an manipulating the grid. This tool can be used to quickly edit or get information about the grid. Please see section gm for further information.
3.5.18
Grid guidelines
The grid resolution should be chosen to obtain a simulation result within an acceptable time frame. In most cases a reasonably good result can be obtained on a coarse grid within less than one hour (in some cases 5 minutes), and high quality results can normally be generated in a few hours (or at least over the night). FLACS v9.0 User’s Manual
3.5 Grid menu
59
Never start a project with a calculation on a grid that will be running for days. If such long simulation times are necessary, always start simulating on a much coarser grid [even if this violates guidelines] to check that the scenario and setup are OK. The user should keep the position of the grid lines in mind while defining the geometry. The geometry details such as walls and decks should be adjusted to the closest grid line when inputted. Thereby the user keeps track of the positioning instead of having the geometry moved in an unwanted direction by the porosity calculation program. In the grid embedding process, it is highly recommended to use Grid→Information in Casd to check different aspects of the grid. Grid sensitivity tests are also recommended. 3.5.18.1
Gas explosion simulations
Attention: The user should always apply cubical grid cells in the combustion region. Deviations from this will give different flame propagation and pressures, and the validation work done is no longer valid. Deviations of the order 10% in aspect ratio is OK, deviations by a factor of 2 in aspect ratio is not OK. If one chooses not to follow this guideline, the results can be somewhat improved by setting a fixed control volume size for the time stepping routine (see section The SETUP namelist, example TIME_STEPPING=" STRICT:L_FIX=1.0" ). Channels and confined vessels and rooms (filled with gas from wall to wall) must always be resolved by a minimum of 5-6 grid cells in smallest direction if flame acceleration shall be modelled. This also applies for pipes where flame acceleration along the pipe is of interest. A pipe connection from one vessel to the next may have less grid cells across the diameter (but preferably more than 1 CV) if only flame transport by pressure difference and not flame acceleration along the pipe shall be modelled. Increase the inner diameter of angles and bends somewhat when modelling pipes with cylinder minus primitives. Remember that one full grid cell is required inside the solid walls around " minus primitive holes" to ensure that the walls will not be leaking. Unconfined gas clouds as well as partially filled clouds should have a minimum of 13 grid cells across the cloud if both sides are unconfined, and a minimum of 10 grid cells in directions where cloud meets confinement on one side (example vertical direction for dense gas cloud in chemical plant). It is not recommended to use non-cubical grids for explosion simulations. As they are often used for dispersion simulations, the dispersion simulation results should be dumped, thereafter the rdfile utility program should be used to transfer the results from the dispersion grid to a grid better suited for explosions, see example below: > run9 rdfile rd111111.n001 rd222222.n001
Here 111111 is the dispersion calculation job number and 222222 is copy of the job, in which the grid has been modified to follow explosion grid embedding guidelines. The grid of job 222222 must be completely inside the grid of 111111. The grid can be stretched outside the combustion region in directions where pressure recordings are not of interest. In directions where pressure wave propagation is of interest, one should not stretch the grid because this will reduce the sharpness and quality of the pressures. A proper distance to external boundaries is important. At least 5-10 grid cells from vent opening to external boundary should be used in situations where the external explosion is not important (small vent area or strong turbulence inside vessel). FLACS v9.0 User’s Manual
60
CASD
In low-congested vessels with significant vent opening, external explosions and reversing of flow may give a strong feedback into the vessel in connection with venting. To pick up this properly, the distance to the external boundaries should be significant (maybe 3-4 times the length of the vessel). EULER boundary will reflect negative pressures, which can destroy results when simulating far field pressure propagation. In this case PLANE_WAVE may be recommended (but then the boundaries must be far away so no products from combustion reaches the boundaries). For unconfined situations, try to have the same distance to boundaries in all directions (use stretching in directions with less interest in results). 3.5.18.2
Blast wave propagation in the far field
Maximum control volume size should be 10% of gas cloud diameter: max CV = 0.1 × ( gascloudvolume)1/3
(3.1)
The grid cells must be approximately cubical in explosion simulations and the cell size should be maintained in directions of interest for blast propagation simulations. For vessel burst the same guidelines applies. If the pressure is much higher than 10 barg, somewhat larger grid cells than this criterion can be acceptable. Remember to use PLANE_WAVE boundary condition and proper distance to the boundaries. 3.5.18.3
Dispersion simulations near field
Calculations of flammable gas requires grid refinement near the leak. The area of the expanded jet (at ambient pressure) must be resolved by one grid cell ( ACV < 2 × A jet ) except for lowmomentum releases of highly buoyant gases such as hydrogen where the guideline ( ACV < 1.25 × A jet ) should be followed. In most cases, a grid refinement near the leak helps in keeping moderate calculation times while getting acceptable results. Grid refinement guidelines for efficient simulation of high velocity jets recommend only refinement across jet direction (not along). CFLC should be increased by refinement factor (CFLV should not be changed). Smoothing from fine grid cells near jet to normal grid cells further away from jet is recommended. If the jet is not along the axes, is impinging or has a low momentum with positive/negative buoyancy, then extending the refined region of the grid in one or more directions may be required. Refinement along the jet may then also be required. If only far field concentrations are of interest, the refinement near the leak may not be needed. Quicker calculations and less stability problems will be seen without the refinement. For dense gas calculations, it may be a good idea to use a finer resolution in the vertical direction near the ground than in the other two directions. Increase CFLC by this refinement factor. Outside the main area of interest, further stretching to the boundaries is recommended to minimize the influence of the boundaries. If in doubt whether the distance to the boundaries influences your results, increase the distance further to check the sensitivity. The general recommended procedure for setting up the grid is: 1. Cover the computational domain with a uniform grid 2. Refine the grid in the near region of a jet (perpendicular to the jet axis) 3. Stretch the grid outside the main region towards the boundaries FLACS v9.0 User’s Manual
3.6 Porosities menu
3.6
61
Porosities menu
The Porosities menu supplies commands for calculating and verifying porosities. Calculate This command starts the porosity calculation program, Porcalc Verify This command starts Flowvis for porosity verification
3.6.1
Calculate
The Calculate command starts the porosity calculation program, Porcalc. The version of Porcalc can be selected on the Preferences dialog in CASD. The default version of Porcalc uses a resolution of 64, i.e. the control volume is divided into 64 parts when calculating the porosities. It can be useful (and necessary) to use the porcalc_16 version if the porosity calculations are very slow, e.g. if the geometry contains slanted primitives. If you are working with a multiblock simulation, only the porosities for the selected block are calculated. Upon the Save and Save As command in the File menu, CASD will warn you about blocks where the grid and/or geometry is changed since the last time you calculated the porosities. Porcalc can also be started from the command line or the Runmanager.
3.6.2
Verify
The Verify command lets you view the calculated porosities. Flowvis is started for porosity verification. A 2D Cut Plane plot for the appropriate job is automatically created with volume and area porosities shown. The Plot Domain dialog box pops up. This dialog box lets you select other planes as wished. By clicking inside a control volume, you can verify the porosity values for that volume. If you are working with a multiblock simulation, only the porosities for the selected block are verified.
3.7
Scenario menu
The purpose of this chapter is to outline how to edit the scenario sections. A short description of each section and the impact on the FLACS simulations will be given. The scenario-file (cs-file) is an ASCII file and it is easy to edit manually as well as using CASD. The items in the Scenario menu are read from a scenario definition file. There is one default scenario definition file, but several other choices can be activated by changing the scenario template. This can be done using the Preferences command in Options menu. For instance, the default+1 template activates several advanced options, especially in the Simulation and Output Control section. The current job should then be saved, closed, and reopened in order for the new options to be available. When a section has been selected, the items in the section are displayed as a list. The method for editing a scenario section depends on the type of section. Sections such as INITIAL_CONDITIONS and IGNITION contain a list of parameters, each with one or more values. A parameter is selected for editing by clicking on it, or by typing the parameter name in the command input field in the main window. FLACS v9.0 User’s Manual
62
CASD
Sections such as SINGLE_FIELD_SCALAR_TIME_OUTPUT and SINGLE_FIELD_3D_OUTPUT contain a list of items, which can be selected. For some sections, each item has a subsection. For SINGLE_FIELD_SCALAR_TIME_OUTPUT the user must select from a list of monitor points for each variable selected. An item is selected by clicking on it, or by typing the item name in the command input field in the main window. To select several items using the mouse apply CTRL or SHIFT keys. If a selected option shall be deselected without selecting an alternative option, deselect by clicking while pressing the CTRL key. Typing an item name in the command input field in the main window selects/deselects the item if previously not selected/selected. Importing the scenario from another job Use the Import command in the File menu to import the scenario from another job. The user will be asked to verify that the current scenario is overwritten by the scenario from the specified file. Saving the scenario The Save and Save As commands in the File menu saves the scenario (together with the rest of the job data).
3.7.1
Monitor points
Monitor points are user defined locations in the simulation domain where one or more variables are to be monitored during the simulation. The maximum number of monitor points allowed is currently 1000. Positions for monitor points are given in the unit selected in Options Preferences (normally meter) When the user has defined all the desired monitor points, he/she may specify a list of variables to be monitored and the relevant monitor numbers for each variable (see the next sub- section). FLACS identifies the 8 surrounding control volume centres and writes an interpolated value of the specified variables to the scalar-time output file (nodes on other side of wall or zero porosity will not be used when interpolating).
Attention: The user should avoid putting monitor points exactly on a grid line or within fractions of a grid cell size from a wall. Usually it is best to enter monitor points according to the grid, not according to the geometry. E.g. if monitor points are placed on each side close to a solid wall they may not necessarily be in two different control volumes (as was the intention). During the porosity calculation the wall will be adjusted to the nearest control volume face (grid line) and might therefore move to the wrong side of a monitor point! It may also be a good idea to ensure that none of the monitor points are inside fully blocked control volumes. FLACS v9.0 User’s Manual
3.7 Scenario menu
63
Figure 3.9: Specification of Monitor Points
FLACS v9.0 User’s Manual
64
CASD
Figure 3.10: How to position the monitor points
3.7.2
Single field scalar time output
As indicated in the previous sub-section, the user must define all the monitor points and panels before he/she can specify the list of output variables. For each output variable, the user may enter one or more numbers indicating the monitor point number(s), or panel number if it is a panel averaged output variable, for which you want this variable to be measured. An example from a scenario-file section is shown below [NP identifies variable P (pressure) whereas NPP identifies variable PP (panel pressure) which gives average pressure load across a surface described by panel]: SINGLE_FIELD_SCALAR_TIME_OUTPUT NP 1 2 3 4 5 NPP 1 2 3
This shows that pressure is reported for monitor points 1-5 and panel pressure is reported for panels 1-3 (The definition of measurement panels is described in the section on Pressure Relief Panels). To select more than 1 monitor point use CTRL+MOUSE+LEFT or CTRL+SHIFT+MOUSE+LEFT, or simply drag the mouse over the list of monitor points while pressing LEFT. If a variable and monitor points are selected by a mistake, these can be deselected by using the CTRL+MOUSE+LEFT to deselect the last monitor point of a variable. FLACS v9.0 User’s Manual
3.7 Scenario menu
65
The most commonly used monitor point variables for explosion simulations are pressure (P), dynamic pressure (DRAG), panel pressure (PP), pressure impulse (PIMP) and sometimes flow velocity (UVW). For dispersion calculations volume gas concentration (FMOLE) and flow velocity (UVW) are among the most commonly used variables. The monitor point results for job 010100 are written to the r1010100.dat3 file which can be read by Flowvis. If ASCII data is required, the r1file-utility program can be used (see the section on FLACS utilities). The first section on the scenario file defines the names and identifiers for all the variables which may be selected for output from FLACS. In order to select alternative units for certain variables (e.g. psi or kPa for pressure or K for temperature) the scenario-file should be manually edited. The variable pressure is described in the top of the scenario-file as follows: NP "P "Pressure"
" 1 "(barg)
"
N
An output in a different unit (psig) can simply be obtained by editing this as follows: NP "P "Pressure"
" 1 "(psig)
"
N
Similar changes can be made for other variables and other units. Please note that the units of time must always be seconds (however, it is possible to change them to ms in Flowvis). A complete list of all variables available can be found in section Output variables in FLACS.
3.7.3
Single field 3D output
This is an output facility in FLACS which enables the user to generate plots of the spatial distribution of the variables (e.g. cut plane plots and volume plots) at different moments in time. The user needs to specify the list of desired variables for SINGLE_FIELD_3D_OUTPUT, an example from scenario-file is shown below (here P, PROD and VVEC are selected in CASD): SINGLE_FIELD_3D_OUTPUT NP NPROD NVVEC NU NV NW
To select more than one variable press the CTRL-key selecting variable 2 and 3 etc. Please observe that when velocity vectors are selected for output (VVEC), directional velocities U, V and W will automatically be selected. These should not be deselected while VVEC remains selected (this is possible e.g. by manual editing of scenario-file), if this is done very strange results will be seen in Flowvis as result file is not consistent with scenario file used by Flowvis to interpret results. The user should be aware that this type of output may give very large files (r3-files). If the user wants to save disk space, the number of output variables and the number of time instants for output must be limited. The r3-files are binary, if ASCII data is wanted the utility program r3file can be applied. The most commonly used variables for explosion modelling are pressure (P), flame (PROD), sometimes gas volume concentration (FMOLE), dynamic pressure (DRAG), maximum pressure (PMAX) FLACS v9.0 User’s Manual
66
CASD
and velocity vectors (VVEC). For dispersion FMOLE and VVEC will be the most common variables to report. In certain situations the variable PMAX may not be written to even if specified (zero results everywhere). This may happen if the simulation job requires more RAM than allowed (e.g. due to self-defined limits in Linux) To specify output times DTPLOT and NPLOT will normally be used (see section Simulation and output control), and sometimes also cc-files (see section Runtime simulation control file). To create animations it is normally recommended to have plots at 100-200 different moments in time. When creating results files to be used for animations a combination of DTPLOT and NPLOT is usually recommended for explosions, for dispersion only DTPLOT should be used. Units of output variables can be changed using the technique described in the previous section. A complete list of all variables available can be found in section Output variables in FLACS.
3.7.4
Simulation and output control
This section describes parameters for general simulation and output control. scenario-file setup which is suitable for gas explosion simulations is listed below: TMAX LAST CFLC CFLV SCALE MODD NPLOT DTPLOT GRID WALLF HEAT_SWITCH
The default
-1 -1 5 0.5 1 1 -1 -1 "CARTESIAN" 1 0
In addition the following entries are available in the default+1 template: TSTART TMIN LOAD STEP KEYS
-1 -1 -1 "" ""
For dispersion simulations, higher value of CFL numbers (20 and 2) are recommended. NPLOT should be -1 (it has no meaning) and a finite value of DTPLOT should be given. A detailed description of each parameter is given below. 3.7.4.1
TMAX
This is the maximum time interval (seconds) that the simulation will last. For explosion simulations, default value in CASD can typically be used. The default value set by CASD is -1, this means there is no maximum time specified, and automatic stop criteria will be applied. Automatic stop criteria will usually work well for explosion calculations. The simulation stops 20 % after either • 90 % of fuel is burnt or pushed out of the domain • 50 % of fuel is burnt or pushed out of the domain (and average pressure becomes negative) FLACS v9.0 User’s Manual
3.7 Scenario menu
67
Automatic stop criteria can not be used for dispersion, and will also be less useful for some special situations. If pressure If far-field blast pressures are of interest, the automatic stop criteria should not be used as it may stop the simulation before blast waves have hit their target. In a gas dispersion simulation TMAX will typically be from a few seconds to a few minutes. Simulation results are not affected by variation of this parameter.
3.7.4.2
LAST
This is the maximum number of iterations allowed for the simulation. The default value set by CASD is -1 that means that there is no limitation to the number of iterations. This value may be changed if additional control of when to stop a simulation is required but this is generally not used.
3.7.4.3
CFLC
This is a Courant-Friedrich-Levy number based on sound velocity. The value of CFLC connects simulation time step length to control volume dimension through signal propagation velocity (in this case the velocity of sound), in the following way: Each time step length is chosen so that sound waves may propagate only a limited distance, which is the average control volume length multiplied by the value of CFLC. The default value set by CASD is 5.0. The time step limit imposed by this criterion is normally dominant in the early phase of an explosion, when flow velocities and combustion rate are still low (see also CFLV). Simulation results may change with this parameter. Therefore, it is not recommended to change this value for explosion simulations as the validation work is nullified. If convergence problems occur (a rare occurrence), CFLC may be reduced by a factor of 2. However, other problems should first be ruled out. Extreme changes of the CFL numbers (i.e. by an order of magnitude) are never recommended for normal simulations. Note that for multi-block simulations a maximum CFLC=0.5 should be used for the BLAST blocks. It is recommended to use CFLC=0.2 in the BLAST blocks in order to ensure numerical stability and good representation of the blast wave. For dispersion simulations, a default value of 20 is normally recommended. This can be increased by the grid refinement factor (if applicable) i.e. if the grid is refined near the leak by a factor of 5, a CFLC number of 20∗5 = 100 may be used (a lower value should be used in case of stability problems). For far-field blast simulations, this should be combined with STEP="KEEP_LOW" in order to keep the time step short even after the explosion is outside the "core" area (more information is given below).
3.7.4.4
CFLV
This is a Courant-Friedrich-Levy number based on fluid flow velocity. The value of CFLV connects simulation time step length to control volume dimension through signal propagation velocity (in this case the fluid flow velocity), in the following way: Each time step length is chosen so that the fluid may propagate only a limited distance, which is the average control volume size multiplied by the Courant number. The default value set by CASD is 0.5. The time step limit imposed by this criterion is normally dominant in the later phase of an explosion, when flow velocities and combustion rate are high (see also CFLC). FLACS v9.0 User’s Manual
68
CASD
Simulation results may change with this parameter. Therefore, it is not recommended to change this value for explosion simulations as the validation work is nullified. If convergence problems occur (rare), CFLV may be reduced by a factor of 2. However, other problems should first be ruled out. Extreme changes of the CFL numbers (i.e. by an order of magnitude) are never recommended for normal simulations. For dispersion simulations, a default value of 2 is normally recommended (a lower value should be used in case of stability problems).
3.7.4.5
SCALE
This parameter is used to scale all linear dimensions in a scenario. This means that a 10 m long explosion vessel is calculated as being 20 m long if SCALE is set to 2.0. Positions and sizes of equipment, gas cloud, ignition region, monitor points, panels etc. are scaled accordingly. The CASD default value is 1.0. This parameter will influence simulation results, typically explosion pressures increase with increasing scale. This is practically never used for realistic geometries.
3.7.4.6
MODD
This is a parameter that may be used to determine how often data for scalar-time plots are written to the results file during a simulation: data are namely stored every MODD timesteps. CASD default is set to 1. This variable does not influence simulation results, only the amount of data stored. This is normally not used in explosion simulations, but a value of MODD=10 (or higher) may be used for long dispersion simulations.
3.7.4.7
NPLOT
This is a parameter that may be used to determine how often data for field plots are written to file during a simulation: data are namely stored at given fuel levels where NPLOT is the number of fuel levels equally spaced between zero and a maximum. Fuel level is defined as the current total mass of fuel divided by the initial total mass of fuel. This output mechanism is not active in the case of a gas dispersion simulation (leaks are specified). This variable does not influence simulation results, only the amount of data stored.
3.7.4.8
DTPLOT
This is the time interval (in seconds) for field output. This is useful in gas dispersion simulations and also in gas explosion simulations when frequent output is required. Note that the field output file will become very large if DTPLOT is set small. This variable does not influence simulation results, only the amount of data which is stored.
3.7.4.9
WALLF
This is a control switch that specifies the use of wall-functions in FLACS. Wall-functions are used to resolve the effect of momentum and thermal boundary layers on the momentum and energy equations in near wall regions. The following choices are available: FLACS v9.0 User’s Manual
3.7 Scenario menu
69
• 0 = OFF • 1 = ON The CASD default value of WALLF is set to 1. This parameter will influence the simulation results. When WALLF equals 1 wall, functions are employed based on theory explained in [Sand and Bakke, 1989]. A slightly modified version of this wall-functions procedure is employed when WALLF equals 2. However, no validation work for WALLF=2 is available. The default value of WALLF=1 should always be used. 3.7.4.10
HEAT_SWITCH
This parameter is meant to control the activation of thermal attributes on objects in Flacs. Default choice is zero (0) as large scale explosions is not much influenced by heat loss. Choosing one (1) will let walls and objects have background temperature, and if gas temperature changes, some heat transfer into or out of gas will take place. This is useful for small-scale confined explosions and dispersion with important heat effects. This can be combined with KEYS="RADIATE=04" in order to activate radiation heat losses (see below). If heat switch is activated all solid surfaces will now be initialized with ambient temperature (in previous versions of FLACS a cs [jobno].HEAT file had to be written). Further heat objects can be specified at different temperatures, see manual. Two models for radiation heat loss can be activated. One simplified model can be activated using a "cs \e [jobno].RAD" file, see description in previous FLACS manual. Alternatively a 6-flux model can be activated with the KEYS-string in the scenario-file or setup-file: KEYS = " RADIATE=04"
This model will calculate gas heat loss (and absorption/scatter) from radiation as well as radiation from hot objects around. Walls absorb 100% of the incoming radiation and emit radiation based on its own temperature. Symmetry planes will reflect 100%. 3.7.4.11
TSTART
Remarks: Only available in the default+1 template. This variable makes it possible to specify a start time for simulation (-1 means not applied => default is zero or time of dump-file). If dump-file exists, one can still adjust TSTART, but the previous history of the simulation can then not be kept using the KEEP_OUTPUT in a setup-file. 3.7.4.12
TMIN
Remarks: Only available in the default+1 template. This variable makes it possible to define a minimum time for simulation. Automatic stop criteria will not activate before TMIN has been reached (-1 means not applied). This can be useful in blast simulations. FLACS v9.0 User’s Manual
70
CASD
3.7.4.13
LOAD
Remarks: Only available in the default+1 template. This makes it possible to load a dump file directly from CASD (instead of using the cc-file). The number of dump file should be specified here. 3.7.4.14
STEP
Remarks: Only available in the default+1 template. It is possibility to give time stepping input (ref. Manual, options for TIME_STEPPING). The options include KEEP_LOW that is recommended for calculations of far-field blast propagation, and effectively prevents time step from growing when explosion calms down. Another option is STRICT:L_FIX=1.00”
that instructs the simulation to use 1m grid size as basis for timestep (and ignore local grid refinement). This can be used instead of increasing the CFLC number as a result of grid refinement. 3.7.4.15
KEYS
Remarks: Only available in the default+1 template. This provides the option for entering setup-file options directly in CASD, e.g. RADIATE=04
for enabling radiation heat loss.
3.7.5
Boundary conditions
In the Boundary condition menu, the user must specific boundary conditions for the outer boundaries of the simulation domain. The lower boundaries in X- Y- and Z-direction are denoted by XLO, YLO and ZLO respectively, and the upper boundaries likewise by XHI, YHI, ZHI. Recommended boundary conditions are as follows: EULER: Euler equations NOZZLE: Nozzle formulation PLANE_WAVE: Plane wave condition WIND: Wind inflow or outflow SYMMETRY boundary condition is generally not used for realistic geometries. In addition, it is also possible to choose EQCHAR, and BERNOULLI, but these boundary conditions are not recommended. FLACS v9.0 User’s Manual
3.7 Scenario menu
71
Remarks: • For most explosion simulations, EULER can and should be used. • For wind and dispersion simulations, NOZZLE boundary condition (similar to EULER) is more robust. • PLANE WAVE boundary condition is recommended for explosion in low confinement and for far field blast propagation. Boundary must be extended far outside the explosion (Flames should not reach boundary). • Different boundaries do not need to have the same condition. • Boundary conditions try to model what happens beyond the boundary. Except for solid walls, this is not straightforward. Sometimes the boundary condition will disturb or even destroy a simulation. Then the user should: 1. Ensure that the chosen boundary conditions are those that fit best to the problem. 2. Consider to increase the Simulation volume and move the boundaries to regions where less steep gradients will cross the boundaries. The details of various boundary conditions are given below:
3.7.5.1
Euler
The inviscid flow equations (Euler equations) are discretized for a boundary element. This means that the momentum and continuity equations are solved on the boundary in the case of outflow. The ambient pressure is used as the pressure outside the boundary. A nozzle formulation is used in the case of inflow or sonic outflow. Warning: EULER boundary condition may give too low explosion pressures in unconfined situations. In such cases, the Simulation volume should be extended and the Plane wave boundary condition should be applied.
3.7.5.2
Nozzle
A nozzle formulation is used for both sub-sonic inflow and outflow and sonic outflow. This condition is suitable for porous areas with small sharp edged holes or grids (e.g. louvres and gratings). A discharge coefficient is calculated from the area porosity and a drag coefficient. NOZZLE condition has shown to give a bit higher explosion pressures than EULER, but it is more robust. Warning: NOZZLE boundary condition may give too low explosion pressures in unconfined situations. In such cases, the Simulation volume be extended and the Plane wave boundary condition should be applied.
3.7.5.3
Plane wave
This boundary condition was designed to reduce the reflection of the pressure waves at open boundaries which occurs when using EULER or NOZZLE. The pressure wave reflection is caused by setting a fixed pressure at the boundary. PLANE_WAVE boundary condition extrapolates the pressure in such a way that reflections are almost eliminated for outgoing waves. FLACS v9.0 User’s Manual
72
CASD
Attention: The pressure might stabilize at a slightly elevated level after an explosion. For low confinement scenarios it is recommended to use Plane wave boundary condition and to extend the domain such that the total volume is about 100 times larger than the volume of the initial gas cloud. Warning: In semi-confined situations where the boundaries are close to the vents, PLANE_WAVE should not be applied. 3.7.5.4
Wind
WIND boundary condition models an external wind field. Velocity and turbulence profiles are specified at the wind boundaries, either by setting some turbulence parameters manually or by choosing one of the atmospheric stability classes, see Pasquill class . WIND boundary conditions are particularly applicable to dispersion scenarios. It is possible to apply WIND on both inflow and outflow boundaries and on boundaries where the flow is parallel to the boundary. Warning: In cases where a generated internal flow has a strong impact on the boundary flow, e.g. gas explosions, WIND should not be used.
Figure 3.11: Specification of Wind boundary condition
Wind speed WIND_SPEED, U0 , is the velocity on the boundary at a given Reference height. It is possible to set WIND_SPEED to positive, zero or negative values, but GexCon recommends to set a postive value and use the Wind direction parameter to specify the direction of the wind. In case of no wind, the user should consider to use another boundary conditions. A uniform velocity FLACS v9.0 User’s Manual
3.7 Scenario menu
73
profile is obtained by setting the Reference height equal to zero. Then the total volumetric flux over the boundary is as follows: V = U0 ∑ An β2n (3.2) n
Wind direction WIND_DIRECTION is a vector and each component may be given a positive, zero or negative value. The sign of this parameter determines the flow direction. A positive value means wind flow in positive direction, that is inflow over the lower boundaries and outflow at the upper boundaries. Wind at an angle different from axis directions may be specified using the WIND_DIRECTION vector. Relative turbulence intensity RELATIVE_TURBULENCE_INTENSITY, IT , is the ratio between the isotropic fluctuating velocity, u0 , and the mean flow velocity U0 : IT =
u0 U0
(3.3)
IT will typically have a value in the range 0.0 to 0.1. This parameter is used to calculate the value for turbulent kinetic energy, k = 3/2u02 , at the boundary. Attention: It is not necessary to set RELATIVE_TURBULENCE_INTENSITY for inflow boundaries when a Pasquill class is specified. When a Pasquill class is set, FLACS will automatically create profiles for velocities and turbulence parameters at the boundary. Turbulence length scale TURBULENCE_LENGTH_SCALE, ` LT is a typical length scale on the boundary. It is used to calculate the rate of dissipation of turbulent kinetic energy, ε at the boundary: ε=
Cµ k3/2 ` LT
(3.4)
For internal flows, the length scale should be about half of the hydraulic diameter. It is not necessary to give a turbulent length scale when a Pasquill class is set. Wind buildup time WIND_BUILDUP_TIME is the time velocities on the boundaries used to rise from zero to WIND_SPEED. A value for WIND_BUILDUP_TIME larger than zero gives a smooth start of the simulation. GexCon recommends to use WIND_BUILDUP_TIME and to start eventual leaks after the wind field has reached steady state. 3.7.5.5
Fluctuating wind
One can also specify fluctuating wind field from the boundaries using a setup-file. Two different frequencies in horizontal directions and one in vertical direction are applied. This can be defined with: VERSION 1.1 $WINDGUST USE=.TRUE. $END FLACS v9.0 User’s Manual
74
CASD
To change the frequencies of the fluctuations, more options must be specified: VERSION 1.1 $WINDGUST USE = .TRUE. AMP = 2.40, 1.90, 1.84 , 2.40, 1.90, 0.00 TAU = 15.0, 10.0, 10.0 , 70.0, 50.0, 0.0 $END
Two different fluctuating periods are assumed along (15s and 70s) and across (10s and 50s) the wind direction. In the vertical direction one period of 10s is used. Fluctuations are done as harmonic periods with average velocity fluctuation equal to 2.4 (along), 1.9 (across) and 1.3 (vertical direction, here constant is multiplied by square root of 2 since only one period is used) times the friction velocity, u∗. 3.7.5.6
Using TRACER mode when simulating dispersion
In flacs2.2.6 a new simulation mode called ’tracer’ has been implemented. In this mode Flacs will only solve a passive transport equation. Below is an example of how to enable the tracer mode: VERSION 1.1 $SETUP KEYS="TRACE:T=100,DT_MUL=Y:5" $END
In the new scenario templates there is an option to define KEYS within the scenario file (no need to use a separate setup-file). The above setup-file (or the KEYS-line defined in the scenario file) will simulate normally until time=100s, thereafter only the fuel transport equation will be solved (flow field will be kept constant). When the flow field equations are switched off at 100s, the time step is at the same time increased by a factor of 5 if using the DT_MUL=Y:5 string. This option can be useful when a dispersion of neutral gas (or small quantities of gas) shall be simulated in an established wind field. This option should be used with care! 3.7.5.7
Effect of temperature gradients
To simulate the effect of an inversion layer, it is possible to define a cold or hot temperature region (layer) in FLACS. One can for instance define a cold valley by using the setup-file as described in previous sections. Boundary conditions can not take a temperature profile. 3.7.5.8
General considerations for boundary conditions
It is generally advantageous to place the outer boundaries of the simulation domain far away from the geometrical extent, but limitations of memory and computing speed may restrict the practical size of the problem, and in most cases one is forced to compromise between quality and cost. Solid wall boundary The solid wall boundary condition is straightforward to model. The velocity vectors are zero at solid walls, both in the tangential and the perpendicular directions. Hence, a zero gradient perpendicular to the boundary, or a fixed value, works well for the scalar variables. Furthermore, wall-functions may improve the modelling of the flow in near-wall regions, both at the outer boundaries and in the interior space. FLACS v9.0 User’s Manual
3.7 Scenario menu
75
External influences In cases where there are obstacles outside the vent openings of a semiconfined volume they should be included in the total simulation volume, because they may have an effect on the explosion. One effect may be that the total venting is reduced due to the external obstacles, especially if they are large and are placed close to the vents. Since the vent flow is changed, also the internal flow past obstacles is modified and the explosion becomes different (higher or lower pressure results). An effect which also may be important is the appearance of an external explosion which will start when the flame reaches any unburnt gas which may have escaped through the vent openings. The pressure waves from the external explosion will propagate back into the semi-confined volume and give rise to higher pressures there. The strength of the external explosion will depend on the local turbulence in the external space, this again depends on the properties of the vent openings and on any obstacles which may be positioned in the external space.
3.7.6
Initial conditions
Initial conditions set values for turbulence fields, temperature and pressure at the beginning of the simulation. Information about the gravity condtions, parameters for the atmospheric boundary layer and the composition of the air is also set here. The default values are as follows: UP-DIRECTION GRAVITY_CONSTANT CHARACTERISTIC_VELOCITY RELATIVE_TURBULENCE_INTENSITY TURBULENCE_LENGTH_SCALE TEMPERATURE AMBIENT PRESSURE AIR GROUND_HEIGHT GROUND ROUGHNESS REFERENCE_HEIGHT LATITUDE SURFACE_HEAT_P1 SURFACE_HEAT_P2 MEAN_SURFACE_HEAT_FLUX PASQUILL_CLASS GROUND_ROUGHNESS_CONDITON
0 0 1 9.8 0.0 0.0 0.0 20.0 100000 "NORMAL" 0 0 0 0 0 0 0 0 0 0 0 "NONE" "RURAL"
Attention: GexCon recommends to use the default values for explosion scenarios. 3.7.6.1
Up-direction
This is a three-component vector defining the upward direction, i.e. opposite of the acceleration due to gravity. The three components denote the spatial components x, y, and z, respectively. It is possible to define any direction, not only directions aligned with the main axes. The vector does not need to be of unit length. 3.7.6.2
Gravity constant
GRAVITY_CONSTANT is the magnitude of the gravitational acceleration, normally 9.8m/s2 . Setting a zero value for this parameter means that there are no gravitational influences on the flow. Note that panels with inertia are also influenced by this parameter if the direction of panel release is in the direction of the Up-direction, see also Pressure relief panels . FLACS v9.0 User’s Manual
76 3.7.6.3
CASD Characteristic velocity
CHARACTERISTIC_VELOCITY, U0 , is used to find values for initial turbulence fields and it should take a positive or a zero value. 3.7.6.4
Relative turbulence intensity
RELATIVE_TURBULENCE_INTENSITY, IT , is the ratio between the isotropic fluctuating velocity, u0 , and the mean flow velocity U0 : u0 IT = (3.5) U0 IT will typically have a value in the range 0.0 to 0.1. This parameter is used together with Characteristic velocity to calculate the value for the turbulent kinetic energy, k = 3/2u02 at the beginning of the simulation. 3.7.6.5
Turbulence length scale
TURBULENCE_LENGTH_SCALE, ` LT , is the length scale of the initial turbulence. A reasonable choice for TURBULENCE_LENGTH_SCALE is a typical value for the grid size ∆ g . ` LT is used to calculate an initial value for dissipation of turbulent kinetic energy, ε:
ε=
3.7.6.6
Cµ k3/2 ` LT
(3.6)
Temperature
The initial temperature, T0 in C ◦ . The default value is T0 = 20.0 C ◦ 3.7.6.7
Ambient pressure
AMBIENT PRESSURE, P0 is the initial pressure in the simulation and the pressure outside the simulation volume. The default ambient pressure is P0 = 100000 Pa = 1 bar. 3.7.6.8
Air
AIR is used to define the composition of the air. "NORMAL" denotes a standard composition of 20.95% oxygen and 71.05% nitrogen in mole fractions. The composition of air can be changed by setting an other oxygen content, either as a mole fraction: AIR
"25%MOLE"
or mass fraction: AIR
"10%MASS"
Changing the air composition will influence the lower flammability limit and upper flammability limit. FLACS v9.0 User’s Manual
3.7 Scenario menu 3.7.6.9
77
Wind boundary layer parameters
In the Wind boundary condition, the reference wind speed, U0 and the direction of the wind are specified. In additions, information about the ground conditions is needed to setup inlet profiles for velocity and turbulence parameters. It is possbile to specify an atmospheric boundary layer stability class, Pasquill class. Then turbulence parameter profiles are generated at the wind boundaries. If no Pasquil class is given, uniform values for k and ε are obtained on the boundary according to the expressions in Relative turbulence intensity and Turbulence length scale. The velocity profile at the wind boundaries is given by the following expression: ´´ ( ∗³ ³ u if z0 > 0 ln z+z0z0 κ u(z) = U0 if z0 = 0
(3.7)
where z is the height abouve the ground, z0 is the atmospheric roughness length and u∗ is a friction velocity. For the neutral and none Pasquill class, u∗ is defined by: u∗ = ln
U0 κ ³z ´
(3.8)
re f
z0
Ground Height GROUND_HEIGHT is height above the ground where the boundary layer actually starts, for instance due to vegetation. Usually GROUND_HEIGHT=0. Ground roughness GROUND_ROUGHNESS, z0 , refers to the aerodynamic roughness length . Typical values for z0 are given in table Typical values for aerodynamic roughness length. z0 should in not be larger than the control volume height close to the surface. aerodynamic roughness length should not be mixed with pipe roughnes etc., but a rule of dump is to relate z0 to the average height ε g of the surface irregularities by z0 = ε g /30. Terrain description z0 (m) Open water, fetch at least 5 km 0.0002 Mud flats, snow; no vegation, no obstacles 0.005 Open flat terrain; grass, few isolated 0.03 obstacles Low crops; occasional large obstacles 0.10 High crops, scattered obstacles 0.25 Parkland, bushes, numerous obstacles 0.5 Regular large obstacle coverage (suburb, 1.0 forest) Table 3.4: Typical values for aerodynamic roughness length
Reference height REFERENCE_HEIGHT, zref is the height relative to the ground where the velocity equals the wind speed. Latitude Due to the rotation of the earth, the height of the atmospheric boundary layer is much larger at equator than at the poles. LATITUDE will only have an effect if a Pasquill class is chosen and the simulation volume is very large (> 200 m) in the z direction.
FLACS v9.0 User’s Manual
78
CASD
Surface heat P1 If the temperature and velocity are known at two different altitudes, for instance from experimental data, it is possible to estimate the surface heat flux. SURFACE HEAT P1 is a vector Z1 U1 T1, where Z1 is the altitude, U1 is the velocity at Z1 and T1 is the temperature at Z1 Surface heat P2 SURFACE HEAT P2 is a vector Z2 U2 T2, where Z2 is the altitude, U2 is the velocity at Z2 and T2 is the temperature at Z2. Z2 must larger than Z1 and U2 must larger than U1. Mean surface heat flux MEAN_SURFACE_HEAT_FLUX is the heat flux in W/m2 from the ground to the flow. This is a parameter in the boundary layer profiles for the unstable Pasquill class A, B, and C. MEAN_SURFACE_HEAT_FLUX will not apply as heat contribution from the ground to the flow in the simulations. Pasquill class Pasquill atmospheric stability classes is a method of categorizing the amount of atmospheric turbulence present. Pasquill (1961) made six stability classes where: A is very unstable B is unstable C is slightly unstable D is neutral E is slightly stable F is stable An overview when to apply the different Pasquill classes is given in table: Pasquill stability classes. Wind speed
Day, strong sun
Day, moderate Night, clouds sun >50% B E B-C E B-C D C-D D D D Table 3.5: The Pasquill stability classes
< 2m/s 2 - 3 m/s 3 - 5 m/s 5 - 6 m/s > 6 m/s
A A-B B C C
Night, clouds <50% F F E D D
Ground roughness conditions At present the only option is: GROUND_ROUGHNESS_CONDITON
3.7.7
"RURAL"
Pressure relief panels
Pressure relief panels are commonly used in the process industry as a mitigating device in the case of an explosion. When the pressure forces on the panel exceed a certain limit, the panel yields and the pressure is relieved. In FLACS, special planar elements inside the simulation volume are FLACS v9.0 User’s Manual
3.7 Scenario menu
79
used to model the effect of explosion relief panels and yielding walls. They are specified as boxes where one, and only one, of the three side lengths must be zero. In addition several parameters such as the yield pressures, area porosities and drag factor may be specified. A panel may be active or passive. An active panel will initially modify the porosity in the region it covers, and again when the pressure difference over the panel exceeds the given limit. If the initial porosity is set to 0 and the final porosity is set to 1 the panel will start as being closed and end up as being open. This is how the behaviour of a real yielding panel is imitated. A passive panel will never modify the porosity, but it may be used to monitor the same data as an active panel. There are 5 different panel types (active panels): Panel type UNSPECIFIED POPOUT HINGED PLASTIC OVERLAY
Description Panel with linear displacement, full parameter set Panel with a linear displacement, limited parameter set Panel with a rotational displacement, limited parameter set Simulates the presence of plastic sheets (commonly used in experiments) Panel properties are ’blended’ with existing equipment Table 3.6: Panel types in FLACS
Note that the PLASTIC and OVERLAY panels can not be chosen directly in CASD. It is possible to edit the scenario file manually to use these panel types, by initially creating an UNSPECIFIED panel and change it to the desired type. It is advised that only experts attempt to do this. The edges of a panel will be automatically adjusted to the closest grid line. It is advised that the user specify panels whose dimensions match the grid, either by adjusting the grid or modifying the panel. One topic which needs special attention is the presence of solid structures close to the panel area. For example structural beams may constitute quite large blockages which must be accounted for, especially since they also occur at the vent openings of typical offshore modules. In such cases the panel area must be defined using several panels with solid beams in between. Smaller beams, often an integral part of the support frame of the pressure relief panels, may be accounted for by adjusting the final panel porosity for the panel itself. In addition to different types of pressure relief panels, an INACTIVE panel used to monitor variables is incorporated in the FLACS code. Use an INACTIVE panel to measure for instance average pressure across surfaces such as decks and walls. The panel will be available under Single field scalar time output in the Scenario menu. Note that if several panels are specified, they should not overlap each other anywhere (in the model it is assumed that the panels are non-overlapping). The full list of parameters from the scenario file is shown below. The type of the panel determines which of these parameters are relevant. INSERT POSITION SIZE MATERIAL PANEL_TYPE OPENING_PRESSURE_DIFFERENCES INITIAL_AND_FINAL_POROSITY
1 0.0, 0.0, 0.0 1.0, 1.0, 0.0 BLUE UNSPECIFIED 0.0, 0.0 0.0, 1.0
(m) (m)
(bar) (-) FLACS v9.0 User’s Manual
80
CASD
WEIGHT DRAG_COEFFICIENT MAXIMUM_TRAVEL_DISTANCE SUB_SIZES
0.0 1.0 0.0 1.0, 1.0
(kg/m2) (-) (m) (m)
The parameters for specifying panel properties are described in the following sections.
3.7.7.1
Insert
Integer identifying the panel.
3.7.7.2
Position
Cartesian coordinates [m] of the corner of the panel (the corner with lowest value of the coordinate in each axis direction).
3.7.7.3
Size
The panel is assumed to be a plane of rectangular shape. The dimension [m] in each of the axis directions is given. One dimension should be zero, this shows how the panel is oriented. If e.g. the dimension in x-direction is zero, the normal of the panel points in x-direction. The other two dimensions should be positive. The panel edges will be adjusted to match the grid lines perfectly in the FLACS code, so it is advised that you only specify panels which coincide with the grid lines, in order to avoid any confusion concerning the geometrical representation of the panel.
3.7.7.4
Material
Colour used when visualizing the panel as part of the geometry considered.
3.7.7.5
Panel type
There are 6 different panel types described here. A panel is either active or passive. They are listed and described in more detail below. Inactive is a passive panel that does not affect the geometry or the flow field. such panels are used to monitor certain variables Unspecified is a panel with linear displacement when the panel yields Popout is a panel with a linear displacement when the panel yields Hinged is a panel with a rotational displacement when the panel yields Plastic is a panel that imitates a plastic sheet (commonly used in experiments) Overlay is a panel superimposed on other equipment FLACS v9.0 User’s Manual
3.7 Scenario menu
81
Inactive panel A passive panel, corresponding to the case when PANEL_TYPE equals INACTIVE, does not affect the numerical simulation in any respect (neither the geometrical description including the porosities, nor the flow field), it is only included in order to monitor variables related to the area which the passive panel occupies (e.g. area-averaged pressure over all the control-volume faces of the panel). Only the parameters INSERT, POSITION, SIZE, and MATERIAL are relevant when the PANEL_TYPE is INACTIVE.
Unspecified, Popout and Hinged panel An active panel (all the panel types except INACTIVE) will in general affect both the porosities representing the geometry and the flow field. For the active panels both an initial and a final area porosity value is given (open area divided by total area). The initial porosity corresponds to the state of the panel before it has started to yield due to external pressure forces. The final porosity corresponds to the state of the panel after it has yielded completely due to pressure forces. Note that each control-volume face covered by the panel is treated separately. Thus it is possible, depending on the scenario, that one part of the panel is open, while at the same time other parts are closed. The way of using the initial and final porosities in the FLACS code depends on the type of panel considered. In the case of the panel types UNSPECIFIED, POPOUT, and HINGED; for each controlvolume face which is covered by the panel, the porosity representing the panel is also the porosity of the control-volume face, regardless of any other equipment or structure included in the geometry.
Overlay and Plastic panel For a panel of type OVERLAY the value of the area porosity determined by the panel is multiplied by the value of the area porosity determined by other geometrical objects (in the absence of any panel) giving the effective area porosity, for each control-volume face covered by the panel. Let us consider an example where the final porosity of the panel is 1 (fully open). After the panel has yielded completely for all the control-volume faces covered by the panel, the area porosities have the same values that they would have if no panel was included in the geometry. So in this case other geometrical objects than the panel itself can not "blow out" together with the yielding panel (cf. the comments above). For a panel of type PLASTIC, and in the case of the initial area-porosity; for each control-volume face which is covered by the panel, the porosity representing the panel is also the porosity of the control-volume face, regardless of any other equipment or structure included in the geometry (as in the case of the panel types UNSPECIFIED, POPOUT, and HINGED). Usually the initial areaporosity of a PLASTIC panel is zero (completely blocked). If the PLASTIC panel has yielded completely, the area porosities have the same values that they would have if no panel was included in the geometry (as in the case of an OVERLAY panel with final area-porosity 1). Note that for a PLASTIC panel the final area-porosity is always assumed to be 1, regardless of the user-given value. Panels of the types PLASTIC or OVERLAY are assumed to be light panels, and inertial forces are neglected by the FLACS code during the dynamical process when the panel yields due to pressure forces. In the case of panels of types UNSPECIFIED, POPOUT, or HINGED, inertial forces are in general taken into account. Attention: The PLASTIC and OVERLAY panels cannot be chosen in the present version of CASD. It is possible to edit the scenario file manually to use these panel types, by initially creating an UNSPECIFIED panel and change it to the desired type. It is advised that only experts attempt to do this. FLACS v9.0 User’s Manual
82
CASD
The +IMP panel types The standard panel types will yield at a give negative or positive pressure. In some situations it can of necessary for panels to also take into account the pressure impulse. To activate this the user can add the suffix +IMP to the end of the panel type string, e.g PANEL_TYPE OPENING_PRESSURE_DIFFERENCES
POPOUT+IMP 0.05 0.005 (bar
bar*s)
When using this the input to OPENING_PRESSURE_DIFFERENCES changes to • Parameter 1 is yield pressure for both negative and positive direction [bar] • Parameter 2 is yield pressure impulse both negative and positive direction [bar∗s] The panel will yield if both the pressure and pressure impulse yield values are exceeded. This is valid for panel types HINGED and POPOUT. Attention: This functionality is not available through CASD. The user must manually edit the cs-file. 3.7.7.6
Opening pressure differences
When the net pressure over the panel (pressure on the negative side relative the coordinate axis minus pressure on the positive side, i.e. the net pressure from the fluid acting on the panel) exceeds limits given by the user, the panel starts to yield. Two opening pressure-differences, in units of [bar], are given by the user, the first corresponds to the case when the net pressure acts in negative direction, the second when the net pressure acts in positive direction. If for example the OPENING_PRESSURE_DIFFERENCES are given by -0.1 and 0.2, the panel starts to yield when the net pressure is less than -0.1bar (negative direction) or greater than 0.2 bar (positive direction). Note that the FLACS code makes the same interpretation of the input data, regardless of the sign of the opening pressure-differences given by the user. Hence, if the OPENING_PRESSURE_DIFFERENCES were given instead by 0.1 and 0.2, or by 0.1 and -0.2, the effect on the numerical simulation would be exactly the same. The opening pressure-differences must be given for the active panels (all the panel types accept INACTIVE). 3.7.7.7
Initial and final porosity
The initial area porosity [-] of the panel (open area divided by total area) corresponds to the state of the panel before it has started to yield due to external pressure forces. Similarly the final porosity corresponds to the state of the panel when it has yielded completely. In the case of a PLASTIC panel, the final area-porosity should be set to 1 by the user in order to correspond to how the FLACS code works (the final area-porosity is always assumed by FLACS to be 1, regardless of the user-given value). Note that the way of using the initial and final porosities in the FLACS code depends on the type of panel considered. See the comments above. The initial and final porosity must be given for the active panels (all the panel types accept INACTIVE). The value of the porosity range from 0 (totally blocked) to 1 (fully open). It is expected that the initial porosity is less than the final porosity (the net pressure opens up the panel). If this is not the case, the FLACS code gives a warning. 3.7.7.8
Weight
The parameter WEIGHT is specified for the panel types UNSPECIFIED, POPOUT, and HINGED (for the other panel types it is not relevant to specify this parameter). This parameter is the mass per FLACS v9.0 User’s Manual
3.7 Scenario menu
83
unit area [kg/m2 ] of the part of the panel that is moved when the net pressure acting on the panel exceeds the opening pressure-difference given by the user. The value of WEIGHT is either zero or positive. If the value of WEIGHT is zero, it is assumed in the FLACS code that the panel is so light that inertial forces can be neglected during the dynamical process when the panel yields due to pressure forces. If the value of WEIGHT is positive, the dynamics when the panel yields and part of the panel is accelerated and move away from its initial position, is modelled in the FLACS code. Note that when you specify the panel type POPOUT, the inertial forces must be there when the panel yields, that is the WEIGHT must be a positive number. Note that it is only the part of the panel that is accelerated and move away from its initial position, that should contribute to the mass per unit area. The mass of a rigid frame that does not yield, should not be included. 3.7.7.9
Drag coefficient
A dimensionless drag-coefficient [-] is given by the user for the panel type UNSPECIFIED (for the other panel types it is not relevant to specify this parameter). The drag coefficient is used when the drag force from the panel on the fluid is modelled (both before, during, and after the panel has yielded). The value of the DRAG_COEFFICIENT is zero or positive. A typical value is 2.0 (the value 2.0 is set by the FLACS code as a fixed preset value for the panel types POPOUT and HINGED).
3.7.7.10
Maximum travel distance
The maximum travel-distance [m] is given by the user for the panel type UNSPECIFIED (for the other panel types it is not relevant to specify this parameter). The use of the maximum travel-distance is based on a rough approximation. The maximum travel- distance is the smallest distance from the initial position of the panel to the position where the yielded panel no longer affects significantly the effective area-porosity at the initial position of the panel. The maximum travel-distance is used to model the effective area-porosity at the initial position of the panel during the dynamical process when the panel yields. A typical value of the MAXIMUM_TRAVEL_DISTANCE is in the order of 1m (for the panel types POPOUT and HINGED there is an in-built model in the FLACS code for the effective area-porosity that does not need a value of the MAXIMUM_TRAVEL_DISTANCE as input from the user). 3.7.7.11
Sub sizes
The parameter SUB_SIZES is given by the user for the panel types POPOUT or HINGED (for the other panel types it is not relevant to specify this parameter). These panel types are assumed to consist of sub-panels mounted on a frame. The width [m] and the height [m] of the sub-panels are given by the user (it is assumed in the numerical model that all the sub-panels are of uniform size). Both the width and the height should be positive. In the case of the panel type HINGED it is assumed that each sub-panel turns on a hinge when it yields. The width of the sub-panel is defined as the dimension in the direction normal to the axis of the hinge. It is important (since it affects the model in FLACS) to specify the width of a HINGED panel first and then the height, in the parameter SUB_SIZES (for a POPOUT panel it is of no importance which dimension that is defined as the width).
3.7.8
Gas composition and volume
This section allows you to define a box shaped cloud region and the gas concentration and composition. The menu is shown below: FLACS v9.0 User’s Manual
84
CASD
POSITION_OF_FUEL_REGION DIMENSION_OF_FUEL_REGION TOXIC_SPECIFICATION VOLUME_FRACTIONS EQUIVALENCE_RATIOS_(ER0_ER9)
0.0, 0.0, 0.0 0.0, 0.0, 0.0 ""