关注phoenics 上海飞熠软件技术有限公司 TUTORIAL- Simulation of a Labyrinth Flow.FLAIR Tutorial 1: Investigating library case I203 FLAIR Tutorial 2: A room with two radiators Tutorial 3 Comfort index in a roomFLAIR Tutorial 4: Fire in a room关注phoenics 上海飞熠软件技术有限公司 TUTORIAL- Simulation of a Labyrinth Flow.This tutorial shows how to set up and run the simulation of steady turbulent flow of air with heat exchange. The geometry is 2-dimensional, of (y-direction) depth 1.0m, with the other dimensions shown in the figurebelow:wherein:•••••••• INLET is an aperture through which air at 20 degrees C enters at a prescribed velocity of 0.5m/s, giving an inlet Reynolds Number of 15000; OUTLET is a fixed-pressure aperture, where the pressure is held close to 0.0 Pa, relative to the external pressure;WALL-W is an adiabatic boundary wall; WALL-E is an adiabatic boundary wall; IN-PLATE is an internal plate, the only function of which is to obstruct the flow; IN-BLOCK is a block of copper; H-BLOCK is a block of aluminium; L-BLOCK is a block of aluminium, supplied uniformly with 100 watts of electrical heat.Preliminary remarks1. The tutorial must therefore explain how to:o introduce objects comprising blocks of different materials;关注phoenics 上海飞熠软件技术有限公司 o provide them with heat sources, if required;o introduce also objects (inlets, outlets, plates) having no volume, which however affect the flow in prescribedways;o prescribe the sizes and positions of each object.2. In addition to physically-meaningful data of the above kind, all CFD codes, and therefore PHOENICS also, requirenon-physical inputs such as:o a grid of computational cells;o instructions about how many iterations to make of the trial-and-error solution process; o how much data to print.3. The non-physical settings suggested below are such as to make the calculations rapid rather than accurate.4. Other tutorials will be provided which will enable you to explore the influences of the settings on accuracy and speed ofcalculation.Starting the tutorialFirst close any windows which are already running PHOENICS modules.Then activate the PHOENICS Satellite module, in its VR-Editor mode, by either:1. clicking on the PHOENICS VR icon on the desktop; or2. entering the Run Modules panel of the PHOENICS Commander, and then clicking the VRE button; or3. typing 'vre' at the Command prompt.Your working directory, i.e. the folder in which will appear all the files which the tutorial creates, will be, according to the method which you have used:1. /PHOENICS/d_priv1 [You can change this, if you wish, by clicking the 'options button' at the top of the satellite panel;or2. that which appears after: set globdata(cwd) in the file:/PHOENICS/d_pc/userpref/pcprefs . [You can change this, ifyou wish, by clicking the 'Customise' button of the PHOENICS Commander Top Page]; or3. that of the prompt itself.In PHOENICS-VR Editor•• When the Editor starts execution, what it shows on the screen at first will depend on what happens to be in your working directory. You can disregard this. [If you see the words: 'Display From Current Q1', click the top-right X so as to close the panel ] Make a fresh start by:o clicking first on 'File';o then on 'Start New Case';o then on 'Core';o then on 'OK' to confirm the resetting.You are now ready to begin. •关注phoenics 上海飞熠软件技术有限公司 Set the domain size••• Click on 'Menu' or 'M' (according to whether the 'C ontrol panel' is visible) [If it is not, you can make it so by clicking on 'View']. Set 'Simulation of a Labyrinth Flow' as the Title. Click on 'Geometry'.Note that clicking on '?' at the top right of any panel and then on a menu button will provide 'help'.•• Change the X-Domain Size to 2.0 m. Click 'OK' to close the Grid mesh settings dialog. Select the working fluid• Click on 'Properties' which will elicit a page like this:• Note that关注phoenics 上海飞熠软件技术有限公司 o the properties are already those of "Air at 20 deg C, 1 atm" by default.o The ambient (external) temperature is 20 degC, and the ambient pressure is 0 Pa relative to 1atm.•• Since these correspond to the present case, there is no setting to make. Click on 'Top menu' to accept.Indicate for which variables solution is required••••• Click on 'Models'. Leave the 'solution for velocities and pressure' ON. Click on'energy equation', and select 'Temperature'. Click on 'OK'. Leave the turbulence models setting as 'KECHEN'. This indicates that the 'KE' and 'EP' (kinetic energy and dissipation rate of turbulence) are to be solved by the CHEN method. Click on 'Top menu', then on 'OK' to exit the Main Menu.Re-sizing the domainThe domain no longer fits the screen; to resize, click on the pull-down arrow next to the 'R' icon on the toolbar; then 'Fit to window'.Introduce the required objects•• Click on the 'object-management' button ('O' on the toolbar or on the Control panel, also sometimes called 'hand set'). This will display the 'object-management' window with a list of objects of which only 'DOMAIN' is currently present.关注phoenics 上海飞熠软件技术有限公司 •• You will now introduce the objects which you require, one-by-one, giving each a position and size, corresponding to those in the above diagram. You will also introduce 'attributes', for example material, velocity or temperature, so completing the description of thesituation to be simulated.1. Introduce the H-BLOCK objecto In the Object-management dialog, click on 'Object', 'New' and 'New Object'. o Change name to H-BLOCK (which stands for 'high block'). o Click on 'Size' and set SIZE of object as:Xsize: tick 'To end'Ysize: tick 'To end'Zsize: 0.05o Click on 'Place' and set Position of object as:Xpos: 0.0Ypos: 0.0Zpos: Tick 'At end'o Click on 'General'.Select Type: Blockage (default).To change the material, click on 'Attributes', click on 'Types', click on SOLIDS then'OK', and then on 'ALUMINIUM at 27 deg C' and 'OK';. Click on 'OK' to exit the Attributes menu, and on 'OK' to close the Object Specification Dialogue Box. o H-BLOCK will now appear in the Object-Management list of objects.2. Introduce the L-BLOCK objecto Click on 'Object', 'New' and 'New Object'. o Change name to L-BLOCK (which stands for 'low block'). o Click on 'Size' and set SIZE of object as:Xsize: tick 'To end'Ysize: tick 'To end'Zsize: 0.05o Note that no position coordinates need to be given because the defaults (0.0, 0.0, 0.0) suffice in this case. o Click on 'General'.Define Type: Blockage (default).Click on 'Attributes', select material as ALUMINIUM as above. Set a HEAT SOURCE of 100 W. Click on 'Adiabatic', and select 'Fixed heat Flux'. Click on Value andenter 100.Click on 'OK' to return to the Object Specification Dialogue Box, and on 'OK' to close the ObjectDialogue Box.3. Introduce the WALL-E objecto Click on 'Object', 'New' and New Object'. o Change name to WALL-E (which stands for 'east wall'). o Click on 'Size' and set SIZE of object as:关注phoenics 上海飞熠软件技术有限公司 Xsize: 0.0Ysize: tick 'To end'Zsize: tick 'To end'o Click on 'Place' and set Position of object as:Xpos: tick 'At end'Ypos: 0.0Zpos: 0.5o Click on 'General'.Define Type: Plate.Note that a Plate differs from a Blockage in having zero thickness. Click on 'OK' to close the Object Specification Dialogue Box.4. Introduce the WALL-W objecto Click on 'Object', 'New' and New Object'.o Change name to WALL-W (which stands for 'west wall').o Click on 'Size' and set SIZE of object as:Xsize: 0.0Ysize: tick 'To end'Zsize: 0.5o Note that no position coordinates need to be given because the defaults (0.0, 0.0, 0.0) suffice in this case. o Click on 'General'.Define Type: Plate.Click on 'OK' to close the Object Specification Dialogue Box.5. Introduce the IN-PLATE objecto Click on 'Object', 'New' and New Object'.o Change name to IN-PLATE (which stands for 'interior plate').o Click on 'Size' and set SIZE of object as:Xsize: 0.0Ysize: tick 'To end'Zsize: 0.65o Click on 'Place' and set Position of object as:Xpos: 0.5Ypos: 0.0Zpos: 0.3o Click on 'General'.Define Type: Plate.Click on 'OK' to close the Object Specification Dialogue Box.6. Introduce the IN-BLOCK objecto Click on 'Object', 'New' and New Object'.o Change name to IN-BLOCK (which stands for 'interior block').o Click on 'Size' and set SIZE of object as:关注phoenics 上海飞熠软件技术有限公司 Xsize: 0.4Ysize: tick 'To end'Zsize: 0.65o Click on 'Place' and set Position of object as:Xpos: 1.0Ypos: 0.0Zpos: 0.05o Click on 'General'.Define Type: Blockage (default).Click on 'Attributes', then on 'Other materials', select 'Solids' then COPPER. Click on 'OK' to return to the Object Dialogue Box, o and on 'OK' to close the Object Dialogue Box.7. Introduce the INLET objecto Click on 'Object', 'New' and New Object'. o Change name to INLET.o Click on 'Size' and set SIZE of object as:Xsize: 0.0Ysize: tick 'To end'Zsize: 0.45o Click on 'Place' and set Position of object as:Xpos: 0.0Ypos: 0.0Zpos: 0.5o Click on 'General'.Define Type: Inlet.Click on 'Attributes' and set the velocity in X-direction to 0.5 m/s. Leave the temperature at 'Ambient', 20 degrees C for the inflowing air. Click on 'OK' to close the Attributes menu, and on 'OK' in the Object Specification Dialogue Box.8. Introduce the OUTLET objecto Click on 'Object', 'New' and New Object'. o Change name to OUTLET.o Click on 'Size' and set SIZE of object as:Xsize: 0.0Ysize: tick 'To end'Zsize: 0.45o Click on 'Place' and set Position of object as:Xpos: tick 'At end'Ypos: 0.0Zpos: 0.05o Click on 'General'.Define Type: Outlet.Leave the default values in the Attributes dialogue box.关注phoenics 上海飞熠软件技术有限公司 Click on 'OK' to exit the Object Dialogue Box.Inspect the 'default' gridThe PHOENICS satellite sets a computational grid which has sufficient lines at least to co-incide with the extremities of all the objects which have been introduced, with some additional lines in-between.This is known as the 'default' grid. It can be displayed by pressing the 'Mesh toggle' button, which in the present case should generate a picture such as the following.• Click anywhere on the image, and the 'Grid Mesh Settings' dialog box will appear.The grid in all three directions is set to 'Auto'. This gives 20 cells in X and Z, and 1 cell in Y. New users are advised not to make any changes at present.The grid is too coarse (i.e. lines are too few) for accuracy but will suffice for the tutorial. Click on 'OK' to close the dialog box.• Click on 'Mesh toggle' again to turn off the mesh display.The remaining solution-control parametersIn later tutorials, you will be invited to click on 'Menu' and then on 'Numerics', in order to inspect and modify other solution-control parameters. However, the defaults which have been set will suffice for the present purposes.关注phoenics 上海飞熠软件技术有限公司 Setting the Probe Location• Before running the solver, it is a good idea to place the probe in a suitable place to monitor the convergence of thesolution. Too close to an inlet, and the value will settle down very quickly before the rest of the solution. Placed in a recirculation zone, it may still show traces of change even though the bulk solution is converged. In this case, to the right of the IN-BLOCK object is fine.• Click on the probe icon on the toolbar or double-click the probe itself, and move the probe to X=1.8, Y=0.5, Z=0.3. Running the SolverYou are now ready to perform the simulation by activating the PHOENICS solver module (known as 'Earth'). In order to do so, click on 'Run', 'Solver'(Earth), and click on 'OK' to confirm.You will then see, in quick succession, indications on the screen of:••• the activation of the module; the reading by the solver of the input data which your tutorial session has prepared for it; the progress of the calculation, represented by the generation of curves which look like this:关注phoenics 上海飞熠软件技术有限公司 •Click the 'END' button to terminate the calculation and to return to the VR Environment. You are now able to launch the VR Viewer, in order to inspect the results of your calculation.Using the VR Viewer.Click on 'Run', 'Post processor', then GUI Post processor (VR Viewer) . Click 'OK' on the file names dialog to accept the default files.• To view:oVelocity vectors - click on the 'Vector toggle' . Vectors are coloured by the current plotting variable, but theirlength is always related to the absolute velocity.oContours - click on the 'Contour toggle'oStreamlines - click on the 'Streamline management' buttonoIso-surfaces - click on the 'Iso-Surface toggle'• To select the plotting variable:关注phoenics 上海飞熠软件技术有限公司 o To select Pressure - click on the 'Select Pressure button' oTo select Velocity - click on the 'Select Velocity button' .oTo select Temperature - click on the 'Select Temperature button' oTo select any other variable - click on the 'Select a Variable button' •• To change the direction of the plotting plane, set the slice direction to X, Y or Z To change the position of the plotting plane, move the probe using the probe position buttons.Alternatively, click on the probe iconProbe Location dialog.• on the toolbar or double-click the probe itself to bring up the A typical plot from this case is:关注phoenics 上海飞熠软件技术有限公司Checking the S ource Balance •• It is very important to know whether the inflows and outflows of mass and energy are in balance. If they are, it is a good sign that the solution is convergent. If they are not, the solution is definitely not converged. For further information on the assessment of Convergence, see the lecture . Open the Object-Management dialog, and right-click on the Domain entry. From the context menu select 'Showresults'.关注phoenics 上海飞熠软件技术有限公司•• This will display the sources and sinks of all variables. The section showing 'Nett source of R1 at...' shows the mass source in kg/s at each inlet and outlet.•• Positive values are inflo ws,negative values are outflows. The 'nett sum' at the end of the section should be close to zero, as all the mass entering must leave. The section showing 'Nett source of TEM1 at...' shows the energy source in Watts for each heat source or sink.关注phoenics 上海飞熠软件技术有限公司•• Positive values are sources, negative values are sinks. The energy source at the inlet is mass flow*specific heat*inlet temperature in Kelvin. The 'nett sum' should also be close to zero. In this case, it is a bit high, showing that convergence is not complete. These balances can also be checked by inspecting the RESULT file. This contains an echo of the inputs, a selection ofthe solution and the source balances. Click on 'File', 'Open file for editing', then 'Result'. Scroll down the file until you reach the section headed 'Sources and sinks'.Saving the results.• In the VR Viewer click on 'File' then 'Save as a case'; create a new folder called'LABYRINTH ' and save as 'CASE1'.These actions will place in the folder the following files:o case1.q1, a copy of the Q1 file created in the tutorial;o case1.res, a copy of the RESULT file printed by the solver at the end of its execution; o case1.phi, a copy of the PHI file which the VIEWER uses as the basis of its displays; o case1.gif, an image representing the approach to convergence, such as was shown above; ando case1.pat and case1.pbc, which there is no need to describe now.These will be discussed in later tutorials.关注phoenics 上海飞熠软件技术有限公司 FLAIR Tutorial 1: Investigating library case I203This tutorial is to illustrate how to load case I203, Ventilation of Hackney Hall from the FLAIR library andinvestigate its settings, run the case and view the results. However, this is a 3D real application and the user does not have to run through the case as it may take a few hours to complete.Starting FLAIR with the default room••• Click the FLAIR icon (a desktop shortcut created by the FLAIR installation program); or Start the VR-Editor by clicking on 'Start', 'Programs', 'FLAIR', then 'FLAIR-VR'. Click on the 'File' button and then select 'Start new case', followed by 'FLAIR' and 'OK'. The FLAIR VR-Environment screen should appear, which shows the default room with the dimensions10mx10mx3m.Loading the library case I203•Click on the 'File - Load from libraries'. This will bring up the 'load from library' dialog as shown below•• If you If you know the number of the case you want to load, type it into the 'case number' data entry box and click OK Otherwise, click on 'Browse' and you will see the following dialog box关注phoenics 上海飞熠软件技术有限公司••Use the '+' next to Special-purpose programs to expand the library listings Click on'Flair library' and 'Large cases'. Case I203 will be on the listings• Select cas e I203 and click 'OK' to confirm the selected case number, and then OK to load it. There willbe a display window which provides a link to the description. Click on the link, you will see the following details:The case was three dimensional, and the flow was steady.The geometry within the domain was created using PHOENICS-VR. As the auditorium was very nearly symmetrical, only half of the domain was modeled.The theatre includes balconies modeled as a series of stepped Cartesian blockages, and a new object was created to model the curved roof and stage.The computational domain is 33m long by 11.4m wide by 12.7m high, a 60*20*50 cells Cartesian grid was used, and PARSOL technique was used to map the non- Cartesian elements of the geometry; e.g. the curved roof and stage.The domain is filled with air, and ideal gas law was used in calculating pressure-temperature-density relationships. The roof, floors and walls were modeled using a material with the following properties; density: 1700kg/m***3, specific heat:800J/kg/degK, thermal conductivity: 0.62W/m/degK. • • • • •关注phoenics 上海飞熠软件技术有限公司 •••••• The effects of gravity were included by way of sources of momentum in the equations for the z component velocity using the built in buoyancy facilities available in PHOENICS. The LVEL model of turbulence is used to close the Reynolds-averaged Nervier-Stokes equations. The air supply under the balcony, on the long side of the theatre, is introduced at 12degC. The incoming air velocity, from this long, slot is fixed at 3m/s, giving a mass flow rate of about 2.5m***3/s. Further air is supplied under the seating, at the rear of the auditorium only, at 18degC. The inflow supply rate is0.7m***3/s, providing by seven 70mm high inlets. The occupant and light gains were modeled using heat sources respectively equivalent to 145W/m**2 and 50W/m**2. The ten outlets at the top of the domain were modeled as 500mm long by 250mm wide fixed pressureopenings.You may close the window after viewing the description. After successful loading, the geometry of caseI203 will appear on the screen as shown below •This case was to simulate the air flow and thermal characteristics in the theartre influenced by theventilation system.Investigating the model settings关注phoenics 上海飞熠软件技术有限公司 •Click on the 'Object management' button,management dialog box as shown below. . You will see all the objects listed on the object• You may double click on each of those objects and investigate their attributes. For example, you doubleclick on object IN37 (reference no. 36). This will bring up the 'Object specification' dialog box as shownbelow.• Click on 'Attributes' button to bring up the 'Object attributes' panel as shown below. The settings showthat the air is coming into the theatre at the speed of 0.16 m/s with the temeparture of 18 degree C.关注phoenics 上海飞熠软件技术有限公司• You can investigate the physical model settings by clicking on 'Main menu' button, and click oneach model panel, for example 'Models' as shown below. As can be seen, the pressure, velocities and temperature are solved. The constant effective turbulence viscosity is used for the turbulence modelling.关注phoenics 上海飞熠软件技术有限公司Running the solver•To obtain a solution for the loaded case, click 'Run - Solver' as shown below. Then click 'OK' to confirm.As the Earth solver starts and the flow calculations commence, two graphs should appear on the screen.The left-hand graph shows the variation of solved variables at the monitoring point that was set during the model definition. The right-hand graph shows the variation of errors as the solution progresses.A 3D real application like this requires 5-6 hours of computing time on a 3Mhz PC. Viewing the results• To run the Viewer, click on 'Run - Post processor - GUI post processor (VR Viewer)'. 关注phoenics 上海飞熠软件技术有限公司 • When the 'File names' dialog appears, click 'OK' to accept the current result files. You may use the'Object management' dialog box to hide some of objects. Click on the 'T' (Select temperature) button,Click on 'C', (select a variable) button,dialog box as shown below on the main control panel. • •on the main control panel to bring up the 'Viewer options'•••• Set the 'Maximum value' to 24 C Click on 'Streamline management' button, . Move the probe position to X = 26; Y = 6.0; Z = 2.46. Click on 'Object' menu and click on'Options' to bring up the 'Stream Options' dialog box and choose thefollowing settings关注phoenics 上海飞熠软件技术有限公司••Click on 'OK' to return to the 'Streamline management' dialog box. Click on 'Object' and select 'New' as shown belowThe following streamlines will be displayed on the screen.关注phoenics 上海飞熠软件技术有限公司 FLAIR Tutorial 2: A room with two radiatorsThis case models a three-dimensional room, consisting of two closed windows, a door. The radiators underneath each window produce a total heat flux of 500W and the walls are fixed at 20 degree C. TheIMMERSOL model is used. The 'Duplicate object' function is used for the creation of the second window and radiator. The material of the radiators is selected from the property database.The following picture shows the geometry.Setting up the modelStart FLAIR with the default roomClick FLAIR icon (a desktop shortcut created by the FLAIR installation program); or • Start the VR-Editor by clicking on 'Start', 'Programs', 'FLAIR', then 'FLAIR-VR'. • Click on the 'File' button and then select 'Start new case', followed by 'FLAIR' and'OK'. •关注phoenics 上海飞熠软件技术有限公司 The FLAIR VR-Environment screen should appear, which shows the default room with the dimensions 10mx10mx3m.Resizing the room• Change the size to 3m in x-direction, 5m in y-direction and 2.7m in z-direction respectively onthe control panel. Alternatively, the user may click on the main 'Menu' and then bring up the'Geometry' panel to change the room size.Click on "Reset" button,'OK'. on the movement panel, and then click on "Fit to window" and •Adding objects to the rooma. add a door• Click on the 'Object' menu and select 'New' (New Object) on the 'Object management' dialogbox to bring up the 'Object specification' dialog boxChange name to 'DOOR'.• Click on the 'Size' button and set 'Size' of the object as: •X: 1.0 mY: 0.0 mZ: 2.0 m• Click on the 'Place' button and se t 'Position' of the object as:X: 2.0 mY: 0.0 mZ: 0.0 mClick on 'General' button and set value of 'Type' to be 'PLATE'.• Click on 'Attributes', and leave the default setting 'Adiabatic'• Click on 'OK' to return to the 'Object specification' dialog box. Click on 'OK' again. •b. add the first window• Click on the 'Object' menu and select 'New' (New Object) on the 'Object management' dialogbox to bring up the 'Object specification' dialog boxChange name to 'WIND1'.• Click on the 'Size' button and set'Size' of the object as: •关注phoenics 上海飞熠软件技术有限公司 X:0.00 mY:0.98 mZ:0.87 m• Click on the 'Place' button and set 'Position' of the object as: X: 0.0 mY: 1.28 mZ: 1.0 m• Click on 'Shape' and a click on 'Geometry' selection will bring up the following dialog boxSelect WINDOW1 from d_object/public/furnture and click on 'Save' • Click on'General' button and select 'PLATE' from the object 'Type' list . • Click on 'Attributes', and then set the 'Energy source' to 'Surface temperature' and value to 15degree C.• Click on 'OK' to return to the 'Object specification' dialog box. Click on 'OK' again to exit from the'Object specification' dialog. •c. add the second window by duplicating the first window••Click on the WIND1 object to highlight it on the 'Object management' dialog. Click on 'Duplicate object' button, . followed by clicking 'OK' to confirm duplicating. Note that the newly duplicated object will be at the same place as the original one.关注phoenics 上海飞熠软件技术有限公司 Click on 'Object' menu and select 'Refresh' on the Object Management dialog box. This will bring the newly duplicated object into the list of objects.• Double click on the newly dup licated object to bring up the "Object specification" dialog box. • Change the name to WIND2• Click on the 'Place' button and set 'Position' of the object as: •X: 0.0 mY: 2.92 m。
