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".
  • What's New in PHOENICS 3.3
    undertaken between mid December 1999 and mid February 2000 to make it again possible for users to read understand and profit from it Particular features of this work have been All coding connected with an individual material property has been collected into a unique file For example all density related coding now appears in GXDENS htm The various inter phase and other source terms have been similarly treated so there is now a GXBUOY HTM for buoyancy sources and a GXIFRIC HTM file for inter phase friction The two phase related files have been removed from the core to D TWOPHS subdirectory The present organization and state of the relevant open source coding can be seen by clicking here The physical significance of the coding has not changed significantly but the improved order and clarity of the open source sector of PHOENICS EARTH has made it possible to resume the enrichment of its contents 6 Input file library improvements For many years PHOENICS has been supplied with a large library of tried and tested Q1 files which could be used as starting points for new simulations for demonstrations and for instructional puprposes However because of file handling limitations which prevailed at the time at which the library was instituted all Q1s were collected into a few very large files and individual cases were made selectable by an ingenious but machine dependent coding sequence in the Satellite With the passage of time the increase of size of the library files became a serious impediment to editing and upgrading with the result that their content remained static and increasingly out of date PHOENICS 3 3 however has a new library system which assigns a distinct file to each case so facilitating access editing loading etc The current contents of the library system can be seen by clicking here The way is now clear for bringing into the library system a suitable selection from the large number of interesting Q1s which have been created by CHAM and by PHOENICS users generally Offers of cases for inclusion will be gratefully accepted and of course the files will carry the names of their originators 7 MUSES the MUltiply SharEd Space procedure The MUSES feature which makes extensive use of PLANT has been further developed during the course of 1999 It is described in the Encyclopaedia article and extensively exploited in the creation of the SAFIR blast furnace special purpose program MUSES is likely to be especially welcomed by PHOENICS user concerned with complex chemical engineering equipment 8 SPP improvements and novelties HOTBOX improvements FLAIR improvements SAFIR 8 HOTBOX improvements HOTBOX the version of PHOENICS which is used for the simulation of electronics cooling problems has been extensively improved in respect of its Virtual Reality user interface its physical modelling capability especially in respect of radiative heat transfer and HOTBOX is available for trial by way of the Remote Computing Service FLAIR improvements What has been said about HOTBOX can be said also about FLAIR the

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

  • Front-end developments for PHOENICS 3.2
    that there exist two SATELLITE executables namely SATEXE EXE and PVREXE EXE SATEXE EXE acts in the same way as the SATELLITE of previous versions and should be used whenever virtual reality features are not required It resides in the sub directory PHOENICS D SATELL D WINDF It is activated at the command line in response to the long established command RUNSAT PVREXE latter resides in the sub directory PHOENICS D MSDEV It is activated at the command line in response to the new command RUNPVR It embodies the SATELLITE which should be used whenever the virtual reality editor or viewer are to be employed and it also plays the part of an environment from which all PHOENICS modules can be run including the just mentioned SATEXE In this role it does everything which the PHOENICS Commander of earlier versions was intended to do and of course it does so in a Windows like manner In respect of the other modules of PHOENICS namely the solver EARTH and the graphics display package PHOTON incorporating the once separate program AUTOPLOT there exists only one of each and they are activated at the command line by the long established commands RUNEAR and RUNPHO 2 Which modules to use 2 1 Which modules to use for VR and non VR When the virtual reality editor or the virtual reality results viewer is required the command to issue is RUNPVR This will however work only when the Q1 file has TALK T in its top line and it is unsuitable for cases in which body fitted coordinates are to be used RUNPVR may also be used when the virtual reality features are not required simply in order that the environmental facilities ie file handling and program execution can be employed For non VR and non environment needs RUNSAT is preferable whether TALK T or TALK F appears in the top line of the Q1 and whether or not body fitted coordinates are required 2 2 Which modules to use for TALK F and TALK T As implied above only RUNSAT should be used for TALK F cases for in these the only function of the the SATELLITE module is to read and interpret the Q1 file and to create the corresponding EARDAT and if VR type objects are contained in the Q1 file FACETDAT files It is often desirable to create a Q1 in TALK T mode using the virtual reality editor to create the geometry and to set the main attributes Thereafter the top line can be changed to TALK F and further changes to the Q1 perhaps for subsequent runs can be made by means of a text 2 3 Which module to use for BFC cases As implied above only RUNSAT should be used for BFC T cases whether TALK T or TALK F The reason is that at present there provision in the Virtual Reality interface for creating BFC grids 2 3 Which module to use for PLANT If PLANT

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

  • The Virtual-Reality user interface: PHOENICS-VR
    as employed for solids The apertures may be of the bursting panel character and indeed the flow rate through them can be any non linear function of pressure difference and time Fan characteristics are introduced in this way They may also consist of sources of gas or liquid which become active after the start of the explosion as a consequence of the passage of the flame or of the pressure wave which precedes it 3 5 Initial conditions Initial condition objects are introduced in the same fashion They may represent the distributions of fuel air ratio within the module or they may define the time location and strength of the ignition source The shapes of such initial condition objects may be box like ellipsoidal cylindrical or wedge shaped and more complicated shapes may be created by combinations of these It should also be stated that the inventory of shapes which are currently supplied can be augmented by any others for which there is sufficient demand to justify the not very large cost of preparing them 3 6 Other input data At data input time users may if they wish indicate their desire for particular pieces of information to be supplied at the end of the not yet started flow simulating computation This is done by introducing output objects into the world This is especially advantageous when either the user has already made previous computations of the same kind and requires only to learn how a particular value for example peak pressure at a particular point will be influenced by a change of input data or the user is availing himself of the remote computing option and desires not to pay for the transmission of excessive amounts of information 4 Viewing results by way of VR 4 1 Probing the domain of study The PHOENICS interface is used for viewing and exploring the results of the simulation as well as for setting up the problem at the start Entry is by way of the VR Viewer option of the PHOENICS Commander choice of which presents the user with a picture of the world which he created but this now has a different set of control buttons Important features to note in the pictures are the probe and its positioning the ability to switch the horizon off and on the ability to represent objects in wire frame i e transparent form and the abilities to change the view point and to zoom in and out in the same manner as at data input time The VR viewer Note the contour scale on the left The horizon can be removed Objects can be viewed as wire frames Streamlines Velocity vectors on planes Temperature contours on a plane Temperature contours on another plane and another and another Successive planes can be plotted and animated Equal temperature surfaces 5 How EARTH responds to VR When the user is satisfied that he has sufficiently defined his world he simply issues a calculate the flow in this command and leaves the rest to PHOENICS The rest means translating the world definition into PHOENICS terms i e writing a Q1 file determining whether the flow is likely to be laminar or turbulent and if the latter selecting a turbulence model creating a computational grid which will represent the geometrical aspects of the world with adequate precision and selecting the solution algorithm and the numerical control parameters in such a manner as will lead to a quick and accurate solution Once the PHOENICS Earth program has performed the computations and dumped the results into a file the contents of this file are processed by a PHOENICS to VR converter The consequence is that as has been explained in section 4 when the user looks at his world again he finds that it is enriched with information concerning temperatures pressures fluid compositions and velocities at all points Of course it is not essential for the user to examine all this information or even any of it Indeed the conventional RESULT file may already contain the two or three pieces on information in which the user is interested at the time for example the answer to the question Is the maximum permissible temperature exceeded anywhere 6 How VR fits into the PHOENICS system When PHOENICS is installed a directory PHOENICS is created and all of the files for PHOENICS are placed into sub directories according to the file type and to which PHOENICS program s it applies The directory structure for files relating to PHOENICS VR is PHOENICS D PHOEVR D WINDF D SATELL D WINDF D SPP D GEOM D LIBS The following panels describe this structure in more detail Directories for the VR Commander D PHOEVR This directory contains files for the VR Commander the executable the help files and run scripts D PHOEVR D WIN This contains the equivalent files for the Windows version of the VR Commander Directories for the VR Editor and Results Viewer D SATELL The VR Editor and Results Viewer are both parts of the SATELLITE program which is in this directory D SATELL D WIN This contains the equivalent files for the Windows version of SATELLITE D SATELL D SPP D GEOM The VR image ie clip art files are stored here divided into sub directories according to the object type D SATELL D SPP D LIBS VR library cases are stored here divided into sub directories according to the SPP name 7 How to run VR 7 1 Using the VR Commander The VR Commander gives the user complete control over the VR Editor EARTH and the VR Results Viewer and allows users to load use and save their cases There are several ways to start the VR Commander namely From the PHOENICS Commander click on the PHOENICS VR button From the DOS prompt type runpvr and press ENTER From the Windows version Desktop if installed double click on the PHOENICS VR icon on the Desktop Once the VR Commander has been

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

  • 2 phase turbulent flow in bubble column reactors RHONE POULENC and turbulent gaseous combustion for the design of furnaces and large utility boilers RADIAN CORPORATION STATOIL Several projects were performed at 50 discount to the Client on the condition that any new features were coded as generally available options in PHOENICS Spin offs from such projects included the 2 phase IPSA enhancements for interfacial momentum transfer and turbulence modelling the P 1 thermal radiation model and the extended SCRS gaseous combustion facility COLLABORATIVE PROJECTS The collaborative ESPRIT project concerned with the creation of PHOENICS CVD was recently completed The partners were ASM CHAM Fraunhofer SIEMANS and TU Delft PHOENICS CVD is an industry specific product which simulates Chemical Vapour Deposition CVD processes which involve laminar flow of multicomponent gas with diffusion processes due to concentration and thermal gradients gas phase and surface chemical reactions convective diffusive and radiative heat transfer and conjugate heat transfer between solids and gases and plasma enhanced chemical reactions PHOENICS CVD can now be purchased as a stand alone package or as an add on to an existing PHOENICS installation CHAM in collaboration with COLT VR Ltd has developed a virtual reality interface to PHOENICS There is a generic PHOENICS VR interface and there are special purpose versions designed for specific industries The VR interface enables users to set up flow simulation tasks for PHOENICS in real life terms placing walls fans apertures etc where they want them and seeing them appear on the screen in a life like manner The first release of PHOENICS VR is now available for PCs or as an Internet linked service During 1995 the European ESPRIT III and PASHA projects concerned with the creation of Parallel PHOENICS were completed Parallel PHOENICS has recently been benchmarked to provide end users with information on its performance and several users have taken delivery of a version with basic PHOENICS functionality DEVELOPMENT ACTIVITIES PHOENICS Development activities which will be covered in more detail in a separate lecture have included the provision of several new features and improvements to many existing ones New features include PHOENICS for WINDOWS 3 1 the extended SCRS gaseous combustion facility the P 1 and Rosseland thermal radiation models turbulence modelling enhancements to IPSA Arbitrary Source or Solid Allocation Procedure ASAP the extension of the CCM multi block option to include IPSA and a sliding mesh feature and an extensive set of higher order convection discretisation schemes for all single phase variables on staggered meshes including BFCs Improvements have been made in a number of areas including the solution algorithm in respect of procuring rapid convergence at very low Reynolds numbers e g capillary flows unification of the PRPS and VPOR treatment for representing completely blocked cells the numerical stability of the k e equations KELIN 3 and the staggered mesh BFC option in that 2d swirling flow may now be simulated with NX 1 Improvements have also been made and continue to be made to the CHEMKIN interface the

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

  • Developments in PHOENICS during 1994
    discretization methods for convection and diffusion 2 5 Parallelization progress 2 1 EARTH generated wall functions Earlier versions of PHOENICS when used for conjugate heat transfer problems generated a burdensome conglomeration of patches These are rendered unnecessary by the the EGWF feature which calculates all the flow resistances at solid fluid interfaces internally EGWF T in the Q1 does it all No more CONPORs See Encylopaedia item EGWF A core library example illustrating its use is case 290 2 2 Extended multi blocking and fine grid embedding PHOENICS 1 6 already possessed grid restructuring capabilities which could be used for economically fitting awkwardly shaped domains and for embedding fine grids within coarser ones See for example the active demonstration series These have been refined and improved for PHOENICS 2 0 but much greater advances are being supplied with PHOENICS 2 1 to be issued at the end of April Specifically whereas the earlier re structuring method involved the use of over lapping cells the new method allows for edge to edge contact Further instead of the long familiar staggered grid and solution for velocity resolutes at cell boundaries the solution for co located velocity resolutes The third departure from previous practice which has rendered the multi block method efficient is reliance on the conjugate gradient solver instead of the long standard Stone type solver This change has been brought about with very little internal change to the logic and storage arrangements of PHOENICS which still thinks it is a structured grid code Some simple passive demonstrations can be viewed as option k of hlp Some multi block grids can be seen by entering Q then F then d grids mb1 d grids mb2 d grids mb3 d grids mb4 2 3 Colocated variables extension 2 4 Alternative discretization methods for convection and diffusion 2 5 Parallelization progress CHAM is participating in EEC supported Projects PASHA EuroPort and SHIPS and it has also been collaborating with Southampton University s Parallel Application Centre As a consequence PHOENICS 2 has now been successfully ported to the Distributed Array Processor of Cambridge Parallel Computing formerly call Active Memory Technology several multi transputer computers manufactured by Parsytech GmbH clusters of IBM RS 6000 and other work stations The ports to 1 and 2 make use of domain decomposition techniques which have been found to give excellent efficiency 3 Utilities and interfaces 3 1 The CHEMKIN interface 3 2 The IGES interface 3 3 New grid generation features 3 4 The Q1 interpreter 3 5 HLP 3 6 The phi file stripper 3 1 The CHEMKIN interface A new interface has been written by JK Worrell of CHAM which connects PHOENICS with the public domain code CHEMKIN created by the Sandia Laboratories It is a PHOENICS add on and is supplied with a special library of Q1 s See About PHOENICS 2 items in description of modules how to run and FORTRAN files 3 2 The IGES interface See About PHOENICS 2 items in description of

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

  • develop versions which can detect which variable has most effect on convergence and to optimise its factors 3 to handle several variables at the same time 4 to extend EXPERT to the making of initial set up choices including those of domain decomposition for parallelization CHAM 7 developments Incorporation of other solvers 1992 1993 16 Incorporation means replacing EARTH by a different equation solving program having special features for example completely open Fortran source code or adaptation to simulation of a particular equipment type or familiarity to an already existing user comunity The arrangement is Satellite solver PHOTON CHAM 8 developments Incorporation of other solvers 1992 1993 16 continued CHAM has already incorporated several other codes in the system including CORA CHAM s pre PHOENICS combustion chamber code LEARN an instructional code from Erlangen University KIVA 1 and KIVA 2 public domain codes for IC engines VVV a series of codes from CHAM MEI Moscow Solvers to be incorporated in 1993 4 include A structured grid code simulating external flow around aircraft and missiles An unstructured grid code used for aircraft design but adaptable also to other external flow problems Subject to not yet completed negotiations a hybrid structured plus unstructured code CHAM 9 developments Interfacing to other packages 1992 1993 16 Interfaces already exist to PATRAN but not yet PDA s latest issue IDEAS FEMGEN FEMVIEW AVS Data visualiser Wavefront VVG CHAM MEI Interfaces to be created during 1993 include CFView Free University Brussels Fieldview Intelligent Light Concept Aries CHAM s policy is to provide or procure interfaces to all packages favoured by PHOENICS users Suggestions are welcome CHAM 10 developments Non staggered grid option 1992 1993 16 Background Some popular codes employ grids for which velocities are stored at cell centres solving equations with variants of the Rhie Chou scheme The cartesian direction velocities are solved for even for polar or BFC grids However all known schemes produce non physical solutions in some circumstances eg ripples in velocity profiles chequer board pressure fields Recent progress CHAM MEI has reviewed all schemes and invented a new one free from the above defects Work for 1993 The new scheme will be tested then attached to PHOENICS CHAM 11 developments Chemical kinetics upgrade 1992 1993 16 Background CHAM possesses and still uses in house the CREK package created by Pratt and Wormeck in the 1970s An interface has been created by Sher and Kravchuk in Israel to the public domain CHEMKIN package The public domain KIVA 2 program has a significant chemical kinetic capability Development tasks for 1993 4 The capabilities strengths and weaknesses of these packages will be studied during the course of funded projects By adaptation incorporation or new development a flexible and user friendly chemical kinetic add on module for PHOENICS will be created CHAM 12 developments Surface to surface radiation 1992 1993 16 Background PHOENICS has the built in six flux model of radiation which is intended for absorbing and emitting media This model can be used

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

  • Early versions of PHOENICS
    come to a new computer system the editor commands of which are unfamiliar to them Movement and searching through the attached file can then be effected by use of regular GUIDE facilities ie n n string string etc Enter H when viewing a file to learn what these facilities are 4 6 2 New lecture material The lecture material concerned with PHOENICS and with CFD generally has been updated and extended A lecture concerned with combustion fundamentals has been added 4 6 3 The CHANGE facility The user is now permitted to change if he wishes the number of lines printed to the screen by entering C when the top menu is on the screen Some preliminary help information is offered at the same time 4 6 4 User s own menu The addition of the user s own menu of items has been made easier by the provision of a skeleton menu which the user can adapt to his purposes simply by editing the item titles and contents in the GDY1 file 4 7 Compatibility between versions 1 4 and 1 5 Contents 4 7 1 Q1 files 4 7 2 GROUND files 4 7 3 EARDAT files 4 7 4 PHI and PHIDA files 4 7 5 XYZ and XYZDA files 4 7 1 Q1 files All Q1 files written for version 1 4 will be read and correctly interpreted by version 1 5 unless they infringe some of the restrictions which the advanced PHOENICS Input Language features entail for example colons in TEXT arguments or the use of commas rather than semi colons as delimiters The reverse is not necessarily true because Q1 files written for version 1 5 may make use of advanced PIL features 4 7 2 GROUND files The principles governing the construction of GROUND coding and its interaction with EARTH are the same for version 1 5 as for version 1 4 GREX3 FTN is the version 1 5 counterpart of the GREX2 FTN issued with version 1 4 It contains everything which its predecessor had but is somewhat richer in content tidier in form and more economical in execution 4 7 3 EARDAT files All EARDAT files written for version 1 4 can be read and correctly interpreted by version 1 5 The reverse is not necessarily true because EARDAT files for version 1 5 may contain data settings which are without significance for version 1 4 There are however very few of these Version 1 5 2 January 1990 has introduced some economies in the reading and writing of EARDAT files which are now shorter than for PHOENICS version 1 4 If for reasons of compatibility users wish to adhere to the 1 4 format they should change the logical variable EAR14 from FALSE to TRUE in SATLIT FTN re compile and relink This feature renders it especially desirable to set NPHI in the SATELLITE to the lowest value consistent with the number of variables which require to be stored 4 7 4 PHI and PHIDA files All PHI and PHIDA files written for version 1 4 can be read and correctly interpreted by version 1 5 The reverse is also true 4 7 5 XYZ and XYZDA files All XYZ and XYZDA files written for version 1 4 can be read and correctly interpreted by version 1 5 The reverse is also true 4 8 Add ons to version 1 5 Contents 4 8 1 Computation speed enhancers 4 8 2 Mathematical formulation variants Introduction The capabilities of PHOENICS can be extended in so many ways that it is impracticable to supply all available add ons with the standard version Instead add ons are being prepared for inclusion when appropriate arrangements are made in a separate sub directory of the PHOENICS system In most cases an add on consists of a specific GROUND subroutine a library of Q1 files and some documentation The add ons which are currently available or can quickly be brought into an issuable state are listed below 4 8 1 Computation speed enhancers 4 8 1 1 The chain solver This add on is a set of Fortran subroutines which can be called from GROUND for the purpose of accelerating the solution of the pressure correction or other scalar equations when the flow domain is labyrinth like ie when it consists of circuitous interconnected passages formed by the blockage of many of the faces of most cells It is useful because the solver built into earth like most solvers which have no place for information about the connectedness of the domain takes a long time to converge The chain solver has been found to converge 50 times as fast as the standard solver in some test problems The add on has been used to enable PHOENICS to represent a flow domain consisting of a three dimensional nuclear reactor vessel coupled to a series of pipes conveying steam and or water to steam generators pressurisers and other vessels The solver was devised by D B Spalding and implemented by M A Serag Eldin in 1987 It was further worked on by A Palacio during 1988 and 1989 4 8 1 2 The cyclic recursion TDMA replacement This add on is a Fortran subroutine which can be supplied to replace the built in tri diagonal matrix solver The latter employs recursive DO loops and so is not vectorisable The add on by contrast is fully vectorisable 4 8 1 3 Multi grid solvers PHOENICS version 1 5 has a user adaptable multi grid solver which has been described above This is more effective that conventional multi grid solvers when blocks of high conductivity material are arbitrarily distributed in a low conductivity fluid Conventional multi grid solvers can be supplied as an add on for acceleration of the convergence of equation sets of which the coefficients vary more smoothly through space 4 8 1 4 Solvers for parallel computers Linear equation solvers are under development by CHAM and its associates which exploit the possibility of using several processors in parallel They are being made available to PHOENICS users as add ons 4 8 2 Mathematical formulation variants 4 8 2 1 Equation formulation schemes This add on consists of a series of FORTRAN subroutines which can be called from GROUND and activated from SATELLITE They have the effect of changing the formulation of the convection and diffusion terms in the finite domain equations to accord with a variety of alternative prescriptions thus replacing the default prescription namely fully implicit upwind Twelve schemes are provided including fully explicit Crank Nicholson and Courant number dependent in respect of time and QUICK and QUICKEST in respect of space This add on has been incorporated in versions 1 6 1003 onwards 4 8 3 1 Particle tracking schemes This add on handles two phase flows by means of a variant of the particle source in cell method of CT Crowe This was first applied to PHOENICS by A Castrejon It allows for full coupling between the particulate and continuous phases which exchange mass momentum and energy This add on has been built into PHOENICS from versions 1 5 4 onwards with the name of GENTRA 4 8 3 2 The scalar equation procedure for free surface tracking This add on is based upon the Imperial College PhD work of Liu Jun It represents a simple and economical method of handling transient surface movement phenomena in two or three dimensions This add on has been incorporated in versions 1 6 1003 onwards 4 8 3 3 The Lagrangian formulation for a distorting medium This add on is based on the work of DB Spalding and HQ Qin It has been applied to the distortion of a solid under impact and in particular to that of a reactive solid which may proceed to detonate 4 8 3 4 The calculation of aerodynamic lift This add on is the work of JP Edwards It assists the calculation of the lift and drag of an airfoil the flow around which has been computed on a body fitted coordinate grid 4 9 Code configuration of version 1 5 Contents 4 9 1 The directory structure 4 9 2 The PREFIX and CONFIG files 4 9 3 Installation procedures 4 9 4 Testing procedures 4 9 5 Benchmarking procedures 4 9 1 The directory structure The PHOENICS 1 4 principle of distinguishing public and private elements of the directory structure has been retained but the names of directories have been modified and the pattern of distribution of files within them has been changed so as to accommodate the present and future enlargements of the PHOENICS family of programs and so as to facilitate delivery and installation The directory structure adopted for this installation can be inspected by returning to the top menu and entering 5 followed by 1 All the directories which do not have priv in their names are to be regarded as being public which means that users have read only access to them 4 9 2 The PREFIX and CONFIG files Each program now has its own CONFIG file residing in its own public directory PREFIX files are individual to users and reside in the private directories d priv1 d priv2 etc 4 9 3 Installation procedures The detailed information on Installing PHOENICS from source code is supplied in TR108 Installation of PHOENICS from binary code on various computer is fully described in TR110 series 4 9 4 Testing procedures PHOENICS has been thoroughly tested before release However the installer should perform the test on each program to confirm that the installation has been successful For the version 1 5 a test directory d test will be created during the installation from source code The installer should set the working directory to d test The testing procedure for each program is described in TR108 For the installation from binary code the checks are performed in the private directory d priv1 through running library case 523 The procedure is described in TR110 4 9 5 Benchmarking procedures The benchmark should be run after installation in order to establish the performance profile of PHOENICS on the computer in question The installer should set the working directory to d phoe15 d test A Q1 file containing eleven test cases from the PHOENICS library is provided These cases can be run from d test by simply typing runben 5 PHOENICS version 1 6 SATELLITE and data input features 5 1 Details of the data input features GSET commands The PROPS file Reading properties from the Q1 file Interpreting the Q1 file library case 20 GREX3 reads Q1 data directly February 3 1991 USTEER reads top of Q1 March 4 1991 general menu improvements April 1991 grid generation menu created and callable the general menu objects introduced into general menu Jan Feb March 1992 Reynolds Stress Model of turbulence activated via Q1 mouse click grid generation method accessed via VIEW PHOENICS data input features introduced during 1990 GSET This new command greatly facilitates grid generation See the help entry and also library cases 510 to 523 527 529 and 536 The PROPS file Material properties set equal to GRND10 are picked up from a file called PRPOS residing in d earth in accordance with the value of the material marker property PRPS Library cases 459 to 467 921 to 950 and many others illustrate its use Reading properties from the Q1 file Properties may also be inserted in the Q1 file Library case 240 illustrates this February 3 1991 GREX3 reads Q1 data directly On this date the new subroutine GXRDQ1 was introduced It was called by GREX3 from Group 1 if LSG3 and USEGRX were set T in SATELLITE to enable EARTH to accept new input data by reading the top of the Q1 file Subsequently LSG3 was replaced by the PIL variable RDQ1 and full information about how to use the feature can now be found in the SATELLITE HELP file by entering RDQ1 The feature not only allows data to be supplied to EARTH without the necessity to run SATELLITE but it also since the extra lines supplied in Q1 are also printed on the RESULT file provides a convenient device for annotating the latter file with descriptive text This feature must be used with caution because some resettings of data may conflict with settings made in the satellite For example whole field solution cannot necessarily be switched ON because necessary storage may not be provided but it may be switched OFF by dividing the relevant ISLN by 5 March 4 1991 USTEER reads top of Q1 The USTEER feature has been extended so that on the occasions requested by the user EARTH re reads the top of Q1 entries which of course the user may have altered by editing since it was last read This greatly increases the number of data items which the user can change during execution New SATELLITE features in version 1 6 2 PHOENICS version 1 6 2 succeeded PHOENICS 1 5 4 and was released in April 1991 This section of GUIDE summarises the main changes General menu improvements The PHOENICS general menu was enhanced in several respects the speed of response was improved in most cases the menu panels were better formulated body fitted coordinates and parabolic flows were made accessible by the menu the choice of turbulence models was increased to three effective viscosity k l and k e each with a simple low Reynolds number option a new grid generation menu see below was created which could be called from the general menu Grid Generation Menu The grid generation and grid handling procedures of PHOENICS underwent major re development The main new features were The possibility of specifying the geometry first and of subordinating the grid to it This facilitates the specification of boundary and internal features and allows the refinement of the grid in specific regions of the computational domain without disturbing the geometry A new grid generation menu facilitating the process of creating Cartesian cylindrical polar and body fitted curvilinear grids Grid planes in body fitted coordinates can be copied translated and rotated to generate new planes Transitions between two different section shapes eg square to circle can be readily generated A new generation of more powerful PIL commands the RSET and GSET suites that simplify the specification and handling of geometry and grids Interpreting Q1 files Library case 20 was supplied for the purpose of enabling Q1 files to be translated into plain language The translation includes information on the grid variables solved variables stored but not solved patches viscosity and solver parameters It is grouped so that different items of information for each variable are kept together Activation is effected by typing NOWIPE T LOAD 20 The user will be asked whether to display the information on the screen or append it to the end of the COPYQ1 file The latter will cause satellite to save the COPYQ1 file not the stack in Q1 and then terminate hence all PIL logic will be lost e g DO loops so this should be used with caution If the translation was displayed on the screen then the user should type LOAD CLEAR to remove the PIL commands that make up library case 20 from the stack once the translation has finished The file which may be inspected by entering SEELIB 20 is written in PIL with much use of its logical capabilities It is supplied as an example which users may wish to augment or replace in accordance with their preferences New SATELLITE features introduced with version 1 6 1003 All three types of grid handled by the menus The general Menu now provides access to the GridMenu for setting up all three grid types available in PHOENICS Cartesian polar and BFC and this is the recommended route for all grid specification The earlier procedure within the general Menu for specifying Cartesian and polar grids region by region is nevertheless still available for users who prefer that method Objects for Cartesian and cylindrical polar grids The main improvements to GridMenu concern the facility for defining objects in Cartesian and cylindrical polar grids Previously defining an object within the solution domain was merely a device for dividing the domain into regions Now however this feature has been extended in two ways the location of boundary conditions subsequently specified through the General Menu can be linked directly to named objects previously defined via the GridMenu the boundaries of an object are stored so that once objects have been defined they may be copied moved renamed or deleted singly or in groups Region boundaries will then be repositioned automatically when the option for matching the grid to the geometry is used again Mouse driven grid generation in GridMenu The VIEW facility for body fitted coordinates has been transformed into a mouse driven grid generation facility which allows the user to drop points and construct lines and frames using a mouse or the arrow keys when a mouse is not available Grids can be matched to frames without leaving the new system and the grid check facility previously available only through the GRDCHK command can be used to inspect the quality of the mesh Grid check uses a colour code to indicate the orthogonality of the grid at each point Settings effected in the VIEW environment are transferred automatically to the GridMenu session and recorded as PIL commands The MENSAV facilities will also record and replay VIEW sessions Conjugate heat transfer in the general menu A new menu option now enables users to simulate conjugate heat transfer problems solving for temperature directly instead of enthalpy As part of this feature a library of material properties can be accessed from the menu to assign properties to fluid and solid parts of the domain users themselves can add materials to this library Appropriate solution control devices viz the block correction feature solving whole field for temperature and harmonic averaging of conductivities are activated automatically Improved specification of boundary conditions in general menu In Cartesian and cylindrical polar coordinates the location of boundary features inlets outlets blockages etc can now be linked to named objects defined during the grid generation procedure This obviates the need to enter the coordinates twice once when defining the grid and again when specifying boundary conditions as was the case in previous versions If an object is subsequently repositioned or re sized then the boundary condition is also changed automatically If an object is deleted any associated boundary conditions will also be deleted without further instructions from the user If a new object is created by copying an existing one the boundary conditions are not automatically copied but a new boundary condition may be linked to the new object Similarly heat sources can be attached to existing block or plate type boundary conditions without re entering their location Advanced user options Several new options have been introduced so as to make the Menu more attractive to experienced PHOENICS users including facilities to assign values to the PIL variables RSG1 30 ISG1 20 LSG1 10 CSG1 10 and the arrays RG IG LG and CG for transmitting data to GROUND specify boundary conditions directly as patch type affected variable coefficient and value enter any PIL command directly which will then be appended to the Q1 file Other convenience improvements The following improvements have also been carried out data specifying the location of and conditions at boundary features are input directly to two menu display panels rather than via interrogation of the user about each data item in turn the scalar dependent variables defaulted to concentration of species C1 C2 etc can be renamed by the user more flexible options for printout and spot value monitoring are available Menu and help panel improvements For all menu options the hot key has been identified by an upper case character in the option name As far as possible this is the first character where this is not possible because two or more options on the panel share the same initial character this first character is lower case and whichever character is chosen as the hot key is upper case The help files for the general Menu and GridMenu have been overhauled to make the information in both more easily understood and in many cases more comprehensive The display panels for both GridMenu and the general Menu have been rationalised so that they all conform to a standard layout for example the option Return to previous panel is chosen by the same hot key on all panels New and modified help entries The following entries in the SATELLITE help dictionary are new or have been modified The entries marked have a corresponding new feature description in this section of GUIDE Entry Notes ABORT Help absent in previous version DELAY Help absent in previous version FREE New entry on free surface flows HOL New entry PROPS File listing updated REYNOLD New entry ROTA Errors corrected SCHEME New entry SEM New entry TURBUL Modified entry GSET Modified GSET V command VIEW Instructions for new version of VIEW included U1AD Modified entry TURMOD Modified entry IURVAL Modified entry SOLIDS There is no longer a need for a SOLIDS patch TSTSWP Entry modified to include graphical monitoring INTRPT New entry MENSAV Modified entry for pausing and adding to MENSAV RSET New RSET G option Changes to FORTRAN files The maximum dimensions for the grid generation arrays are now set in PARAMETER statements in the file SATLIT PHOENICS features introduced during 1992 Items introduced during the first three months SATELLITE The new menu The features added by the new menu system are direct solution for temperature defining objects through grid generation objects can be copied moved or deleted specifying boundaries linked with objects more options such as assigning Real Integer or Logicals specifying source terms through patches help entries extensions the Reynolds Stress Model of turbulence can be activated through Q1 the mouse click grid generation method is now accessed via VIEW mensav files for grid generation have been updated grid generation examples have been included in section 12 5 of the Input Library Flow simulation features in version 1 6 6 arranged alphabetically CASE 20 of the input file library interpreting the Q1 file This library file has been supplied in file lib6 for the purpose of enabling Q1 files to be translated into plain language The translation includes information on the grid variables solved variables stored but not solved patches viscosity and solver parameters It is grouped so that different items of information for each variable are kept together Activation is effected by typing NOWIPE T LOAD 20 The user will be asked whether to display the information on the screen or append it to the end of the COPYQ1 file The latter will cause satellite to save the COPYQ1 file not the stack in Q1 and then terminate hence all PIL logic will be lost e g DO loops so this should be used with caution If the translation was displayed on the screen then the user should type LOAD CLEAR to remove the PIL commands that make up library case 20 from the stack once the translation has finished The file which may be inspected by entering SEELIB 20 is written in PIL with much use of its logical capabilities It is supplied as an example which users may wish to augment or replace in accordance with their preferences CLDA A new GXCLDA FTN subroutine was added It performs the operations associated with the conservative low dispersion algorithm more efficiently than when this algorithm is implemented by way of the neighbour cell technique CONJUGATE heat transfer The conjugate heat transfer capability of PHOENICS 1 5 has been improved in several respects a An ASCII file houses now the library of materials which can therefore be easily customised expanded by the user Material properties can now be constant or dependent on other properties They can also be read from Q1 directly by EARTH b Boundary conditions including thermal links at the fluid solid interface are set up automatically c Heat transfer coefficients can be optionally calculated and printed out by Earth CONTACT resistances The Group 12 feature allowing diffusion coefficients to be modified patchwise has been extended to allow the addition of resistances to heat transfer brought about by the inclusion of thin sub grid scale sheets of poorly conducting material CORNER coordinates The following statement in GREX3 arranges for the print out of corner coordinates when STORE XCEN YCEN ZCEN appears in the Q1 file PRNCEN LBNAME XCEN NE 0 OR LBNAME YCEN NE 0 OR 1 LBNAME ZCEN NE 0 CORIOLIS forces for horizontal flows The variable CORIOL when set to a value other than zero causes momentum sources per unit mass which are equal to CORIOL V1 for U1 CORIOL U1 for V1 and to corresponding expressions for the second phase velocities This facility is useful when atmospheric and hydrospheric simulations are being performed An example can be found in the two phase section of the active demo set The DENPCO feature This new logical introduces densities into the calculation of pressure correction coefficients in order to promote convergence when strong density variations are present EXPERT system A set of auto pilot devices have been introduced which make in flight adjustments to the numerical parameters such as relaxation factors in order to speed up the convergence of the solution procedure FNxxxsubroutines FNDSDX Y SCAL d scal dx FNDSDY Y SCAL d scal dy FNDSDZ Y SCAL d scal dz FNIFDV Y X A LOGIC conditioned divide FNIFML Y X A LOGIC conditioned multiply FNMAX Y X maximum of X and Y FNVSLP Y A Slip velocity calculation FN71 FN72 FULLY DEVeloped flows It is now possible to solve for velocities without solving for pressures or indeed activating convection terms This is useful when it is desired to compute a fully developed pipe plane walled channel or Couette flow situation without as was formerly necessary simulating the flow in a long duct either elliptically or parabolically The solution is then one dimensional It is necessary to prescribe a longitudinal pressure gradient in order to create finite longitudinal velocities Library case 951 957 illustrates the use of this feature If the duct wall is neither circular nor an infinite plane the flow will be two dimensional Then pressure must be solved for and cross duct convection must be taken into account Otherwise the situation is as above An example will be found in the one phase section of the active demo series GROUP 12 features The following patch names in which stands for any character cause the terms indicated for the variable in question in the region occupied by the patch to be multiplied by the third CO argument of the corresponding COVAL GP12CON all convection terms GP12SOR all built in sources GP12CNE the east face convections GP12CNN the north face convections GP12CNH the high face convections GP12DFE the east face diffusions GP12DFN the north face diffusions GP12DFH the high face diffusions The fourth VAL argument has no significance assigned to it so far These facilities are especially useful when PHOENICS is to be used for fluid flow analysis in one part of a domain and for stresses in solids analysis in another Examples can be found in the active demo series under the heading stress analysis in solids The INTERPHASE mass transfer and friction FN98 and FN99 have been re written and made directly accessible to users They are to be found in file gxsor ftn IPARAB 4 A new hyperbolic option IPARAB 4 has been added for the efficient solution of wholly supersonic flows This feature allows the lateral pressure variations to influence also the longitudinal velocity when PARAB T in regions where the flow is supersonic It can be developed for mixed supersonic and subsonic flows and will be when specific needs arise Library case 156 exemplifies its use KINETIC heating sources for TEM1 TEM2 If TERMS TEM1 Y or TERMS TEM2 Y the kinetic heating source terms become active in the same way as for enthalpy Warning This may not always be desirable for example when stress analysis is being simulated and the velocity variables have the significance of displacements LINEAR links Highly conducting sub grid scale links can now be introduced by special PATCH and COVAL commands They are useful for representing for example thin metal plates embedded in a poorly conducting medium The LINK feature See special topics grid restructuring LOW REynolds Number turbulence A simple low Reynolds number feature is provided in GXENUT It is activated by setting PRNDTL 1 and PRT 1 which have no other significance to non default values The effective viscosity formula used is vist vist nom amin1 1 0 A vist nom visl B where vist turbulent kinematic viscosity to be used vist nom nominal high Reynolds number value of vist visl laminar viscosity A PRNDTL 1 B PRT 1 The term vist nom visl represents the local Reynolds number of turbulence Suitable values for A and B may be 0 01 and 4 0 MONITORING of spot values etc A graphical display of spot values and residuals has been provided for use on personal computers and work stations It was first activated by setting TSTSWP 12345 Later TSTSWP n has been caused to elicit an updating of these values only every n sweeps Where the PHOENICS graphical driver allows an interruption the EARTH run can be stopped to alter the relaxation factors or to end the run At the end of the EARTH run the monitoring graphs can be plotted on any hard copy device present in the PHOENICS graphics system NEIGHBOUR patch improvement The capabilities of the subroutine GXNEPA n file GXMODS FTN have been extended so that neighbours can be displaced in both INDVAR space AND time or distance The use of the new features is illustrated in a Q1 file which tests the conservative low dispersion algorithm It appears in the library as case number 491 NFUSER The user is now allowed to reserve for himself a part of the F array In the satellite he must set NFUSER to the number of locations which he needs Then in GROUND coding he can refer to the I th of his reserved locations as F L0USER I PLOTting of residuals and corrections When a solved for variable is given the name xxx and the statement STORE xxxR appears in the Q1 file the field of xxxR contains the values of the residuals of variable xxx before entry to the linear equation solver These values can be printed in the result file or viewed via PHOTON or AUTOPLOT in the same was as other stored variables When a solved for variable is given the name xxx and the statement STORE xxxC appears in the Q1 file the field of xxxR contains the values of the corrections to variable xxx after exit from the linear equation solver These values can be printed in the result file or viewed via PHOTON or AUTOPLOT in the same was as other stored variables Note this works at present only for nz 1 The PROPS file Material properties set equal to GRND10 are picked up from a file called PRPOS residing in d earth in accordance with the value of the material marker property PRPS Library cases 459 to 467 921 to 950 and many others illustrate its use P2 the second phase pressure STORE P2 has always introduced storage for and print out of the second phase pressure but no use has been made of it Now coding within EARTH ensures that if P2 is stored the pressure acting on the second phase in the momentum balance is P1 P2 so P2 is the EXTRA pressure experienced by phase 2 It is the responsibility of the user to provide coding which ascribes values to P2 but an example has been supplied in GREX3 group 19 where P2 is set proportional to R2 This is useful for the simulation of then floating layers for example oil slicks An example can be found in the two phase section of the active demo set PROPERTIES read from the Q1 file Properties may be inserted in the Q1 file Library case 240 illustrates this SETBFC and MOVBFC Setting SETBFC T causes GREX3 to call a new subroutine GXBFGR which exemplifies the use of GROUND coding to create a body fitted coordinate grid Setting MOVBFC T further causes GXBFGR to be called at each time step in order that the grid can be changed and if storage is provided for CONI CONJ and CONK GREX3 will calculate the grid movement contributions to the convection fluxes Library case 966 illustrates the use of this feature SLIP velocity feature for two phase flow Storing the variable SLPU SLPV and or SLPW now ensures that the differences between first and second phase velocities are computed and stored at the end of iterations on a slab They can be used for the plotting of slip velocity vectors in PHOTON The relevant sequence in GREX3 is SLIPVL NOT ONEPHS AND LBNAME SLPU NE 0 OR 1 LBNAME SLPV NE 0 OR LBNAME SLPW NE 0 Examples can be found in the two phase flow section of the active demo set SOLIDS It is no longer necessary to create a SOLIDS patch in order to activate the SOLIDS feature because EARTH test for the values of PRPs which have been set See SATELLITE HELP STAR NAME PATCHES If a PATCH name begins name where name is the four character name of a solved or stored variable a corresponding COVAL for variable PHI will introduce a source of PHI equal to CO PHI NAME VAL where CO is the third argument of the COVAL and NAME stands for the local value of the variable having the name which has been referred to Since such power law patches have hitherto had to be introduced via GROUND coding it is thought that their easy creation via the Q1 file will often be found convenient If the referenced variable has a 2 character name the patch name must include the blanks because the internal to EARTH test is made on the first 5 characters of the PATCH name Thus PATCH P1 A followed by COVAL P1 A H1 1 E3 2 0 will create a source of first phase enthalpy which is proportional to pressure squared but PATCH P1A followed by COVAL P1A H1 1 E3 2 0 will not An example is SOLVE VOSQ KE REAL CONST1 CONST1 4 0 PATCH VOSQ1 1 NX 1 1 1 1 1 1 COVAL VOSQ1 KE CONST1 0 5 This provides a source equal to 4 0 KE VOSQ 0 5 The feature is illustrated by library cases 413 414 and 415 STOPJOB If during execution the user creates and saves a file called STOPJOB execution stops after completion of final sweep print out The relevant coding is in GREX3 Group 19 This device may not work on all computer systems SUPPRESSing the making of pressure correction The PIL setting RELAX P1 LINRLX 0 0 now prevents pressure corrections from being added to pressure which therefore cannot change The full velocity adjustments are still made however TURBULENCE energy generation rates printed This is now effected when a variable named GEN1 is stored and the default print out provisions are used TWO EQUATION turbulence models augmented Provision has been made for both the Saffman Spalding k W model and the Kolmogorov model k omega model VIST and LEN1 formulae for these models have also been provided and also wall function entries Testing is in progress UNDER RELAXation for pressure correction The PIL setting RELAX P2 LINRLX factor where factor is a real number multiplies the continuity errors by factor before they enter the

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

  • Science Service Centre, Moscow
    CHAM Home Page Please choose the font Cyrillic Alphabet Windows before clicking below Science Service Centre Moscow

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

web-archive-uk.com, 2016-10-26