Flow123D - Technical University of Liberec
Contents
1. Existence of elements with the same elm number Reference to non existing node Reference to non existing material see below Difference between number of nodes and actual number of lines in nodes section Difference between number of elements and actual number of lines in elements section e By default 1 2 3 FLOW uses meshes with number of tags 2 tag1 is number of region in which the element lies tag2 is number of material reference to MTR file in the element 9 e Currently line type 1 triangle type 2 and tetrahedron type 4 are the only supported types of elements Existence of an element of different type stops the calculation e Wherever possible we use the file extension MSH It is not required but highly re comended 10 Material properties file format version 1 0 The file is divided in two sections header and data The extension MTR is highly recomended for files of this type MaterialFormat 1 0 file type data size EndMaterialFormat Materials number of materials material number material type lt material type specific data gt text EndMaterials Storativity material number lt storativity coefficient gt text EndStorativity Geometry material number geometry type lt geometry type specific coefficient gt text EndGeometry Sorption material number substance id sorption type lt sorption type specific data gt text EndSorption S2. of neighbouring and their specific data type 10 Edge with common nodes i e sides of elements with common nodes Pos sible many elements type 11 Edge with specified sides i e sides of the edge are explicitly defined Possible many elements type 20 Compatible i e volume of an element with a side of another element Only two elements type 30 Non compatible i e volume od an element with volume of another element Only two elements type type specific data Description 10 n_elm eid1 eid2 number of elements and their ids 11 n_sid eid1 sid1 eid2 sid2 number of sides their elements and local ids 20 eid1 eid2 sid2 coef Elm 1 has to have lower dimension 30 eid1 eid2 coef Elm 1 has to have lower dimension coef is of the double type other variables are ints 17 Comments concerning 1 2 3 FLOW e Every inconsistency or error in the NGH file causes stopping the calculation These are especially Multiple usage of the same nezghbour number Difference between number of neighbours and actual number of data lines Reference to nonexisting element Nonsence number of side e The variables sid must be nonegative and lower than the number of sides of the particular element 18 Sources file format version 1 0 The file is divided in two sections header and data The extension SRC is highly recomended for fi
3. type of water in our jargon text char arbitrary text description of the fracture notes etc up to 256 chars This is an optional parameter 15 Types of boundary conditions and their data type 1 Boundary condition of the Dirichlet s type type 2 Boundary condition of the Neumann s type type 3 Boundary condition of the Newton s type type type specific data Description 1 scalar Prescribed value of pressure or piez head 2 flux Prescribed value of flux through the boundary 3 scalar sigma Scalar value and the o coefficient scalar flux and sigma are of the double type Ways of defining the place for the boundary condition where where 1 Condition on a node 2 Condition on a generalized side where 3 Condition on side for element with only one external side where lt where data gt Description 1 node id Node id number according to MSH file 2 elm id sid id Elm id number local number of side 3 elm id Elm id number The variables node id elm id sid id are of the int type Comments concerning 1 2 3 FLOW We assume homegemous Neumman s condition as the default one Therefore we do not need to prescribe conditions on the whole boundary If the condition is given on the inner edge it is treated as an error and stops calculation Any inconsistence in the file stops calculation Bad num
4. 1 node3s time 1 comp 1 node1 time2 comp 1 node2 time2 comp 1 node3 time2 comp2 nodel timel comp2 node2 time1 comp2 nodes time1 comp2 nodel time2 comp2 node2 time2 comp2 nodes time2 text2d is a list of 4 double precision numbers coord1 coord style index where coord and coord2 give the coordinates of the leftmost element of the 2D string in screen coordinates index gives the starting index of the string in tert2d chars and style is currently unused comp3 node1 time1 comp3 node2 time1 comp3 node3 time1 comp3 node1 time2 comp3 node2 time2 comp3 node3s time2 21 tert2d chars is a list of nb text2d chars characters Substrings are separated with the character which is a forbidden character in regular strings textd is a list of 5 double precision numbers coord1 coord2 coord3 style index where coord1 coord2 and coord3 give the coordinates of the leftmost element of the 3D string in model real world coordinates index gives the starting index of the string in tert3d chars and style is currently unused A text3d chars is a list of nb tert3d chars chars Substrings are separated with the character More information about GMSH can be found at its homepage http www geuz org gmsh Comments concerning FFLOW20 e FFLOW20 generates POS file with four views Elements pressure edges pressure in terelement fluxes and complex view First three views shows raw data results ob tained by t
5. Technical university of Liberec Faculty of mechatronics informatics and interdisciplinary studies Flow123D Numerical simulation software for flow and solute transport problems in combination of fracture network and continuum Documentation of file formats and brief user manual O Severyn M Hokr J Kralovcova J Brezina J Kopal M Tauchman Liberec 20 11 2008 Flow123D Flow123D is simulating software based on Borland C Builder 6 0 It enables to solve the task of underground water flow in heterogenous rock solute transport and their inter action with rock Considered interaction with rock are non equilibrium mobile immobile pore exchange and non linear adsorption with independent parameters in each zone mo bile immobile and each area fracture continuum rock The flow is based on mixed hybrid FEM The supported task of flow are steady state flow unsteady state flow and variable density flow Calculation is supported on compatible or incompatible multidimenzional meshes Solute transport is solved with the operator splitting Convection is solved with the FVM Mobile immobile pore exchange is solved with using analytic solution and non linear adsorption is solved numerically Principle for calculation are files of mesh msh material mtr neighbours ngh boundary conditions of flow bcd eventually are needed files of boundary conditions of transport tbc initial conditions of transport tic or initial condition
6. _porosity is turned on in the ini file 12 material type material type specific data Description 11 k K k 11 a A K a k 0 2 fk Ke a kz 0 a hier aes K a Kg dki gi a A K gt 6 7 f 0 22 Ar Qy A K 0 A 23 Ar Qy Gin Ae A zy Qy k 0 0 31 k K 0k0 00k ks 0 0 33 ks ky kz K 0 k 0 0 0 k ky kry kez 36 ka ky kz key kez kys K key ky be kzz kyz kz a00 31 a A K 0a0 i 0 a a 0 0 33 r Qy az A K 0 a 0 0 0 a Ay Ary Qgrz 36 Qi Qy Oe zy Arz Tye ASK E Cig Mai Gays Qrz Qyz Qz Note all variables k kv ky kz kay Kez kyz Q Qx Qy Qz Ary Asz Qyz are of the double type sorption type sorption type specific data Description 1 kp 1 s kpe 2 kell MAA all s kre 3 1 mazjr 3 1 Kris c 3 Kyi hee MH sL MY s ETE Note all variables kp kr a Kg 8 are of the double type text char is a text description of the material up to 256 chars optional 13 This parameter is Comments concerning 1 2 3 FLOW e If number of materials differs from actual number of material lines in the file it stops the calculation 14 Boundary conditions file format version 1 0 The file is divided in two sections header and data BoundaryFormat 1 0 file type data size EndBoundaryFormat BoundaryConditions number of conditions condition number t
7. ber of conditions multiple definition of condition reference to non existing node etc At least one of the conditions has to be of the Dirichlet s or Newton s type This is well known fact from the theory of the PDE s Local numbers of sides for where 2 must be lower than the number of sides of the particular element and greater then or equal to zero The element specified for where 3 must have only one external side otherwise the program stops 16 Neighbouring file format version 1 0 The file is divided in two sections header and data The extension NGH is highly recomended for files of this type NeighbourFormat 1 0 file type data size EndNeighbourFormat Neighbours number of neighbours neighbour number type lt type specific data gt EndNeighbours where file type int is equal 0 for the ASCII file format data size int the size of the floating point numbers used in the file Usually data size sizeof double number of neighbours int Number of neighbouring defined in the file neighbour number int is the number index of the n th neighbouring These numbers do not have to be given in a consecutive or even an ordered way Each number has to be given only onece multiple definition are treated as inconsistency of the file and cause stopping the calculation type int is type of the neighbouring lt type specific data gt format of this list depends on the type Types
8. ces for the density doubles flow calculation Section Constants KEY TYPE DEFAULT DESCRIPTION g double 1 0 Gravity acceleration rho double 1 0 Density of fluid Section Run KEY TYPE DEFAULT DESCRIPTION Log_file string mixhyb log Name of log file Screen_verbosity int 8 Amount of messages printed on the screen 0 no messages 7 all messages Log_verbosity int 8 Amount of messages printed to the log file 0 no messages 7 all messages Pause_after_run YES NO NO If set to YES the program waits for a key press before it finishes Section Solver KEY TYPE DEFAULT DESCRIPTION Use_last_solution YES NO NO If set to YES uses last known solution for chosen solver Solver_name string matlab Command for calling external solver Supported solvers are petsc isol and matlab Solver_params string NULL Optional parameters for the external solver passed on the command line or PETSc options if the PETSC solver is chosen see doc petsc_help Keep_solver_files YES NO NO If set to YES files for solver are not deleted after the run of the solver Manual_solver_run YES NO NO If set to YES programm stops after writing input files for solver and lets user to run it Use_control_file YES NO NO If set to YES programm do not create control file for solver it uses gi
9. ement Line 2 nodes Triangle 3 nodes Quadrangle 4 nodes Tetrahedron 4 nodes Hexahedron 8 nodes Prism 6 nodes Pyramid 5 nodes Second order line 3 nodes Second order triangle 6 nodes Second order tetrahedron 10 nodes Point 1 node OON DOK WD re o number of tags gives the number of tags for the n th element By default Gmsh generates meshes with two tags and reads files with an arbitrary number of tags see below tag is an integer tag associated with the n th element By default the first tag is the number of the physical entity to which the element belongs the second is the number of the elementary geometrical entity to which the element belongs the third is the number of a mesh partition to which the element belongs node number list is the list of the node numbers of the n th element separated by white space without commas The ordering of the nodes is given in Gmsh node ordering for second order elements the first order nodes are given first followed by the nodes associated with the edges followed by the nodes associated with the faces if any The ordering of these additional nodes follows the ordering of the edges faces given in Gmsh node ordering More information about GMSH can be found at its homepage http www geuz org gmsh Comments concerning 1 2 3 FLOW e Every inconsistency of the file stops the calculation These are Existence of nodes with the same node number
10. files Now only the value 1 GMSH like files is accepted Mesh string NULL Name of file containig definition of the mesh for the problem Material string NULL Name of file with hydraulical properties of the elements Boundary string NULL Name of file with boundary condition data Neighbouring string NULL Name of file describing topology of the mesh Sources string NULL Name of file with definition of fluid sources This is optional file if this key is not defined calculation goes on without sources Section Transport KEY TYPE DEFAULT DESCRIPTION Transport_on YES NO NO If set YES program compute transport too Sorption YES NO NO If set YES program include sorption too Dual_porosity YES NO NO If set YES program include dual porosity too Concentration string NULL Name of file with initial concentration Transport_BCD string NULL Name of file with boundary condition for transport Transport_out string NULL Name of transport output file Transport_out_im string NULL Name of transport immobile output file Transport_out_sorp string NULL Name of transport sorbed output file Transport_out_im_sorp string NULL Name of transport sorbed immobile output file N_substances int 1 Number of substances Subst_names string undefined Names of the substances separated by commas Substances_density_scales list of 1 0 Scales of substan
11. gt 20 lt vector pyramid values gt lt tensor pyramid values gt lt text2d gt lt text2d chars gt lt text8d gt lt text3d chars gt EndView where file type is an integer equal to 0 in the ASCII file format data size is an integer equal to the size of the floating point numbers used in the file usually data size sizeof double view name is a string containing the name of the view max 256 characters nb time steps is an integer giving the number of time steps in the view nb scalar points nb vector points points in the view are integers giving the number of scalar points vector nb tert2d nb text3d are integers giving the number of 2D and 3D text strings in the view nb tert2d chars nb text3d chars are integers giving the total number of characters in the 2D and 3D strings time step values is a list of nb time steps double precision numbers giving the value of the time or any other variable for which an evolution was saved scalar point value vector point value are lists of double precision numbers giving the node coordinates and the values associated with the nodes of the nb scalar points scalar points nb vector points vector points for each of the time step values For example vector triangle value is defined as coord1 node1 coordi1 node2 coord1 node3 coord2 node1 coord2 node2 coord2 node3 coord3 nodel coord3 node2 coord3 node3 comp 1 nodel time1 comp 1 node2 time 1 comp
12. he solver without any interpolations smoothing etc The fourth view contains data processed in this way Elements pressure Contains only scalar triangle values Triangles are the same as the elements of the original mesh We prescribe constant value of the pressure on the element as it was calculated by the solver as the unknown p Therefore the three values on every triangle are the same Edge pressure Contains only scalar line values The lines are the same as the edges of the elements of the original mesh We prescribe constant value of the pressure on the edge as it was calculated by the solver as the unknown A Therefore the two values on every edge are the same Interelement flux Contains vector point values and scalar triangle values The scalar triangle values carry no information all values are set to 0 these are in the file only to define a shape of the elements The points for the vector point values are midpoints of the sides of the elements The vectors are calculated as un where u is value of the flux calculated by the solver and n is normalized vector of outer normal of the element s side Complex view Contains scalar triangle values and vector point values The scalar triangle values shows the shape of the pressure field The triangles are the the same as the elements of the original mesh Values of pressure in nodes are interpolated from ps and As The vector point values shows the velocity of the flow in the cent
13. ial point xe ye ze end point delta cylinder radius Specify_elm_type YES NO NO If set to YES next param specify type of prefered elements If set to NO each element is included Output_elm_type int 1 Spefify type of element dimension 1 1D line 2 2D triangle 3 3D tetrahedron BTC_elms list of ints undefined List of the breakthrough curve elements ints this concentrations are written to seperate file with extension btc FCs_params double 4 zero Params of flow cross section x y z 1 plane of cut general equation output values are written by coordinate of axis x 0 y 1 z 2 Pos_format string ASCII Format of the POS output file ASCII BIN opening a binary file in the GMSH is much faster Description Options controling output file of the programm Section Density KEY TYPE DEFAULT DESCRIPTION Density_implicit YES NO NO NO explicit iteration simple flow update YES implicit iteration more accurate flow update Density_max_iter int 20 Maximum number of iterations for implicit density calcultation Eps_iter double le 5 Stopping criterium for iterations maximum norm of pressure difference Write_iterations YES NO NO Write conc values during iterations to POS file Mesh file format version 2 0 The mesh file format comes from the GMSH system Following text is copied from the GMSH documentation Version 2 0 of the MSH file format
14. is Gmsh s new native mesh file format It is very similar to the old one Version 1 0 but is more general it contains information about itself and allows to associate an arbitrary number of integer tags with each element The MSH file format version 2 0 is divided in three sections defining the file format MeshFormat EndMeshFormat the nodes Nodes EndNodes and the elements Elements EndElements in the mesh MeshFormat 2 0 file type data size EndMeshFormat Nodes number of nodes node number x coord y coord z coord EndNodes Elements number of elements elm number elm type number of tags lt tags gt node number list EndElements where file type is an integer equal to 0 in the ASCII file format data size is an integer equal to the size of the floating point numbers used in the file usually data size sizeof double number of nodes is the number of nodes in the mesh node number is the number index of the n th node in the mesh Note that the node numbers do not have to be given in a consecutive or even an ordered way x coord y coord z coord are the floating point values giving the X Y and Z coordinates of the n th node number of elements is the number of elements in the mesh elm number is the number index of the n th element in the mesh Note that the elm numbers do not have to be given in a consecutive or even an ordered way elm type defines the geometrical type of the n th el
15. les of this type SourceFormat 1 0 file type data size EndSourceFormat Sources number of sources source number type eid density EndSources where file type int is equal 0 for the ASCII file format data size int the size of the floating point numbers used in the file Usually data size sizeof double number of sources int Number of sources defined in the file source number int is the number index of the n th source These numbers do not have to be given in a consecutive or even an ordered way Each number has to be given only onece multiple definition are treated as inconsistency of the file and cause stopping the calculation type int is type of the source This variable is still unused eid int is id number of the element where the source lies density double is the density of the source in volume of fluid per time unit Possitive values are sources negative are sinks Comments concerning 1 2 3 FLOW e Every inconsistency or error in the SRC file causes stopping the calculation These are especially Multiple usage of the same source number Difference between number of sources and actual number of data lines Reference to nonexisting element 19 ASCII post processing file format version 1 2 File format of this file comes from the GMSH system Following text is copied from the GMSH documentation The ASCII post processing file is divided in several section
16. of flow fic Number and type of required input files are depended on the type of the problem File of mesh is generated by using software GMSH which is distributed under the terms of the GNU GPL www geuz org File of neigbours is generated with using program NGH Structure of all input files are defined in the files description in detail Required inputs Problem dependent inputs mtr file msh file ngh file FLOW123D fbc file Figure 1 Scheme of calculation Output of the program generates pos files supported by the GMSH Eventualy it is possible using text output files for whole area specified area or elements Flow123D ini file format Flow123D version 03 10 08 Note All string values have maximal length MAXBUFF 1 1023 Section Global KEY TYPE DEFAULT DESCRIPTION Problem_type int NULL Type of solved problem Currently supported 1 steady saturated flow 3 variable density saturated flow Description string undefined Short description of solved problem any text Stop_time double 1 0 Time interval of the whole problem time units Save_step double 1 0 The output with transport is written every Save_step time units Density step double 1 0 Time interval of one density iteration in the varible density calculation type 3 time units Section Input KEY TYPE DEFAULT DESCRIPTION File_type int 1 Type of the input
17. orptionFraction material number lt sorption fraction coefficient gt text EndSorptionFraction DualPorosity material number lt mobile porosity coefficient gt lt immobile porosity coefficient gt lt nonequillibrium coefficient substance 0 gt lt nonequilibrium coefficient substance n 1 gt text EndDualPorosity Reaction material number lt mobile porosity coefficient gt lt immobile porosity coefficient gt lt nonequillibrium coefficient substance 0 gt lt nonequilibrium coefficient substance n 1 gt text EndReaction where file type int is equal 0 for the ASCII file format data size int the size of the floating point numbers used in the file Usually data size sizeof double number of materials int Number of materials defined in the file 11 material number int is the number index of the n th material These numbers do not have to be given in a consecutive or even an ordered way Each number has to be given only onece multiple definition are treated as inconsistency of the file and cause stopping the calculation exception Sorption section material type int is type of the material see table lt material type specific data gt format of this list depends on the material type lt storativity coefficient gt double coefficient of storativity geometry type int type of complement dimension parameter only for 1D and 2D mate rial for 1D element is s
18. res of the elements 22
19. s one format section enclosed between PostFormat EndPostFormat tags and one or more post processing views enclosed between View EndView tags PostFormat 1 2 file type data size EndPostFormat View view name nb time steps nb scalar points nb vector points nb tensor points nb scalar lines nb vector lines nb tensor lines nb scalar triangles nb vector triangles nb tensor triangles nb scalar quadrangles nb vector quadrangles nb tensor quadrangles nb scalar tetrahedra nb vector tetrahedra nb tensor tetrahedra nb scalar herahedra nb vector hexahedra nb tensor hexahedra nb scalar prisms nb vector prisms nb tensor prisms nb scalar pyramids nb vector pyramids nb tensor pyramids nb tert2d nb tert2d chars nb tert3d nb text3d chars lt time step values gt lt scalar point values gt lt vector point values gt lt tensor point values gt lt scalar line values gt lt vector line values gt lt tensor line values gt lt scalar triangle values gt lt vector triangle values gt lt tensor triangle values gt lt scalar quadrangle values gt lt vector quadrangle values gt lt tensor quadrangle values gt lt scalar tetrahedron values gt lt vector tetrahedron values gt lt tensor tetrahedron values gt lt scalar hexahedron values gt lt vector herahedron values gt lt tensor hexahedron values gt lt scalar prism values gt lt vector prism values gt lt tensor prism values gt lt scalar pyramid values
20. upported type 1 cross section area for 2D element is supported type 2 thickness lt geometry type specific coefficient gt double cross section for 1D element or thickness for 2D element substance id int refers to number of transported substance numbering starts on 0 sorption type int type 1 linear sorption isotherm type 2 Freundlich sorption isotherm type 3 Langmuir sorption isotherm lt sorption type specific data gt format of this list depends on the sorption type see table Note Section Sorption is needed for calculation only if Sorption is turned on in the int file lt sorption fraction coefficient gt double ratio of the mobile solid surface in the contact with mobile water to the total solid surface this parameter section is needed for calculation only if Dual_porosity and Sorption is together turned on in the ini file lt mobile porosity coefficient gt double ratio of the mobile pore volume to the total volume this parameter is needed only if Transport_on is turned on in the ini file lt immobile porosity coefficient gt double ratio of the immobile pore volu me to the total pore volume this parameter is needed only if Dual_porosity is turned on in the ini file lt nonequilibrium coefficient substance i gt double nonequilibrium coefficient for substance i Vi 0 n 1 where n is number of transported substances this parameter is needed only if Dual
21. ven file Control_file string NULL Name of control file for situation when Use_control_file YES NSchurs int 2 Number of Schur complements to use Valid values are 0 1 2 The last one should be the fastest Section Solver parameters KEY TYPE DEFAULT DESCRIPTION Solver_accuracy double le 6 When to stop solver run value of residum of matrix Useful values from le 4 to le 10 Bigger number faster run less accuracy Note For aditional documentation see manual of the solver i isol manual Section Output KEY TYPE DEFAULT DESCRIPTION Write_output file YES NO NO If set to YES writes output file Output_file string NULL Name of the output file type 1 Output_file_2 string NULL Name of the output file type 2 Output digits int 6 Number of digits used for floating point numbers in output file Output_file type int 1 Type of output file 1 GMSH like format 2 Flow data file 3 both files two separate names POS_set_view YES NO NO Write a header setting the view in GMSH to POS POS_view_params double 8 000 x y z angle of rotation RotationX 111 x y z scaling ScaleX 00 x y screen position shift TranslationX Write ftrans_out YES NO NO If set to YES writes output file for ftrans Cross_section YES NO NO If set to YES uses cross section output Cs_params double 7 zero Params for cross section x0 y0 z0 init
22. ype lt type specific data gt where lt where data gt number of tags lt tags gt text EndBoundaryConditions where file type int is equal 0 for the ASCII file format data size int the size of the floating point numbers used in the file Usually data size sizeof double number of conditions int Number of boundary conditions defined in the file condition number int is the number index of the n th boundary condition These num bers do not have to be given in a consecutive or even an ordered way Each number has to be given only onece multiple definition are treated as inconsistency of the file and cause stopping the calculation type int is type of the boundary condition See below for definitions of the types lt type specific data gt format of this list depends on the type See below for specification of the type specific data for particular types of the boundary conditions where int defines the way how the place for the contidion is prescribed See below for details lt where data gt format of this list depends on where and actually defines the place for the condition See below for details number of tags int number of integer tags of the boundary condition It can be zero lt tags gt number of tags int list of tags of the boundary condition Values are separated by spaces or tabs By default we set number of tags 1 where tag1 defines group of boundary conditions
Download Pdf Manuals
Related Search