web-archive-uk.com


Web directory, archive
Search web-archive-uk.com:


Find domain in archive system:
web-archive-uk.com » UK » C » CHAM.CO.UK

Total: 682

Choose link from "Titles, links and description words view":

Or switch to "Titles and links view".
  • MOVING.HTM
    face areas cell volumes angles and direction cosines which PHOENICS needs The user does not have to attend to these matters If as is usual the fluid can penetrate the cell faces the convection fluxes which enter the balance equations of PHOENICS depend on the grid velocity as well as on the fluid velocity PHOENICS calculates the grid movement contribution by way of the swept volume of the cell face This depends on the following twenty four 24 quantities the 3 cartesian coordinates of each of the four corners of the cell face at the start of the time step and the corresponding coordinates of the corners at the end of the time step The swept volume calculation is effected by a call to the EARTH sub routine VOLUM which is to be seen in Group 8 of GREX3 and it leads to values of the arrays CONI CONJ and CONK for which three dimensional storage should be provided by appropriate settings in the Q1 file A difficulty about this practice is that VOLUM always returns a positive quantity it is therefore necessary to determine independently whether the swept volume tends to increase the volume of the cell in question or to decrease it In the former case the relevant CON eg CONI for an east face must be given a negative sign otherwise the sign should be positive The explanation is that the CONs have the same signific ance as those based upon the fluid velocity in respect of tendency to increase the mass of material in a cell A positive u velocity at an east face gives a positive CON and a positive cell wall velocity has the opposite tendency Examples of how the sign may be fixed are to be found in GREX3 where the determination is

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_enc/moving.htm (2016-02-15)
    Open archived version from archive


  • VRSRCINF.HTM
    be set as SOURCE of TEM1 at PIPE1 is 1 e 5 ZG with VOLU SOURCE of TEM1 at PIPE2 is 1 e 5 ZG with VOLU etc The formula for Nusselt number is Nu 0 3 Re 0 6 Pr 1 3 where the Prandtl number equals 0 7 For a data of this problem the Nusselt number equals 1 9 Also the In Form is used for calculation of average Nusselt number without the creation special Ground code similarly to the previous example The results are dumped in a special INFOROUT file as SOLID TEMP SUM 385547 SOLID NUMBER SU 9370 00 AVERAG SOL TEMP 41 1470 TOTAL SOLID VOL 8 886132E 04 AVERAGE Nusselt 3 38313 Calculated average Nusselt number equals 3 38 Graphical display of results The temperature field on the XY plane is here As in the previous example the results have shown insignificant change of the temperature field in the Z direction The complete Q1 file is library case v147 2 3 The simulation of a heat exchange by heat conductivity about filament Description This case modeling a heat exchange by heat conductivity about a filament The filament is simulated by the cylinder shape with properties of steel It is surrounded by air The air temperature has constant significance 18 degree at all domain boundaries The filament is heated up by a electrical current which occurs in the conductor at presence the filament votage The emitted heat is calculated as the filament votage divided by thermal resistance which is function of the conductor temperature The VR editor can not set any form of dependence of a source from solved or stored variables In this case it should use In Form In Form statements in the Q1 file The heating source per unit volume of the HEATING VR object can be set as INFORM13BEGIN SOURCE of TEM1 at HEATING is 220 2 0 1 TEM1 with VOLU STORED of PRPS at HEATING is 111 with TSTSTR INFORM13END In Form can be used for the setting or change of the material number of VR object In this case in the beginning the material of conductor is described as the steel In other words PRPS variable assigns significance 111 in cells located inside HEATING VR object Graphical display of results The temperature field is here The Q1 file of this problem data library case v148 2 4 The simulation of a heat exchange about hot pipe Description In this example modeling of a heat exchange by convection is considered at a cross flow of cold air about a single hot pipe The cold air with temperature of 18 degrees moves into the domain in a direction of a X axis The flow is laminar The Reynolds number equals 26 The pipe with a hot heat transport medium is simulated by the cylinder shape with properties of steel The direction of a symmetry axis of a pipe coincides a direction of Z axis The total heat at a cylinder object is calculated from formulas for the heat tranfer from hot air to cold air through a pipe wall as QFlux U Area Outside Tem Hot Tem Wall where the heat transfer coefficient from hot air to a pipe wall is U 1 Area Outside Alfa Area Inside 1 Cond Steel Thick Wall the heat exchange coefficient from hot air to a pipe wall is Alfa Nu Cond Air 2 Radius Inside and Nu 4 364 for a constant heat flux from hot air to the pipe wall The heating source per whole PIPE object can be set by In Form as SOURCE of TEM1 at PIPE is QFlux with WHOLOB LINE where LINE flag is used for setting a linearised source For calculation of Nusselt number can use the formula of which is fair for a external cross flow of a round pipe located in depth of a pipe array Nu 0 3 Re 0 6 Pr 1 3 where the Prandtl number equals 0 7 For a data of this problem the Nusselt number equals 1 9 The diameter of a cylinder can be used the characteristic size Nu Alfa Diam Inside Cond Air It is possible to calculate Nusselt number by In Form without programming of a special Ground code by next formula Alfa Q total Area pipe T Wall T cold where Q total is a total heat flux Area pipe is a total pipe area T Wall is a average pipe temperature T cold equals 18 In Form statements in the Q1 file The appropriate In Form statements are as follows Echo InForm settings for Group 13 INFORM13BEGIN Pi Pi number REAL Pi Pi 3 14159 PLEN Pipe length REAL PLEN PLEN 0 1 RadOu The outside radius of the pipe REAL RadOu RadOu 0 02 RadIn The inside radius of the pipe REAL RadIn RadIn 0 015 ArOu The area of outside pipe surface REAL ArOu ArOu 2 Pi RadOu PLEN ArIn The area of inside pipe surface REAL ArIn ArIn 2 Pi RadIn PLEN Thick The thickness of a pipe wall REAL Thick Thick RadOu RadIn CondA The thermal conductivity of the air REAL CondA CondA 0 0258 CondS The thermal conductivity of the steel REAL CondS CondS 43 Nu 4 364 for a constant heat flux from hot air to the pipe wall REAL Nu Nu 4 364 Alfa The heat exchange coefficient from hot air to a pipe wall REAL Alfa Alfa Nu CondA 2 RadIn U The heat transfer coefficient from hot air to a pipe wall REAL U U 1 ArOu Alfa ArIn 1 CondS Thick Thot The average temperature of hot air inside a pipe REAL Thot Thot 50 QFlux The heat flux from hot air to cold air through a pipe wall CHAR QFlux QFlux U ArOu Thot TEM1 where TEM1 is the temperature of a pipe wall QFlux SOURCE of TEM1 at PIPE is QFlux with WHOLOB LINE REYNO Reynolds number REAL REYNO WIN DIAM

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_enc/vrsrcinf.htm (2016-02-15)
    Open archived version from archive

  • The PHOENICS FacetFix Utility
    1 stl building 0 2 stl etc If a name has been given the sequence number will be appended to that In addition a file named outputfile 0 stl is also created This contains a single facet which sets the overall size of all the sub objects It can be used by the Editor to place all the sub objects in their correct relative positions Once the files have been selected the Run FacetFix button at the bottom of the dialog launches facetfix exe with arguments derived from the contents of the dialog window 4 Filtering the data Optionally the user might want a filter file to define multiple sub regions for extracting parts of the model individual buildings etc Rather than running repeatedly with different inputs in the filter box inputs the sub regions can be specified in a filter file Specify the filterfilename with the name and path if necessary of your filter file if no filter file is specified and the filter box entries are blank then all the facets in the model file will be used and mended The filter file allows sections of the model to be output from the original data perhaps a single building from a town Example 1 runfacet with input file building stl will read the whole building stl file and correct facet direction and mend holes and output a single dat file A large case with 1 400 000 facets takes about 1 minute on a 3GHz P4 machine Example 2 runfacet with filter file filters tolerance 0 0001 and input file building stl will extract buildings in areas specified in the file filters Any file name without spaces in is satisfactory it is suggested that a name meaningful to the user with the extension fil should be used In the example if the domain of an object file is about 3000 metres long then the tolerance option will consider vertices closer than 0 3 metres to be identical 5 Format of the filters file The file is ascii Comment lines start with a Other lines define filter volumes there is no limit to the number of filter volumes which may be specified other than disk space memory size and cpu time Each filter defines a box consisting of 6 numbers minimum X max X min Y max Y min Z max Z and some options There are currently 2 options exclude if the word exclude appears in the filter then only facets with all vertices within the filter box will be used facets with any one or more vertices outside the box will be excluded The default is to include any facet whose extent includes the filter box Some facets that are much larger than the volume do not have any vertices inside the volume but could pass through the filter volume hence the use of the extent test output the word output followed by a name without spaces in it will output the filtered and corrected facets to

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_docs/facetfix/facetfix.htm (2016-02-15)
    Open archived version from archive


  • TR326: Appendix A - The PHOENICS-VR Colour Palette
    Encyclopaedia Index Contents list Appendix A The PHOENICS VR Colour Palette The colour palette used in PHOENICS VR is as shown below Contents list

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_docs/tr326/appendxa.htm (2016-02-15)
    Open archived version from archive

  • Tr326: VR Viewer
    with reduced pressure range Averaged when ticked contours are drawn using interpolated and averaged values The values at the cell centres are averaged to the cell corners The corner values are then averaged to the position of the plotting plane When not ticked the raw data value at the cell centre is used to block fill the entire cell as shown below The probe then also shows the raw data Contours with Averaging turned off If the contours do not change dramatically when switching between averaging on and off it is an indication that the grid is reasonable If they do change dramatically the grid is almost certainly not fine enough to resolve the details of the flow Opaqueness By default contours are opaque Setting the opaqueness value to less than 100 makes the contours transparent Reverse colours so that red equates with the minimum value and blue maximum Update palette changes the colours used to draw the contours and colour the vectors A maximum of 24 colours bands can be defined The palette is automatically adjusted to keep blue at the low end and red at the high end Selecting four colours results in Red Yellow green and Blue contours The individual colours can then be further adjusted by clicking on the colour sample and using the colour chooser dialog to fine tune the colour The modified palette can be saved to an ASCII text file or a previously saved palette can be loaded Saved palettes can be loaded from Macros Greyscale when ticked the contours are drawn using a greyscale instead of the normal colours Contour drawn in Greyscale Vector Toggle Vector Options A left click turns the display of vectors on and off The display plane is selected by the Slice Direction buttons The variable used to colour the vectors is selected with the Select pressure temperature velocity variable buttons A right click brings up the Vector Options page of the Viewer Options dialog This dialog can be left open Any changes made are implemented in the image immediately The same dialog is opened from Settings Vector options When entering values in the data entry boxes press the tab key to signal end of input Show vectors turns the display of vectors on and off This is the same as clicking the vector toggle icon Vectors scaled by scale factor or reference velocity By default vectors are scaled relative to the largest velocity in the domain This can make it difficult to compare results from two cases where the maximum velocity is different because the vectors will be different lengths for the same speed This can be avoided by setting a user defined reference velocity and using the same value for both cases Note that whilst making the vector scale larger makes the vectors longer making the reference velocity larger makes the vectors shorter Vector size limit when ticked the value in the input box represents the longest vector Vectors representing higher velocities will be drawn as if they were this size The aim is to enable low velocity vectors to be seen without making high velocity vectors ridiculously long Unlimited vectors Limited vectors for same case In IPSA cases this option has a further effect When ticked the vector is multiplied by the volume fraction of the currently selected phase This has the effect of removing vectors from regions where the current phase is absent Vector Type switches between Total vectors using all available vector components In plane vectors which omit the velocity component normal to the current plotting plane This can be used to bring out secondary flow features which are often masked by the dominant flow direction which is much stronger Vector phase switches the vector and streamline plotting between Phases 1 and 2 for IPSA or The User option allows up to three variables to be used to construct the vectors Leaving a component set to None will omit it from the vector calculation An extra variable Vector magnitude is added to the list of variables available for plotting The Select Velocity button will select Phase 1 or Phase 2 velocity or Vector magnitude according to the vector phase setting Similarly the Select Temperature button will select Phase 1 or Phase 2 temperature Vector intervals controls the plotting of vectors When set to say 3 only every third vector will be drawn This can prevent vector plots on very fine grids from appearing like contour plots Vector appearance This provides the alternatives Coloured by displays a list of available variables The vector will be coloured by the local values of the selected variable Note This will also change the variable used for colouring the contours and iso surfaces Fixed Colour when not ticked the default the vectors are coloured by the current plotting variable unless contours are displayed In that case vectors are drawn in black if the image background is light e g white or in white if the image background is dark e g black When ticked all vectors are drawn in the chosen colour regardless of the background colour or whether contours are displayed or not The colour is selected from a colour chooser dialog Vector line width sets the width in pixels of the lines used to draw the vectors Increasing this to 2 or 3 can improve the appearance of saved images especially at reduced resolution Draw vectors as 3D arrows when not ticked the default the vectors are drawn as lines with a 2D arrow head When ticked the vectors are drawn as solid 3D bars with cones as vector heads Line arrows with line width 3 3D vectors with line width 3 Iso surface Toggle Iso surface Options A left click turns the display of iso surfaces on and off The variable used to colour the display is selected with the Select pressure temperature velocity variable buttons The iso surface value is usually the value of the selected variable at the probe location An exact value for the iso surface can also be set from the Iso surface options dialog box An iso surface will be drawn for each stored slice at the value of the current variable at the location the probe was when the slice was stored This becomes especially apparent in a BFC multi block case when several slices may be stored to build up a larger image If the iso surface value is changed from probe value to User set only a single surface will be drawn at the set value A right click brings up the Iso surface Options page of the Viewer Options dialog This dialog can be left open Any changes made are implemented in the image immediately The same dialog is opened from Settings Iso surface options Show iso surface contours turns the display of iso surfaces on and off This is the same as clicking the iso surface toggle icon Iso Surface value switches the iso surface between the value at the probe location and a user specified value Opaqueness setting the opaqueness to less than 100 will make the iso surface transparent The iso surface opaqueness is shared with the contour opaqueness changes made to one will affect the other Streamline Management Stream Options Control over the generation of streamlines is exercised through the Streamline Management Panel shown below It is a listbox similar to the Object Management Panel although the contents of the columns are different From the left they are index type X Y Z location for the origin of the streamline and visibility flag The Streamline Management Panel has three menus Object Action and Animate which are described below Object Action Animate Streamline Object Menu Menu item New will generate streamlines according to the Options currently selected The second Menu item Options will bring up a dialog where various options for streamline generation may be selected The streamline options menu may also be started by right clicking on the streamline icon on the handset The options on this dialog are Stream line mode Lines Arrows or Ribbons Ribbons can show rotation in the flow but have the potential to look messy if the two lines making up the ribbon move too far apart Stream line width This sets the width of the streamline in pixels The default is one pixel Stream line direction Upstream Downstream Both Flight time limits This sets the minimum and maximum flight time in seconds for which the track will be drawn Time zero is at the track start point Tracks drawn upstream will display negative time values Stream line coloured by Total time Flight time Current plotting variable Time along track When coloured by Track time Flight time or Total time the contour scale will show time unless Contours Vectors or Iso surfaces are turned on When coloured by Total time the time scale will reflect the minimum maximum times for all the tracks When coloured by Flight time the time scale will reflect the minimum maximum flight times When coloured by current plotting variable the variable used to colour streamlines is selected with the Select pressure temperature velocity variable buttons When coloured by Track time each track will go from blue to red and the scale will show the times for the last track drawn Stream line start From probe position Along a line Around a circle The dialogs for starting along a line and around a circle are shown below Start along a line Start around a circle When starting along a line the icon will reset the start or end point to the current probe location When starting around a circle the centre of the circle is always the current probe position and the circle lies in the current plotting plane The line mode for each streamline is stored Changing the Coloured by mode or the Flight time limits will change the way all existing streamlines are drawn Select All will highlight all the streamlines for subsequent action Refresh will repopulate the listbox and clear any current selections Close closes the management panel The line mode for each streamline is stored Changing the width Flight time or Coloured by mode will change the way all existing streamlines are drawn Streamline Action Menu Actions will apply to all currently selected streamlines When creating streamlines or groups of streamlines there are occasions when not all are desired in the final plot The first three actions therefore give the ability to turn individual streamline visibility on or off or to delete the streamline altogether The Streamline location item brings up a dialog which enables the starting location for individual streamlines to be modified The last three menu items enable the streamline type to be changed Streamline Animate Menu Once generated the streamlines can be animated to give a better impression of the flow pattern The default animation display balls radio button active is to send a single grey sphere along the length of each streamline The length of the animation will be 100 frames and will repeat until the animation is terminated The animation control dialog is shown here Increasing the number of frames makes the animation smoother but will slow down the refresh rate Increasing the number of balls per streamline gives a clearer impression of the path but also slows the animation down The best result is a compromise between the number of frames and the number of balls for the particular graphics card in question When the Show streamline checkbox is unticked only the balls or vectors will be drawn The streamline balls are usually grey When the Colour by value checkbox is ticked they will be coloured in the same way as the streamline itself The size of the balls is controlled by the Animation ball size input box When the vector radio button is active the balls are replaced by moving vectors which show the local flow direction and speed The Animation ball size input box is replaced by Vector scale Vectors look best with the streamline display turned off When the line segments radio button is active the balls are replaced by line segments denoting a set time period The default time period is 10 of the total flight time The Animation ball size input box is replaced by Segment length The Save button opens a dialog from which an animated GIF file or an AVI of the animation can be saved The file will contain all the frames chosen for the animation these can be a sub set of the frames making up the full animation For example 1000 frames may be needed for smooth movement but only 100 frames need be saved to AVI making a movie lasting 10 seconds at the default 10 frames per second Streamlines from GENTRA Tracks GENTRA particle tracks can be displayed as streamlines in the VR Viewer The track information is written to a particle history file called by default ghis To display the tracks click the Macro button on the handset and enter the name of the history file as the name of the RUN macro file Click OK to close the Macro functions dialog All the track information in the file will be read and a streamline will be generated for each track in the file If only selected tracks are required create a file with any convenient name and insert into it the lines history read ghis m n where ghis is the name of the history file m is the number of the first track to read and n if present is the number of the last track to read The file can contain any number of such pairs of lines This file can now be used as a VR Viewer macro file The new streamlines will appear in the Streamline management panel and can be operated on as normal streamlines except Cannot change start location Cannot change to ribbon In a transient case particle tracks will only be displayed for the time steps which fall within their flight time Slice Toggle This button toggles the display of the saved slices on and off It is not possible to store slices or delete slices when the slice toggle is off Slice Management The Slice Management Panel offers a limited subset of facilities described for the Streamline Management Panel described above it has no Options or Animation items This makes a copy of the contours and vectors on the slice plane When the slice plane is moved the copied image will remain until deleted In this way it is possible to build up a composite image with many planes The saved slices have an orange border to differentiate them from the current slice which has a white or black border The image below shows two X slices a Y slice and a Z slice The contour and vector display in the saved slices always follows that of the current slice Animation Toggle Animation Options A left click will cause the Viewer to re create the current image from each of the saved intermediate files if any are available The frequency with which intermediate files are saved by the solver is set from the Output Dump Settings panel of the Main Menu A right click on the Animation Toggle brings up a dialog from which the animation can be controlled The dialog controls the first and last sweeps or steps to be plotted and the sweep or step frequency of plotting Save animation saves an animated GIF or AVI of the animation The file name and type will be prompted for each time the animation button is pressed Save animation as MACRO opens the MACRO dialog in order to save the current screen set up as a macro file The saved macro will end with the ANIMATE keyword Units for time display allows the choice of seconds milliseconds minutes hours or days for the time shown in the top right corner of the screen Once an animation is started it can be stopped by clicking the Animation Toggle again or by clicking in the graphics window and pressing any key on the keyboard Slice Direction X Y Z This selects the plane to be used for displaying contours and vectors The plane always passes through the probe location so moving the probe moves the plotting plane This plane is also used when starting streamlines around a circle In BFC these buttons refer to the I J and K grid co ordinate directions and not to the Cartesian space directions As the probe moves by one cell at a time the slice also moves by a cell at a time In a multi block case the contours and vectors are drawn for the block that contains the probe If plots are required from more than one block at a time use the Slice Management button to store the current slice then move the probe to the next block create another slice and so on until the entire picture has been built up Select Pressure This selects the pressure variable P1 as the current display variable It will be used to colour contours vectors streamlines and iso surfaces Select Velocity This selects the absolute velocity as the current display variable It will be used to colour contours vectors streamlines and iso surfaces For two phase cases the velocity will be that of the current vector phase Select Temperature This selects the temperature variable TEM1 or TMP1 if it is stored as the current display variable It will be used to colour contours vectors streamlines and iso surfaces For two phase cases the temperature will be that of the current vector phase Select a Variable This brings up the Contour Options page of the Viewer Options dialog box Any variable can be selected as the current display variable Probe position hand set controls These down up buttons move the probe in the X Y and Z directions The distance moved each time a button is pressed is controlled by the Snap size set in the Editor Parameters dialog box The bigger the snap size the faster the probe moves In BFC the probe can only be moved from cell centre to cell centre The probe location is always in IX IY IZ It is possible to double click on the data entry box next to the buttons and type in an exact position The new position is read when another data entry box is selected To see which cell the probe is in right click on the Axis Toggle View Options button and turn Cell posn on The probe cell location will be displayed in the bottom right corner of the graphics screen Show Probe Location The probe is also a pickable item Double clicking on the probe or clicking the icon on the toolbar causes the following dialog to appear It contains a quick summary of the data at the current probe location The current variable can be changed by choosing from the pull down list The probe location can be moved through the model either in physical space X Y Z position or by cell location When moving in physical space the Cell location boxes will show the cell centre nearest to the probe When moving by Cell location the X Y Z boxes will show the physical location of the probe The Set view centre button centres the view on the probe Display of Minimum and Maximum Value Locations The other pages of the above dialog show the low and high spots for the current variable If the Reveal checkbox is ticked then the Low and High Spots will be revealed as blue and red spheres respectively These indicate the position of the minimum and maximum values respectively of the current plotting variable The Set view centre button centres the view on the Low or High spot and the Move probe here button moves the probe to the position of the High or Low Spot bringing with it the contours or vectors if they are displayed In order to control the size of the spheres open the Viewer options dialog by clicking on the Select variable icon then moving to the final page Options Viewer Options Options Dialog This dialog is reached by clicking the Select variable icon or from Settings Contour options then selecting the Options tab This page contains options to Show or hide the grid equivalent to button Show or hide slices equivalent to button Show or hide the minimum and maximum locations Change the size of the ball used to display the minimum and maximum locations Show or hide the probe equivalent to button Show or hide cell position equivalent to right click on Axis toggle button Show hide move axes equivalent to right click on Axis toggle button Show hide move plot title equivalent to right click on Axis toggle button Show hide move contour key equivalent to right click on Axis toggle button Boxed key when ticked the contour key is represented by a series of coloured boxes with contour bounds displayed Otherwise contours values are represented by a colour coded numerical value Use coarse contours on rotation Useful if the model is very large and rotating the domain appears sluggish Plot Variable Profile Clicking on the icon leads to the Graph Options dialog The aim is to create a line plot of a selected variable along a straight line between two points The variable to plot is first chosen from the pull down selection this may be any of the saved variables in the PHI file Next the start and end point for the line should be selected By default these will be the low and high spot for the current plotting plane The icon will reset the start or end point to the current probe location The number of intervals used for the profile is restricted to between 2 and 5000 The Show points check box toggles the display of the profile points in the main VR Viewer window as shown above The probe will be moved along the line joining the start and end points and the value of the selected variable will be saved at the set number of intervals If contour averaging is turned on the data will be taken from the averaged field If it is turned off it will be taken from the raw data and sharp steps may be seen at cell boundaries Distance is measured along the line with zero denoting the start point The Plot button will display the data as a graph using the Title and X and Y axis labels from the dialog To close the plot use the left mouse button to click on this window A hard copy of the profile plot may be obtained using the Save window as option from the File menu The print option for a profile plot is not currently available The data points are automatically saved to the file specified in the Filename input box By default it is varname profile csv where varname is the name of the selected variable This file can be easily imported to Excel It is also compatible with AUTOPLOT The Main Menu dialog The Main Menu in the Viewer is much simpler than in the Editor The Title can be temporarily reset but changes are lost on return to the Editor The Type and Domain size items are for information only they cannot be changed The Probe Position and scaling factors can be changed but are also available elsewhere The Display unit system can be used to change the units on the display between SI FPS and CGS The following three images use SI FPS and CGS units systems respectively Note that it is only the contour scale probe value and average value which are changed the data is left untouched The probe location is also shown in the selected unit set The Object Dialog Box Just as in the VR Editor double clicking on an object will display the object dialog box Unlike in the VR Editor the values in the size and position boxes cannot be changed they display the size and position of the current object Similarly the Geometry button on the Shape page is inactive although textures can still be modified The Go To and Hide buttons work as in the VR Editor The Colour button allows the colour and or transparency of an object to be temporarily changed Permanent changes can only be made in the VR Editor and become permanent when the working files are saved The Show nett sources button on the Options page has several functions depending on the object selected The possibilities are To show the sources and sinks associated with this object The values are read from the RESULT file and are displayed in a scrolling window To view the entire RESULT file click File Open file for editing RESULT When the display window is closed the contents are automatically written to a text file in the local working directory The name of the file is the same as the title of the display dialog with a txt extension If Show nett sources is selected for the Domain in the Object Management Dialog all the sources for all objects and all force and moment information will be displayed in the dialog and subsequently saved to the local text file To show the integrated pressure force on the object This only works for BLOCKAGE objects made of a solid material if Output of forces and moments is On on the Main Menu Output panel This data is also saved to the same file as the sums of sources In a transient case the Results button will plot a time history of the selected variable if a Point history object has been selected The data points are automatically saved to the file with the default name objname varname history csv where objname is the name of the object and varname is the name of the selected variable This file can be easily imported to Excel It is also compatible with AUTOPLOT The time plot can be saved to a macro by pressing F4 overwrite current or F5 append to current when the variable selection dialog is open Show nett sources can also be found on the Object Management Dialog action menu and on the right click context menu Contours and Vectors on the Surfaces of Objects To colour an object by the current plotting variable select the object right click on it and select Surface contour from the context menu Alternatively select the object or objects from the Object Management Dialog and then select Surface contours from the Action Menu To plot contours or vectors on an arbitrary surface create a Plotting Surface object This can be done in the Viewer or Editor from Settings New Object Plotting surface or from the Object management Dialog Object New Plotting Surface Double click the new object to open its Specification dialog set any desired shape and place it where required In the Viewer select the new object right click and select Surface contours Now select Surface vectors Note that surface contours and surface vectors cannot be displayed simultaneously on the same plotting surface To do so copy the surface then plot contours on one and vectors on the other VR Viewer Object context menu When an object is selected the right mouse button will bring up a context menu for the current object Most options which allow for the modification of an object have been disabled except Modify Colour However there are five items which are only available in the viewer Surface contour contour values of the current variable on the surface of an object Surface vector plot vectors on the surface of the current object Strictly the vector tails are at the centres of cells affected by the object Dump surface values if Surface contour ticked writes a file objname varname csv which contains four columns of data x y z value which were used to generate the surface contours So if the current object was named FENCE and the current variable is Pressure the filename would be FENCE Pressure csv Dump object profile if Surface contour ticked writes a file similar to the above although it will only dump the values along the current slice The filename is composed as objname varname plane planeloc csv where objname is the object name varname the current variable plane the currently select plane X Y or Z and planeloc the current position of that plane This latter will be moved to the cell centre in the file Both csv files can be easily imported into Excel for further processing They are also compatible with AUTOPLOT Show nett sources See description in section above VR Viewer Scripting Macro Facility The function of the VRV Macro A typical use of a VRV macro is to record a particular view and magnification and then restore it whenever required for example when preparing images for a report Using a macro in this way will guarantee that a sequence of images from different phi files perhaps generated on separate occasions can all use the identical view and magnification settings A macro file can contain a single image or a sequence of images separated by PAUSE commands Macro command files are ASCII text files and can be edited or created using any convenient text file editor Saving VR Viewer Macros Macro files can be created by the following means Press the MACRO button on the Hand set This will bring up the MACRO Functions dialog Select Save as new Yes enter a filename then click OK The current view settings will be written to the selected file The default filename is vrvlog When the setting Full or partial is set to Full all settings will be written for each frame This gives a long macro file but guarantees that the views will be reproduced exactly regardless of the starting condition When the setting Full or partial is set to Partial only non default settings will be written for the first frame For subsequent frames only settings which have changed from the previous fame will be written This results in a much shorter macro file but the final effect may depend on the state of the Viewer view direction object visibility etc when the macro is run Press the MACRO button on the Hand set Select Append to old Yes then click OK The current view settings will be added to the end of the selected file Individual views will be separated by a PAUSE command Pressing F4 or clicking the F4 icon will overwrite the currently selected Save as new file Pressing F5 or clicking the F5 icon will add to the currently selected Append file Individual views will be separated by a PAUSE command Hand editing a file with any text editor Running VR Viewer Macros The macro commands can be copied into the Q1 input file by placing them between the statements VRV USE and ENDUSE These two lines and all the macro lines must start in column 3 or more to ensure that they are treated as comments by the VR Editor Macros can be run by as follows Press the MACRO button on the Hand set then Run Macro Yes Click Ok and the selected file will be read and any macro commands in it will be executed The default file name is Q1 Pressing F3 or clicking the F3 icon will run the currently selected macro file VR Viewer Macro Commands The commands set the state of the VR Viewer settings The screen image is updated when a PAUSE command is encountered or the end of the file is reached The image is the outcome of the final states of all the settings The commands making up the macro language can be divided into a number of groups Only those commands that change a default or change an existing setting need be placed in the macro file The individual commands can be shortened as long as the remaining part is unique Short cuts to commands Setting the File Name Setting the View Setting the Position of Screen Items Setting the Variable to be Plotted Controlling the Plot Elements Setting the Colour Palette Line Plots Time History Plots X Cycle Options Controlling Slices Probe Position Low and High Spots General Display Toggles Controlling the Display of Objects Clipping Plane Options Saving Images Progress Controls GENTRA Particle Tracks PHOTON USE Files Extracting multiple data points Setting the File Name FILE name xyz name for BFC case Sets the name of the PHI and XYZ file to plot FILE Read the next saved PHI file equivalent to F8 FILE Read the previous saved PHI file equivalent to F7 USE file Read commands from another file Use files can be nested to a depth of 5 Setting the View VIEW x y x Sets the View direction UP x y z Sets the Up direction VIEW CENTRE x y z Sets the Cartesian co ordinates of the view centre the point about which the image rotates VIEW DEPTH depth For the Windows Viewer sets the View Depth The default is 3 0 A value of 100000 makes the view isometric VIEW TILT angle For the DOS Unix Viewer sets the perspective angle The default is 0 8 A value of 0 0 makes the view isometric SCALE scalex scaley scalez Sets the overall domain scaling factors NEARPLANE depth Sets the near plane position The View Centre Scale factors and View Depth Tilt settings can be seen by clicking on Reset Setting the Position of Screen Items The primary keyword POSITION is followed by a secondary keyword to identify which item is under consideration The two integers represent the normalised location in the range 0 0 1 0 for the first character of the item The origin is at the top left hand corner of the client area Each item also has an associated ON OFF switch to indicate whether it is to be drawn or not POSITION CELL x y Shows the current probe position in terms of cells CELPOS ON OFF POSITION CONTOURKEY x y Positions the contour colour scale CONTOUR SCALE ON OFF POSITION TITLE x y Positions the run title TEXT ON OFF POSITION PROBE x y Positions the echo of the probe value PROBE ON OFF POSITION AXIS x y Positions the location of the axes AXIS ON OFF Labeling the Plot TEXT CLEAR Deletes all previous text items TEXT size colour Sets size 1 large 4 small and colour immediately followed by string Sets the actual text to be placed on the plot Xpos Ypos Sets the normalised position of the first character as above Setting the Variable to be Plotted VARIABLE name Sets the name of the current plotting variable VARIABLE RANGE min max Sets the minimum and maximum values for the plot If the minimum maximum values are changed from the default i e the current field values the set values will be used for all subsequent plots for this variable Setting specific minimum and maximum values ensures that contour plots from different PHI files are scaled consistently Controlling the Plot Elements Contours CONTOUR ON OFF CLEAR Clear also deletes all saved slices CONTOUR FILL ON OFF Controls the display of the filled contour bands CONTOUR LINE ON OFF Controls the display of the contour lines between bands CONTOUR LINE COLOUR index MULTI Controls the colour index of the contour lines between bands MULTI indicates using contour band colours CONTOUR LINE WIDTH ncwid Controls the width in pixels of the contour lines between bands CONTOUR SCALE ON OFF Controls the display of the contour scale CONTOUR OPAQUENESS iopaq Sets the contour opaqueness to iopaq CONTOUR BLANK ON OFF Sets the out of range transparency CONTOUR AVERAGE ON OFF Sets contour averaging on or

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_docs/tr326/vr-view.htm (2016-02-15)
    Open archived version from archive

  • ENC_CVA.HTM
    full in a lecture on the mathematical basis of MFM Here it suffices to say that the equations are the same as those which are encounterd in single fluid computational fuid dynamics except that in all places the dependent variable i e the CVA is multiplied by the mass fraction of the particular fluid in question and additional terms are present which represent fluid to fluid transfer 3 General implementation in PHOENICS 3 1 How to set up a CVA solving simulation In order to be recognised as a CVA when a flow simulation involving MFM is in progress a solved for variable must be given a name which has a letter other than F as the first character because F is reserved for the mass fractions of the various fluids has digits 1 2 3 etcetera up to 99 as the next one or two characters which denote to which fluid the particular CVA applies and has the letter C as the final character For example the names H1C H2C H3C H4C and H5C might represent the enthalpies of fluids 1 2 3 4 and 5 however the choice of H as the first character does not of itself determine this Its significance as enthalpy must be given by way of the boundary conditions and source terms that are provided for each variable The source terms appropriate to CVAs are of two kinds namely what might be called conventional sources such as for enthalpy the energy gained or lost per unit mass by absorption or emission of radiation or for species concentration the nett gain of species mass in unit time of unit mass of fluid and MFM related sources which express the loss and gain of the attribute in question by reason of micro mixing The PHOENICS user is responsible for specifying the former but the latter are computed and applied automatically by PHOENICS 3 2 A library example a Description and purpose The example to be considered is L231 which concerns the generation of smoke in an idealised indeed one dimensional combustor The smoke generating reaction is also idealised its rate being supposed to depend only on the unburned fuel mass fraction and the temperature as explained here Case L231 is a modification of case L230 in which the flame was supposed to be adiabatic in case L231 by contrast heat loss is allowed The absence of heat loss from case L230 made it unnecessary to employ any CVA at all because the unburned fuel mass fraction and the temperature of each fluid remained constant therefore the rate of smoke production within each fluid could be computed consequently the total rate of production could be obtained by summing the individual rates In case L231 the heat loss reduces the temperature of each fluid and that reduction lowers also but more than proportionately the production rate of smoke Therefore at least the temperature must appear as a CVA and if it influences the heat loss rate so must

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_enc/enc_cva.htm (2016-02-15)
    Open archived version from archive

  • MICA Software Delivery: Version 1.2 (December 1997)
    Object Dialogue Box COFFUS has been extended so that the PLATE object automatically produces the required radiative and convective boundary conditions for a fixed temperature bounding wall The user must specify the temperature in o C and the wall emissivity as indicated in Figure 4 At present the only other acceptable boundary condition in COFFUS for PLATE objects is zero heat flux i e adiabatic b BLOCKAGE objects Walls of a finite thickness or wedges which are 3d objects should be constructed by use of the BLOCKAGE object Similar to the PLATE object it is possible to specify a solid blockage with a hole in it by creating the whole wall and then superimposing the hole later This is quicker than splitting the wall into several solid regions so as to define the hole it also means that the hole can be of any shape for which clip art exists e g cylindrical However note that when viewing the geometry the blockage will appear entirely solid c BURNER objects BURNER type objects are 2d and are comprised of 2 co axial cylinders The shapes cylpipe clip art type should be selected for this type of object As shown in Figures 5 and 6 the attributes panel provides access to a dialogue box from which the user can define a the location of the burner centre b the burner axis direction with toggles between ve and ve and c the inner and outer diameters of the burner object Figure 5 Burner Object Dialogue Box Primary Inlet Figure 6 Burner Object Dialogue Box Secondary Inlet The user can bring up the dialogue box for the Primary Inlet Properties by clicking on the Inner Inlet Properties panel and that for the Secondary Inlet Properties by clicking on the Outer Inlet Properties panel As shown in Figure 5 the Primary Inlet Properties are the gas and coal mass flow rates in kg s the gas composition by way of the species mass fractions of O 2 N 2 CO 2 H 2 O and NO the Swirl number as defined by the ratio of a constant tangential velocity to the axial injection velocity the swirl direction which toggles between clockwise and anti clockwise the gas and coal inlet temperatures in o C and the inlet turbulent intensity in As shown in Figure 6 the Secondary Inlet Properties are the gas mass flow rate in kg s the gas composition by way of the species mass fractions of O 2 N 2 CO 2 H 2 O and NO the Swirl number as defined by the ratio of a constant tangential velocity to the axial injection velocity the swirl direction which toggles between clockwise and anti clockwise the gas inlet temperature in o C and the inlet turbulent intensity The swirl velocity as deduced from the Swirl number and direction is ve for clockwise swirl into the solution domain when using a right handed coordinate system Note that VR currently sets the inlet turbulence intensity to 5 independent of the value specified in the dialogue box d OUTLET objects OUTLETs are 2d objects COFFUS VR currently assumes that both phases exit through an outflow boundary The dialogue box permits user specification of the pressure coefficient for the gas phase and the external pressure value relative to the furnace operating pressure which is used by both phases The pressure coefficient is usually set to unity and COFFUS VR arranges that the solid phase pressure coefficient is scaled with the solid phase density and also that in the event of inflow the incoming phasic enthalpies are equal to those in the cell For single phase and non reacting flow the default CORE VR practice is retained which requires user specification of the external temperature in o C e NULL objects In this version the computational grid is based on the number of cells in each direction specified by the user and the objects defined There is no provision for a fine grid region to be specified If the automatically calculated grid does not provide enough resolution in some areas it may be possible to remedy the defect by creating small null objects these have no effect on the problem specification but introduce additional regions into which at least one cell will be placed f INLET objects INLETs are 2d objects which can be used for rectangular or square inlets of gas and or coal At present this facility is restricted to specifying an orthogonal injection in terms of the mass flow rate and as such it does not allow the possibility of an inclined injection stream of coal and gas As shown in Figure 7 the user must define the following inlet properties the gas and coal mass flow rates in kg s the gas composition by way of the species mass fractions of O 2 N 2 CO 2 H 2 O and NO the Swirl number as defined by the ratio of a constant tangential velocity to the axial injection velocity the gas and coal inlet temperatures in o C and the inlet turbulent intensity in Figure 7 Inlet Object Dialogue Box 2 4 The use of PIL fragments It is possible for the user to introduce additional functionality by using PIL coding with which PHOENICS users have long been familiar Examples of this approach are provided for most of the MICA application sectors these can act as a guide for other users PIL fragment location There are three different types of PIL fragment a application specific COFFUS for the coal furnace b object specific for a particular application c case specific The first type is used for global settings related to an application these are settings that can be expected to be required in all such cases although they can be changed if necessary by means of case specific PIL code The second type is used for settings that are expected to be required whenever a particular type of object e g INLET is used for a particular application The third type is used for settings global or object related that are specific to the current case COFFUS VR uses PIL fragments only of the first and third types The first type is in the file COFFUS2 which is located in the directory PHOENICS D SATELL D MEN MICA COFFUS Case specific PIL coding is located at the end of the Q1 file for the case in question examples can be found in the COFFUS library Q1 files located in directory PHOENICS D SATELL D VRLIBS MICA COFFUS PIL coding of both types will be included automatically on exit from the VR Editor the settings will be included in the Q1 and EARDAT files that are created at that time If MONITR T is inserted at the end of the Q1 file then for inspection purposes only COFFUS VR will write a file Q1EAR on leaving the VR Editor this will contain the input settings in conventional PIL format Note that changes to PIL fragments of the first and second type will only influence simulations run locally for remote runs the standard PIL fragments on the remote machine will be accessed Thus changes for remote runs should always be made at the end of the Q1 file PIL fragment descriptions The PIL coding in COFFUS2 which defines the DOMAIN settings for COFFUS is as follows PIL settings for DOMAIN of COFFUS BOOLEAN BURN INERT NOXCAL RADCAL SIZECH 2 Inert 1 Combustion IF ATTRIB 31 EQ 2 THEN BURN F ELSE BURN T ENDIF SPEDAT SET COFFUS BURN L BURN 2 6 flux 1 no radiation IF ATTRIB 17 EQ 2 THEN RADCAL T ELSE RADCAL F ENDIF SPEDAT SET COFFUS RADCAL L RADCAL 1 particle size change 2 no size change IF ATTRIB 19 EQ 1 THEN SIZECH T ELSE SIZECH F ENDIF Coal characterisation INTEGER JCOAL JCOAL ATTRIB 12 IF JCOAL EQ 1 THEN Lignite ENDIF IF JCOAL EQ 2 THEN Bituminous ENDIF IF JCOAL EQ 3 THEN Anthracite ENDIF IF JCOAL EQ 4 THEN Subbituminous ENDIF IF JCOAL EQ 5 THEN Semianthracite ENDIF Coal proximate analysis mass fractions as fired YWATM water YASHM ash n b YRAWC 1 0 YWATM YASHM YCHAM char if YCHAM 0 volatiles if YCHAM 0 REAL YWATM YASHM YCHAM YRAWC YWATM 1 E 2 ATTRIB 23 YASHM 1 E 2 ATTRIB 24 YCHAM 1 E 2 ATTRIB 25 YRAWC 1 E 2 ATTRIB 26 SPEDAT SET COFFUS YWATM R YWATM SPEDAT SET COFFUS YASHM R YASHM SPEDAT SET COFFUS YCHAM R YCHAM SPEDAT SET COFFUS INCOAL R YRAWC Dry coal ultimate analysis mass fractions YCDRY carbon YODRY oxygen YSDRY sulphur YHDRY hydrogen YNDRY nitrogen REAL YHDRY YCDRY YODRY YNDRY YSDRY YHDRY 1 E 2 ATTRIB 27 YCDRY 1 E 2 ATTRIB 28 YODRY 1 E 2 ATTRIB 29 YNDRY 1 E 2 ATTRIB 30 YSDRY 1 YASHM YHDRY YCDRY YODRY YNDRY YSDRY 0 0132 SPEDAT SET COFFUS YCDRY R YCDRY SPEDAT SET COFFUS YODRY R YODRY SPEDAT SET COFFUS YSDRY R YSDRY SPEDAT SET COFFUS YHDRY R YHDRY SPEDAT SET COFFUS YNDRY R YNDRY Group 7 Variables STOREd SOLVEd NAMEd the following STORE SOLUTN coding is inactive as currently coded in s5 cof for STORE TMP1 IF NOT ONEPHS THEN SOLUTN U2 Y Y N N N Y SOLUTN V2 Y Y N N N Y SOLUTN W2 Y Y N N N Y SOLUTN R1 Y Y N N N Y SOLUTN R2 Y Y N N N Y SOLUTN H2 Y Y N N N Y SOLUTN H1 P P N P P P STORE TMP2 Gas phase mass fractions SOLVE YCHX YCO2 YCO YH2O YO2 STORE YN2 Particle phase mass fractions SOLVE COL2 CHA2 WAT2 STORE ASH2 ENDIF Particle Size Change REAL SMDIAM SMDIAM ATTRIB 21 1 E 6 IF SIZECH THEN SOLVE PHIS STORE SIZE FIINIT PHIS 1 0 TERMS PHIS P P P P N P FIINIT SIZE SMDIAM ENDIF Thermal Radiation REAL ABSORB SCAT INTEGER IHRADL IF RADCAL THEN SPEDAT SET COFFUS RADCAL L RADCAL ABSORB RADIA SCAT RADIB RADIAT FLUX ABSORB SCAT H1 SOLUTN H1 P P N P P P PATCH RADISO FREEVL 1 NX 1 NY 1 NZ 1 LSTEP if IHRADL 0 then set CP1 in MJ kg for use in GXRADI IHRADL 1 SPEDAT SET RADI IHRADL I IHRADL relaxation REAL DTRAD DTRAD 0 1 RELAX RADX FALSDT DTRAD RELAX RADY FALSDT DTRAD RELAX RADZ FALSDT DTRAD reset RADZ to slabwise to aid convergence JRH 230797 SOLUTN RADZ P P N P P P ENDIF Group 8 Terms Devices TERMS H1 N P P P P P TERMS H2 N P P P P P Gas phase mass fractions assign variables to phases TERMS YO2 P P P P Y P TERMS YCHX P P P P Y N TERMS YCO P P P P Y P TERMS YCO2 P P P P Y P TERMS YH2O P P P P Y P Particle phase mass fractions assign variables to phases TERMS COL2 P P P P N P TERMS CHA2 P P P P N P TERMS WAT2 P P P P N P CSG3 CNGR Group 9 Properties REAL RG56 RG57 TMP1 GRND2 TMP1A TEMP0 TMP1B 1 CP1 ENUL 1 8E 5 RHO2 2590 TMP1 GRND2 TMP1A TEMP0 TMP1B 1 CP1 TMP2 GRND2 CP2 1 8E3 TMP2A TEMP0 TMP2B 1 CP2 PRNDTL H2 1 E10 RHO1 GRND TMP1 GRND RG56 GRND RG57 2 0 SPEDAT SET COFFUS RG56 R RG56 SPEDAT SET COFFUS RG57 R RG57 TMP1B 1000 TMP2 GRND2 CP2 1 8 1 E3 TMP2B 1 CP2 TMP1A 0 0 TMP2A 0 0 TEMP0 0 0 Enthalpies solved in MJ kg HUNIT 1 E 6 CP1 CP1 HUNIT PRNDTL WAT2 1 E10 PRNDTL COL2 1 E10 PRNDTL CHA2 1 E10 Group 10 Inter Phase Transfer Processes INTEGER ISLIP ISLIP ATTRIB 18 STORE CFIP SPEDAT SET COFFUS SMDIAM R SMDIAM IF ISLIP EQ 1 THEN No slip ISLIP 1 for now high CFIPS eventually use EQUVEL T CFIPS 1 E6 ELSE Slip islip 2 RLOLIM 1 E 5 CFIPS GRND STORE SLIP REYN CINT YCHX 0 0 CINT YO2 0 0 CINT YH2O 0 0 CINT YCO 0 0 CINT YCO2 0 0 CINT WAT2 0 0 CINT CHA2 0 0 CINT COL2 0 0 Gas particles heat transfer CINT H1 0 0 CINT H2 0 0 PATCH LHEATRA FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL LHEATRA H1 GRND3 GRND3 COVAL LHEATRA H2 GRND3 GRND3 storage for relaxation of heat source STORE QDOT FIINIT QDOT 0 0 special relaxation factors for heat source RELAX QDOT LINRLX 1 0 RESREF QDOT 0 3 ENDIT QDOT 0 3 PRT QDOT 1 0 PRNDTL QDOT 1 0 REAL THCON THCON 0 0458 SPEDAT SET COFFUS GASCON R THCON ENDIF Group 11 Initialise Var Porosity Fields REAL HCALC TREFE TINI TCALC FIINIT R2 1 E 5 FIINIT R1 1 0 FIINIT R2 FIINIT YCHX 0 FIINIT YCO 0 IF BURN THEN FIINIT YO2 0 0 FIINIT YH2O 0 12 FIINIT YN2 0 73 FIINIT YCO2 0 15 ELSE FIINIT YO2 0 232 FIINIT YH2O 0 FIINIT YN2 0 768 FIINIT YCO2 0 ENDIF FIINIT ASH2 1 0 FIINIT CHA2 0 0 FIINIT WAT2 0 0 FIINIT COL2 0 0 TREFE 273 0 TINI 340 IF BURN THEN TCALC TREFE 600 0 ELSE TCALC TREFE TINI ENDIF compute the enthalpy of the air stream RG57 2 0 HCALC 1 0E 3 0 77 0 97035 1 493E 4 TCALC RG57 TCALC HCALC HCALC 1 0E 3 0 23 1 0802 3 265E 5 TCALC RG57 TCALC FIINIT H1 HCALC FIINIT H2 CP2 TCALC FIINIT TEMP1 TCALC FIINIT TEMP2 TCALC FIINIT DEN1 PRESS0 FIINIT TEMP1 287 398 GROUP 13 Boundary conditions and special sources COAL DRYING PATCH VAPORIS0 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL VAPORIS0 P1 FIXFLU GRND5 COVAL VAPORIS0 P2 FIXFLU GRND5 COVAL VAPORIS0 YH2O ONLYMS 1 0 COVAL VAPORIS0 H1 ONLYMS GRND5 PATCH VAPORIS2 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL VAPORIS2 H2 FIXFLU GRND5 COVAL VAPORIS2 WAT2 FIXFLU GRND5 Patch to counter the transfer of species due to mass transfer PATCH VAPORIS5 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL VAPORIS5 COL2 GRND9 GRND9 COVAL VAPORIS5 CHA2 GRND9 GRND9 COVAL VAPORIS5 WAT2 GRND9 GRND9 IF SIZECH THEN COVAL VAPORIS5 PHIS GRND9 GRND9 ENDIF STORE VAPO FIINIT VAPO 0 0 RELAX VAPO LINRLX 1 0 RESREF VAPO 0 3 ENDIT VAPO 0 3 PRT VAPO 1 0 PRNDTL VAPO 1 0 REAL ADEVOL EDEVOL RG55 C1EBU C2EBU INTEGER IG14 IG16 MODHET IORDER IKDMEA RAW COAL VOLATILISATION Raw coal Y Volatiles 1 Y Char IF BURN THEN PATCH DEVOLAT0 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL DEVOLAT0 P1 FIXFLU GRND1 COVAL DEVOLAT0 P2 FIXFLU GRND1 COVAL DEVOLAT0 YCHX ONLYMS GRND1 COVAL DEVOLAT0 YO2 ONLYMS GRND1 COVAL DEVOLAT0 H1 ONLYMS GRND1 Patch for volatile phase 2 species PATCH DEVOLAT2 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL DEVOLAT2 COL2 FIXFLU 0 0 COVAL DEVOLAT2 CHA2 FIXFLU GRND1 COVAL DEVOLAT2 COL2 GRND1 0 0 ENDIF ADEVOL constant A EDEVOL constant E R in volat model ADEVOL 2000 0 EDEVOL 2 3E4 8 130 SPEDAT SET COFFUS ADEVOL R ADEVOL SPEDAT SET COFFUS EDEVOL R EDEVOL store volatilization rate special relaxation STORE VRAT FIINIT VRAT 0 0 RELAX VRAT LINRLX 1 0 RESREF VRAT 0 3 ENDIT VRAT 0 3 PRT VRAT 1 0 PRNDTL VRAT 1 0 Use this line to de activate devolat skip rg 21 0 0 Patch to counter the transfer of non volatile species IF BURN THEN PATCH DEVOLAT5 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL DEVOLAT5 COL2 GRND9 GRND9 COVAL DEVOLAT5 CHA2 GRND9 GRND9 COVAL DEVOLAT5 WAT2 GRND9 GRND9 IF SIZECH THEN COVAL DEVOLAT5 PHIS GRND9 GRND9 ENDIF ENDIF Ask litec what this is suggest better name than IG16 IG16 1 indicates CO as product 2 indicates CO2 IG16 1 SPEDAT SET COFFUS IG16 I IG16 Combustion Ask litec what these are CHx O2 CO2 H20 Not active patch combustc phasem 1 nx 1 ny 1 nz 1 1000 coval combustc ychx fixflu grnd2 coval combustc yo2 fixflu grnd2 coval combustc yh2o fixflu grnd2 coval combustc yco2 fixflu grnd2 ig16 2 coval combustc h1 fixflu grnd2 TWO STEP HOMOGENEOUS COMBUSTION OF VOLATILES YCHX Ask litec what these are Store combustion rate for relaxation STORE COM1 FIINIT COM1 0 0 LITER COM1 1 RELAX COM1 LINRLX 1 0 RESREF COM1 0 3 ENDIT COM1 0 3 PRT COM1 1 0 PRNDTL COM1 1 0 CHx O2 CO H20 seguida por CO O2 CO2 IF BURN THEN PATCH COMBUSTA PHASEM 1 NX 1 NY 1 NZ 1 LSTEP COVAL COMBUSTA YCHX FIXFLU GRND2 COVAL COMBUSTA YO2 FIXFLU GRND2 COVAL COMBUSTA YH2O FIXFLU GRND2 COVAL COMBUSTA YCO FIXFLU GRND2 COVAL COMBUSTA H1 FIXFLU GRND2 PATCH COMBUSTB PHASEM 1 NX 1 NY 1 NZ 1 LSTEP COVAL COMBUSTB YO2 FIXFLU GRND2 COVAL COMBUSTB YCO2 FIXFLU GRND2 COVAL COMBUSTB YCO FIXFLU GRND2 COVAL COMBUSTB H1 FIXFLU GRND2 ENDIF Ask litec what these are Store combustion rate for relaxation and products for printout STORE COM2 FIINIT COM2 0 0 RELAX COM2 LINRLX 1 0 RESREF COM2 0 3 ENDIT COM2 0 3 PRT COM2 1 0 PRNDTL COM2 1 0 Liter 1 selects mixing rate ep k as combustion rate LITER COM2 1 C1EBU C2EBU EBU constants in combustion model C1EBU 4 0 C2EBU 0 0 SPEDAT SET COFFUS C1EBU R C1EBU SPEDAT SET COFFUS C2EBU R C2EBU HETEROGENEOUS COMBUSTION OF CHAR C S O2 CO2 IF BURN THEN PATCH BURNOUT0 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL BURNOUT0 P1 FIXFLU GRND5 COVAL BURNOUT0 P2 FIXFLU GRND5 COVAL BURNOUT0 H1 ONLYMS GRND1 PATCH BURNOUT2 FREEVL 1 NX 1 NY 1 NZ 1 LSTEP COVAL BURNOUT2 H1 FIXFLU GRND5 COVAL BURNOUT2 YO2 FIXFLU GRND5 COVAL BURNOUT2 CHA2 FIXFLU GRND5 COVAL BURNOUT2 YCO GRND9 1 0 Ask litec what this is

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_chron/mica/coffus.htm (2016-02-15)
    Open archived version from archive

  • MICA models for CFD
    project description Click below for abridged versions of some of the Final Reports on the MICA project by CHAM Final overview UK Building Research Establishment Fire Research Station UK Building Research Establishment flows around buildings Marine flows the Swedish Meteorological

    Original URL path: http://www.cham.co.uk/phoenics/d_polis/d_chron/mica/mica.htm (2016-02-15)
    Open archived version from archive



  •  


web-archive-uk.com, 2017-12-18