WO2016072515A1 - 流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム - Google Patents

流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム Download PDF

Info

Publication number
WO2016072515A1
WO2016072515A1 PCT/JP2015/081402 JP2015081402W WO2016072515A1 WO 2016072515 A1 WO2016072515 A1 WO 2016072515A1 JP 2015081402 W JP2015081402 W JP 2015081402W WO 2016072515 A1 WO2016072515 A1 WO 2016072515A1
Authority
WO
WIPO (PCT)
Prior art keywords
regular expression
expression
flow pattern
word
pattern
Prior art date
Application number
PCT/JP2015/081402
Other languages
English (en)
French (fr)
Inventor
貴之 坂上
知郎 横山
Original Assignee
国立研究開発法人科学技術振興機構
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 国立研究開発法人科学技術振興機構 filed Critical 国立研究開発法人科学技術振興機構
Priority to JP2016557834A priority Critical patent/JP6401288B2/ja
Priority to US15/524,897 priority patent/US10540460B2/en
Priority to EP15858028.2A priority patent/EP3217303A4/en
Publication of WO2016072515A1 publication Critical patent/WO2016072515A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Definitions

  • the present invention relates to a flow pattern regular expression creating method, a regular expression creating apparatus, and a computer-executable program.
  • an optimization technique has been developed in which a fluid simulation is repeatedly performed on a structure while changing the design variable of the structure using an annealing method (simulated annealing method), a genetic algorithm method, or the like.
  • the search range is inevitably limited due to these restrictions.
  • the possibility that the optimal structure is locally optimal cannot be excluded. That is, there is a problem that the optimization range of the derived structure depends on where the search range is set, depending on the experience and intuition of the engineer as to where the search range is set.
  • the applicant of the present application is obtained by numerical calculation or experiment in the “flow pattern word expression method, word expression device, and program (hereinafter also referred to as“ word expression theory ”)” of Patent Document 1.
  • word expression theory assigns this character string to the flow and expresses the quantitative characteristics (eg, lift-drag ratio) of the flow that appears as a qualitative partial character string.
  • the present invention has been made in view of the above, and a flow pattern regular expression creating method, regular expression creating device, and a new expression method capable of providing a one-to-one correspondence with a flow pattern, And it aims at providing a program.
  • the present invention provides a regular expression creation method for creating a regular expression of a flow pattern in a multi-connected external region having topologically N (where N is an integer of 1 or more) holes.
  • the graph representation is a tree T with a unique root, a label, and a direction with respect to the structurally stable Hamilton vector field H defined by the flow pattern.
  • H (V, E) is assigned (where V is a set of points called vertices, E is a set of edges connecting the vertices), and is preferably visualized as a planar graph.
  • the graph representation is such that the parent vertex is v, the child vertex is w, the label assigned to the parent vertex v is l (v), and the child vertex is w. If the assigned labels are l (w) and v's child vertex set ⁇ (v), then v's child vertex set ⁇ (v) is rearranged according to a predetermined ordering rule, and w ⁇ (v) , L (v) to l (w) are desirably drawn side by side from left to right.
  • the flow pattern is (1) a suction spring out of two types of flow patterns that can be taken topologically in a single connected outer region having one hole.
  • Pattern I with a pair and two ss-ad-saddle connections
  • Pattern II with a pair of saddle points, a homoclinic saddle connection and two ss-saddle connections that connect it
  • a pattern O that does not have a suction spring pair in a double connected external region with two holes It is desirable that the number is one or more.
  • the graph creating step includes a counterclockwise ss-orbit that is closest to the suction and spring-out pair when the flow pattern has a suction and spring-out pair.
  • the region is converted so that it becomes the outermost region, and the trajectory that appears in the (ss-) saddle connection diagram of the converted flow pattern is extracted from the entire region, and the trajectory that appears in the (ss-) saddle connection diagram from the entire region.
  • a vertex is set for the connected component obtained by excluding all, the connected component on the outermost side is set as the root, the current component is set as the root, and the connected components that are in contact with the current component are children of the current component.
  • the flow pattern defines five types of operations that can be taken topologically when one hole is added to the flow pattern, starting from one or more of the pattern words. It is desirable that the flow pattern diagram is created by repeating the operation of assigning any one of the operated words until the number of holes becomes N.
  • the word expression conversion step it is possible to perform forward substitution on a regular expression of a basic type among regular expressions corresponding to a tree with a root, a label, and a direction. It is desirable to convert the allowable regular expression into the word expression.
  • the word expression includes two types of flow patterns that have two holes in addition to two kinds of flow patterns that can be taken topologically in a single connected outer region having one hole.
  • a pattern word that defines a total of three types of flow patterns, with the addition of a pattern that does not have a pair of suction and outflow in the externally connected region, it can be used topologically when adding one hole to the flow pattern. It is desirable that it is a symbol word formed by assigning any one of the operation words defining five types of operations to the number of added holes.
  • an upper limit and a lower limit of the design parameter are set, and the upper limit and the lower limit of the design parameter are defined.
  • a plurality of parameters are selected from the parameter region, and a flow experiment and / or numerical calculation is performed for each of the selected plurality of parameters, and the word expression and / or the numerical calculation is performed on the result of the experiment and / or numerical calculation.
  • a design parameter selection step of assigning a regular expression and selecting a design parameter having the word expression and / or the regular expression indicating an optimum state among the assigned word expression and / or the regular expression as a candidate for the design parameter It is desirable to include.
  • the design parameter selection step further includes a step of performing an optimization design using the selected design parameter.
  • a transition is made from the regular expression of the current flow pattern to the regular expression of the target flow pattern. It is desirable to include a second design parameter selection step of changing the design parameter so as to select the design parameter that becomes the regular expression of the target flow pattern.
  • the regular expression and word expression of the current flow pattern are changed to the regular expression and word expression of the target flow pattern. It is desirable to change the design parameters to select the design parameters that are the regular expressions and word expressions of the target flow pattern.
  • the second design parameter selection step transitions from the regular expression and word expression of the current flow pattern to the regular expression and word expression of the target flow pattern.
  • the present invention provides a regular expression for creating a regular expression of a flow pattern in a multi-connected outer region having topologically N (where N is an integer of 1 or more) holes.
  • the present invention provides a regular expression for creating a regular expression of a flow pattern in a multi-connected outer region having topologically N (where N is an integer of 1 or more) holes.
  • a program installed in a creation device a graph expression creation step for creating a graph expression corresponding to the flow pattern on a one-to-one basis, and a regular expression for creating a regular expression from the graph expression created in the graph expression creation step
  • An expression creating step is executed by a computer.
  • a flow pattern regular expression creation method, regular expression creation device, and program capable of providing a new expression method that can be associated with a flow pattern on a one-to-one basis are provided. There is an effect that becomes possible.
  • FIG. 1 is an explanatory diagram for explaining the outline of the present invention.
  • FIG. 2-A is a flowchart for explaining a flow pattern regular expression creation method according to the present embodiment.
  • FIG. 2-B is a flowchart for explaining an example of the graph expression creation step of FIG. 2-A.
  • FIG. 3 is a flowchart shown in the outline of the present embodiment.
  • FIG. 4 is a diagram schematically showing a structure-stable flow pattern.
  • FIG. 5 is a diagram describing all the characteristic trajectories (streamlines) that perform topological classification of structurally stable flow patterns in a region.
  • FIG. 6 is a diagram schematically showing three types of flow patterns serving as an initial structure.
  • FIG. 7 is a diagram schematically showing five types of operations for forming a structurally stable flow by adding one hole.
  • FIG. 8 is a diagram showing all classifications of flow patterns when there are two structures and a uniform flow.
  • FIG. 9 is a diagram showing all the flow patterns generated by the operations B 0 , B 2 , and C and the parent-child relationships excited correspondingly to the connected subset (vertex set).
  • FIG. 10A is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10-B is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10C is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10A is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10-B is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10D is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 10E is a flowchart for explaining a conversion process of a saddle connection diagram tree in the O series.
  • FIG. 12 is a diagram showing the local structure of the ss-saddle connection diagram generated by the operations A 0 , A 2 , and C and the visualization of the graph expression obtained correspondingly.
  • FIG. 13-A is a flowchart for explaining the conversion process to the ss-saddle connection diagram tree in the I and II series.
  • FIG. 13-A is a flowchart for explaining the conversion process to the ss-saddle connection diagram tree in the I and II series.
  • FIG. 13B is a flowchart for explaining the conversion process to the tree of the ss-saddle connection diagram in the I and II series.
  • FIG. 13-C is a flowchart for explaining the conversion process to the tree of the ss-saddle connection diagram in the I and II series.
  • FIG. 13-D is a flowchart for explaining the conversion process of the ss-saddle connection diagram tree to the I and II series.
  • FIG. 14A is a diagram showing a graph expression of a structurally stable streamline pattern having the word expression IA 0 C and its regular expression.
  • FIG. 14B is a flowchart for explaining a process of converting a tree into a regular expression.
  • FIG. 15 is a diagram showing five fundamental types TN and regular expressions N corresponding to the five fundamental types TN .
  • FIG. 16 is a diagram showing a list of 14 forward replacements.
  • FIG. 17 is a diagram showing changes in forward replacement and streamline structure when operations A 0 , A 2 , and C are performed.
  • FIG. 18 shows excitation along with the change in streamline structure obtained when operations B 0 , B 2 , and C are applied to the closed orbit included in the simple connected components ⁇ 0 and ⁇ 2 in the (ss-) saddle connection diagram. It is a figure for demonstrating forward replacement of the regular expression to do.
  • FIG. 16 is a diagram showing a list of 14 forward replacements.
  • FIG. 17 is a diagram showing changes in forward replacement and streamline structure when operations A 0 , A 2 , and C are performed.
  • FIG. 18 shows excitation along with the change in streamline structure obtained when operations B 0 , B 2 , and C are applied to the closed orbit included in the simple connected components ⁇
  • FIG. 19 is a diagram illustrating a list of backward replacement [B 2 ], [B 0 ], and [C] of regular expressions induced by the reverse operation of the operations B 0 , B 2 , and C.
  • FIG. 20 is a diagram illustrating a backward replacement of a regular expression corresponding to a change in the local streamline structure in the (ss-) saddle connection diagram.
  • Figure 21 is a regular expression o ⁇ (o 2 (o 2 (o 0 (- 0 (- 0, - 0))), + 2 (+ 2)) of the structural stability Hamiltonian vector field with ss-saddle It is a figure which shows connection diagram.
  • FIG. 22-A is a flowchart for explaining a process of converting a regular expression into a word expression.
  • FIG. 22-B is a flowchart for explaining a process of converting a regular expression into a word expression.
  • FIG. 22-C is a flowchart for explaining a process of converting a regular expression into a word expression.
  • FIG. 23 is a block diagram illustrating an example of a regular expression creation device to which the exemplary embodiment is applied.
  • FIG. 24 is a flowchart for explaining a method of designing an object in a fluid using word expressions and regular expressions.
  • FIG. 25 is a diagram for explaining an example of a specific example for explaining a method for designing an object in a fluid using word expressions and regular expressions.
  • FIG. 26 is a diagram for explaining an example of a specific example for explaining a method for designing an object in a fluid using word expressions and regular expressions.
  • FIG. 27 is an explanatory diagram for explaining a method for optimizing a flow around an object in a fluid using word expressions and regular expressions.
  • FIG. 28 is an explanatory diagram for explaining a method for optimizing a flow around an object in a fluid using word expressions and regular expressions.
  • FIG. 29 is a flowchart for explaining an example of a method for optimizing a flow around an object in a fluid using word expressions and regular expressions.
  • FIG. 1 is an explanatory diagram for explaining the outline of the present invention.
  • Patent Document 1 International Publication No. 2014/041917, Japanese Patent Application No. 2012-203601, PCT / JP2013 / 070939
  • this character string is given to the flow pattern obtained by numerical calculation and experiment, and the word expression that expresses the quantitative characteristic (for example, lift-drag ratio) of the flow that appears there as a qualitative partial character string.
  • a theory was proposed.
  • Patent Document 1 has already filed a method for creating a word expression from a flow pattern (P1).
  • a method for creating a fluid transition path from a word expression (P2) has been filed.
  • a unique word expression can be given to the flow pattern, but a plurality of different flow patterns correspond to the character string given in the word expression. Therefore, in the flow pattern regular expression creating method of the present embodiment, a graph expression corresponding to the flow pattern is created one-to-one (T1), and a regular expression is created from the graph expression (T2).
  • T3 regular expression is created from the regular expression (T3).
  • the fluid transition path can be created directly from the word expression. Furthermore, in the flow pattern regular expression creation method of the present embodiment, the fluid transition path can also be created directly from the regular expression (T4).
  • FIG. 2A is a flowchart for explaining the outline of the flow pattern regular expression creation method according to the present embodiment.
  • the regular expression creation method of the flow pattern according to the present embodiment is based on topologically N (where N is an integer of 1 or more) holes in a multi-connected outer region.
  • a regular expression of a flow pattern is created, and is roughly divided into a graph expression creation process (step S1) for creating a graph expression corresponding to a flow pattern on a one-to-one basis, and a graph expression creation process of step S1.
  • a regular expression creating step (step S2) for creating a regular expression from the graph expression; and a word expression converting step (step S3) for converting the regular expression created in the regular expression creating step of step S2 into a word expression.
  • V is called a vertex
  • a set of points, E is a set of edges connecting vertices) and may be visualized as a planar graph.
  • the graph representation is such that the parent vertex is v, the child vertex is w, the label assigned to the parent vertex v is l (v), and the label assigned to the child vertex w is l (w), v.
  • the child vertex set ⁇ (v) is used, the child vertex set ⁇ (v) of v is rearranged according to a rule of a predetermined order relation, and from w ( ⁇ ) to l (w) for w ⁇ (v). You may include what was drawn with the arrow arranged from left to right.
  • the flow patterns are (1) two ss- ⁇ -saddle connections, which have a suction-and-out pair of two types of flow patterns that can be taken topologically in a single connected outer region having one hole.
  • One or more pattern words consisting of a pattern II having a homoclinic saddle connection and two ss-saddle connections, and (3) a pattern O having no suction spring-out pair in a double-connected external region having two holes. It may be.
  • the word representation does not have a suction-out pair in a double connected outer region with two holes.
  • the operation words that define the five types of operations that can be taken topologically when adding one hole to the flow pattern with respect to the pattern words that define a total of three types of flow patterns with added patterns It is good also as a symbol word formed by giving any one word by the number of the added holes.
  • FIG. 2-B is a flowchart for explaining an example of the graph expression creation process of FIG. 2-A.
  • the graph expression creation process first, when the flow pattern has a suction-and-out pair, the region including the counterclockwise ss-orbit closest to the suction-and-out pair is the first. Conversion is performed so as to be in the outer region (step S31).
  • the flow pattern starts from one or more of the above-described pattern words, and is one of the operation words that define five types of operations that can be taken topologically when a hole is added to the flow pattern. It is good also as a flow pattern figure with the number of the holes created by repeating the operation which provides a word until the number of holes becomes N pieces.
  • the route of the ss-saddle connection diagram (flow pattern diagram) of I-word and II-word is converted as shown in FIG.
  • the trajectory that appears in the (ss-) saddle connection diagram of the converted flow pattern is extracted from the entire region (step S32).
  • a vertex is set to the connected component obtained by excluding all trajectories appearing in the (ss-) saddle connection diagram from the entire region, and the outermost connected component is set as the root (step S33).
  • the current component is set as the root (step S34).
  • a connected component that touches the boundary with the current component is a child of the current component, labels are assigned according to the trajectory corresponding to the boundary, and they are arranged in accordance with a predetermined label order relationship (step S35).
  • the child of the current component is set as the current component, and step S35 is repeatedly executed until there are no more children (step S36).
  • the flow pattern regular expression creation method of the present embodiment may execute a design parameter selection step of selecting design parameter candidates for an object in a fluid.
  • a design parameter selection step for example, when selecting a design parameter candidate for an object in a fluid, an upper limit and a lower limit of the design parameter are set, and a parameter region defined by the upper limit and the lower limit of the design parameter is set. A plurality of parameters are selected, and a flow experiment and / or numerical calculation is performed for each of the selected parameters, and the word expression and / or the regular expression is used for the result of the experiment and / or numerical calculation.
  • a design parameter having the word expression and / or the regular expression indicating an optimal state may be selected as the design parameter candidate.
  • the design parameter selection step may further include a step of performing optimization design using the selected design parameter.
  • the flow pattern regular expression creation method of the present embodiment may execute a second design parameter selection step of selecting a design parameter for an object in a fluid.
  • the second design parameter selection step for example, when the optimum design parameter is selected for an object in the fluid, the current flow pattern regular expression is changed to the target flow pattern regular expression.
  • the design parameter may be changed to a design parameter that is a regular expression of the target flow pattern.
  • the word expression is used in addition to the regular expression so that the current flow pattern regular expression and word expression are changed to the target flow pattern regular expression and word expression.
  • the design parameter may be changed, and a design parameter that becomes a regular expression and a word expression of the target flow pattern may be selected.
  • one or a plurality of paths for transitioning from the current flow pattern regular expression and word expression to the target flow pattern regular expression and word expression are acquired, and the acquired 1 Or, select a path that can be transitioned from the current flow pattern by changing the design parameter among multiple paths, and change the design parameter around the changed design parameter for the selected path. You may decide to select the design parameter used as the regular expression and word expression of a flow pattern.
  • FIG. 3 is a flowchart showing an outline of the word expression algorithm.
  • step SA-1 in order to form a word representation of the flow pattern in a connected external region having N holes topologically, topologically in a single connected external region having one hole.
  • a pattern word that defines a total of three types of flow patterns for example, I, II, O is assigned (step SA-1).
  • an operation word defining five types of operations that can be taken topologically when adding one hole to the flow pattern with respect to the word given in step SA-1 (for example, , A 0 , A 2 , B 0 , B 2 , C) (Step SA-2), and repeats the process of Step SA-2 until the number of holes becomes N
  • step SA-3 a word expression corresponding to the connected external region having N holes is formed.
  • the pattern word assignment (step SA-1) is performed first and the operation word assignment (step SA-2) is performed later.
  • the present invention is not limited to this, and the operation word assignment is performed first.
  • the pattern word may be given later.
  • the “connected external region” is an expression including a single connected external region and a multiple connected external region
  • the “single connected external region” is in two dimensions (plane).
  • a region is a region with a single hole
  • a “multiple connected external region” is a region in two dimensions (plane) with a plurality of holes.
  • the expression “hole” here is a mathematical abstract expression, but various expression forms may be applied. For example, when focusing on a uniform flow that occurs relatively when a vehicle moves, if there is one or more structures on the cross-section of the vehicle along the uniform flow, the region is simply It can be treated as a connected external region or a multiple connected external region.
  • one or more holes and a single connection or multiple connection Can be treated as a connected external region.
  • a flow such as “one or more obstacles in the flow” is treated as a flow in the connected external region in this embodiment.
  • an isolated vortex structure or a flow structure (such as an elliptical stopping point) having a periodic orbit around it can be regarded as a “hole”.
  • Topological is a technical term in mathematics and refers to a field of geometry, also called topology (topology).
  • topology topology
  • triangles and quadrilaterals are considered to be different figures due to the difference in the number of their corners, but from the viewpoint of topological geometry, triangles and quadrilaterals can move to each other by deforming rubber bands, for example.
  • topology topology
  • the connected external region is expressed as D ⁇ (M) with respect to the number of holes M + 1. For example, when there is only one hole, it is a single connected external region D ⁇ (0), and when there are two holes, it is a double connected external region D ⁇ (1).
  • topological structure a specific structure that characterizes the flow
  • Those that cannot be deformed by continuous (ie, not cut or stuck) deformation of a particular structure are considered different flows.
  • An uncompressed flow may be handled as the “flow” in the present embodiment.
  • the incompressibility of a fluid refers to a property that hardly changes its volume even when force is applied. When we consider normal water and air flow on the scale of daily life, it is generally safe to consider this flow framework. Note that the present invention is not limited to this, and a compressible flow may be handled in calculation.
  • the fluid in the present embodiment may be handled as a non-viscous fluid.
  • An inviscid fluid can treat the boundary condition as a “slip condition”.
  • the boundary condition is the “zero boundary condition”.
  • the word conversion algorithm for non-viscous fluids may be described in particular, but this embodiment can also be applied to viscous fluids by appropriate boundary expansion.
  • FIG. 4 is a diagram schematically showing a flow pattern.
  • FIG. 4C schematically shows a boundary having four ⁇ 4-saddle points (the definition will be described later).
  • the vortex is an element that creates a flow rotating around it.
  • the uniform flow is a basic flow in terms of river flow, and is a flow that crosses the entire region.
  • the uniform flow is a flow such as a relative air flow or a water flow as viewed from an observer riding on the vehicle. That is, in the coordinate system of the moving object, even if air or water is actually stationary, it is possible to imagine the relative flow that flows from the infinity point.
  • the element that constitutes a uniform flow is called a 1-source-sink point, and the correctness of the rule can be proved mathematically.
  • a plane may be projected onto a spherical surface by a projection method called stereographic projection (stereo projection / stereo projection) in mathematics.
  • stereographic projection stereo projection / stereo projection
  • the uniform flow has a flow structure such as a flow source and a suction pair in the spherical north pole, and it can be mathematically shown to correspond to the flow field as shown in FIG.
  • the position of the north pole and the south pole can be appropriately shifted using the fact that the spherical surface has a highly symmetric shape, so the point at infinity is the south pole,
  • a flow as shown in FIG. 4B can be obtained near the origin corresponding to the south pole.
  • the flow field of the entire plane can be expressed as a bounded area having a shape as shown in FIG.
  • D ⁇ (M) is moved to the unit circle outer region by selecting the pair of suction and spring at infinity point as the origin and one of M + 1 boundaries by this conversion method.
  • a connected region inside a unit circle including M boundaries inside is sometimes written as D z (M). All subsequent descriptions are made using this D z (M). Therefore, the expression as shown in FIG. 4B is equivalent to a flow in which a uniform flow is contained in the entire plane through an appropriate projection method. In the description of the present embodiment, it is convenient to schematically show the flow, and therefore, it may be expressed in the drawing using such a projection method.
  • FIG. 5 is a diagram describing all the characteristic trajectories (streamlines) that perform topological classification of the structure-stable flow in such a connection region D z (M).
  • ss-orbit the trajectory returning from the pair of suction and outflow and returning to itself.
  • Fig. 5 (b) the trajectory that exits the pair of suction and outflow and connects to the boundary is called ss- ⁇ -saddle connection, and the trajectory is connected as shown in Fig. 5 (c).
  • the point on the border is called ss- ⁇ -saddle.
  • FIG. 5 (e) the trajectory leading out from a point on a certain boundary and leading to a point on the same boundary, not from a pair of suction and outflow, is called ⁇ -saddle connection, and FIG. 5 (d) The points on the boundary connected by this are called ⁇ -saddle.
  • a hyperbolic stop point not on the boundary as shown in FIG. 5 (h) is called a saddle point (saddle point).
  • the trajectory leading to this saddle point is called ss-saddle connection.
  • a closed curve orbit that creates a boundary or a vortex is called a closed orbit
  • Fig. 5 (i) a trajectory that leaves the saddle point and returns to itself. It is called homoclinic saddle connection. It can be mathematically proved that the target structurally stable flow can only be expressed by a combination of these trajectories.
  • step SA-2 in the above-described step SA-2, one hole and the flow structure associated therewith are added to the flow of the connected external region D z (M ⁇ 1) having M holes.
  • the flow field of many multiple connected external regions D z (M) is constructed recursively. Therefore, the simplest hole is a single connected external region D z (0) or a double connected external region D z (1), which is the initial structure of these inductive structures given in step SA-1 ing.
  • FIG. 6 is a diagram schematically showing three types of structure-stable flow patterns as an initial structure.
  • FIGS. 6A and 6B there are two types of flows, a pattern I and a pattern II, in which a hole is in one single connected external region D z (0). Both of these patterns have a pair of suction and spring, and it can be proved mathematically that there are only these two types.
  • the double-coupled outer region D z (1) is composed of these, but does not have a suction-and-out pair. Is not constructed from here, the initial flow necessary for constructing the flow is the pattern O schematically shown in FIG. 6 (c).
  • step SA-2 the five types of operations that can be taken topologically are: 1) A 0 operation that replaces one ss-orbit with one saddle point, a homoclinic saddle connection that connects the ss-orbit and a hole inside, and two ss-saddle connections.
  • FIG. 7 is a diagram schematically showing five types of operations for forming a structurally stable flow by adding one hole.
  • operation A 0 and A 2 is performed for a single ss-orbit. Further, as shown in FIG. 7 (b), the operation B 0 and B 2 is performed for a single closed orbit. Further, as shown in FIG. 7C, the operation C is performed on the boundary already having ⁇ -saddle. It can be mathematically proved that there are only five types of operations that enable such a process while maintaining structural stability.
  • a flow pattern diagram obtained by repeating this step is referred to as ss-saddle connection diagram when starting from patterns I and II, and as a saddle connection diagram when starting from pattern O.
  • FIG. 8 is a diagram showing the entire classification of the structure-stable flow pattern when there are two structures and a uniform flow. As shown in FIG. 8, all the flow patterns in the double connected region D z (1) are given to the patterns I and II of the initial structure in the single connected external region D z (0) by assigning operation words. Can be described. However, all the flow patterns shown in FIG.
  • the constraints will be described as follows. That is, as described above with reference to FIG. 7, since the operations A 0 and A 2 are performed on one ss-orbit, the presence of one ss-orbit is indispensable as a premise for performing this operation. It becomes. Further, since the operations B 0 and B 2 are performed on one closed orbit, the presence of one closed orbit is indispensable as a premise for performing this operation. Further, since the operation C is performed on a boundary having two or more ⁇ -saddles, the presence of a boundary having ⁇ -saddles is indispensable for performing this operation. Therefore, the rules for arranging the pattern words are different depending on where the pattern words are started from I, II, and O. The rules for arranging the character strings starting from the pattern words I, II, and O, which are derived based on the above constraint conditions, will be described below.
  • Such a character string is called an O-series word (O-word), and the correctness of the rule can be mathematically proved.
  • O-word O-series word
  • I-1 Operations that can be performed are all of A 0 , A 2 , B 0 , B 2 , and C.
  • I-2 the word expression starting from I lists these five types of operation words.
  • I-2) In order to include the word B 0 or B 2 in the word expression of the operation sequence, C or A 0 must always exist before that.
  • I-word Such a character string is called an I-series word (I-word), and the correctness of the rule can be mathematically proved.
  • Such a character string is called a II-series word (II-word), and the correctness of the rule can be mathematically proved.
  • the graph T (V, E) is a set given as a pair of a set of points (vertex set) called “vertex” V and a set E called “edge” connecting the vertices. It is.
  • a graph can have various structures, but in the graph representation theory of the present embodiment, the following set of graphs having a specific structure in the entire set of graphs is considered.
  • Tree and tree refer to a graph in which two arbitrary vertices are connected only by one edge.
  • a routed graph refers to a graph in which a specific vertex (hereinafter referred to as a root) exists.
  • a root a specific vertex
  • “Directed” means that all edges have a parent-child order.
  • the edge from the vertex v ⁇ V to w ⁇ V in the directed graph is expressed as v ⁇ w ⁇ E.
  • v is called the parent of w and w is called the child of v.
  • Write ⁇ (v) to represent the set of all children of vertex v ⁇ V. That is, ⁇ (v): ⁇ w ⁇ V
  • the number of child vertices # ⁇ (v) included in the set is called v out-degree, and conversely, the number of edges entering v is called v in-degree.
  • a graph is labeled as being labeled with a specific label assigned to all vertices.
  • a tree with root, label, and orientation is considered. Further, the root is considered to be a vertex where in-degree is zero, that is, there is no incoming edge.
  • the vertex label any one of ⁇ o ⁇ , o 0 , o 2 , + ⁇ , + 0 , + 2 , ⁇ , ⁇ 0 , ⁇ 2 ⁇ is assigned.
  • the label is o ⁇ (o 0, o 2 , + ⁇ , + 0, + 2, - ⁇ , - 0, or - 2) a subset of the vertex set has become Vo ⁇ (Vo 0, Vo 2 , V + ⁇ , V + 0 , V + 2 , V ⁇ ⁇ , V ⁇ 0 , or V ⁇ 2 ).
  • V the label assigned to the vertex.
  • the following set is defined.
  • Vo Vo ⁇ ⁇ Vo 0 ⁇ Vo 2
  • V + V + ⁇ ⁇ V + 0 ⁇ V + 2
  • V- V- ⁇ ⁇ V- 0 ⁇ V- 2
  • V 0 Vo 0 ⁇ V + 0 ⁇ V- 2
  • V 2 Vo 2 ⁇ V + 2 ⁇ V- 2
  • ⁇ symbols represent a disjointed union.
  • ⁇ (v) can be divided into a child vertex set of vertex v as follows.
  • H a structurally stable Hamiltonian vector field expressed by O-word on the two-dimensional region D z (M), and let D be its saddle connection diagram.
  • a set C H D z (M) ⁇ D consisting of several connected components (including infinite periodic orbits) is an open set. it can.
  • Each connected component is set as a vertex set V. Apex of the root is a connected component in the outermost of these connected components, when periodic orbits included as the label is counterclockwise + phi, when a clockwise - assign phi.
  • FIG. 9 shows all streamline patterns generated by the operations B 0 , B 2 , and C, and the parent-child relationships excited correspondingly to the connected subset (vertex set).
  • (A) is B 0
  • (b) is B 0
  • (c) is B 0 B 0
  • (d) is B 2 C k ⁇ 1 , k ⁇ 1
  • (e) is B 2 C 1 , l ⁇ 1
  • (F) is B 2 C k + l ⁇ 1 , k, l ⁇ 1
  • (g) is B 0 B 2 C k , k ⁇ 1
  • (h) is B 0 B 2 C l , l ⁇ 1.
  • v represents a parent connected subset, and its child connected subset is expressed as w, y, z, and the like. Since the broken line represents a set of parent connected parts, only one periodic trajectory included therein is drawn with a direction.
  • FIG. 9 (a) can be constructed character pattern outward 8 made by the operation B 0.
  • the number of elements included in the set of child connected components is # ⁇ + 0 (v), # ⁇ -0 (v) ⁇ 2, and any non-negative integer can be selected for the child vertex sets # ⁇ + 2 (v) and # ⁇ - 2 (v).
  • the order is cyclically arranged with respect to the connected portion written as j in FIG. That is, a specific connected portion is selected as y 1 from these, and thereafter the connected components are arranged counterclockwise.
  • the child connected components written as j j can be simply arranged counterclockwise.
  • N a set of non-negative integers
  • N ⁇ a set of finite sequences of non-negative integers
  • G a tree
  • s, s ′, t, u vertex ids
  • T a subset of vertex ids.
  • FIGS. 10A to 10E are flowcharts for explaining the conversion process of a saddle connection diagram to a tree in the O series.
  • the conversion process to the tree of the saddle connection diagram in the O series shown in FIGS. 10-A to 10-E can be executed by a device such as a computer.
  • step S101 input is set as saddle connection diagram D
  • step S101 saddle connection diagram D is converted so that there is root in outermost
  • step S106 it is determined whether or not root 0 has the shape (pattern) shown in FIG. If root 0 has the form shown in FIG. 9A (“Yes” in step 106), the vertex ids of w 1 and w 2 are set to 00, 01, and the tree G is set as shown in FIG. 9A.
  • step S106 if root 0 does not have the shape of FIG. 9A ("No" in step 106), the process proceeds to step S108.
  • step S108 it is determined whether or not root 0 has the left shape of FIG. If root 0 has the left shape of FIG. 9B (“Yes” in step 108), the vertex id of w is set to 00, the tree G is the label of the tree of FIG. 9B. Is replaced with a set of id and label (ie, G is ((0, ⁇ ⁇ ) ⁇ (00, ⁇ 0 ))), and ⁇ 00 ⁇ is added to the set T (step S109). That is, T ⁇ TU ⁇ 00 ⁇ . Thereafter, the process proceeds to step S114.
  • step S108 if root 0 does not have the left shape of FIG. 9B ("No" in step S108), the process proceeds to step S110.
  • step S110 it is determined whether or not root 0 has the shape shown in FIG. If root 0 has the shape of FIG. 9D (“Yes” in step S110), y 1 ,. . . , Y k with vertex ids 00,. . . , 0k ⁇ 1, and the tree G is defined as a tree in which the label of the tree in FIG. 9D is replaced with a set of id and label, and ⁇ 00,. . . , 0k-1 ⁇ is added (step S111). That is, T ⁇ TU ⁇ 00,. . . , 0k-1 ⁇ . Thereafter, the process proceeds to step S114.
  • step S110 if root 0 does not have the shape of FIG. 9D (“No” in step S11), tree G is defined as ⁇ ⁇ (step S112), and this tree G is output. (Step S113), the process ends.
  • step S114 it is determined whether or not an element t greater than s exists in T.
  • s is the largest element in T (“No” in step S114)
  • all the vertex ids of the tree G are removed, the remaining tree is set as G (step S115), and the tree G is output (step S113). ), The process is terminated.
  • step S117 it is determined whether or not the outer boundary of the vertex s is composed of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the shape shown in FIG.
  • the outer boundary of the vertex s is composed of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the shape of FIG. 9A (“Yes” in step 117)
  • w 1 , w 2 Is defined as s0, s1
  • the tree G is defined as a tree obtained by replacing the label of the tree in FIG. 9A with a set of id and label, and ⁇ s0, s1 ⁇ is added to T (step S118). . Thereafter, the process returns to step S114.
  • step S117 when the outer boundary of the vertex s is not one homoclinic saddle connection and saddle in step S117, or the inner boundary of the vertex s does not have the shape of FIG. 9A (“No” in step S117). The process proceeds to step S119.
  • step S119 it is determined whether or not the outer boundary of the vertex s is composed of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the left shape of FIG. 9B. If the outer boundary of the vertex s consists of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the left shape of FIG. 9B (“Yes” in step S119), the vertex id of w Is defined as s0, the tree G is defined as a tree in which the label of the left tree in FIG. 9B is replaced with a set of id and label, and ⁇ s0 ⁇ is added to T (step S120). Thereafter, the process returns to step S114.
  • step S119 if the outer boundary of the vertex s is not one homoclinic saddle connection and saddle in step S119, or if the inner boundary of the vertex s does not have the left shape of FIG. 9B (“No” in step S119) '), The process proceeds to step S121.
  • step S121 it is determined whether or not the outer boundary of the vertex s is composed of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the shape shown in FIG. 9D.
  • the outer boundary of the vertex s is composed of one homoclinic saddle connection and saddle, and the inner boundary of the vertex s has the shape of FIG. 9D (“Yes” in step 121)
  • the tree G is defined as a tree in which the label of the tree in FIG. 9D is replaced with a set of id and label, and T is set to ⁇ s0,. . . , Sk-1 ⁇ is added (step S122). Thereafter, the process returns to step S114.
  • step S121 if the outer boundary of the vertex s is not one homoclinic saddle connection and saddle in step S121, or the inner boundary of the vertex s does not have the shape shown in FIG. 9D (“No” in step S121). The process proceeds to step S123.
  • step S123 it is determined whether or not the vertex s has the left shape of FIG.
  • the process proceeds to step S124.
  • step S124 the one sigma label sigma 0 vertex s is + (i.e., flow in the s is either a counter-clockwise) to determine.
  • ⁇ of the label ⁇ 0 of the vertex s is + (that is, the flow in s is counterclockwise) (“Yes” in step S124)
  • the vertex ids of w 1 , w 2 , and w 3 are set.
  • s0, s1, s2 are defined, and the tree G is defined as a tree in which the label of the left tree in FIG. 9C is replaced with a set of id and label, and ⁇ s0, s1, s2 ⁇ is added to T (step S125). Thereafter, the process returns to step S114.
  • step S124 when ⁇ of label ⁇ 0 of vertex s is not +, that is, when ⁇ of label ⁇ 0 of vertex s is ⁇ (that is, the flow in s is clockwise) (
  • step S124, “No”) the vertices id of w 1 , w 2 , w 3 are defined as s1, s2, s0, and the tree G is set to the left tree label in FIG. The replaced tree is determined, and ⁇ s0, s1, s2 ⁇ is added to T (step S126). Thereafter, the process returns to step S114.
  • step S123 when the vertex s does not have the left shape of FIG. 9C (“No” in step S123), the process proceeds to step S127.
  • step S127 it is determined whether or not the vertex s has the right shape of FIG.
  • the process proceeds to step S128.
  • step S1208 it is determined whether or not ⁇ of label ⁇ 0 of vertex s is + (that is, the flow in s is counterclockwise).
  • ⁇ of label ⁇ 0 of vertex s is + (ie, the flow in s is counterclockwise) (“Yes” in step S128)
  • the vertex ids of w 1 and w 2 are set to s0 and s1.
  • the tree G is defined as a tree in which the label of the right tree in FIG. 9C is replaced with a set of id and label, and ⁇ s0, s1 ⁇ is added to T (step S129). Thereafter, the process returns to step S114.
  • step S128 when ⁇ of label ⁇ 0 of vertex s is not +, that is, when ⁇ of label ⁇ 0 of vertex s is ⁇ (that is, the flow in s is clockwise) ( “No” in step S128), the vertex ids of w 1 and w 2 are defined as s1 and s0, and the tree G is defined as a tree in which the label of the right tree in FIG. , T is added to ⁇ s0, s1 ⁇ (step S130). Thereafter, the process returns to step S114.
  • step S127 when the vertex s does not have the right shape of FIG. 9C (“No” in step S127), the process proceeds to step S131.
  • step S131 it is determined whether the vertex s has the shape shown in FIG.
  • y 1 ,. . . , Y l , w are the vertex ids s1,. . . , Sl, s0
  • the tree G is defined as a tree in which the label of the tree in FIG. 9H is replaced with a set of id and label, and ⁇ s0,. . . , Sl ⁇ are added (step S132). Thereafter, the process returns to step S114.
  • step S131 when the vertex s does not have the shape of FIG. 9H (“No” in step S131), the process proceeds to step S133.
  • step S133 it is determined whether or not the vertex s has the left shape of FIG.
  • the tree G is defined as a tree in which the label of the left tree in FIG. 9G is replaced with a set of id and label, and ⁇ s0,. . . , Sk + 1 ⁇ are added (step S134). Thereafter, the process returns to step S114.
  • step S133 when the vertex s does not have the left shape of FIG. 9G (“No” in step S133), the process proceeds to step S135.
  • step S135 it is determined whether or not the vertex s has the right shape in FIG.
  • w, z 1 ,. . . , Z k vertices id are s0, s1,. . . , Sk
  • the tree G is defined as a tree in which the label of the right tree in FIG. 9G is replaced with a set of id and label, and ⁇ s0,. . . , Sk ⁇ is added (step S136). Thereafter, the process returns to step S114.
  • step S135 when the vertex s does not have the right shape of FIG. 9G (“No” in step S135), the process proceeds to step S138 (in this case, as shown in step S137, the vertex s is in the form of FIG. 9 (f)).
  • step S138 it is determined whether or not ⁇ of label ⁇ 0 of vertex s is + (that is, the flow in s is counterclockwise).
  • ⁇ of label ⁇ 0 of vertex s is + (ie, the flow in s is counterclockwise) (“Yes” in step S138)
  • the tree G is defined as a tree obtained by replacing the label of the tree in FIG. 9F with a set of id and label, and ⁇ s0,. . . , Sl + k ⁇ 1 ⁇ are added (step S139). Thereafter, the process returns to step S114.
  • step S138 when ⁇ of label ⁇ 0 of vertex s is not +, that is, when ⁇ of label ⁇ 0 of vertex s is ⁇ (that is, the flow in s is clockwise) (step S138).
  • “No”) y 1 ,. . . , Y k , z 1 ,. . . , Z l with the vertex ids sl,. . . , Sl + k-1, s0,. . . , Sl ⁇ 1, and the tree G is defined as a tree obtained by replacing the label of the tree in FIG. 9F with a set of id and label, and ⁇ s0,. . . , Sl + k ⁇ 1 ⁇ is added (step S140). Thereafter, the process returns to step S114.
  • the above operation is the conversion process to the tree of the saddle connection diagram in the O series.
  • D is an ss-saddle connection diagram that can be expressed in I-word or II-word.
  • the broken line is one ss-orbit included in the connected component that is the root, along with its direction.
  • (a) is a ss-saddle connection diagram expressed by the word expression IA 0 A 0 , and the flow direction of the ss- ⁇ -saddle connection is from right to left.
  • (B) is a ss-saddle connection diagram having the word expression IA 0 A 0, and the direction of the ss- ⁇ -saddle connection is from left to right.
  • (C) is a ss-saddle connection diagram having the word expression IIA 0 A 0, and the direction of the periodic orbit included in the outermost connected component is clockwise.
  • the word expression IIA 0 A 0 is the same as (c), but the direction of the periodic orbit is opposite.
  • each ss-saddle connection diagram an image of a ss-saddle connection diagram is shown which is formed by a conformal mapping that infinitely maps a point in the connected component that is the root.
  • this mapped image all ss-orbits included in the route are counterclockwise.
  • the connected component (the curve drawn with a broken line is directly above 1-source-sink point). Included area) becomes the root. If the same ss-saddle connection diagram has the same shape, but the direction of the flow is reversed as shown in FIG. 11B, the route becomes a connected component immediately below 1-source-sink point from the definition. .
  • the ss-saddle connection diagram with the word representation IIA 0 A 0 in FIG. 11 (c) when the connected component directly above the 1-source-sink point reverses the direction (FIG. 11 (d)). The connected component directly below is the root.
  • the connected component of the root in the image is the outermost component. Note that the ss-orbit is always counterclockwise (see the lower part of each figure in FIG. 11).
  • the orientation edge between the vertices is defined as shown in FIG. 11 by using a conformal map obtained by mapping the root connected component to the outermost side.
  • FIG. 12 shows the visualization of the local structure of the ss-saddle connection diagram generated by the operations A 0 , A 2 , and C and the corresponding graph representation.
  • the parent connected component v is expressed as o ⁇ (o 0 or o 2 ).
  • (a) is A 0
  • (b) is A 0
  • (c) is A 0 A 0
  • (d) is A 2 C k (k ⁇ 0)
  • (e) is A 2 C l ( l ⁇ 1)
  • (f) are A 2 2 C k + 1 (k, l ⁇ 1)
  • (g) is A 0 A 2 C k (k ⁇ 0)
  • (h) is A 0 A 2 C l
  • the local ss-saddle connection diagram constructed by l ⁇ 1) and the parent-child relationship of the connected components generated corresponding thereto, and the counterclockwise ss-orbit included in the parent connected component v are indicated by broken lines. .
  • the label of the parent connected component is expressed as o * (where * is either ⁇ , 0, or 2). Since the route is determined so that the direction of the ss-orbit is always counterclockwise, the direction in which the trajectory included in the connected component as a child flows is automatically determined.
  • the parent connected component v has two child connected components w ⁇ Vo 0 and y ⁇ V + 0 , so from o * to draw the two arrows to the o 0 and + 0.
  • the depending on choosing to which connecting component the parent, o 0 or - can be an arrow pointing to 0.
  • two types of patterns are created as shown in FIG. 12 (c), so the three arrows to o 0 , + 0 , and ⁇ 0 are changed accordingly. or o 0 and draw - two of the edge of the 0 can be drawn.
  • the l pieces of arrow - can be drawn to 2.
  • FIGS. 13-A to 13-D are flowcharts for explaining the conversion process of the ss-saddle connection diagram in the I and II series into a tree.
  • the conversion process of the ss-saddle connection diagram in the I and II series shown in FIGS. 13-A to 13-D into a tree can be executed by a device such as a computer.
  • step S143 if root 0 does not have the shape of FIG. 12A (“No” in step S143), the process proceeds to step S145.
  • step S145 it is determined whether or not root 0 has the left shape of FIG. If root 0 has the left shape of FIG. 12B (“Yes” in step S145), the vertex id of w is set to 1, the tree G is the label of the tree of FIG. The tree replaced with the set of id and label (that is, G is defined as ((0, o ⁇ ) ⁇ (1, o 0 )), and ⁇ 1 ⁇ is added to T (step S146). The process proceeds to step S149.
  • step S145 if root 0 does not have the left shape of FIG. 12B (“No” in step S145), the process proceeds to step S148 (in this case, as shown in step S147, root). 0 has the shape of FIG. 12 (d)).
  • step S148 w, z 1 ,. . . , Z k vertices id are set to 1, 00,. . . , 0k ⁇ 1, and the tree G is defined as a tree in which the label of the tree in FIG. 12D is replaced with a set of id and label, and ⁇ 1,0,. . . , 0k ⁇ 1 ⁇ is added, and the process proceeds to step S149.
  • step S149 it is determined whether or not an element t larger than u exists in T. If u is the largest element in T (“No” in step S149), all vertex ids of the tree G are removed, the remaining tree is set as G (step S150), and the tree G is output (step S150). S151), the process ends.
  • step S153 it is determined whether u is a natural number (that is, corresponding to the ids of o ⁇ , o 0 , o 2 ). If u is not a natural number (“No” in step S153), the process proceeds to S114 in FIG. 10-B, which is an O-word process.
  • step S153 if u is a natural number in step S153 (“Yes” in step S153), the process proceeds to step S155.
  • step S155 the outer boundary of the vertex u is composed of an outward homoclinic ss-saddle connection and a circle S (see FIG. 6), and the inner boundary of u has the shape of FIG. 12 (a). Determine whether or not.
  • the outer boundary of the vertex u is composed of one homoclinic ss-saddle connection and a circle S, and the inner boundary of u has the shape shown in FIG. 12A (“Yes” in step S155)
  • w, y Vertices id is defined as u + 1, u0, and the tree G is replaced with a set of id and label in FIG.
  • step S155 if the outer boundary of the vertex u does not consist of an outward homoclinic ss-saddle connection and a circle S, or if the inner boundary of u does not have the shape of FIG. (“No” in S155), the process proceeds to step S157.
  • step S157 it is determined whether or not the outer boundary of the vertex u is composed of the outward homoclinic ss-saddle connection and the circle S, and the inner boundary of u has the left shape of FIG.
  • the outer boundary of the vertex u is composed of one homoclinic ss-saddle connection and a circle S, and the inner boundary of u has the left shape of FIG. 12B (“Yes” in step S157)
  • the tree G is determined as a tree in which the label of the left tree in FIG. 12B is replaced with a set of id and label, and ⁇ u + 1 ⁇ is added to T (step S158). Thereafter, the process returns to step S149.
  • step S157 if the outer boundary of the vertex u is not composed of the outward homoclinic ss-saddle connection and the circle S, or the inner boundary of u does not have the left shape of FIG. (“No” in step S157), the process proceeds to step S159.
  • step S159 it is determined whether or not the outer boundary of the vertex u is composed of an outward homoclinic ss-saddle connection and a circle S, and the inner boundary of u has the shape shown in FIG.
  • the outer boundary of the vertex u is composed of one homoclinic ss-saddle connection and a circle S, and u has the shape of FIG. 12D (“Yes” in step S159)
  • the tree G is defined as a tree in which the label of the tree in FIG. 12D is replaced with a set of id and label, and ⁇ u + 1, u0,. . . , Uk-1 ⁇ is added (step S160). Thereafter, the process returns to step S149.
  • step S159 if the outer boundary of the vertex u does not consist of the outward homoclinic ss-saddle connection and the circle S in step S159, or if the inner boundary of u does not have the shape of FIG. ("No" in S159), the process proceeds to step S161.
  • step S161 it is determined whether or not the vertex u has the left shape of FIG. If the vertex u has the left shape of FIG. 12C (“Yes” in step S161), the vertex ids of w, y 1 , y 2 are defined as u + 1, u0, u1, and the tree G is A tree in which the label of the left tree in FIG. 12C is replaced with a set of id and label is determined, and ⁇ u + 1, u0, u1 ⁇ is added to T (step S162). Thereafter, the process returns to step S149.
  • step S161 when the vertex u does not have the left shape in FIG. 12C (“No” in step S161), the process proceeds to step S163.
  • step S163 it is determined whether or not the vertex u has the right shape in FIG.
  • the vertex ids of w 1 and w 2 are determined as u + 1, u0, and the tree G is changed to FIG. )
  • the process proceeds to step S149.
  • step S163 when the vertex u does not have the right shape of FIG. 12C (“No” in step S163), the process proceeds to step S165.
  • step S165 it is determined whether or not the vertex u has the shape shown in FIG.
  • the vertex u has the shape shown in FIG. 12G (“Yes” in step S165)
  • the tree G is defined as a tree in which the label of the tree in FIG. 12G is replaced with a set of id and label, and ⁇ u + 1, u0,. . . , Uk ⁇ is added (step S166). Thereafter, the process returns to step S149.
  • step S165 when the vertex u does not have the shape of FIG. 12G (“No” in step S165), the process proceeds to step S167.
  • step S167 it is determined whether or not the vertex u has the left shape in FIG.
  • the vertex id of y l u + 1, u0, u1 ,. . . , Ul, and the tree G is defined as a tree in which the label of the left tree in FIG. 12H is replaced with a set of id and label, and ⁇ u + 1, u0, u1,. . . , Ul ⁇ is added (step S168).
  • the process returns to step S149.
  • step S167 when the vertex u does not have the left shape of FIG. 12H (“No” in step S167), the process proceeds to step S169.
  • step S169 it is determined whether or not the vertex u has the right shape in FIG.
  • the vertex u has the right shape in FIG. 12H (“Yes” in step S169)
  • the tree G is defined as a tree in which the label of the right tree in FIG. 12H is replaced with a set of id and label, and T is set to ⁇ u + 1, u0,. . . , Ul-1 ⁇ is added (step S170). Thereafter, the process returns to step S149.
  • step S169 when the vertex u does not have the right shape in FIG. 12H (“No” in step S169), the process proceeds to step S171.
  • step S171 it is determined whether or not the vertex u has the shape shown in FIG.
  • the vertex u has the shape shown in FIG. 12F (“Yes” in step S171)
  • the vertex ids are u + 1, u0,. . . , Uk-1, uk,. . . , Ul + k ⁇ 1
  • the tree G is defined as a tree in which the label of the tree in FIG. 12F is replaced with a set of id and label, and T is set to ⁇ u + 1, u0,. . . , Ul + k-1 ⁇ is added (step S172). Thereafter, the process returns to step S149.
  • step S171 when the vertex u does not have the shape of FIG. 12F (“No” in step S171), the process proceeds to step S173.
  • step S173 it is determined whether or not the vertex u has the right shape of FIG.
  • the vertex id of y is defined as u0
  • the tree G is defined as the right tree in FIG. 12B.
  • the tree is determined by replacing the label with a set of id and label, and ⁇ u0 ⁇ is added to T (step S174). Thereafter, the process returns to step S149.
  • step S173 when the vertex u does not have the right shape of FIG. 12B (“No” in step S173), the process proceeds to step S175.
  • step S175 it is determined whether or not the vertex u has the shape shown in FIG.
  • the tree G is defined as a tree in which the label of the tree in FIG. 12E is replaced with a set of id and label, and ⁇ u0,. . . , Ul-1 ⁇ is added (step 176). Thereafter, the process returns to step S149.
  • step S175 when the vertex u does not have the shape of FIG. 12E in step S175 (“No” in step S175), the process returns to step S149 (in this case, as shown in step S177, the vertex u is (It has the shape of FIG. 12 (a)).
  • the above operation is the conversion process from ss-saddle connection diagram in I / II series to tree.
  • Proposition 3.1 A streamlined phase structure of a stable Hamilton vector field with each 1-source-sink point corresponds to a tree expression with a root, a label, and an orientation.
  • the route selection method is not unique. Because, for any connected component that includes a circular boundary that does not have ⁇ -saddle inside, there is a continuous mapping, so that the connected component can always be the outermost connected component. is there. Therefore, for the saddle connection diagram, there are different graph representations as many as the number of connected components. In order to remove such ambiguity, the graph representations that can be transferred by these continuous maps are different.
  • V 1 and V 2 For the two structural stability Hamiltonian vector field V 1 and V 2, V a 1 ⁇ V 2 becomes equivalence relation, V 1 exist homeomorphism on certain homeomorphism D z (M) and V 2 Are defined by the fact that the outer boundaries of V 1 and V 2 correspond to each other without changing their orientation.
  • defines an equivalence relation, and a tree with a unique root, a label, and a direction can be given to this equivalence class.
  • Proposition 3.2 Equivalent class streamline phase structure of ⁇ in a structure-stable Hamiltonian vector field corresponds to tree representation with unique root, labeled, and orientation.
  • FIG. 14-A shows visualization of a graph expression of a structure-stable streamline pattern given by the same word expression IA 0 C and its regular expression.
  • IA 0 C Graph representation of flow expressed by I-word and II-word
  • the connected component of the route is selected directly below or directly above 1-source-sink point. Yes. As you can see, all regular expressions are different.
  • FIG. 14-A shows a graph representation of a structurally stable streamline pattern having the word representation IA 0 C and its regular expression.
  • the route corresponds to the connected component immediately below 1-source-sink point from its definition.
  • the root connected component is copied to the outermost connected component on the right side.
  • FIG. 14B is a flowchart for explaining a process of converting a tree into a regular expression.
  • the conversion process of a tree into a regular expression shown in FIG. 14-B can be executed by a device such as a computer.
  • id is given to all points of the tree as preprocessing. However, ids of different points are different.
  • input is set as tree G (step S178).
  • V is set as the vertex set of tree G
  • X is set as empty set ⁇
  • regular expression N (s, ⁇ ⁇ ) () (step S179).
  • step S180 the result of reading the children of s from the left is expressed as (s (1) , ⁇ (1) ),. . . , (S (h) , ⁇ (h) ).
  • (s, ⁇ * ) () in the regular expression N is changed to ⁇ * (s (1) , ⁇ (1) ) (),. . . , (S (h) , ⁇ (h) ) ()), and replace T with ⁇ s (1) ,. . . , S (h) ⁇ and ⁇ s ⁇ are added to X.
  • the above operation is the process of converting the tree into a regular expression.
  • an allowable regular expression refers to all regular expressions that can be obtained by applying forward replacement of a regular expression to a regular expression of a fundamental type among regular expressions corresponding to trees with roots, labels, and directions. .
  • fundamental type is, + ⁇ , - ⁇ , is a regular expression that contains the o ⁇ .
  • FIG. 15 shows five fundamental type TN and a regular expression N corresponding to the five fundamental type TN .
  • These regular expressions correspond to a (ss-) saddle connection diagram including a route in the graph expression. That, + ⁇ (+ 0, + 0) and, o ⁇ (o 2, + 2, ...., + 2) , such as regular expression of ss-saddle connection diagram with a word representation consisting OB 0 and IC k It is a regular expression.
  • FIG. 16 is a list of 14 forward replacements.
  • the forward permutation A 0+ is given as o 0 ⁇ o 0 (o 0 , + 0 ).
  • a connected component that does not have a streamline phase structure therein will be referred to as a simple connected component.
  • Simple connected components correspond to vertices that do not have child connected components in the graph representation.
  • forward replacement shown in FIG. 16 is performed when the operations A 0 , A 2 , B 0 , B 2 , and C are applied to ss-orbit and closed orbit in the simple connected component. It describes how the streamline structure and its regular expression change.
  • FIG. 17 shows the change of forward replacement and streamline structure when the operations A 0 , A 2 , and C are performed. More specifically, FIG. 17 shows changes in the local structure of streamlines obtained by applying the operations A 0 , A 2 , and C to the ss-orbit in the simple connected component o 0 or o 2 in the ss-saddle connection diagram. Is a forward replacement of a regular expression corresponding to.
  • forward replacement as shown in FIG. 18 can be defined when the operations B 0 , B 2 , and C are applied to the closed orbit included in the simple connected component.
  • FIG. 18 shows excitation along with the change in streamline structure obtained when operations B 0 , B 2 , and C are applied to the closed orbit included in the simple connected components ⁇ 0 and ⁇ 2 in the (ss-) saddle connection diagram. It is a figure for demonstrating forward replacement of the regular expression to do.
  • a set of regular expressions obtained by applying forward replacement to this fundamental type recursively is defined as follows.
  • the regular expression is O-type (or I, II-type) when the regular expression is composed of the above rules from the basic form O, OB 0 , OB 2 C k (or IC k or II). Say. Since all local streamline structures obtained by the operation are covered in FIGS. 17 and 18, the following can be understood from the configuration method.
  • the opposite of () is an operation of replacing a homoclinic saddle connection pair (or a circle boundary having exactly two ⁇ -saddles) with a circle boundary.
  • the reverse operation of the operation C is to delete a simple connected component surrounded by a ⁇ -saddle connection that is attached to a circle boundary having two or more ⁇ -saddles.
  • the backward replacement can be expressed as a change in the arrangement of character strings corresponding to the structure to be deleted of the allowable regular expression.
  • FIG. 19 is a diagram illustrating an example of a backward replacement of a regular expression corresponding to a change in the local streamline structure in the (ss-) saddle connection diagram.
  • the connected component ⁇ 0 is an outward homoclinic saddle connection and is surrounded as shown in FIG. 19 (a)
  • the simple connected component surrounded by ⁇ -saddle connection is deleted from the (ss-) saddle connection diagram.
  • backing replacement ⁇ 0 [[ ⁇ 2 ]]
  • ⁇ ⁇ 0 is obtained.
  • the parent connected component ⁇ 0 is surrounded by the inward homoclinic saddle connection (FIG. 19B)
  • the backward replacement of the forward replacement C 1 B 2 is given by one of the followings according to the sign of ⁇ 2 from the rule (1) of the order relation with respect to the label.
  • the backward replacement is ⁇ 0 ([[ ⁇ 0 , ⁇ 0 ]]) ⁇ ⁇ 0 .
  • FIG. 19 (d) As backward replacement is shown in FIG. 19 (d) from order relation rules (1) when the connecting component sigma 0 is surrounded inwardly homoclinic saddle connection, + 0 ([ [+ 0, + 0]] , ⁇ 0 ) ⁇ + 0 ( ⁇ 0 ) or ⁇ 0 ( +0 , [[ ⁇ 0 , ⁇ 0 ]]) ⁇ ⁇ 0 ( +0 ).
  • forward replacement B 0- is compared with three back replacements ⁇ 0 ([[ ⁇ 0 ( ⁇ 0 )])) ⁇ ⁇ 0 , + 0 ([[ +0 ( ⁇ 0 )]], ⁇ 0 ) ⁇ + 0 ( ⁇ 0 ), ⁇ 0 ( +0 , [[ ⁇ 0 ( +0 )]]) ⁇ ⁇ 0 ( +0 ).
  • Forward replacement C l B 0+ in Figure 18 (d) and FIG. 18 (e) induces the following backward replacement.
  • Backward replacement corresponding to deleting a simple connected component surrounded by ⁇ -saddle connection is fundamental type IC k (k ⁇ 0) and forward replacement A 2 C k , B 2 C k (k ⁇ 0), C l A 0+ , C l A 0 ⁇ , C l B 0+ , C l B 0 ⁇ (l ⁇ 0), C O l , C ⁇ l , C l A 2 C k , (l + k ⁇ 0), C l B 2 This corresponds to the reverse operation of C k .
  • the backward replacement is a sequence (+ 2 ,..., + 2 ) or ( ⁇ 2 ) including simple connected components surrounded by ⁇ -saddle connections attached to the circle boundary. , ...,- 2 ) represents that one of the simple connected components is deleted.
  • Backward replacement of Fundamental type IC k (k ⁇ 1) is as follows.
  • the connected component o 0 is outward or inward homoclinic saddle connection as shown in FIGS. 19E and 19F.
  • backward replacement obtained by erasing the ⁇ -saddle connection surrounded by a concatenation component + 2 is as follows.
  • forward replacement C O l , C ⁇ l , C l A 0+ , C l A 0- , C l B 0+ , C l B 0- (l ⁇ 0) are all defined for the connected component o 2 or ⁇ 2 surrounded by ⁇ -saddle connection, the back replacement for it is naturally defined as it can.
  • FIG. 20 shows a list of back replacement [B 2 ], [B 0 ] and [C] of regular expressions derived by the reverse operation of the operations B 0 , B 2 and C.
  • the word expression of the regular expression finally obtained when the backward replacement in FIG. 20 is applied as much as possible to the last given allowable regular expression is as follows. This has proven to be mathematically correct.
  • Lemma 4.2 When the given allowable regular expression N can no longer be applied to all the backward replacements included in [B 0 ], [B 2 ], [C], the word expression for the regular expression is one of the following: It is.
  • N O 2
  • the word representation of N is either IA 0 p + m A 2 c-2-p-m or IIA 0 c-2 .
  • c, p, m is, o * contained in N, + 0, - is either 0.
  • a partial word string CB 0 B 2 is obtained according to the extracted structure.
  • FIGS. 22A to 22C are flowcharts for explaining the process of converting the regular expression into the word expression.
  • the conversion processing of regular expressions into word expressions shown in FIGS. 22A to 22C can be executed by an apparatus such as a computer.
  • the word W is set as an empty set ⁇ (step S185).
  • step S186 it is determined whether sigma containing innermost of ⁇ 2 ⁇ ( ⁇ 2) or innermost of sigma 2 sigma 0 containing (sigma 2) are included in the regular expression N. If innermost of sigma 2 to comprise sigma phi (sigma 2) or innermost of sigma 2 sigma 0 containing (sigma 2) are included in the regular expression N is ( "Yes" in step S186), the process proceeds to step S187.
  • step S187 the innermost ⁇ 2 included in ⁇ ⁇ ( ⁇ 2 ) or ⁇ 0 ( ⁇ 2 ) is extracted from the regular expression N, W is replaced with B 2 W (step S188), and the process returns to step S186.
  • step S186 if the sigma containing innermost of ⁇ 2 ⁇ ( ⁇ 2) or innermost of sigma 2 sigma 0 containing (sigma 2) is not included in the regular expression N (step S186 "No"), The process proceeds to step S189.
  • step S189 + 2 including innermost of + 2 in the regular expression N (+ 2, - 2, ..., - 2) determines whether there.
  • step S189 + 2 including innermost of + 2 in the regular expression N (+ 2, - 2, ..., - 2) If there is no ( "No" in step S189), step S191 Migrate to
  • step S191 a regular expression of innermost in the N - containing 2 - 2 (+ 2, ..., + 2, - 2) determines whether there.
  • step S191 innermost of in a regular expression N - containing 2 - 2 (+ 2, ..., + 2, - 2) does not exist ( "No" in step S191), step S193 Migrate to
  • step S193 it is determined whether or not two consecutive innermost ⁇ 0 and ⁇ 0 are included in the regular expression N.
  • two consecutive innermost ⁇ 0 and ⁇ 0 are included in the regular expression N (“Yes” in step S193), such ⁇ 0 and ⁇ 0 are extracted from the regular expression N (step S194), and W is Replace with B 0 W (step S195), and return to step S193.
  • step S193 when two consecutive innermost ⁇ 0 and ⁇ 0 are not included in the regular expression N in step S193 (“No” in step S193), the process proceeds to step S196.
  • step S196 it is determined whether or not the innermost ⁇ 0 ( ⁇ 0 ) is included in the regular expression N.
  • the innermost ⁇ 0 ( ⁇ 0 ) is included in the regular expression N (“Yes” in step S196)
  • such ⁇ 0 ( ⁇ 0 ) is extracted from the regular expression N (step S197), and W is set to B 0. Replace with W (step S195), and proceed to step S193.
  • step S196 if the innermost ⁇ 0 ( ⁇ 0 ) is not included in the regular expression N (“No” in step S196), the process proceeds to step S198.
  • step S198 the regular expression N is added to two or more consecutive ⁇ 2 ,. . . , ⁇ 2 , whether or not the innermost ⁇ 2 is included is determined.
  • the regular expression N includes two or more consecutive ⁇ 2 ,. . . , ⁇ 2 includes an innermost ⁇ 2 (“Yes” in step S198), such ⁇ 2 is extracted from the regular expression N (step S199), and W is replaced with CW ( Step S200) and return to Step S198.
  • step S198 the regular expression N is added to two or more consecutive ⁇ 2 ,. . . , ⁇ 2 does not include the innermost ⁇ 2 (“No” in step S198), the process proceeds to step S201.
  • step S201 it is determined whether or not ⁇ 2 that is a child of o * exists in the regular expression N. If the o * child and innermost ⁇ 2 exists in the regular expression N (“Yes” in step S201), such ⁇ 2 is removed from the regular expression N (step S199), and W is replaced with CW. (Step S200), the process returns to Step S198.
  • step S201 if the o * innermost of sigma 2 child does not exist in the regular expression N ( "No" in step S201), and proceeds to step S202.
  • step S202 innermost a mu 2 in children of sigma 2 to determine whether there regular expression N. If ⁇ 2 and innermost ⁇ 2 are present in the regular expression N (“Yes” in step S202), such ⁇ 2 is extracted from the regular expression N (step S203), and W is replaced with CW ( Step S200) and return to Step S198.
  • step S202 if there is no innermost ⁇ 2 as a child of ⁇ 2 in step S202 (“No” in step S202), the process proceeds to step S204.
  • step S204 it is determined whether or not ⁇ 2 exists in the regular expression N. If ⁇ 2 exists in the regular expression N (“Yes” in step S204), the process returns to step S186.
  • step S204 if ⁇ 2 does not exist in the regular expression N in step S204 (“No” in step S204), the process proceeds to step S205.
  • step S205 it is determined whether or not ⁇ 0 ( ⁇ 0 ) exists in the regular expression N.
  • ⁇ 0 ( ⁇ 0 ) exists in the regular expression N (“Yes” in step S205)
  • the process returns to step S186.
  • step S205 if ⁇ 0 ( ⁇ 0 ) does not exist in the regular expression N in step S205 (“No” in step S205), the process proceeds to step S206.
  • step S206 it is determined whether or not there are two ⁇ 0 and ⁇ 0 consecutive in the regular expression N. If there are two consecutive ⁇ 0 and ⁇ 0 in the regular expression N (“Yes” in step S206), the process returns to step S186.
  • step S206 if there are no two ⁇ 0 and ⁇ 0 consecutive in the regular expression N in step S206 (“No” in step S206), the process proceeds to step S207.
  • step S207 the regular expression N is equal to or a sigma phi. If the regular expression N is sigma phi ( "Yes" in step S207), the W was replaced by OW (step S208), the word W is output (step S212), the process ends.
  • step S207 if the regular expression N is not sigma phi ( "No" in step S207), and proceeds to step S209.
  • step S209 it is determined whether the regular expression N comprises o 2.
  • the regular expression N includes o 2 (“Yes” in step S209), W is replaced with IA 0 p + m A 0 c-2-p ⁇ m W (step S210), and the word W is output (step S210). S212), the process ends.
  • c, p, m is o * contained in N, + 0, - the number of 0.
  • step S208 if the regular expression N does not contain o 2 ( "No" in step S208), W is substituted with IIA 0 c-2 W a (step S211), the word W is output (step S212), the process ends.
  • the above processing is conversion processing of regular expressions into word expressions.
  • FIG. 23 is a block diagram illustrating an example of the regular expression creation device 100 to which the present exemplary embodiment is applied, and conceptually illustrates only the portion related to the present exemplary embodiment of the configuration.
  • the regular expression creation device 100 in this embodiment schematically includes at least a control unit 102 and a storage unit 106.
  • the regular expression creation device 100 further includes an input / output control interface unit 108 and a communication control interface.
  • the unit 104 is provided.
  • the control unit 102 is a CPU or the like that comprehensively controls the entire regular expression creating apparatus 100.
  • the communication control interface unit 104 is an interface connected to a communication device (not shown) such as a router connected to a communication line or the like, and the input / output control interface unit 108 is connected to the input device 112 or the output device 114.
  • the storage unit 106 also stores various data (for example, FIG. 4 to FIG. 9, FIG. 11 to FIG. 12, FIG.
  • each unit of the regular expression creating apparatus 100 is connected to be communicable via an arbitrary communication path. Further, the regular expression creation device 100 is communicably connected to the network 130 via a communication device such as a router and a wired or wireless communication line such as a dedicated line.
  • the various data, database, table (simulation result file 106a, streamline diagram file 106b, regular expression file 106c, etc.) stored in the storage unit 106 are storage means such as a fixed disk device.
  • the storage unit 106 stores various programs, tables, files, databases, web pages, and the like used for various processes.
  • the simulation result file 106a is a simulation result storage unit that stores data indicating a simulation result mathematically simulated by the simulation unit 102a.
  • the simulation result file 106a includes values of design variables indicating the shape of a structure, and hydrodynamic simulation results of a predetermined fluid (such as an ocean current or an air current) with respect to the structure (fluid pressure and flow at each spatial coordinate). Data indicating an orientation or the like).
  • the simulation result file 106a may store, as simulation results, data input in advance via the input device 112 through model measurement in a laboratory such as a wind tunnel experiment.
  • the streamline diagram file 106b is streamline data storage means for storing data indicating streamlines such as a streamline diagram.
  • the streamline data stored in the streamline diagram file 106b may be data indicating streamlines analyzed by the streamline analysis unit 102b based on data indicating simulation results.
  • the regular expression file 106c is regular expression storage means for storing regular expression data.
  • the regular expression data stored in the regular expression file 106c is a character string or the like.
  • the input / output control interface unit 108 controls the input device 112 and the output device 114.
  • the output device 114 in addition to a monitor (including a home television), a speaker can be used (hereinafter, the output device 114 may be described as a monitor).
  • the input device 112 a keyboard, a mouse, a microphone, and the like can be used.
  • control unit 102 has an internal memory for storing a control program such as an OS (Operating System), a program defining various processing procedures, and necessary data. And the control part 102 performs the information processing for performing various processes by these programs.
  • the control unit 102 includes a simulation unit 102a, a streamline analysis unit 102b, a graph expression creation unit 102c, a regular expression creation unit 102d, and a word expression conversion unit 102e in terms of functional concept.
  • the simulation unit 102a is a simulation means for simulating a fluid with respect to a structure.
  • the above-described design parameter selection step (see FIG. 24), the second design parameter selection step (see FIG. 29, etc.) Execute.
  • the simulation unit 102a is not limited to a simulation in a two-dimensional plane, and may perform a fluid simulation in a three-dimensional space.
  • the simulation part 102a may optimize a structure using a well-known optimization method.
  • the simulation unit 102a performs a fluid simulation on the structure while repeatedly changing the design variable for determining the shape of the structure using an annealing method, a genetic algorithm method, or the like, thereby obtaining an appropriate shape of the structure.
  • the simulation unit 102a stores data indicating the simulation result in the simulation result file 106a.
  • the simulation unit 102a determines the value of a design variable indicating the shape of the structure, and the hydrodynamic simulation result of a predetermined fluid (such as an ocean current or an air current) with respect to the structure (the pressure or flow direction of the fluid in each spatial coordinate). Or data indicating resistance) may be stored.
  • the streamline analysis unit 102b is streamline analysis means for performing streamline analysis.
  • the streamline analysis unit 102b may derive a streamline diagram by performing streamline analysis on the simulation result of the simulation unit 102a.
  • the streamline analysis unit 102b creates a streamline diagram from a numerical simulation or experiment data stored in the simulation result file 106a using a known method. Specifically, the streamline analysis unit 102b calculates all the saddle points, 1-source-sink points, etc. from the numerical simulation results, and then calculates the contours of the stream function having the same value as the stream function value at that point.
  • a flow diagram can be created by drawing all the contours of the flow function having the same value as the value of the flow function on the boundary.
  • the streamline analysis unit 102b may perform streamline analysis after converting it to two-dimensional data in a cross section of the structure.
  • the plane used as a cross section is arbitrary, but preferably, the streamline analysis unit 102b may convert it into two-dimensional data in a cross section along the direction of the fluid flow (uniform flow). For example, in vehicles such as trains, automobiles, and aircraft, a cross section may be generated along the traveling direction.
  • the streamline analysis part 102b may extract the characteristic structure which satisfy
  • the streamline analysis unit 102b stores the created streamline diagram data in the streamline diagram file 106b.
  • the graph expression creating unit 102c is a graph expression creating means for executing the above-described graph expression creating step (step S1). For example, the sad connection in the O series shown in FIGS. 10-A to 10-E. By executing the conversion process of the diagram to the tree and the conversion process of the ss-saddle connection diagram in the I and II series shown in FIGS. 13-A to 13-D, etc. (Where N is an integer greater than or equal to 1) A graph representation corresponding to the flow pattern in the multi-connected external region having a number of holes is created in a one-to-one correspondence.
  • the child vertex set ⁇ (v) When the child vertex set ⁇ (v) is used, the child vertex set ⁇ (v) of v is rearranged according to a rule of a predetermined order relation, and from w ( ⁇ ) to l (w) for w ⁇ (v). You may include what was drawn by arranging the arrows from left to right.
  • the flow pattern includes (1) two ss- ⁇ -saddle connections, which have a suction-out flow pair of two types of flow patterns that can be taken topologically in a single connected outer region having one hole.
  • the flow pattern is converted so that the region containing the counterclockwise ss-orbit closest to the suction-and-out pair is the outermost region, and the converted flow
  • the trajectory that appears in the (ss-) saddle connection diagram of the pattern is extracted from the entire region, and the vertex is set to the connected component obtained by excluding all the trajectories that appear in the (ss-) saddle connection diagram from the entire region.
  • the connected component in is set as the root
  • the current component is set as the root
  • the connected components that are in contact with the current component are the children of the current component
  • labels are assigned according to the trajectory corresponding to the boundary, and the labels are assigned in a predetermined order.
  • the flow pattern starts from one or more of pattern words, and is one of operation words defining five types of operations that can be taken topologically when adding one hole to the flow pattern. May be a flow pattern diagram created by repeating the operation to give N until the number of holes reaches N.
  • the regular expression creating unit 102d is a regular expression creating means for realizing the above-described regular expression creating step (step S2).
  • the regular expression creating unit 102d performs, for example, a tree conversion process shown in FIG.
  • a regular expression is created from the graph expression created by the graph expression creation unit 102c.
  • the word expression conversion unit 102e is a word expression creating means for realizing the above-described word expression conversion step (step S3).
  • the word expression conversion unit 102e converts the word expressions into the regular expressions shown in FIGS. 22A to 22C.
  • the regular expression created by the regular expression creation unit 102d is converted into a word expression.
  • the word expression represents a pattern that does not have a suction-and-out pair in a double connected outer region having two holes. Any of the operation words that define five types of operations that can be taken topologically when adding one hole to the flow pattern for the added pattern words that define three types of flow patterns
  • a single word may be a symbol word formed by adding the number of holes added.
  • the regular expression creation device 100 may be connected to the external system 120 via the network 130.
  • the communication control interface unit 104 performs communication control between the regular expression creation device 100 and the network 300 (or a communication device such as a router). That is, the communication control interface unit 104 has a function of communicating data with other terminals via a communication line.
  • the network 130 has a function of connecting the regular expression creation device 100 and the external system 200 to each other, such as the Internet.
  • the external system 120 is connected to the regular expression creation device 100 via the network 130, and word representation is performed on an external database related to various data such as simulation result data and streamline diagram data, and connected information processing devices.
  • the external system 120 may be configured as a WEB server, an ASP server, or the like.
  • the hardware configuration of the external system 120 may be configured by an information processing device such as a commercially available workstation or personal computer, and an accessory device thereof.
  • Each function of the external system 120 is realized by a CPU, a disk device, a memory device, an input device, an output device, a communication control device, and the like in the hardware configuration of the external system 120 and a program for controlling them.
  • design method is used to determine the optimum parameters for a design object through the shape and arrangement of objects (hereinafter referred to as design objects) such as bridges and piers, and the control of the flow around them.
  • design objects such as bridges and piers
  • Assumption 1 There are design parameters that can be changed in the design object (shape, arrangement, flow control device, etc.).
  • Assumption 2 In designing a design object, an “optimal state” corresponding to the problem is set, and the optimal state can be described as a feature of the streamline structure. For example, when the state where the separation vortex is confined is set to the optimum state, it is expected that the wings maximize the lift, and the piers minimize the drag.
  • Assumption 3 In these designs, experiment and numerical calculation (and word expression / regular expression of streamline pattern obtained in it) can be performed by changing design parameters. Hereinafter, it is limited to a two-dimensional flow. In the case of 3D, it may be possible to design by taking a cross-section and making it 2D, etc. (Note that not all cases are possible in 3D).
  • an ideal situation is set as a word expression or regular expression of a streamline pattern from the beginning for such a situation, and a parameter area is searched so that the situation is realized. Will be able to. In this way, it becomes possible to search for as many and as quickly as possible design parameters “candidates” for realizing optimum settings. Starting with those candidates, performing a known optimization method increases the possibility that many and highly feasible design parameters can be obtained.
  • FIG. 24 is a flowchart for explaining a design method of an object in a fluid using word expressions and regular expressions, that is, the above-described design parameter selection process.
  • the flowchart shown in FIG. 24 is executed by the simulation unit 102a.
  • an upper limit and a lower limit of possible design parameters are determined (step S311).
  • Word expressions and regular expressions are assigned to the obtained results (step S314). It is determined whether or not the word expression has an IC partial word representing the optimum state or is a structural expression having a regular expression corresponding thereto (step S315). If the word expression has an IC partial word representing the optimum state or a structural expression having a regular expression corresponding to the IC partial word (“Yes” in step S315), the word expression is adopted as a design parameter candidate (step S316). On the other hand, if the word expression is not a structural expression having an IC partial word representing the optimum state or a corresponding regular expression (“No” in step S315), it is not adopted as a design parameter candidate.
  • step S317 it is determined whether there is a design parameter candidate (step S317). If there is a design candidate parameter (“Yes” in step S317), an existing optimization is performed on the design parameter candidate. A method is applied (step S319). In addition to applying the existing optimization method, it is ideal not only to quantitatively optimize the flow, but also to calculate and monitor word expressions and regular expressions at each stage of the optimization process. It is possible to optimize the streamline structure at the same time.
  • step S317 if there is no design parameter candidate (“No” in step S317), the parameter area is further divided (step S318), and the processes of steps S313 to S317 are executed.
  • FIG. 25 and FIG. 26 are diagrams for explaining an example of a specific example for explaining a method of designing an object in a fluid using the above-described word expression and regular expression.
  • FIG. 25 it is assumed that there are two vortex structures (vortex 1 and vortex 2) behind the object 200 placed in a uniform flow. It is assumed that the uniform flow speed (or the traveling speed of the object) is U, and the object 200 can rotate at the angular velocity G.
  • U can vary from 0 to 1.1
  • G can vary from ⁇ 1.6 ⁇ 2 ⁇ to 1.6 ⁇ 2 ⁇ (the upper and lower limits of the design parameters are determined in step S311). .
  • These parameter sections are divided into five equal parts to divide the design parameter area (step S312 above).
  • the results of executing steps S313 to S319 are (a) to (f) in FIG.
  • the division of these parameter areas is further repeated to find the best parameter area that maximizes the quantitative quantity such as lift, thereby confining the vortex and maximizing the lift. It is possible to know the speed of the object 200 and the rotation speeds.
  • FIG. 27 and 28 are explanatory diagrams for explaining a method for optimizing a flow around an object in a fluid using word expressions and regular expressions.
  • FIG. 29 is a flowchart for explaining a method for optimizing a flow around an object in a fluid using a word expression and a regular expression, that is, the above-described second design parameter selection step. The flowchart shown in FIG. 29 is executed by the simulation unit 102a.
  • FIG. 27 shows a case where the flat blade 300 is placed at a constant angle with respect to the uniform flow U.
  • FIG. 28 shows streamline patterns formed around the flat blade 300 in FIG. 27 and word expressions and regular expressions of the streamline patterns.
  • FIG. 28A shows a streamline pattern extracted from a flow created around a flat blade 300 placed in a uniform flow U (“flow” can be calculated by a fluid experiment or numerical calculation). (Referred to as “pattern A”).
  • FIG. 28C shows a pattern (referred to as “pattern B”) for confining the rotating flow structure on the flat blade 300.
  • the design parameters of the flat blade 300 are suitably changed to find design parameters that realize the pattern B (target pattern) from the pattern A (current pattern), and increase the lift applied to the flat blade 300. The method to make it explain is demonstrated.
  • step A the word expression (Japanese Patent Application No. 2012-203601 (PCT / JP2013 / 070939)) and the algorithm of the present embodiment are applied to pattern A and pattern B.
  • the pattern A includes a word expression IA 0 CCB 0 and a regular expression N A (o ⁇ (o 0 (o 2 , ⁇ 0 ( ⁇ 0 , ⁇ 0 ), + 2 , + 2 )).
  • the pattern B the word representation IA 0 CCB 0 and o ⁇ (o 2 (o 0 , + 0, - 2 (- 0, - 0, + 2))) to impart a regular expression N B of (Step S400, see FIGS. 28A and 28C). Since the word expressions of pattern A and pattern B are exactly the same, they cannot be distinguished in the prior art, but the regular expressions are different, so using these makes it possible to distinguish between the two flows. It becomes possible to specify the target shape of the optimum design.
  • step S401 it is determined whether or not the regular expression N A of the pattern A and the regular expression N B of the pattern B are the same (step S401).
  • the path P ⁇ is output (step S402). If the regular expression N A of the pattern A and the regular expression N B of the pattern B are different (“No” in step S401), the process proceeds to step S403.
  • step B a transition path between pattern A and pattern B is determined (step S403).
  • step C the design parameters are changed so as to make a transition along the route R i (steps S404 to S415).
  • the route is determined, and the changed design parameter that matches the regular expression of the pattern B is output as the optimum design parameter.
  • step S411 If i ⁇ N (“Yes” in step S411), the process returns to step S405 to search for the next route, while if i ⁇ N is not satisfied (“No” in step S411), all Since it is impossible to optimize the route, "impossible” is output (step S412).
  • the design parameter is changed from the pattern A in FIG. 28A to search for a word expression and a regular expression in FIG. If the pattern shown in FIG. 28B is not found even after searching all the design parameters that can be changed, it is determined that the optimum design along this route is impossible, and the pattern shown in FIG. ) Repeat the same for routes other than the shortest route.
  • step S415) If a design parameter that achieves the desired pattern B in FIG. 28C is suddenly found, this is the optimum design parameter, so this design parameter is output and the flow is terminated (step S415).
  • this design parameter is selected, a range of design parameters that can be changed around this design parameter is specified, and it is divided. Then, an experiment or numerical calculation is performed on each of them to search whether there is a pattern B in FIG. When the design parameter that becomes the pattern B in FIG. 28C is found, this is the optimum design parameter, so this design parameter is output and the flow is terminated (step S415).
  • step S412 If it is not found even after searching the range of design parameters that can be changed, it is determined that optimization on this route is impossible, and another route obtained in step B is searched. Even if all the routes acquired in step B are searched, if the design parameter that becomes the pattern B in FIG. 28C is not found, it is determined that the optimization is impossible (step S412).
  • the regular expression creation device 100 performs processing in a stand-alone form
  • the regular expression creation device 100 performs processing in response to a request from a client terminal, and the processing result is transmitted to the client terminal. You may make it return to.
  • all or a part of the processes described as being automatically performed can be manually performed, or all of the processes described as being manually performed can be performed.
  • a part can be automatically performed by a known method.
  • each illustrated component is functionally conceptual and does not necessarily need to be physically configured as illustrated.
  • each device of the regular expression creation device 100 in particular, the processing functions performed by the control unit 102, all or any part thereof are interpreted by a CPU (Central Processing Unit) and the CPU. It may be realized by a program to be executed, or may be realized as hardware by wired logic.
  • the program is recorded on a recording medium to be described later, and is mechanically read by the regular expression creating apparatus 100 as necessary. That is, in the storage unit 106 such as a ROM or HD, a computer program for performing various processes by giving instructions to the CPU in cooperation as an OS (Operating System) is recorded. This computer program is executed by being loaded into the RAM, and constitutes the control unit 102 in cooperation with the CPU.
  • OS Operating System
  • the computer program may be stored in an application program server connected to the regular expression creation device 100 via an arbitrary network 300, and may be downloaded in whole or in part as necessary. Is possible.
  • the program according to the present invention may be stored in a computer-readable recording medium, or may be configured as a program product.
  • the “recording medium” includes a memory card, USB memory, SD card, flexible disk, magneto-optical disk, ROM, EPROM, EEPROM, CD-ROM, MO, DVD, and Blu-ray (registered trademark). It includes any “portable physical medium” such as Disc.
  • program is a data processing method described in an arbitrary language or description method, and may be in any form such as source code or binary code. Note that the “program” is not necessarily limited to a single configuration, but is distributed in the form of a plurality of modules and libraries, or in cooperation with a separate program typified by an OS (Operating System). Including those that achieve the function. In addition, a well-known structure and procedure can be used about the specific structure for reading a recording medium in each apparatus shown in embodiment, a reading procedure, or the installation procedure after reading.
  • Various databases and the like (simulation result file 106a, streamline diagram file 106b, word expression file 106c, etc.) stored in the storage unit 106 are memory devices such as RAM and ROM, fixed disk devices such as hard disks, flexible disks, and the like.
  • the storage means such as an optical disk stores various programs, tables, databases, web page files, and the like used for various processes and website provision.
  • the regular expression creating apparatus 100 may be configured as an information processing apparatus such as a known personal computer or workstation, or may be configured by connecting any peripheral device to the information processing apparatus. Further, the word expression creation device 100 may be realized by installing software (including programs, data, and the like) that causes the information processing device to realize the method of the present invention.
  • the specific form of distribution / integration of the devices is not limited to that shown in the figure, and all or a part of them may be functional or physical in arbitrary units according to various additions or according to functional loads. Can be distributed and integrated. That is, the above-described embodiments may be arbitrarily combined and may be selectively implemented.
  • a regular expression for creating a regular expression of a flow pattern in a multi-connected external region having topologically N (where N is an integer of 1 or more) holes is a creation method for creating a graph expression corresponding to the flow pattern in a one-to-one correspondence, a regular expression creating step for creating a regular expression from the graph expression created in the graph expression creation step, Therefore, there is an effect that it is possible to provide a new expression method that can be associated with the flow pattern on a one-to-one basis.
  • V I a set of points called vertices
  • E is a set of edges connecting vertices
  • the graph representation is such that the parent vertex is v, the child vertex is w, the label assigned to the parent vertex v is l (v), and the label assigned to the child vertex w is Is the l (w), v child vertex set ⁇ (v), the child vertex set ⁇ (v) of v is rearranged according to a rule of a predetermined order relation, and l (v ) To l (w) are drawn side by side from left to right, so that it is possible to provide a graph expression in which the parent-child connection relationship can be easily recognized.
  • the flow pattern has (1) a suction / outflow pair of two types of flow patterns that can be taken topologically in a single connected outer region having one hole.
  • Pattern I having two ss- ⁇ -saddle connections
  • (2) having two pairs of flow patterns that can be taken topologically in the single connected outer region having one hole.
  • One saddle point a pattern II having a homoclinic saddle connection and two ss-saddle connections connecting the saddle points
  • one or the other of the pattern O having no suction spring-out pair in the double connected external region having two holes Since there were multiple, all the basics Re can be given a regular expression for a pattern, an effect that will handle the regular expression of a specific flow pattern.
  • the regular expression created in the regular expression creating step is provided with a word expression conversion step that converts the regular expression into a word representation, it is possible to convert the regular expression into a word representation. There is an effect.
  • the word expression is a double connection having two holes in addition to two kinds of flow patterns that can be topologically taken in a single connection outer region having one hole.
  • Any of the operation words that specify the topologically five types of operations for the pattern words that define a total of three types of flow patterns with the addition of a pattern that does not have a pair of suction and outflow in the external area Is a symbolic word formed by assigning the same number of holes as the number of added holes, so that the regular expression can be taken from all the basic flow patterns and topologically There is an effect that the operation can be converted into a defined word expression.
  • the regular expression creating method when designing a structure in a flow field, a flow pattern that can be easily handled without relying on experience or intuition for the flow pattern that can be taken for the structure.
  • the regular expression creating method, regular expression creating apparatus, and program can be provided.
  • word expressions and regular expressions are structural, such as bridge pier design, breakwater arrangement, port contamination removal, wind power blade design, train pantograph structure, optimal arrangement of oil fences, etc. It is extremely useful in various fields involving object design and arrangement. Also, it can be applied to fields such as sports mechanics such as structural design of sports equipment.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

 流れパターンと1対1に対応させることが可能な新たな表現方法を提供することが可能な流れパターンの正規表現作成方法、正規表現作成装置、および、プログラムを提供することを目的とする。位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成方法であって、前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、を備えている。

Description

流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム
 本発明は、流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラムに関する。
 従来、気流や水流などの流れ場における最適な構造物の設計を行うために、風洞実験のほか、大規模な数値計算による流体シミュレーションが用いられている。
 例えば、焼きなまし法(Simulated Annealing法)や遺伝的アルゴリズム法等を用いて構造物の設計変数を変化させながら、当該構造物に対する流体シミュレーションを繰り返し行う最適化技術が開発されている。
 また、近年、流体等の数理モデルを構築して、流れのパターンのトポロジーを数理的に扱うことのできるアルゴリズムやプログラムが開発されている。
 しかしながら、従来の構造物設計の最適化手法においては、繰り返し行われる大規模計算により計算時間と設計コストの増大化を招くので、それらの制約上、探索範囲を限定せざるを得ず、導き出される最適な構造物が局所最適なものである可能性を排除できないという問題点を有していた。すなわち、探索範囲をどこにするのかは、技術者の経験と直感に頼らざるを得ず、どこに探索範囲を設定したかによって、導き出される構造物の最適化結果が左右されるという問題があった。
 これに対して、本願出願人は、特許文献1の「流れパターンの語表現方法、語表現装置、および、プログラム(以下、「語表現理論」とも称する)」では、数値計算や実験で得られた流れに対して、この文字列を付与し、そこに現れる流れの定量的な特徴(例えば揚抗比)を定性的な部分文字列として表現する語表現理論を提案した。これにより、流れ場における構造物設計を行うにあたって、構造物に対して採り得る流れパターンを経験や直感に頼ることなく容易に扱うことが可能となった。
国際公開第2014/041917号
Tomasz Kaczynski, Konstantin Mischaikow, Marian Mrozek著 "Computational Homology" Springer,2000年
 しかしながら、流れの制御問題などを考えるときには,前記語表現理論は有用なものであるが,与えられたハミルトンベクトル場に対して固有な語表現を与えるアルゴリズムは存在するものの、逆に与えられた文字列に対しては複数の異なるハミルトンベクトル場が対応している。すなわち、語表現は流れの特定に使えるものの、語表現だけでは特定の流れを限定できない。
 このように、語表現のみによる流れの最適制御の理論は困難であることから、語表現と合わせて、別の新しい、流れパターンと1対1に対応する何らかの表現方法が必要とされている。
 本発明は、前記に鑑みてなされたもので、流れパターンと1対1に対応させることが可能な新たな表現方法を提供することが可能な流れパターンの正規表現作成方法、正規表現作成装置、および、プログラムを提供することを目的とする。
 上述した課題を解決するため、本発明は、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成方法であって、前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、を含むことを特徴とする。
 また、本実施の形態の好ましい形態によれば、前記グラフ表現は、前記流れパターンで規定される構造安定なハミルトンベクトル場Hに対して、固有のルート付き、ラベル付き、及び向き有りのツリーT=(V,E)を割り当て(但し、Vは頂点と呼ばれる点の集合、Eは、頂点の間を結ぶエッジの集合である)、平面グラフとして可視化したものであることが望ましい。
 また、本実施の形態の好ましい形態によれば、前記グラフ表現は、親の頂点をv、その子の頂点をw,親の頂点vに割り当てられたラベルをl(v)、子の頂点wに割り当てられたラベルをl(w),vの子頂点集合Γ(v)とした場合、vの子頂点集合Γ(v)を所定の順序関係のルールに従って並び替え、w∈Γ(v)について、l(v)からl(w)への矢印を左から右に並べて描画することを含むことが望ましい。
 また、本実施の形態の好ましい形態によれば、前記流れパターンは、(1)一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、(2)前記一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、(3)二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンO、の1又は複数であることが望ましい。
 また、本発明の好ましい態様によれば、前記グラフ作成工程は、前記流れパターンが吸い込み湧き出し対を持つ場合は、前記流れパターンを吸い込み湧き出し対に最も近い反時計回りのss-orbitを含む領域が一番外側領域になるように変換し、前記変換した流れパターンの(ss-)saddle connection diagramに現れる軌道を領域全体から抽出し、領域全体から(ss-)saddle connection diagramに現れる軌道をすべて除外して得られる連結成分に頂点を設定し、一番外側にある連結成分をルートとし、カレント成分をルートに設定し、前記カレント成分と互いに境界を接するような連結成分をカレント成分の子として、境界にあたる軌道に応じてラベルを割り当て、当該ラベルを所定の順序関係に従って並べ、前記カレント成分の子をカレント成分に設定して、子がなくなるまで繰り返すこと、を含むことが望ましい。
 また、本発明の好ましい態様によれば、前記流れパターンは、前記パターン語の1又は複数から開始し、流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を付与する操作を、穴の数がN個となるまで繰り返して作成された流れパターン図であることが望ましい。
 また、本実施の形態の好ましい形態によれば、さらに、前記正規表現作成工程で作成されたグラフ表現を語表現に変換する語表現変換工程と、を含むことが望ましい。
 また、本発明の好ましい態様によれば、前記語表現変換工程では、ルートつき,ラベル付き,及び方向つきのツリーに対応する正規表現のうちの基本タイプの正規表現に前方置換を施すことが可能な許容正規表現を前記語表現に変換することが望ましい。
 また、本実施の形態の好ましい形態によれば、前記語表現は、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語に対して、前記流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を、追加された穴の数だけ付与することにより形成された記号語であることが望ましい。
 また、本発明の好ましい態様によれば、流体中の物体に対して設計パラメータの候補を選択する場合に、前記設計パラメータの上限及び下限を設定し、当該設計パラメータの上限と下限で規定されるパラメータ領域から複数のパラメータを選択し、選択した複数のパラメータに対して、それぞれ流れの実験及び/又は数値計算を行い、実験及び/又は数値計算の結果に対して、前記語表現及び/又は前記正規表現を割り当て、割り当てた前記語表現及び/又は前記正規表現のうち、最適状態を示す前記語表現及び/又は前記正規表現を有する設計パラメータを、前記設計パラメータの候補として選択する設計パラメータ選択工程を含むことが望ましい。
 また、本発明の好ましい態様によれば、前記設計パラメータ選択工程は、さらに、前記選択した設計パラメータを使用して最適化設計を行う工程を含むことが望ましい。
 また、本発明の好ましい態様によれば、さらに、流体中の物体に対して最適な設計パラメータを選択する場合に、現在の流れパターンの前記正規表現から目的とする流れパターンの前記正規表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの前記正規表現となる設計パラメータを選択する第2の設計パラメータ選択工程を含むことが望ましい。
 また、本実施の好ましい態様によれば、前記第2の設計パラメータ選択工程では、前記現在の流れパターンの前記正規表現及び語表現から目的とする流れパターンの前記正規表現及び語表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの前記正規表現及び語表現となる設計パラメータを選択することが望ましい。
 また、本発明の好ましい態様によれば、前記第2の設計パラメータ選択工程は、前記現在の流れパターンの前記正規表現及び語表現から前記目的とする流れパターンの前記正規表現及び語表現へ遷移するための1又は複数の経路を取得し、前記取得した1又は複数の経路のうち、設計パラメータを変更して前記現在の流れパターンから遷移可能な経路を選択し、前記選択した経路に対して、前記変更した設計パラメータを中心にして設計パラメータを変更し、前記目的とする流れパターンの前記正規表現及び語表現となる設計パラメータを選択することが望ましい。
 また、上述した課題を解決するため、本発明は、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成装置であって、前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成手段と、前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成手段と、を備えたことを特徴とする。
 また、上述した課題を解決するため、本発明は、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成装置に搭載されるプログラムであって、前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、をコンピュータに実行させることを特徴とする。
 本発明によれば、流れパターンと1対1に対応させることが可能な新たな表現方法を提供することが可能な流れパターンの正規表現作成方法、正規表現作成装置、および、プログラムを提供することが可能となるという効果を奏する。
図1は、本発明の概略を説明するための説明図である。 図2-Aは、本実施の形態に係る流れパターンの正規表現作成方法を説明するためのフローチャートである。 図2-Bは、図2-Aのグラフ表現作成工程の一例を説明するためのフローチャートである。 図3は、本実施形態の概要に示すフローチャートである。 図4は、構造安定な流れのパターンを模式的に示した図である。 図5は、領域における構造安定な流れパターンの位相的分類を行う特徴的な軌道(流線)を全て記述した図である。 図6は、初期構造となる3種類の流れパターンを模式的に示す図である。 図7は、穴を一つ付け加えて構造安定な流れを構成する5種類の操作を模式的に示した図である。 図8は、二つの構造物と一様流がある場合の流れパターンの全分類を示す図である。 図9は、操作B,B,Cによって生成されるすべての流れパターンとそれに対応して連結部分集合(頂点集合)に励起される親子関係を示す図である。 図10-Aは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図10-Bは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図10-Cは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図10-Dは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図10-Eは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図11は、ルートとなる連結成分をC=D(M)\Dから選ぶ方法を説明するための図である。 図12は、操作A,A,Cによって生成されるss-saddle connection diagramの局所構造とそれに対応して得られるグラフ表現の可視化を示す図である。 図13-Aは、I,II系列におけるss-saddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図13-Bは、I,II系列におけるss-saddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図13-Cは、I,II系列におけるss-saddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図13-Dは、I,II系列におけるss-saddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。 図14-Aは、IACなる語表現をもつ構造安定な流線パターンのグラフ表現とその正規表現を示す図である。 図14-Bは、ツリーの正規表現への変換処理を説明するためのフローチャートである。 図15は、5つのfundamental type Tとそれに対応する正規表現Nを示す図である。 図16は、14個のforward replacementのリストを示す図である。 図17は、操作A,A,Cを施したときのforward replacementと流線構造の変化を表した図である。 図18は、(ss-)saddle connection diagramにおける単純連結成分σとσに含まれるclosed orbitに操作B,B,Cを適用したときに得られる流線構造の変化に伴って励起する正規表現のforward replacementを説明するための図である。 図19は、操作B,B,Cの逆操作により誘導される正規表現のbackward replacement [B],[B],[C]のリストを示す図である。 図20は、(ss-)saddle connection diagramにおける局所流線構造の変化に対応する正規表現のbackward replacementを示す図である。 図21は、正規表現oφ(o(o(o(-(-,-)))),+(+))を持つ構造安定なハミルトンベクトル場のss-saddle connection diagramを示す図である。 図22-Aは、正規表現の語表現への変換処理を説明するためのフローチャートである。 図22-Bは、正規表現の語表現への変換処理を説明するためのフローチャートである。 図22-Cは、正規表現の語表現への変換処理を説明するためのフローチャートである。 図23は、本実施の形態が適用される正規表現作成装置の一例を示すブロック図である。 図24は、語表現と正規表現を使用した流体中の物体の設計方法を説明するためのフローチャートである。 図25は、語表現と正規表現を使用した流体中の物体の設計方法を説明するための具体例の一例を説明するための図である。 図26は、語表現と正規表現を使用した流体中の物体の設計方法を説明するための具体例の一例を説明するための図である。 図27は、語表現と正規表現を使用した流体中の物体まわりの流れの最適化手法を説明するための説明図である。 図28は、語表現と正規表現を使用した流体中の物体まわりの流れの最適化手法を説明するための説明図である。 図29は、語表現と正規表現を使用した流体中の物体のまわりの流れの最適化方法の一例を説明するためのフローチャートである。
 以下に、本発明に係る構造安定な流れパターン(「流線パターン」とも称する)の正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラムの実施形態を図面に基づいて詳細に説明する。本発明の構成要素は、本明細書の図面に一般に示してあるが、様々な構成で広く多様に配置し設計してもよいことは容易に理解できる。したがって、本発明の装置、方法、及びプログラムの実施の形態についての以下のより詳細な説明は、特許請求の範囲に示す本発明の範囲を限定するものではなく、単に本発明の選択した実施形態の一例を示すものであって、特許請求の範囲に示す本発明と矛盾無く装置、方法、及びプログラムについての選択した実施形態を単に示すものである。当業者は、特定の細目の1つ以上が無くても、又は他の方法、部品、材料でも本発明を実現できることが理解できる。また、フローチャートの各ステップを実行する順番はフローチャートに示す順番に限定されるものではなく、本発明の趣旨を逸脱しない範囲で順番を変更してもよい。なお、この実施の形態によりこの発明が限定されるものではない。また、下記実施の形態における構成要素には、当業者が容易に想定できるもの又は実質的に同一のものが含まれる。本願出願人が既に出願した国際公開第2014/041917号(PCT/JP2013/070939号)及び特願2013-230678号の明細書及び図面の内容は援用により本明細書に組み込まれる。
 特に以下の実施形態においては、本発明を二次元流体のシミュレーションに適用した例について説明することがあるが、この場合に限られず、三次元流体のシミュレーションにおける任意の断面(構造物の断面等)についても、同様に本発明を適用することができる。
[本発明の概略]
 図1は、本発明の概略を説明するための説明図である。本願出願人が既に出願した特許文献1(国際公開第2014/041917号、特願2012-203601号、PCT/JP2013/070939号)の「流れパターンの語表現方法、語表現装置、および、プログラム」では、数値計算や実験で得られた流れパターンに対して,この文字列を付与し,そこに現れる流れの定量的な特徴(例えば揚抗比)を定性的な部分文字列として表現する語表現理論を提案した。
 また、本願出願人が既に出願した特願2013-230678(本願出願時未公開)の「流体遷移経路取得装置、流体遷移経路取得方法、及びプログラム」では、語表現から流体遷移経路を取得しているが、語表現を使って遷移の可能性をすべて網羅した上で,あとは文字列ではなく流れのパターンのパターンマッチングを行う必要があった。
 図1において、上述したように、特許文献1では、流れパターンから語表現を作成する(P1)方法について既に出願している。また、特願2013-230678において、語表現から流体遷移経路を作成する(P2)方法を出願している。語表現では、流れパターンに対して固有な語表現を与えることができるが、語表現で与えられた文字列に対しては複数の異なる流れパターンが対応している。そこで、本実施の形態の流れパターンの正規表現作成方法では、流れパターンに1対1に対応するグラフ表現を作成し(T1)、グラフ表現から正規表現を作成する(T2)。また、本実施の形態の流れパターンの正規表現作成方法では、正規表現から語表現を作成する(T3)。正規表現から作成される語表現は、流れパターンと1対1に対応したものとなるので、語表現から直接、流体遷移経路を作成することが可能となる。さらに、本実施の形態の流れパターンの正規表現作成方法では、正規表現から直接、流体遷移経路を作成することもできる(T4)。
[流れパターンの正規表現作成方法]
 図2-A~図22を参照し、本実施の形態に係る流れパターンの正規表現作成方法について説明する。本実施の形態に係る流れパターンの正規表現作成方法はコンピュータ等の装置により実行可能である。図2-Aは、本実施の形態に係る流れパターンの正規表現作成方法の概略を説明するためのフローチャートである。
 図2-Aに示すように、本実施の形態に係る流れパターンの正規表現作成方法は、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成ものであり、大別すると、流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程(ステップS1)と、ステップS1のグラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程(ステップS2)と、ステップS2の正規表現作成工程で作成された正規表現を語表現に変換する語表現変換工程(ステップS3)と、を備えている。
 グラフ表現は、流れパターンで規定される構造安定なハミルトンベクトル場Hに対して、固有のルート付き、ラベル付き、向き有りツリーT=(V,E)を割り当て(但し、Vは頂点と呼ばれる点の集合、Eは、頂点の間を結ぶエッジの集合である)、平面グラフとして可視化したものであってもよい。
 また、グラフ表現は、親の頂点をv、その子の頂点をw,親の頂点vに割り当てられたラベルをl(v)、子の頂点wに割り当てられたラベルをl(w),vの子頂点集合Γ(v)とした場合、vの子頂点集合Γ(v)を所定の順序関係のルールに従って並び替え、w∈Γ(v)について、l(v)からl(w)への矢印を左から右に並べて描画されたものを含むことにしてもよい。
 また、流れパターンは、(1)一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、(2)前記一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、(3)二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンO、からなるパターン語の1又は複数であることにしてもよい。
 また、語表現は、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語に対して、前記流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を、追加された穴の数だけ付与することにより形成された記号語であることとしてもよい。
 図2-Bは、図2-Aのグラフ表現作成工程の一例を説明するためのフローチャートである。図2-Bにおいて、グラフ表現作成工程では、まず、流れパターンが吸い込み湧き出し対を持つ場合は、流れパターンをその吸い込み湧き出し対に最も近い反時計回りのss-orbitを含む領域が一番外側領域になるように変換する(ステップS31)。流れパターンは、例えば、上述のパターン語の1又は複数から開始し、流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を付与する操作を、穴の数がN個となるまで繰り返して作成された穴の数がN個の流れパターン図としてもよい。具体的には、後述するように、I-word,II-wordのss-saddle connection diagram(流れパターン図)のルートが、図11に示すようになるように変換する。次に、変換した流れパターンの(ss-)saddle connection diagramに現れる軌道を領域全体から抽出する(ステップS32)。領域全体から(ss-)saddle connection diagramに現れる軌道をすべて除外して得られる連結成分に頂点を設定し、一番外側にある連結成分をルートとする(ステップS33)。続いて、カレント成分をルートに設定する(ステップS34)。カレント成分と互いに境界を接するような連結成分をカレント成分の子とし、境界にあたる軌道に応じてラベルを割り当て,それらを所定のラベルの順序関係に従って並べる(ステップS35)。カレント成分の子をカレント成分に設定して、子がなくなるまでステップS35を繰り返し実行する(ステップS36)。
 また、本実施の形態の流れパターンの正規表現作成方法は、流体中の物体に対して設計パラメータの候補を選択する設計パラメータ選択工程を実行してもよい。設計パラメータ選択工程では、例えば、流体中の物体に対して設計パラメータの候補を選択する場合に、前記設計パラメータの上限及び下限を設定し、当該設計パラメータの上限と下限で規定されるパラメータ領域から複数のパラメータを選択し、選択した複数のパラメータに対して、それぞれ流れの実験及び/又は数値計算を行い、実験及び/又は数値計算の結果に対して、前記語表現及び/又は前記正規表現を割り当て、割り当てた前記語表現及び/又は前記正規表現のうち、最適状態を示す前記語表現及び/又は前記正規表現を有する設計パラメータを、前記設計パラメータの候補として選択してもよい。
 また、設計パラメータ選択工程は、さらに、選択した設計パラメータを使用して最適化設計を行う工程を含むことにしてもよい。
 また、本実施の形態の流れパターンの正規表現作成方法は、流体中の物体に対して設計パラメータを選択する第2の設計パラメータ選択工程を実行してもよい。例えば、第2の設計パラメータ選択工程は、例えば、流体中の物体に対して最適な設計パラメータを選択する場合に、現在の流れパターンの正規表現から目的とする流れパターンの正規表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの正規表現となる設計パラメータを選択することにしてもよい。
 また、第2の設計パラメータ選択工程では、正規表現に加えて語表現を使用して、現在の流れパターンの正規表現及び語表現から目的とする流れパターンの正規表現及び語表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの正規表現及び語表現となる設計パラメータを選択することにしてもよい。
 また、第2の設計パラメータ選択工程では、現在の流れパターンの正規表現及び語表現から目的とする流れパターンの正規表現及び語表現へ遷移するための1又は複数の経路を取得し、取得した1又は複数の経路のうち、設計パラメータを変更して現在の流れパターンから遷移可能な経路を選択し、選択した経路に対して、変更した設計パラメータを中心にして設計パラメータを変更し、目的とする流れパターンの正規表現及び語表現となる設計パラメータを選択することにしてもよい。
 以下、本実施の形態の流れパターンの正規表現作成方法をその原理と共に、詳細に説明する。
(1.語表現理論)
 まず、本実施の形態の流れパターンの正規表現作成方法では、上述の語表現理論の一部を使用しているので、語表現理論について簡単に説明する。
(1-1.語表現アルゴリズム)
 ここで、上述した語表現を形成させる語変換アルゴリズムの概要について以下に説明する。ここで、図3は、語表現アルゴリズムの概要を示すフローチャートである。
 図3に示すように、まず、位相幾何学的にN個の穴を有する連結外部領域における流れパターンの語表現を形成させるために、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語(例えば、I,II,O)のうちのいずれか一語を付与する(ステップSA-1)。
 そして、本実施の形態は、ステップSA-1にて付与された語に対して、流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語(例えば、A,A,B,B,C)のうちのいずれか一語を付与し(ステップSA-2)、当該ステップSA-2の処理を穴の数がN個となるまで繰り返し行うことにより(ステップSA-3)、N個の穴を有する連結外部領域に対応する語表現を形成させる。なお、この例では、パターン語の付与(ステップSA-1)が先に行われ、操作語の付与(ステップSA-2)が後に行われるが、これに限られず、操作語の付与を先に行い、後にパターン語の付与を行ってもよいものである。
 ここで、本実施の形態において、「連結外部領域」とは、単連結外部領域と多重連結外部領域を含む表現であり、「単連結外部領域」とは、二次元(平面)の中にある領域で、一つの穴があいたものをいい、「多重連結外部領域」とは、二次元(平面)の中にある領域で、複数の穴があいたものをいう。ここでいう「穴」という表現は、数学的な抽象表現であるが、応用上はさまざまな表現形態であってもよい。例えば、乗り物が移動した際に相対的に発生する一様な流れに注目した場合、一様流に沿った乗り物の断面に構造物が一つ又は複数あった状況があれば、その領域は単連結外部領域または多重連結外部領域として扱える。例えば、空気中を進む航空機等の飛行物体や、空気中を進む電車や自動車等の地上走行車両、海上または海中を進む船等の乗り物についても、一または複数の穴と、単連結または多重連結の連結外部領域として扱える。換言すれば、「流れの中に一または複数の障害物がある」ような流れを、本実施形態では、連結外部領域における流れとして扱う。また、孤立した渦構造や周囲に周期軌道を持つような流れ構造(楕円型の停留点など)も「穴」とみなすことができる。
 上述のように、流れパターンに対して語表現を割り当てることによって、連結外部領域における流れの位相的分類が可能となる。位相的というのは数学の専門用語であり、トポロジー(位相幾何学)とも呼ばれる幾何学の一分野を指す。古典的な幾何学では、三角形や四角形はその角の数の違いにより異なる図形と見なすが、位相幾何学の立場ではそういった細かい情報を見ず、三角形と四角形を例えば輪ゴムを変形させて互いに移りあえるという立場にたって同じ図形と見なす。すなわちすべての多角形は、円と同じと見なす。一方で、円の中にもう一つ円がくりぬかれているような円環領域があったときに円と円環領域は一つの輪ゴムの変形では互いに変形できないので、異なる図形と見なす。連結外部領域では、空いている穴の数が異なれば、それらの図形は位相的に異なると見なす一方で、穴の形さえ同じであれば、その穴の形が円であろうが三角形であろうが線分であろうが同じものとみなす。このような理由から流れの領域を特徴づけるのは穴の数のみであるため、本実施形態では、穴の数M+1に対して連結外部領域をDζ(M)と表現する。例えば、穴が1つのみの場合は、単連結外部領域Dζ(0)であり、穴が2つ在る場合は、二重連結外部領域Dζ(1)となる。
 流れの位相的分類を扱うときは、流れを特徴づけるある特定の構造(「位相構造」と呼ぶ。)を捕まえて、その特定の構造をもった流れが二つあったときに、その双方を特定の構造の連続的な(すなわち切ったり貼ったりしない)変形によって互いに変形できないものは異なる流れと見なす。
 ここで、本実施の形態において扱う流れの構成要素について図4~図6を参照して説明する。本実施の形態における「流れ」として非圧縮流れを扱ってもよい。流体の非圧縮性とは、力を加えてもほとんどその体積を変えないような性質をいう。通常の水や大気の流れは、日常生活のスケールで考える場合は、こうした流れの枠組みで考えても概ね差し支えはない。なお、本発明は、これに限られず、圧縮性のある流れを計算上取り扱ってもよい。
 また、本実施の形態における流体は、非粘性流体として扱ってもよい。非粘性流体は、境界条件を「すべり条件」として扱うことができる。なお、粘性流体では、境界条件が「ゼロ境界条件」となるが、粘性流体でも、粘性が大きくないとき、境界付近は渦が層状の流れとなっており、その流れの外側境界を、仮想的に拡張することで「すべり条件」として扱うことができる。したがって、本実施形態では、特に、非粘性流体における語変換アルゴリズムにおいて説明することがあるが、適切な境界の拡張により、粘性流体についても本実施形態を適用することができる。
 本実施の形態で扱う流れの構成要素は、以下の3つである。一つは障害物、一つは渦(vortex)、一つは一様流(uniform flow)である。障害物は連結外部領域における穴のことであり、この形状は位相的分類を考えている上で円としておいても数学的に与える結果に影響はないことが一般的な流体の数学理論から導くことができる。ここで、図4は、流れのパターンを模式的に示した図である。なお、図4(c)は、4つの∂-saddle点(定義は後述)を持つ境界を模式的に示している。
 図4(a)に示すように、渦はそのまわりで回転する流れを作る要素である。一様流とは、川の流れでいえば基本的な流れのことであり、領域全体を横切るような流れである。また、乗り物等の移動物体の場合、一様流とは、その乗り物に乗っている観察者から見た、相対的な気流や水流等の流れである。すなわち、移動物体の座標系においては、実際には空気や水が静止していても、無限遠点から流れてくる相対的な流れを観念することができる。
 一様流を構成する要素は、吸い込み湧き出し対(1-source-sink point)と呼び、そのルールの正しさは数学的に証明可能である。
 ここで、一様流ではなく吸い込み湧き出し対と呼ぶことには理由がある。その理由の説明のため、いくつかの数学的な解説を以下に行う。一様流が存在するとき、考えている流れの領域は無限に拡がる平面であり、その中に一つまたは複数の穴(障害物)が埋め込まれているような連結外部領域であるが、このような流れを模式的(Schematic)に表現する上では扱いにくい。そのため、数学におけるstereographic projection(ステレオ射影・立体射影)と呼ばれる射影法によって、平面を球面の上に投影して扱ってもよい。この場合、平面における無限遠点は球面の北極に、平面の原点は南極に対応させることができる。
 このようにしておくと、一様流は球面の北極における流れの湧き出しと吸い込み対といった流れ構造となり、図4(b)のような流れ場に対応していることが数学的に示せる。さらに模式的に流れ場を表現するためには、球面は対称性の高い形状であるということを利用して、北極と南極の位置を適当にずらすことができるので、無限遠点を南極に、ある円形の穴(障害物)の中心を北極にあわせてから、再びstereographic projectionを使って平面に投影すると、南極に対応する原点付近の近くでは図4(b)のような流れができる。さらに無限遠点に中心を持つ円境界は平面の外側円境界に投影されるので、結果として平面全体の流れ場が、例えば図5のような形の有界な領域として表現できる。なお、連結外部領域Dζ(M)をこの変換方法により,無限遠点の吸い込み湧き出し対を原点に、また、M+1個の境界の1つを選んで、それを単位円外部領域に移したときに得られる、M個の境界をその内部に含む単位円内部の連結領域をD(M)と書く。以後の記載はすべてこのD(M)を用いて行う。したがって、図4(b)のような表現は適切な射影法を通じて、平面全体の中に一様流が入っている流れと等価なものである。本実施形態の説明では、流れを模式的に示すために便利であることから、このような射影法を用いて図に表現する場合がある。
 図5は、このような連結領域D(M)における構造安定な流れの位相的分類を行う特徴的な軌道(流線)を全て記述した図である。図5(a)に示すように、まず吸い込み湧き出し対から出て自分自身に戻ってくる軌道をss-orbitと呼ぶ。この線一本一本が連結外部領域での一様流の流線を表す。次に、図5(b)に示すように、吸い込み湧き出し対から出て境界の上につながる軌道をss-∂-saddle connectionと呼び、図5(c)に示すように、その軌道がつながっている境界上の点をss-∂-saddleと呼ぶ。
 また、図5(e)に示すように、吸い込み湧き出し対からではなく、ある境界の上の点から出て同じ境界上の点につながる軌道を∂-saddle connectionと呼び、図5(d)に示すように、これによってつながれている境界上の点を∂-saddleと呼ぶ。また、図5(h)に示すような、境界上にない双曲型の停留点をsaddle point(サドル点)と呼ぶが、図5(f)に示すように、吸い込み湧き出し対から出る軌道で、このsaddle pointにつながる軌道をss-saddle connectionと呼ぶ。また、図5(g)に示すように、境界や渦の回りを作る閉曲線軌道をclosed orbitと呼び、図5(i)に示すように、saddle pointから出てそれ自身に戻るような軌道をhomoclinic saddle connectionと呼ぶ。対象とする構造安定な流れは、これらの軌道の組み合わせによってしか表現されないことが数学的に証明できる。
 本実施形態では、上述したステップSA-2において、穴の数がM個ある連結外部領域D(M-1)の流れに、一つの穴とそれに伴う流れの構造を付け足すことで、一つ穴の多い多重連結外部領域D(M)の流れ場を帰納的に構成していく。そのため、もっとも簡単な穴が一つの単連結外部領域D(0)や、二重連結外部領域D(1)で、これらの帰納的構成の初期構造になるものをステップSA-1で与えている。
 具体的には、上述した基本となる流れパターンは、3種類あり、
 1)吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、
 2)吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、および
 3)吸い込み湧き出し対をもたないパターンO、
 である。ここで、図6は、初期構造となる3種類の構造安定な流れパターンを模式的に示す図である。
 すなわち、図6(a)と(b)に示すように、穴が一個の単連結外部領域D(0)にある流れは、パターンIとパターンIIの2種類存在する。これらのパターンは、ともに吸い込み湧き出し対を持ち、数学的にはこれら2種類しかないことが証明できる。一様流を仮定した吸い込み湧き出し対を持つような流れに対して、原則的には二重連結外部領域D(1)は、これらから構成されるが、吸い込み湧き出し対を持たない流れはここから構成されないので、その流れを構成するために必要な初期流れが図6(c)に模式的に示されるパターンOである。なお、これらの位相構造は表現の簡便さのため、吸い込み湧き出し対を丸囲みSと図示し、ss-orbitやclosed orbitは、無限に存在するので表現せず、以後、模式的に図6(d)や(e)のように簡略的に表す。また、図6(c)に示すように、吸い込み湧き出し対を持たない二重連結外部Dζ(1)における流れのパターンclosed orbitも、すべて書かないで図6(f)のように簡略に記す。
(1-2.操作語の説明)
 帰納的に流れを構成していくために、穴を一つとそれに伴う流れの構造を追加するという「操作」について、図7および図8を参照して説明する。すなわち、穴の数がM個ある連結外部領域D(M-1)の流れに、一つの穴を加えて多重連結外部領域D(M)の流れを求める操作について説明する。
 上述したステップSA-2において、位相幾何学的に採り得る5種類の操作は、
 1)一本のss-orbitを、一つのsaddle point、それを結び内部に穴をもつhomoclinic saddle connectionと二つのss-saddle connectionに置き換えるA操作、
 2)一本のss-orbitを、二つのss-∂-saddle connectionと新たに追加した境界上の二つの∂-saddleに置き換えるA操作、
 3)一本のclosed orbitを、一つの穴とsaddle pointを追加して8の字をした2本のhomoclinic軌道に置き換えるB操作、
 4)一本のclosed orbitを、新たに追加した穴の境界上に二つ∂-saddleをつけて一本の∂-saddle connectionでつなぐような軌道に置き換えるB操作、および、
 5)既に2k個(k>0)の∂-saddleをもつ境界に、新たに二つの∂-saddleを付け加えて一本の∂-saddle connectionでつなぎ内部に新たに付け加えた穴を置くC操作、
 である。ここで、図7は、穴を一つ付け加えて構造安定な流れを構成する5種類の操作を模式的に示した図である。
 図7(a)に示すように、操作AとAは一本のss-orbitに対して行われる。また、図7(b)に示すように、操作BとBは、一本のclosed orbitに対して行われる。また、図7(c)に示すように、操作Cは、既に∂-saddleを持つ境界に対して行われる。なお、構造安定性を維持しながらそのようなことを可能にする操作が、上記の5種類しかないことが数学的に証明可能である。このステップを繰り返してできる流れパターン図を、パターンI,IIから始めたときはss-saddle connection diagramと呼び、パターンOから始めたときは、saddle connection diagramと呼ぶ。
 ステップSA-1で与えた初期構造の3種類の流れパターンI,II,Oから、これらの操作を行って一つずつ穴を増やす(M→M+1)ことで、多くの穴を持つ領域における流れが帰納的に構成されるので、本実施形態では、その操作を表す操作語の列を文字列と見なして列挙することで、流れ場の語表現(Word representation)を得ることができる。ここで、図8は、二つの構造物と一様流がある場合の構造安定な流れパターンの全分類を示す図である。図8に示すように、単連結外部領域D(0)における初期構造のパターンI,IIに対して、操作語を付与することで、二重連結領域D(1)における全ての流れパターンを記述することができる。ただし、図8に示す全ての流れパターンは、2種類(I,II)×5種類(A,A,B,B,C)の計10種類とはなっていない。すなわち、操作語は、5種類の操作語を任意に並べてできるわけではなく、数学的な理由から様々な制約がつく。
 ここで、制約について説明すれば以下のようになる。すなわち、図7を用いて上述したように、操作AとAは、一本のss-orbitに対して行われるので、この操作を行う前提として、一本のss-orbitの存在が不可欠となる。また、操作BとBは、一本のclosed orbitに対して行われるので、この操作を行う前提として、一本のclosed orbitの存在が不可欠である。また、操作Cは、2個以上の∂-saddleを持つ境界に対して行われるので、この操作を行う前提として、∂-saddleを持つ境界の存在が不可欠となる。そのため、パターン語をI,II,Oのどこからはじめるかによって並べ方のルールは異なる。上記の制約条件に基づいて導出される、各パターン語I,II,Oからはじまる文字列の並べ方のルールについて、以下に説明する。
 まず、湧き出し吸い込み対を持たないOのパターン語からはじめる場合には、次のようなルールがある。Oから始まる語表現に対して、それが構造安定な流れを表すために以下が必要十分である。
O-1)実際に施すことができる操作は、B,B,Cのみであり、その結果、Oから始まる語表現はこれら三つの語を列挙したものとなる。
O-2)操作列の語表現においてCの語が含まれるためには、それ以前に必ずBが存在しなければならない。
 このような文字列をO系列の語(O-word)と呼び、そのルールの正しさは数学的に証明可能である。
 次に、パターン語Iから始まる語表現については、以下のルールが成り立つ必要がある。
I-1)実施可能な操作はA,A,B,B,Cの全てであり、その結果、Iから始まる語表現はこれら5種類の操作語を列挙したものである。
I-2)操作列の語表現において、BあるいはBの語が含まれるためには、それ以前に必ずCかAが存在しなければならない。
 このような文字列をI系列の語(I-word)と呼び、そのルールの正しさは数学的に証明可能である。
 最後にパターン語IIから始まる語表現については、以下のルールが成り立つ必要がある。
 II-1)実施可能な操作はA,B,B,Cであり、その結果、IIから始まる語表現はこれら四つの語を列挙したものである。
 II-2)操作列の語表現においてCの語が含まれるために
は、それ以前に必ずBが存在しなければならない。
 このような文字列をII系列の語(II-word)と呼び、そのルールの正しさは数学的に証明可能である。
(2.グラフ表現作成工程)
上記図2のグラフ表現作成工程(ステップS1)について詳細に説明する。まず、本実施の形態で使用しているグラフ理論について説明する。
(2-1.グラフ理論)
 グラフT=(V,E)とは、「vertex(頂点)」と呼ばれる点の集合(頂点集合)Vとその頂点の間を結ぶ「edge(エッジ)」と呼ばれる集合Eのペアとして与えられる集合である。一般にグラフは多様な構造を持ちうるが、本実施の形態のグラフ表現理論では、グラフ全体の集合におけるある特定の構造をもった以下のグラフの集合を考える。
 1)ツリー,木(Tree)とは、任意の二つ頂点が一本のedgeによってのみ結合されているようなグラフを指す。
 2)ルート付き(Rooted)グラフとは、ある特定の頂点(以下、ルート(root)と呼ぶ)が存在しているグラフを指す。ルート付きグラフに対しては、このルートから各頂点v∈Vへのエッジの連結による最短経路を考えることができるので、これをvの高さ(height)と呼び、記号をht(v)と書く。これによりルート付きツリーTに対してはht(T):=maxv∈Vht(v)によって、ツリーT自体の高さを考えることができる。
 3)グラフが向き有り(directed)であるとは、すべてのエッジに親子の順序が入っているようなものを指す。向き有りグラフにおける、頂点v∈Vからw∈Vへのエッジは、v⇒w∈Eと表わす。このとき、vはwの親、wはvの子と呼ぶ。Г(v)と書いて、頂点v∈Vの子供全体の集合を表す。すなわち、Г(v):={w∈V│v⇒w∈E}。また、その集合に含まれる子供頂点の数#Г(v)を持って,vのout-degreeとよび、逆にvに入ってくるエッジの数をvのin-degreeと呼ぶ。
 4)グラフがラベル付き(labelled)であるとは、全ての頂点に特定のラベルが割り当てられているようなものを指す。
 本実施の形態においては、ルート付き、ラベル付き、向き有りのツリーを考える。またそのルートはin-degreeがゼロ、すなわちそこに入ってくるエッジが存在しないような頂点を考える。頂点のラベルとしては、{oφ,o,o,+φ,+,+,-φ,-,-}のいずれかを割り当てるものとする。このとき,ラベルがoφ(o,o,+φ,+,+,-φ,-,あるいは-)となっている頂点集合の部分集合をVoφ(Vo、Vo、V+φ、V+、V+、V-φ,V-,あるいはV-2)と表す。頂点vに対して、l(v)と書けばそれは頂点に割り当てられたラベルを表すことにする。以後の便宜のため、以下のような集合を定義しておく。Vo=Voφ∪Vo∪Vo、V+=V+φ∪V+∪V+、V-=V-φ∪V-∪V-、V=Vo∪V+∪V-、V=Vo∪V+∪V-、ただし、∪の記号は、disjoint unionを表している。加えて、頂点vの子頂点集合は、Г(v)は以下のように分割することができる。
 Гo(v)=Г(v)∩Vo,Гo(v)=Г(v)∩Vo,Г+(v)=Г(v)∩V+,Г-(v)=Г(v)∩V-,Г+(v)=Г(v)∩V+,Г-(v)=Г(v)∩V-,Г(v)=Г(v)∩V,Г(v)=Г(v)∩V,Г(v)=Г(v)∩Vである。
 さらに、頂点v∈Vに対して「符合」をsgn(v)と書けば、v∈V+∪Vならばsgn(v)=+、v∈Vならばsgn(v)=-を与えるものとする。以後ラベルを並べて流れの表現を与えるので、これらのラベルの出現順序には以下のような順序関係のルール(1)を定めておく。
       o>o>+>->+>-・・・(1)
(2-2.二次元構造安定なハミルトンベクトル場のツリー表現)
 本実施の形態のグラフ表現作成工程に係る二次元構造安定なハミルトンベクトル場のツリー表現について説明する。以下、O-wordで表現される流れパターンのグラフ表現と、I,II-wordで表現される流れパターンのグラフ表現について説明する。
(2-2-1.O-wordで表現される流れパターンのグラフ表現)
 Hを二次元領域D(M)上のO-wordで表現される構造安定なハミルトンベクトル場とし、Dをそのsaddle connection diagramとする。このハミルトンベクトル場Hに対して固有のルート付き、ラベル付き、及び向き有りのツリーT=(V,E)を割り当てる方法とその平面グラフとしての可視化アルゴリズムを以下に説明する。
 まず、領域全体D(M)からDをすべて取り除くと、いくつかの(無限の周期軌道を含む)円環開集合である連結成分からなる集合C=D(M)\Dが構成できる。この各連結成分を頂点集合Vとする。ルートの頂点はこれらの連結成分のうちもっとも外側にある連結成分とし、そのラベルとして含まれている周期軌道が反時計回りであるときは+φ、時計回りであるときは-φを割り当てる。
 次にエッジは以下のように作成する。二つの頂点v、w∈Vに対して、これらの頂点に対応する連結成分の閉包の共通部分(つまり成分の共通境界)が1次元dim(cl(v)∩cl(w))=1、かつ、vがwの外側の成分となっている場合に、その二つの間に向き有りエッジv⇒w∈Eを構成する。さらに、この共通部分(境界)集合cl(v)∩cl(w)がsaddle connectionであり、子連結成分wが反時計回り(あるいは、時計回り)の周期軌道を持つとき、wのラベルとして+(あるいは-)を割り当てるものとする。同様に、その共通部分(境界)集合cl(v)∩cl(w)が∂-saddle connectionであり、w内の周期軌道の向きが反時計回り(あるいは時計回り)であれば、wのラベルを+(あるいは-)とする。
 Cの内部にはss-orbitは一つも存在しないので、連結成分vの子頂点集合Г(v)は、Г(v)=Г+(v)∪Г-(v)∪Г+(v)∪Г-(v)となることに注意する。Vo=Vo=φでもあるので、頂点集合Vの部分集合V=V+∪V-およびV=V+∪V-となっている。これらの部分頂点集合に入っている元の数は全ての構造安定な流れパターンが操作B,B,Cによって構成されていることから見積もることができる。
 図9は、操作B,B,Cによって生成されるすべての流線パターンとそれに対応して連結部分集合(頂点集合)に励起される親子関係を示している。(a)はB、(b)はB、(c)はB、(d)はBk-1,k≧1、(e)はB,l≧1、(f)は、B k+l-1,k,l≧1、(g)はB,k≧1、(h)はB,l≧1を示している。図9において、vは親連結部分集合を表し、その子連結部分集合はw,y,zなどのように表現されている。破線は親連結部分の集合を表すため、そこに含まれている周期軌道を一本だけ向きとともに描いている。
 図9(a)に示すように、操作Bによって作られた外向き8の字パターンが構成できる。vとhomoclinic saddle connectionを境界として共有する二つの連結成分w∈Vおよびw∈Vで、その符合が同じ、つまりsgn(v)=sgn(w)=sgn(w)となるものが存在する。他方、Bで作られたパターンが図9(b)に示すように内向き8の字パターンであったとき、親成分vをどの連結成分にとるかによって、二種類のvとw∈Vの二つのエッジ(一つはsgn(v)=sgn(w)、もう一つはsgn(v)≠sgn(w)となっている)が構成できる。
 図9(c)に示すように、操作Bによっては二種類の親子エッジ関係が形成される。左側のパターンでは、親連結成分vから、sgn(v)=sgn(w)=sgn(w)かつsgn(v)≠sgn(w)となる三つの子連結成分w,w,w∈Vへのエッジが構成できる。
 右側のパターンでは、親連結成分vからsgn(v)=sgn(w)≠sgn(w)となる二つの子連結成分w,w∈Vへのエッジが構成できる。以下同様に、図9(d)~(f)に示すように、各操作Bk-1,B,B k+l-1(k≧1,l≧1)によってk個のvと同じ向きの周期軌道を含む連結成分と逆向きの周期軌道を含むl個の連結成分間の親子関係のエッジが順次構成される。
 図9(g)は、操作Bによって構成された流線パターンであり、これによりvからw,w∈V、あるいは、vからw∈Vへの同じ符合を持つ領域へのエッジと、k個の反対の符合を持つ子連結成分z∈V(j=1,...,k)へのエッジが構成される。
 最後に、操作B(l≧1)によって作られる図9(h)の流線パターンに対応しては、sgn(v)≠sgn(w)となるw∈Vへの一本のエッジと、同じ符合を持つl本のy∈V(j=1,...,l)へのエッジができる。
 これ以外に、操作B,B,Cによって新しい親子関係を生成するものは存在しないので、以上から、子連結成分の集合に含まれる元の数は#Γ+(v),#Γ-(v)≦2であり、子頂点集合#Γ+(v)、#Γ-(v)については任意の非負整数を選ぶことができる。
 次に、与えられたsaddle connection diagramから得られるグラフ表現をどのように平面グラフとして描画するかを説明する。まず、全ての連結部分v∈Vに対して、その子連結部分集合Γ(v)に含まれる元をラベルに関する順序関係のルール(1)に応じて並び替えておく。ただし、複数の同じラベル+(あるいは-)が含まれているときはサイクリックな順序で並べる(このことは、#Γ+(v)、#Γ-(v)≦2ということから一義的に実現できる。)。一方、ラベル+および-に対応する子連結部分集合Γ+(v)およびΓ-(v)については、その並べ方は以下のようなルールに従うものとする。
 図9にあるyと書かれた連結部分に対してその順序をサイクリックに並べる。すなわち、この中から特定の連結部分をyとして選んで、後は反時計回りに連結成分を並べる。同じ図中zと書かれた子連結成分については単に反時計回りに並べることができる。
 以下で説明する(2-2-1-1.O系列におけるsaddle connection diagramのツリーへの変換処理)及び(2-2-2-1.I,II系列におけるss-saddle connection diagramのツリーへの変換処理)について、以下のように定義する。Nで非負整数のなす集合とし、Nで非負整数のなす有限列の集合とする(すなわち、N=N∪N∪N・・・)。N上に以下のように半順序<を定義する。Nの元s=s・・・s,t=t・・・tに対して、「n<m」または「n=mかつs=t,...,sk-1=tk-1かつs<tとなる自然数k>0が存在する事である」が成り立つ時、s<tと定義する。以下では、頂点集合VからNへの単射関数S:V→Nを与える。このとき、頂点vに対して、S(v)を頂点idと呼ぶ。ここで、単射性から、頂点とその頂点idを同一視する。さらに、rootの頂点idは0と置く(頂点とidを同一視するため)。また、以下の構成から、頂点idは、ツリーの高さ優先探索の順番に対応していることが分かる。以下、σ∈{+,-}を表し、μ=-σとし、Gでツリー、s,s’,t,uで頂点id、Tで頂点idの部分集合を表す。また、以下で、図9,12を参照するが、流れの向きは無視するものとする。
(2-2-1-1.O系列におけるsaddle connection diagramのツリーへの変換処理)
 図10-A~図10-Eは、O系列におけるsaddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。コンピュータ等の装置により、図10-A~図10-Eに示す、O系列におけるsaddle connection diagramのツリーへの変換処理を実行可能である。
 図10-A~図10-Eにおいて、inputをsaddle connection diagram Dとし、saddle connection diagram Dをoutermostにrootがあるよう変換し(ステップS101)、s=0,T={0}と置く(ステップS102)。root 0の外側境界が反時計回りであるか否かを判断する(ステップS103)、root 0の外側境界が反時計回りである場合は(ステップS103の「Yes」)、σ=+,ツリーG=(0,+φ)とし(ステップS104)、root 0の外側境界が反時計回りでない場合、すなわち、時計回りである場合は(ステップS103の「No」)、σ=-,ツリーG=(0,-φ)として(ステップ105)、ステップS106に移行する。
 ステップS106では、root 0が図9(a)の形(パターン)をしているか否かを判断する。root 0が図9(a)の形をしている場合は(ステップ106の「Yes」)、w,wの頂点idを00,01と定め、ツリーGを、図9(a)のツリーのラベルをidとラベルの組に置き換えたツリーG(すなわち、Gを(00,σ)←(0,σφ)→(01,σ))と定め、集合Tに{00,01}を追加する(ステップS107)。すなわち、T={0,00,01}と置き換える。この置換操作を、T←TU{00,01}と表す。この後、処理はステップS114に移行する。
 他方、ステップS106において、root 0が図9(a)の形をしていない場合は(ステップ106の「No」)、ステップS108に移行する。
 ステップS108では、root 0が図9(b)の左の形をしているか否かを判断する。root 0が図9(b)の左の形をしている場合には(ステップ108の「Yes」)、wの頂点idを00と定め、ツリーGを、図9(b)のツリーのラベルをidとラベルの組に置き換えたツリー(すなわち、Gを((0,σφ)→(00,σ)))と定め、集合Tに{00}を追加する(ステップS109)。すなわち、T←TU{00}とする。この後、処理はステップS114に移行する。
 他方、ステップS108において、root 0が図9(b)の左の形をしていない場合は(ステップS108の「No」)、ステップS110に移行する。
 ステップS110において、root 0が図9(d)の形をしているか否かを判断する。root 0が図9(d)の形をしている場合は(ステップS110の「Yes」)、y,...,yの頂点idを00,...,0k-1と定め、ツリーGを、図9(d)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、集合Tに{00,...,0k-1}を追加する(ステップS111)。すなわち、T←TU{00,...,0k-1}とする。この後、処理はステップS114に移行する。
 他方、ステップS110において、root 0が図9(d)の形をしていない場合は(ステップS110の「No」)、ツリーGをσφと定め(ステップS112)、このツリーGをoutputして(ステップS113)、処理を終了する。
 ステップS114において、Tの中にsより大きい元tが存在するか否かを判断する。sがTの中で最大元の場合は(ステップS114の「No」)、ツリーGの頂点idを全て取り除いて、残りのツリーをGとし(ステップS115)、ツリーGをoutputして(ステップS113)、処理を終了する。
 他方、ステップS114において、sがTの中で最大元でない場合は(ステップS114の「Yes」)、s’=min{t∈T|s<t}(すなわち、s’をTの中でsの次に大きい元)として、sをs’と置換して、(s,σ)←(s’,σ’)とする(ステップS116)。この後、処理はステップS117に移行する。
 ステップS117において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sの内部境界が図9(a)の形をしているか否かを判断する。頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sがの内部境界が図9(a)の形をしている場合は(ステップ117の「Yes」)、w,wの頂点idをs0,s1と定め、ツリーGを、図9(a)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,s1}を追加する(ステップS118)。この後、処理はステップS114に戻る。
 他方、ステップS117において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleでない、もしくは、頂点sの内部境界が図9(a)の形をしていない場合は(ステップS117の「No」)、ステップS119に移行する。
 ステップS119において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sの内部境界が図9(b)の左の形をしているか否かを判断する。頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sの内部境界が図9(b)の左の形をしている場合は(ステップS119の「Yes」)、wの頂点idをs0と定め、ツリーGを、図9(b)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0}を追加する(ステップS120)。この後、処理はステップS114に戻る。
 他方、ステップS119において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleでない、もしくは、頂点sの内部境界が図9(b)の左の形をしていない場合は(ステップS119の「No」)、ステップS121に移行する。
 ステップS121において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sの内部境界が図9(d)の形をしているか否かを判断する。頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleからなり、頂点sの内部境界が図9(d)の形をしている場合は(ステップ121の「Yes」)、y,...,yの頂点idを{s0,...,sk-1}と定め、ツリーGを、図9(d)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sk-1}を追加する(ステップS122)。この後、処理はステップS114に戻る。
 他方、ステップS121において、頂点sの外側境界が1つのhomoclinic saddle connectionとsaddleでない、もしくは、頂点sの内部境界が図9(d)の形をしていない場合は(ステップS121の「No」)、ステップS123に移行する。
 ステップS123において、頂点sが図9(c)の左の形をしているか否かを判断する。頂点sが図9(c)の左の形をしている場合は(ステップS123の「Yes」)、ステップS124に移行する。
 ステップS124において、頂点sのラベルσのσが+であるか(すなわち、sの中の流れは反時計回りであるか)否かを判断する。頂点sのラベルσのσが+である(すなわち、sの中の流れは反時計回りである)場合は(ステップS124の「Yes」)、w,w,wの頂点idをs0,s1,s2と定め、ツリーGを、図9(c)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,s1,s2}を追加する(ステップS125)。この後、処理は、ステップS114に戻る。
 他方、ステップS124において、頂点sのラベルσのσが+でない場合、すなわち、頂点sのラベルσのσが-である(すなわち、sの中の流れは時計回りである)場合は(ステップS124の「No」)、w,w,wの頂点idをs1,s2,s0と定め、ツリーGを、図9(c)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,s1,s2}を追加する(ステップS126)。この後、処理は、ステップS114に戻る。
 他方、ステップS123において、頂点sが図9(c)の左の形をしていない場合は(ステップS123の「No」)、ステップS127に移行する。
 ステップS127において、頂点sが図9(c)の右の形をしているか否かを判断する。頂点sが図9(c)の右の形をしている場合は(ステップS127の「Yes」)、ステップS128に移行する。
 ステップS128において、頂点sのラベルσのσが+である(すなわち、sの中の流れは反時計回りである)か否かを判断する。頂点sのラベルσのσが+である(すなわち、sの中の流れは反時計回りである)場合は(ステップS128の「Yes」)、w,wの頂点idをs0,s1と定め、ツリーGを、図9(c)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,s1}を追加する(ステップS129)。この後、処理はステップS114に戻る。
 他方、ステップS128において、頂点sのラベルσのσが+でない場合、すなわち、頂点sのラベルσのσが-である(すなわち、sの中の流れは時計回りである)場合は(ステップS128の「No」)、w,wの頂点idをs1,s0と定め、ツリーGを、図9(c)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,s1}を追加する(ステップS130)。この後、処理はステップS114に戻る。
 他方、ステップS127において、頂点sが図9(c)の右の形をしていない場合は(ステップS127の「No」)、ステップS131に移行する。
 ステップS131において、頂点sが図9(h)の形をしているか否かを判断する。頂点sが図9(h)の形をしている場合は(ステップS131の「Yes」)、y,...,y,wの頂点idをs1,...,sl,s0と定め、ツリーGを、図9(h)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sl}を追加する(ステップS132)。この後、処理はステップS114に戻る。
 他方、ステップS131において、頂点sが図9(h)の形をしていない場合は(ステップS131の「No」)、ステップS133に移行する。
 ステップS133において、頂点sが図9(g)の左の形をしているか否かを判断する。頂点sが図9(g)の左の形をしている場合(ただし、k=0の場合も含む)は(ステップ133の「Yes」)、w,w,z,...,zの頂点idをs0,s1,s2,...,sk+1と定め、ツリーGを、図9(g)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sk+1}を追加する(ステップS134)。この後、処理はステップS114に戻る。
 他方、ステップS133において、頂点sが図9(g)の左の形をしていない場合は(ステップS133の「No」)、ステップS135に移行する。
 ステップS135において、頂点sが図9(g)の右の形をしているか否かを判断する。頂点sが図9(g)の右の形をしている場合(ただし、k=0の場合も含む)は(ステップS135の「Yes」)、w,z,...,zの頂点idをs0,s1,...,skと定め、ツリーGを、図9(g)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sk}を追加する(ステップS136)。この後、処理はステップS114に戻る。
 他方、ステップS135において、頂点sが図9(g)の右の形をしていない場合は(ステップS135の「No」)、ステップS138に移行する(この場合、ステップS137に示すように、頂点sは図9(f)の形をしている)。
 ステップS138において、頂点sのラベルσのσが+である(すなわち、sの中の流れは反時計回りである)か否かを判断する。頂点sのラベルσのσが+である(すなわち、sの中の流れは反時計回りである)場合は(ステップS138の「Yes」)、y,...,y,z,...,zの頂点idをs0,...,sk-1,sk,...,sl+k-1と定め、ツリーGを、図9(f)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sl+k-1}を追加する(ステップS139)。この後、処理はステップS114に戻る。
 ステップS138において、頂点sのラベルσのσが+でない場合、すなわち、頂点sのラベルσのσが-である(すなわち、sの中の流れは時計回りである)場合は(ステップS138の「No」)、y,...,y,z,...,zの頂点idをsl,...,sl+k-1,s0,...,sl-1と定め、ツリーGを、図9(f)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{s0,...,sl+k-1}を追加する(ステップS140)。この後、処理はステップS114に戻る。
 以上の操作が、O系列におけるsaddle connection diagramのツリーへの変換処理である。
(2-2-2.I-wordおよびII-wordによって表現される流れのグラフ表現)
 HをD(M)内に1-source-sink pointを持つ構造安定なハミルトンベクトル場、Dをそのss-saddle connection diagramとする。このとき、C=D(M)\Dは、closed orbitを含む円環開領域またはss-orbitを含む開円板領域となる連結成分からなる。これに対して、ルート付き、ラベル付き向き有りツリーを以下のように構成する。O-wordの時と同様に、各Cの連結成分を頂点集合に対応づける。ルートとなる連結成分はCの中で1-source-sink pointに最も近い開円板領域で、その閉包は1-source-sink pointを含みかつ内部に時計回りのss-orbitを含むものを選ぶ。この連結成分にラベルoφを割り当てる。このような決め方によりルートは一義的に決定することができる。
 図11は、ルートとなる連結成分をC=D(M)\Dから選ぶ方法を説明するための図である。ここで、DはI-wordあるいはII-wordで表現できるss-saddle connection diagramとする。破線はルートとなる連結成分に含まれるss-orbitを一本、その向きとともに書いたものである。
 図11において、(a)は、語表現IAで表現されたss-saddle connection diagramであり、ss-∂-saddle connectionの流れの向きは右から左である。(b)は、語表現IAを持つss-saddle connection diagramで、ss-∂-saddle connectionの向きは左から右である。(c)は、語表現IIAを持つss-saddle connection diagramで、その最も外側にある連結成分に含まれた周期軌道の向きが時計回り。(d)は、語表現IIAは(c)と同じであるが周期軌道の向きが反対である。各ss-saddle connection diagramの下には、ルートとなっている連結成分内にある点を無限大に写すような等角写像によってできるss-saddle connection diagramの像を示している。この写像による像では、ルートに含まれるss-orbitはすべて反時計回りになっている。
 例えば、図11(a)に示すような、IAなる語表現を持つss-saddle connection diagramを考えると、1-source-sink pointの真上にある連結成分(破線で書いた曲線が含まれる領域)がルートになる。もし、同じss-saddle connection diagramが同じ形でもその流れの向きが、図11(b)のように反対になれば、その定義からルートは1-source-sink pointの真下にある連結成分となる。図11(c)の語表現IIAを持つss-saddle connection diagramに対しては、1-source-sink pointの真上の連結成分が、向きを逆転すると(図11(d))、その真下の連結成分がルートとなる。このルートとなっている連結成分にある点を一つ選び、それを無限遠点に写すような等角写像でss-saddle connection diagramを写すと、その像において、ルートの連結成分は最も外側の連結成分となり、さらにその中でss-orbitは常に反時計回りになることに注意する(図11の各図の下の部分参照)。
 こうして得られ連結成分からなる頂点集合に対して、頂点間の向き付けエッジを、ルート連結成分が最も外側になるよう等角写像で写したものを使って図11に示したように定義する。二つの連結成分v,w∈Vの閉包の共通部分の次元が1、すなわちdim(cl(v)∩cl(w))=1かつvがwの外側にあるとき、vからwへの向き付けエッジを構成する。ss-orbitを持つ連結成分w∈Vに対して、もしcl(v)∩cl(w)がss-saddle connection(あるいはss-∂-saddle connection)を含む場合、その子頂点wにはo(あるいはo)のラベルを割り当てる。内部に反時計回りのclosed orbitを持つ連結成分w∈Vに対して、もしcl(v)∩cl(w)がsaddle connection(あるいは∂-saddle connection)となっている場合、子の連結成分wにはラベル+(あるいは+)を付与する。同様にして、ラベル-および-を時計回りのclosed orbitを持つ連結成分に付与する。
 この結果、ハミルトンベクトル場Hのss-saddle connection diagramに対して固有のルート付き、ラベル付き、向き有りツリーT=(V,E)を対応させることができる。I-wordあるいはII-wordの語表現を持つ構造安定なハミルトンベクトル場Wのグラフ表現T=(V,E)を平面グラフとして可視化する方法は以下の通りである。
 まず、子連結成分の集合に関して、Г(v)=Г(v)∪Г(v)∪Г(v)となっていることに注意する。いま、Г(v)に入っている元は予め順序関係のルール(1)に従って並びかえられており、また同じラベルを持つГ+(v)∪およびГ-(v)の元については、O-wordの可視化の時と同様にサイクリックに反時計回りに並べられているとしておく。VやVに入っているclosed orbitを持つ連結成分vに対する平面グラフへの可視化は、フローチャート(図10-Aから図10-E)を用いればよいので、後はVに入っている連結成分に対するグラフの描画について考えればよい。
 図12は、操作A,A,Cによって生成されるss-saddle connection diagramの局所構造とそれに対応して得られるグラフ表現の可視化を示したものである。図12において、親連結成分vはoφ(oまたはo)と表現されている。
 図12において、(a)はA、(b)はA、(c)はA、(d)はA(k≧0)、(e)はA(l≧1)、(f)はA k+l(k,l≧1)、(g)は、A(k≧0)、(h)はA(l≧1)によって構成される局所的なss-saddle connection diagramとそれに対応して生成される連結成分の親子関係、親連結成分vに含まれる反時計回りのss-orbitは破線で示している。親連結成分のラベルはo(ただし、*はφ、0あるいは2のいずれか)と表現されている。ss-orbitの方向は常に反時計回りになるようにルートを定めているので、その子になる連結成分に含まれる軌道の流れる方向は自動的に決まる。
 より具体的に説明すると、図12(a)にある外向きhomoclinic saddle connectionに対して、親連結成分vは二つの子連結成分w∈Voとy∈V+を持つので、oからoと+への二本の矢印を描画する。図12(b)にあるような内向きhomoclinic saddle connectionに対しては、親をどの連結成分に選ぶかに応じて、oあるいは-への矢印を描くことができる。操作Aで作られる構造に対しては、図12(c)に示すように二種類のパターンが作られるので、それに応じてo、+、-への三本の矢印を描くかoと-への二本のエッジを描くことができる。
 操作A(k≧0)で作られる図12(d)のようなk個のclosed orbitを含む連結成分{z,...,z}を持つ局所的なss-saddle connection diagramの局所構造に対応して、oおよびk個の+の矢印を描く。操作A(l≧1)で作られるl個の時計回りのclosed orbitを含む子連結成分{y,...,y}を持つ局所流線構造(図12(e))に対しては、l個の矢印を-へ描けばよい。図3(f)には、A k+l(k,l≧1)で作られる一個の子連結成分w∈Voとl個の子の連結成分{y,...,y}とk個の子連結成分{z,...,z}を持つ局所構造と、それに対応してk個の+とl個の-への矢印が描けることが示されている。最後に、操作A(k≧0)で与えられる図12(g)(あるいは図12(h)のようなA(l≧1))の構造に対しては、o,-(あるいはo,+)の矢印を描いてその右側にk個の+(あるいはl個の-)への矢印を描く。
(2-2-2-1.I,II系列におけるss-saddle connection diagramのツリーへの変換処理)
 図13-A~図13-Dは、I,II系列におけるss-saddle connection diagramのツリーへの変換処理を説明するためのフローチャートである。図13-A~図13-Dに示す、I,II系列におけるss-saddle connection diagramのツリーへの変換処理は、コンピュータ等の装置によって実行可能である。
 図13-A~図13-Dにおいて、inputをss-saddle connection diagram Dとして、所定の座標変換をして、図11に示すような、outermostにrootがあるよう変換し(ステップS141)、u=0,T={0},ツリーG=(0,oφ)と置く(ステップS142)。
 ステップS143において、root 0が図12(a)の形をしているか否かを判断する。root 0が図12(a)の形をしている場合は(ステップS143の「Yes」)、w、yの頂点idを1,00と定め、ツリーGを、図12(a)のツリーのラベルをidとラベルの組に置き換えたツリー(すなわち、Gを(1,o)←(0,oφ)→(00,+))と定め、Tに{1,00}を追加する(すなわち、T={0}をT={0,1,00}と置換する)(ステップS144)。この後、処理はステップS149に移行する。
 他方、ステップS143において、root 0が図12(a)の形をしていない場合は(ステップS143の「No」)、ステップS145に移行する。
 ステップS145において、root 0が図12(b)の左の形をしているか否かを判断する。root 0が図12(b)の左の形をしている場合は(ステップS145の「Yes」)、wの頂点idを1と定め、ツリーGを、図12(b)のツリーのラベルをidとラベルの組に置き換えたツリー(すなわち、Gを((0,oφ)→(1,o))と定め、Tに{1}を追加する(ステップS146)。この後、処理はステップS149に移行する。
 他方、ステップS145において、root 0が図12(b)の左の形をしていない場合は(ステップS145の「No」)、ステップS148に移行する(この場合、ステップS147に示すように、root 0が図12(d)の形をしている)。
 ステップS148において、w,z,...,zの頂点idを1,00,...,0k-1と定め、ツリーGを、図12(d)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{1,00,...,0k-1}を追加し、ステップS149に移行する。
 ステップS149において、Tの中にuより大きい元tが存在するか否かを判断する。uがTの中で最大元の場合は(ステップS149の「No」)、ツリーGの頂点idを全て取り除いて、残りのツリーをGとおき(ステップS150)、ツリーGをoutputして(ステップS151)、処理を終了する。
 他方、ステップS149において、uがTの中で最大元でない場合には(ステップS149の「Yes」)、u’=min{t∈T|u<t}(すなわち、u’をTの中でuの次に大きい元)として、(u,σ)を(u’,σ’)と置換する(ステップS152)。この後、処理はステップS153に移行する。
 ステップS153において、uが自然数であるか(すなわち、oφ,o,oのidに対応するか)否かを判断する。uが自然数でない場合は(ステップS153の「No」)、O-wordの処理である図10-BのS114に移行する。
 他方、ステップS153において、uが自然数の場合は(ステップS153の「Yes」)、ステップS155に移行する。
 ステップS155において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと、丸囲みS(吸い込み湧き出し対:図6参照)からなり、uの内側境界が図12(a)の形をしているか否かを判断する。頂点uの外側境界が1つのhomoclinic ss-saddle connectionと丸囲みSからなり、uの内側境界が図12(a)の形をしている場合は(ステップS155の「Yes」)、w、yの頂点idをu+1,u0と定め、ツリーGを図12(a)のツリーのラベルをidとラベルの組に置き換えたツリー(すなわち、Gを(u+1,o)←(u,o)→(u0,+))と定め、TをTU{u+1,u0}と置換する(ステップS156)。この後、処理はステップS149に戻る。
 他方、ステップS155において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと丸囲みSからなっていない、または,uの内側境界が図12(a)の形をしていない場合は(ステップS155の「No」)、ステップS157に移行する。
 ステップS157において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと丸囲みSからなり、uの内側境界が図12(b)の左の形をしているか否かを判断する。頂点uの外側境界が1つのhomoclinic ss-saddle connectionと丸囲みSからなり、uの内側境界が図12(b)の左の形をしている場合は(ステップS157の「Yes」)、wの頂点idをu+1と定め、ツリーGを、図12(b)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1}を追加する(ステップS158)。この後、処理はステップS149に戻る。
 他方、ステップS157において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと丸囲みSからなっていない、または,uの内側境界が図12(b)の左の形をしていない場合は(ステップS157の「No」)、ステップS159に移行する。
 ステップS159において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと丸囲みSからなり、uの内側境界が図12(d)の形をしているか否かを判断する。頂点uの外側境界が1つのhomoclinic ss-saddle connectionと丸囲みSからなり、uが図12(d)の形をしている場合は(ステップS159の「Yes」)、w,z,...,zの頂点idをu+1,u0,...,uk-1と定め、ツリーGを、図12(d)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,...,uk-1}を追加する(ステップS160)。この後、処理はステップS149に戻る。
 他方、ステップS159において、頂点uの外側境界が外向きhomoclinic ss-saddle connectionと丸囲みSからなっていない、または,uの内側境界が図12(d)の形をしていない場合は(ステップS159の「No」)、ステップS161に移行する。
 ステップS161において、頂点uが図12(c)の左の形をしているか否かを判断する。頂点uが図12(c)の左の形をしている場合は(ステップS161の「Yes」),w,y,yの頂点idをu+1,u0,u1と定め、ツリーGを、図12(c)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,u1}を追加する(ステップS162)。この後、処理はステップS149に戻る。
 ステップS161において、頂点uが図12(c)の左の形をしていない場合は(ステップS161の「No」)、ステップS163に移行する。
 ステップS163において、頂点uが図12(c)の右の形をしているか否かを判断する。頂点uが図12(c)の右の形をしている場合は(ステップS163の「Yes」),w,wの頂点idをu+1,u0と定め、ツリーGを、図12(c)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0}を追加する(ステップS164)。この後、処理はステップS149に移行する。
 他方、ステップS163において、頂点uが図12(c)の右の形をしていない場合は(ステップS163の「No」)、ステップS165に移行する。
 ステップS165において、頂点uが図12(g)の形をしているか否かを判断する。頂点uが図12(g)の形をしている場合は(ステップS165の「Yes」)、w、y、z,...,zの頂点idをu+1,u0,u1,...,ukと定め、ツリーGを、図12(g)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,...,uk}を追加する(ステップS166)。この後、処理はステップS149に戻る。
 他方、ステップS165において、頂点uが図12(g)の形をしていない場合は(ステップS165の「No」)、ステップS167に移行する。
 ステップS167において、頂点uが図12(h)の左の形をしているか否かを判断する。頂点uが図12(h)の左の形をしている場合は(ステップS167の「Yes」)、w,z,y,...,yの頂点idをu+1,u0,u1,...,ulと定め、ツリーGを、図12(h)の左のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,u1,...,ul}を追加する(ステップS168)。ステップS149に戻る。
 ステップS167において、頂点uが図12(h)の左の形をしていない場合は(ステップS167の「No」)、ステップS169に移行する。
 ステップS169において、頂点uが図12(h)の右の形をしているか否かを判断する。頂点uが図12(h)の右の形をしている場合は(ステップS169の「Yes」)、w,y,...,yの頂点idをu+1,u0,...,ul-1と定め、ツリーGを、図12(h)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,...,ul-1}を追加する(ステップS170)。この後、処理はステップS149に戻る。
 他方、ステップS169において、頂点uが図12(h)の右の形をしていない場合は(ステップS169の「No」)、ステップS171に移行する。
 ステップS171において、頂点uが図12(f)の形をしているか否かを判断する。頂点uが図12(f)の形をしている場合は(ステップS171の「Yes」)、w,z,...,z,y,...,y,の頂点idをu+1,u0,...,uk-1,uk,...,ul+k-1と定め、ツリーGを、図12(f)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u+1,u0,...,ul+k-1}を追加する(ステップS172)。この後、処理はステップS149に戻る。
 ステップS171において、頂点uが図12(f)の形をしていない場合は(ステップS171の「No」)、ステップS173に移行する。
 ステップS173において、頂点uが図12(b)の右の形をしているか否かを判断する。頂点uが図12(b)の右の形をしている場合は(ステップS173の「Yes」)、yの頂点idをu0と定め、ツリーGを、図12(b)の右のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u0}を追加する(ステップS174)。この後、処理はステップS149に戻る。
 他方、ステップS173において、頂点uが図12(b)の右の形をしていない場合は(ステップS173の「No」)、ステップS175に移行する。
 ステップS175において、頂点uが図12(e)の形をしているか否かを判断する。頂点uが図12(e)の形をしている場合は(ステップS175の「Yes」)、y,...,yの頂点idをu0,...,ul-1,と定め、ツリーGを、図12(e)のツリーのラベルをidとラベルの組に置き換えたツリーと定め、Tに{u0,...,ul-1}を追加する(ステップ176)。この後、処理はステップS149に戻る。
 他方、ステップS175において、頂点uが図12(e)の形をしていない場合は(ステップS175の「No」)、ステップS149に戻る(この場合、ステップS177に示すように、頂点uは、図12(a)の形をしている)。
 以上の操作が、I,II系列におけるss-saddle connection diagramからのツリーへの変換処理である。
(2-2-3.グラフ表現とその正規表現)
 以上のことから、1-source-sink pointを持つss-saddle connection diagramに対して、そのルートを一義的に決定でき、また、そこから操作A,A,B,B,Cによって構成されるss-saddle connection diagramの局所構造と、そこから導入される連結成分の間の親子関係は図9と図12の中で全て表現されているので、以下のことが示されたことになる。
 Proposition 3.1:各1-source-sink pointを持つ構造安定なハミルトンベクトル場の流線位相構造には固有のルート付き、ラベル付き、向き有りツリー表現が対応する。
 1-source-sink pointを持たない構造安定なハミルトンベクトル場については、ルートの選び方は一義的ではない。なぜなら、内部に∂-saddleをもたない円境界を含む任意の連結成分に対して、ある連続写像が存在して、それによってその連結成分を常に最も外側の連結成分にすることができるからである。したがって、saddle connection diagramについては、その連結成分の数だけ異なるグラフ表現が存在することになる。このような曖昧さを取り除くために、こうした連続写像によって移りあってできるグラフ表現はそれぞれ別ということにする。すなわち、二つの構造安定なハミルトンベクトル場VとVに対して、V~Vなる同値関係を、ある同相写像D(M)上の同相写像が存在してVとVの各軌道が向きを変えることなく互いに写りあいVとVの外側境界が対応していることによって定義する。このとき、~は同値関係を定め、この同値類に対して固有のルート付き、ラベル付き、向き有りツリーを与えることができる。
 Proposition 3.2:構造安定なハミルトニアンベクトル場の~による同値類の流線位相構造には固有のルート付き、ラベル付き、向き有りツリー表現が対応する。
(3.正規表現作成工程)
 上記図2の正規表現作成工程(ステップS2)について詳細に説明する。さて、グラフ理論でよく知られた事実として、任意のルート付き、ラベル付き、向き有りツリーには固有の「正規表現(regular expression)」なるものを考えることができる。T=(V,E)を上記方法で与えられた構造安定なハミルトンベクトル場Hに対して与えられたグラフ表現とする。このグラフ表現に対して、その正規表現は以下のようにして帰納的に与えられる。まず、ルート以外のすべての頂点のin-degreeはすべて1であることに注意する。もし、ht(T)=0ならグラフはルートのみ、すなわちV={v}であり、その正規表現はl(v)である。いま、高さht(T)=n-1の正規表現Nがあった時に、そこにある頂点集合Tn-1={v,v,...,v}とすると、各頂点vの子頂点集合Г(v)={vi1,vi2,...,vimi}とそれに対応するラベルl=l(v)(i=1,...,m)を使うと、Nに対してlをl(li1,li2,...,limi)に置き換えることによって新しい正規表現が構成できる。ただし、lik=l(vik)(k=1,...,m)である。このようにして構成されたTの正規表現をNTHとすると、その構成手法から以下の命題を得る。
 Proposition 3.3:ルート付き、ラベル付き、向き有りツリーTから正規表現Nへの写像N=f(T)は全単射である。
 Proposition 3.1とProposition 3.2の主張していることは、言い換えれば、(ss-)saddle connection diagram Hからそのグラフ表現Tへの写像T=f(H)は単射(1対1対応)となるので、上のProposition 3.3と合わせれば、合成写像N=f(f(H))も単射になることがわかる。すなわち、同じ語表現を持つ構造安定なハミルトンベクトル場に対して異なる正規表現が与えられ、語表現で区別できないパターンも正規表現を使えば全て区別が可能となる。その実施例として、図14-Aに同じ語表現IACで与えられる構造安定な流線パターンのグラフ表現の可視化とその正規表現を示す。(2-2-2.I-wordおよびII-wordによって表現される流れのグラフ表現)で説明したようにルートの連結成分は1-source-sink pointの真下あるいは真上にあるものを選んでいる。これからもわかるように全ての正規表現は異なっている。
 図14-Aは、IACなる語表現をもつ構造安定な流線パターンのグラフ表現とその正規表現を示している。ルートはその定義から1-source-sink pointの真下の連結成分に対応している。グラフ表現を得るために、ルート連結成分を最も外側の連結成分へと写したものをその右側に描いている。
(3-1.ツリーの正規表現への変換処理)
 図14-Bは、ツリーの正規表現への変換処理を説明するためのフローチャートである。図14-Bに示す、ツリーの正規表現への変換処理は、コンピュータ等の装置により実行可能である。以下の処理では、前処理として,ツリーの全ての点にidが与える。ただし,異なる点のidは異なるものとする。図14-Bにおいて、inputをツリーGとする(ステップS178)。VをツリーGの頂点集合,s=0,T={0},Xを空集合φ,正規表現N=(s,σφ)()と置く(ステップS179)。
 ステップS180において、sの子を左から読んだ結果を(s(1),σ(1)),...,(s(h),σ(h))とする。このとき、正規表現Nの中の(s,σ)()をσ(s(1),σ(1))(),...,(s(h),σ(h))())と置換し、Tに{s(1),...,s(h)}、Xに{s}を追加する。
 ステップS181において、♯V=♯Xであるか否かを判断する。♯V=♯Xである場合は(ステップS181の「Yes」)、正規表現Nをoutputして、処理を終了する。
 他方、ステップS181において、♯V=♯Xでない場合には(ステップS181の「No」)、s’∈T-Xを一つ選び(注:選び方に依らない)、(s,σ)を(s’,σ’)と置換して(ステップS182)、ステップS180に戻る。
 以上の操作が、ツリーの正規表現への変換処理である。
(4.語表現変換工程)
 上記図2の語表現変換工程(ステップS3)について詳細に説明する。
(4-1.正規表現から語表現への変換アルゴリズム)
(4-1-1.許容正規表現(admissible regular expression))
 以上説明したように、正規表現とss-saddle connection digramは1対1に対応する(単射)ことがわかるので、正規表現からその語表現へのアルゴリズムを与えることができる。ただし、ss-saddle connection diagramから正規表現を与える合成写像f(f(H))は全射ではないので、すべての正規表現にこのアルゴリズムを適用することはできず、可能なのは合成写像f(f(H))の像となっている正規表現のみである。このような正規表現を許容正規表現(admissible regular expression)と称する。許容正規表現を特徴づけるために、以下のような正規表現のfundamental type(基本形)とforward replacement(前方置換)を定義する。すなわち、許容正規表現とは、ルートつき,ラベル付き,および方向つきのツリーに対応する正規表現のうち,fundamental typeの正規表現に対して,正規表現のforward replacementを施してできる全ての正規表現をいう。
 まず、fundamental typeは、+φ、-φ、oφを含む正規表現である。図15は、5つのfundamental type Tとそれに対応する正規表現Nを示している。これらの正規表現はグラフ表現におけるルートを含む(ss-)saddle connection diagramに対応している。すなわち、+φ(+、+)や、oφ(o、+,....,+)といった正規表現はOBやICなる語表現を持つss-saddle connection diagramの正規表現となっている。
 次に、正規表現の形式的な文字列の変換ルールを定め、これらをforward replacement(前方置換)と呼ぶことにする。これらはラベル列b,...,bを正規表現中とラベルlの下に付け加える変換規則であり、以後、l→l(b,...,b)のような形で表現する。図16は、14個のforward replacementのリストである。例えば、A0+なる前方置換は、o→o(o,+)のように与えられる。
 各前方置換の内容を説明するために、その内部に流線の位相構造を持たないような連結成分を単純連結成分(simple connected component)と呼ぶことにする。単純連結成分はグラフ表現においては子連結成分を持たないような頂点に対応している。このとき、図16にあるforward replacementは操作A,A,B,B,Cを単純連結成分の中にあるss-orbitやclosed orbitに施したときにss-saddle connection diagramの局所流線構造とその正規表現がどのように変化するかを記述したものである。
 図17は、操作A,A,Cを施したときのforward replacementと流線構造の変化を表したものである。より詳細には、図17は、ss-saddle connection diagramにおける単純連結成分oあるいはoにあるss-orbitに操作A,A,Cを適用して得られる流線の局所構造の変化に対応する正規表現のforward replacementである。
 操作Aを単純連結成分oにあるss-orbitに施したとき、図17(a)に示すような、外向きhomoclinic saddle connectionと、図17(b)に示すような、内向きhomoclinic saddle connectionが得られる。この変化にともなって派生する正規表現のforward replacementを、それぞれA0+:o→o(o,+)およびA0-:o→o(o(-))のように与える。Forward replacement Aは、図17(c)に示すような操作Aを単純連結成分oにあるss-orbitに施したときに得られる正規表現の変換規則である。操作C,C,C(k,l≧0)を単純連結成分oのss-orbitに施したときに得られる4つのforward replacement C0+,C0-,C ,Cである(図17(d)-(g)参照)。
 さらに、操作B,B,Cを単純連結成分に含まれるclosed orbitに適用したときに図18に示すようなforward replacementが定義できる。図18は、(ss-)saddle connection diagramにおける単純連結成分σとσに含まれるclosed orbitに操作B,B,Cを適用したときに得られる流線構造の変化に伴って励起する正規表現のforward replacementを説明するための図である。
 操作Bを適用して、図18(a)と図18(b)に示すような、外向きおよび内向きhomoclinic saddle connectionが得られたときのforward replacementは、それぞれB0+:σ→σ(σ,σ)およびB0-:σ→σ(σ(μ))である。操作B,C,C,C,Cによって誘導されるforward replacementは、それぞれ図18(c)~図18(g)に示すようなB,C0+,C0-,Cσ ,Cである。最後に、forward replacement Cについては、sgn(σ)の符合に応じて、ラベル間の順序関係のルール(1)を考慮すると、二種類のforward replacementが得られることに注意する。
 このfundamental typeに対して帰納的にforward replacmentを適用して得られる正規表現の集合を以下のように定義する。
Figure JPOXMLDOC01-appb-M000001
 
 正規表現がO-type(あるいはI,II-type)であるとは、その正規表現が基本形O,OB,OB(あるいはIC or II)から上のルールで構成されるときをいう。操作によって得られる全ての局所流線構造が図17及び図18で網羅されているので、その構成方法から以下のことがわかる。
Figure JPOXMLDOC01-appb-M000002
 
 すべての構造安定なハミルトンベクトル場は、5個の操作A,A,B,B,Cをfundamental pattern 3個のO,I,IIに施して作られることに留意すれば以下を得る。このことは、数学的に正しいことが証明されている。この事実から次に与える語表現変換アルゴリズムに入力可能な許容正規表現は上の帰納的な方法によって定義されたものであることが示された。
Figure JPOXMLDOC01-appb-M000003
 
(4-1-2.語表現への変換アルゴリズム)
 以上から分かるように、これから与える語表現への変換アルゴリズムはfundamental typeに各操作A,A,B,B,Cを施した時に得られるforward replacementを順次繰り返して得られる正規表現だけを入力として受け取るこができる。
 したがって、与えられた許容正規表現について、各操作の作る特徴的流線構造を抽出したときにおこる正規表現の変化であるbackward replacement(後方置換)を定義することができれば、そのbackward replacementを使って正規表現を簡単にすると同時に抜き出された構造に対応する語表現の文字列を割り当てることによって語表現への変換が可能になる。
 すなわち、一つのss-saddle connectionを持つhomoclinic saddle point(あるいはss-∂-saddle connectionを持つ円境界)を消去する操作A(あるいはA)の逆操作であり、操作B(あるいはB)の逆はhomoclinic saddle connectionのペア(あるいは、ちょうど二つの∂-saddleを持つ円境界)を円境界に置き換える操作である。二つ以上の∂-saddleを持つ円境界についている∂-saddle connectionによって囲まれた単純連結成分を消去することが操作Cの逆操作である。ここで、backward replacementは、許容正規表現の消去されるべき構造に対応した文字列の並びの変化として表現できることに注意する。
 まずは、Bの逆操作によって定義されるbackward replacementのリストを与える。これらはfundamental type OB(k=0)およびforward replacement B,C(k=0,l≧0)の逆操作に対応する。Fundamental type OB:σφ→σφ(σ)に対するbackward replaementは単純にσφ(σ)→σφと書ける。見やすさのため、backward replaementによって消去される正規表現の文字列に、[数4],[数5],[数11],図19,及び図20では下線を付し、文章中では、[[ ]]を付す。図18(c)のforward replacement Bに対するbackward replacementは、その親連結成分σのまわりの構造に応じて異なっている。
 図19は、(ss-)saddle connection diagramにおける局所流線構造の変化に対応する正規表現のbackward replacementの一例を示す図である。
 もし、連結成分σが外向きhomoclinic saddle connectionで、図19(a)のように囲まれているとき、∂-saddle connectionに囲まれた単純連結成分を(ss-)saddle connection diagramから消去すれば、backward replacement σ([[σ]])→σを得る。一方、親連結成分σが内向きhomoclinic saddle connectionに囲まれているとき(図19(b))、そのbackward replacement σ(μ,[[σ]])→σ(μ)である。Forward replacement Cのbackward replacementは、ラベルに対する順序関係のルール(1)からσの符合に応じて以下のいずれかで与えられる。
Figure JPOXMLDOC01-appb-M000004
 
 Fundamental type OBとforward replacement B0+,B0-,C0+,C0-の逆操作は、(ss-)saddle connection diagramからhomoclinic saddle connectionのペアを円境界に置き換えることに対応する。Fundamental type OBに対するbackward replacementは、外向きhomoclinic saddle connectionに囲まれているときはσφ([[σ,σ]])→σφ、内向きhomoclinic saddle connectionに囲まれているときには、σφ([[σ(μ)]])→σφで与えられる。図18(a)のforward replacement B0+の中にある連結成分σが外向きhomoclinic saddle connectionに囲まれているとき(図19(c)),そのbackward replacementは、σ([[σ,σ]])→σで定義される。連結成分σが内向きhomoclinic saddle connectionに囲まれているときのbackward replacementは順序関係のルール(1)より図19(d)に示すように、+([[+,+]],-)→+(-)あるいは-(+,[[-,-]])→-(+)のいずれかである。
 同様の議論から、forward replacement B0-に対して三つのbackward replacement σ([[σ(μ)]])→σ,+([[+(-)]],-)→+(-),-(+,[[-(+)]])→-(+)が得られる。図18(d)と図18(e)のforward replacement C0+,C0-(l≧0)は以下のbackward replacementを誘導する。
Figure JPOXMLDOC01-appb-M000005
 
 ∂-saddle connectionに囲まれた単純連結成分を消去することに対応するbackward replacementはfundamental type IC(k≠0)およびforward replacement A,B(k≠0),C0+,C0-,C0+,C0-(l≠0),C ,Cσ ,C,(l+k≠0),Cの逆操作に対応している。操作Cは円境界にいくつでも適用できることから、そのbackward replacementは、円境界についている∂-saddle connectionに囲まれた単純連結成分を含む列(+,...,+)や(-,...,-)から、単純連結成分のどれか一つを消去することを表現する。Fundamental type IC(k≧1)のbackward replacementは以下のとおりである。
Figure JPOXMLDOC01-appb-M000006
 
 図17(c)に示したforward replacement A(k≠0)に対して、連結成分oが図19(e)や図19(f)のように外向きあるいは内向きhomoclinic saddle connectionに囲まれている場合、∂-saddle connectionに囲まれた単純連結成分+を消去して得られるbackward replacementはそれぞれ以下のようになる。
Figure JPOXMLDOC01-appb-M000007
 
 同様に、図18(c)のforward replacement B(k≧0)において、連結成分σが外向き、内向きのhomoclinic saddle connectionに囲まれているとき、そのbackward replacementは以下のようになる。
Figure JPOXMLDOC01-appb-M000008
 
 図17(d)~図17(f)及び図18(d)~図18(f)に示すforward replacement C ,Cσ ,C0+,C0-,C0+,C0-(l≧0)はすべて、∂-saddle connectionに囲まれた連結成分oやσに対して定義されているので、それに対するbackward replacementは自然に以下のように定義できる。
Figure JPOXMLDOC01-appb-M000009
 
 図17(g)と図18(g)のforward replacement C,C(l+k≠0)については、単純連結成分+,-のいずれを消去するかに応じて、また順序関係のルール(1)から、それぞれ二つの異なるbackward replacementが得られる。
Figure JPOXMLDOC01-appb-M000010
 
 図20は、操作B,B,Cの逆操作により誘導される正規表現のbackward replacement [B],[B]および[C]のリストを示している。最後に与えられた許容正規表現に対して、図20中のbackward replacementを可能な限り施していった場合に最終的に得られる正規表現の語表現は以下のようになる。このことは、数学的に正しいことが証明されている。
 Lemma 4.2:与えられた許容正規表現Nが[B],[B],[C]に含まれるすべてのbackward replacementがもはや施せなくなったとき、その正規表現に対する語表現は以下のいずれかである。
 1)Nの語表現はO
 2)Nの語表現はIA p+m c-2-p-mかIIA c-2のいずれかである。ただし、c,p,mは、Nに含まれるo,+,-のいずれかである。
 最後に、この変換アルゴリズムの実施例を示す。図21は、正規表現oφ(o(o(o(-(-,-)))),+(+))を持つ構造安定なハミルトンベクトル場のss-saddle connection diagramである。これにbackward replacement C(l=0),B,ICを順次適用すれば、正規表現は以下のようになる。
Figure JPOXMLDOC01-appb-M000011
 
 この結果、取り出した構造に応じて部分語列CBを得る。残った正規表現oφ(o(o(o(-))))の中にあるo+0-0の数を数えると、c=4,p=0,m=1を得るので、oが存在することと、p+m=1,c-2-p-m=1から、結果として語表現IACBを得ることができる。
(4-1-3.正規表現の語表現への変換処理)
 図22-A~図22-Cは、正規表現の語表現への変換処理を説明するためのフローチャートである。図22-A~図22-Cに示す正規表現の語表現への変換処理は、コンピュータ等の装置で実行可能である。以下、σ∈{+,-}を表し、μ=-σとする。図22-A~図22-Cにおいて、まず、inputを正規表現Nとする(ステップS184)。語Wを空集合φと置く(ステップS185)。
 ステップS186において、innermostなσを含むσφ(σ)またはinnermostなσを含むσ(σ)が正規表現Nに含まれるか否かを判断する。innermostなσを含むσφ(σ)またはinnermostなσを含むσ(σ)が正規表現Nに含まれる場合は(ステップS186の「Yes」)、ステップS187に移行する。
 ステップS187において、σφ(σ)またはσ(σ)に含まれるinnermostなσを正規表現Nから抜き、WをBWと置換して(ステップS188)、ステップS186に戻る。
 他方、ステップS186において、innermostなσを含むσφ(σ)またはinnermostなσを含むσ(σ)が正規表現Nに含まれない場合は(ステップS186の「No」)、ステップS189に移行する。
 ステップS189において、正規表現Nの中にinnermostな+を含む+(+,-,...,-)が存在するか否かを判断する。正規表現Nの中にinnermostな+を含む+(+,-,...,-)が存在する場合は(ステップS189の「Yes」)、このinnermostな+を正規表現Nから抜き出し(ステップS190)、WをBWと置換して(ステップS188)、ステップS186に戻る。
 他方、ステップS189において、正規表現Nの中にinnermostな+を含む+(+,-,...,-)が存在しない場合は(ステップS189の「No」)、ステップS191に移行する。
 ステップS191において、正規表現Nの中にinnermostな-を含む-(+,...,+,-)が存在するか否かを判断する。正規表現Nの中にinnermostな-を含む-(+,...,+,-)が存在する場合は(ステップS191の「Yes」)、このinnermostな-を正規表現Nから抜き出し(ステップS192)、WをBWと置換し(ステップS188)、ステップS186に戻る。
 他方、ステップS191において、正規表現Nの中にinnermostな-を含む-(+,...,+,-)が存在しない場合は(ステップS191の「No」)、ステップS193に移行する。
 ステップS193において、連続する2つのinnermostなσ,σが正規表現Nに含まれるか否かを判断する。連続する2つのinnermostなσ,σが正規表現Nに含まれる場合は(ステップS193の「Yes」)、このようなσ,σを正規表現Nから抜き出し(ステップS194)、WをBWと置換し(ステップS195)、ステップS193に戻る。
 他方、ステップS193において、連続する2つのinnermostなσ0,σが正規表現Nに含まれない場合は(ステップS193の「No」)、ステップS196に移行する。
 ステップS196において、innermostなσ(μ)が正規表現Nに含まれるか否かを判断する。innermostなσ(μ)が正規表現Nに含まれる場合は(ステップS196の「Yes」)、このようなσ(μ)を正規表現Nから抜き(ステップS197)、WをBWと置換し(ステップS195)、ステップS193に移行する。
 ステップS196において、innermostなσ(μ)が正規表現Nに含まない場合は(ステップS196の「No」)、ステップS198に移行する。
 ステップS198において、正規表現Nに、連続する2つ以上のσ,...,σの中にinnermostなσが含まれているか否かを判断する。正規表現Nに、連続する2つ以上のσ,...,σの中にinnermostなσが含まれている場合は(ステップS198の「Yes」)、このようなσを正規表現Nから抜き(ステップS199)、WをCWと置換して(ステップS200)、ステップS198に戻る。
 他方、ステップS198において、正規表現Nに、連続する2つ以上のσ,...,σの中にinnermostなσが含まない場合は(ステップS198の「No」)、ステップS201に移行する。
 ステップS201において、oの子でinnermostなσが正規表現Nに存在するか否かを判断する。oの子でinnermostなσが正規表現Nに存在する場合は(ステップS201の「Yes」)、このようなσを正規表現Nから抜き(ステップS199)、WをCWと置換して(ステップS200)、ステップS198に戻る。
 他方、ステップS201において、oの子でinnermostなσが正規表現Nに存在しない場合は(ステップS201の「No」)、ステップS202に移行する。
 ステップS202において、σの子でinnermostなμが正規表現Nに存在するか否かを判断する。σの子でinnermostなμが正規表現Nに存在する場合は(ステップS202の「Yes」)、このようなμを正規表現Nから抜き(ステップS203)、WをCWと置換し(ステップS200)、ステップS198に戻る。
 他方、ステップS202において、σの子でinnermostなμが存在しない場合は(ステップS202の「No)、ステップS204に移行する。
 ステップS204において、正規表現Nの中にσが存在するか否かを判断する。正規表現Nの中にσが存在する場合は(ステップS204の「Yes」)、ステップS186に戻る。
 他方、ステップS204において、正規表現Nの中にσが存在しない場合は(ステップS204の「No」)、ステップS205に移行する。
 ステップS205において、正規表現Nの中にσ(μ)が存在するか否かを判断する。正規表現Nの中にσ(μ)が存在する場合は(ステップS205の「Yes」)、ステップS186に戻る。
 他方、ステップS205において、正規表現Nの中にσ(μ)が存在しない場合は(ステップS205の「No」)、ステップS206に移行する。
 ステップS206において、正規表現Nに連続する2つのσ,σが存在するか否かを判断する。正規表現Nに連続する2つのσ,σが存在する場合は(ステップS206の「Yes」)、ステップS186に戻る。
 他方、ステップS206において、正規表現Nに連続する2つのσ,σが存在しない場合は(ステップS206の「No」)、ステップS207に移行する。
 ステップS207において、正規表現Nがσφであるか否かを判断する。正規表現Nがσφの場合は(ステップS207の「Yes」)、WをOWと置換し(ステップS208)、語Wをoutputして(ステップS212)、処理を終了する。
 他方、ステップS207において、正規表現Nがσφでない場合は(ステップS207の「No」)、ステップS209に移行する。
 ステップS209において、正規表現Nがoを含むか否かを判断する。正規表現Nがoを含む場合は(ステップS209の「Yes」)、WをIA p+m c-2-p-mWと置換し(ステップS210)、語Wをoutputして(ステップS212)、処理を終了する。ただし,c,p,mはNに含まれるo,+,-の数とする。
 他方、ステップS208において、正規表現Nがoを含まない場合は(ステップS208の「No」)、WをIIA c-2Wと置換し(ステップS211)、語Wをoutputして(ステップS212)、処理を終了する。
 以上の処理が、正規表現の語表現への変換処理である。
 つづいて、上述した本実施の形態の流れパターンの正規表現作成方法をコンピュータにより実施するための装置構成について説明する。なお、以上の本実施の形態による方法を、人またはコンピュータにより実施してもよく、以下の実施形態による処理等を人により実施する場合に用いてもよいものである。
[正規表現作成装置の構成]
 次に、本実施の形態に係る正規表現作成装置の構成について図23を参照して説明する。図23は、本実施の形態が適用される正規表現作成装置100の一例を示すブロック図であり、該構成のうち本実施形態に関係する部分のみを概念的に示している。
 図23に示すように、本実施形態における正規表現作成装置100は、概略的に、制御部102と記憶部106を少なくとも備え、本実施形態において、更に、入出力制御インターフェース部108と通信制御インターフェース部104を備える。ここで、制御部102は、正規表現作成装置100の全体を統括的に制御するCPU等である。また、通信制御インターフェース部104は、通信回線等に接続されるルータ等の通信装置(図示せず)に接続されるインターフェースであり、入出力制御インターフェース部108は、入力装置112や出力装置114に接続されるインターフェースである。また、記憶部106は、正規表現作成装置100で上述の正規表現作成方法を実行するために必要な各種のデータ(例えば、図4~図9、図11~図12,図14-A,図15~図21に示すものをデータ化したもの等)、データベース、およびテーブル等を格納する装置である。これらのデータ、データベース、およびテーブル等は、制御部102が正規表現作成方法を実行する際に必要により参照される。これら正規表現作成装置100の各部は任意の通信路を介して通信可能に接続されている。更に、この正規表現作成装置100は、ルータ等の通信装置および専用線等の有線または無線の通信回線を介して、ネットワーク130に通信可能に接続されている。
 記憶部106に格納される各種のデータ、データベース、およびテーブル(シミュレーション結果ファイル106a、流線図ファイル106b、および、正規表現ファイル106c等)等は、固定ディスク装置等のストレージ手段である。例えば、記憶部106は、各種処理に用いる各種のプログラム、テーブル、ファイル、データベース、および、ウェブページ等を格納する。
 これら記憶部106の各構成要素のうち、シミュレーション結果ファイル106aは、シミュレーション部102aにより数理的にシミュレーションされた、シミュレーション結果を示すデータを記憶するシミュレーション結果記憶手段である。例えば、シミュレーション結果ファイル106aは、構造物の形状を示す設計変数の値や、その構造物に対する所定の流体(海流や気流等)の流体力学的シミュレーション結果(各空間座標における流体の圧力や流れの向き等)を示すデータであってもよい。なお、シミュレーション結果ファイル106aは、風洞実験等の実験室内でのモデル計測などを通じて事前に入力装置112を介して入力されたデータをシミュレーション結果として記憶してもよい。
 また、流線図ファイル106bは、流線図等の流線を示すデータを記憶する流線データ記憶手段である。例えば、流線図ファイル106bに記憶される流線データは、シミュレーション結果を示すデータに基づいて流線解析部102bにより解析された流線を示すデータであってもよい。
 また、正規表現ファイル106cは、正規表現データを記憶する正規表現記憶手段である。例えば、正規表現ファイル106cに記憶される正規表現データは、文字列等である。
 図23に戻り、入出力制御インターフェース部108は、入力装置112や出力装置114の制御を行う。ここで、出力装置114としては、モニタ(家庭用テレビを含む。)の他、スピーカを用いることができる(なお、以下においては出力装置114をモニタとして記載する場合がある)。また、入力装置112としては、キーボード、マウス、およびマイク等を用いることができる。
 また、図23において、制御部102は、OS(Operating System)等の制御プログラムや、各種の処理手順等を規定したプログラム、および、所要データを格納するための内部メモリを有する。そして、制御部102は、これらのプログラム等により、種々の処理を実行するための情報処理を行う。制御部102は、機能概念的に、シミュレーション部102a、流線解析部102b、グラフ表現作成部102c、正規表現作成部102d、および、語表現変換部102eを備える。
 このうち、シミュレーション部102aは、構造物に対する流体のシミュレーションを行うシミュレーション手段であり、例えば、上述の設計パラメータ選択工程(図24参照)および第2の設計パラメータ選択工程(図29等参照)等を実行する。ここで、シミュレーション部102aは、2次元平面におけるシミュレーションに限られず、3次元空間における流体のシミュレーションを行ってもよいものである。また、シミュレーション部102aは、公知の最適化手法を用いて、構造物の最適化を行ってもよい。例えば、シミュレーション部102aは、構造物の形状を決定する設計変数を、焼きなまし法や遺伝的アルゴリズム法等を用いて繰り返し変化させながら、当該構造物に対する流体シミュレーションを行って、適切な構造物の形状(例えば、水流に対する抵抗の少ない橋脚の形状等)を求めてもよい。なお、本実施形態において、シミュレーション部102aは、シミュレーション結果を示すデータをシミュレーション結果ファイル106aに格納する。例えば、シミュレーション部102aは、構造物の形状を示す設計変数の値や、その構造物に対する所定の流体(海流や気流等)の流体力学的シミュレーション結果(各空間座標における流体の圧力や流れの向きや抵抗等)を示すデータを格納してもよい。
 また、流線解析部102bは、流線解析を行う流線解析手段である。ここで、流線解析部102bは、シミュレーション部102aによるシミュレーション結果に対し流線解析を行って流線図を導出してもよい。例えば、流線解析部102bは、シミュレーション結果ファイル106aに記憶された、数値シミュレーションや実験のデータから、公知の手法を用いて流線図を作成する。具体的には、流線解析部102bは、数値シミュレーション結果から、saddle pointや1-source-sink pointなどをすべて計算した後、その点における流れ関数の値と同じ値を持つ流れ関数の等高線をすべて描画し、また、境界(boundary)上の流れ関数の値と同じ値を持つ流れ関数の等高線をすべて描画することにより流線図の作成が可能となる。なお、3次元のシミュレーション結果の場合、流線解析部102bは、構造物における断面における2次元のデータに変換してから、流線解析を行ってもよい。断面とする平面は任意であるが好適には、流線解析部102bは、流体の流れ方向(一様流)の方向に沿った断面で2次元データに変換してもよい。例えば、列車や自動車や航空機等の乗り物においては、進行方向に沿って断面を生成してもよい。また、流線解析部102bは、Computational Homology(非特許文献1)に記載の技術等を用いて、流れ場から条件を満たす特徴的な構造を抽出してもよい。なお、本実施形態において、流線解析部102bは、作成した流線図データを流線図ファイル106bに格納する。
 また、グラフ表現作成部102cは、上述のグラフ表現作成工程(ステップS1)を実行するためのグラフ表現作成手段であり、例えば、図10-A~図10-Eに示す、O系列におけるsaddle connection diagramのツリーへの変換処理、および図13-A~図13-Dに示す、I,II系列におけるss-saddle connection diagramのツリーへの変換処理等を実行することにより、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンに1対1に対応するグラフ表現を作成する。ここで、グラフ表現は、流れパターンで規定される構造安定なハミルトンベクトル場Hに対して、固有のルート付き、ラベル付き、向き有りツリーT=(V,E)を割り当て(但し、Vは頂点と呼ばれる点の集合、Eは、頂点の間を結ぶエッジの集合である)、平面グラフとして可視化したものであってもよい。また、グラフ表現は、親の頂点をv、その子の頂点をw,親の頂点vに割り当てられたラベルをl(v)、子の頂点wに割り当てられたラベルをl(w),vの子頂点集合Γ(v)とした場合、vの子頂点集合Γ(v)を所定の順序関係のルールに従って並び替え、w∈Γ(v)について、l(v)からl(w)への矢印を左から右に並べて描画されたものを含んでいてもよい。また、流れパターンは、(1)一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、(2)前記一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、(3)二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンO、の1又は複数であってもよい。また、流れパターンが吸い込み湧き出し対を持つ場合は、流れパターンをその吸い込み湧き出し対に最も近い反時計回りのss-orbitを含む領域が一番外側領域になるように変換し、変換した流れパターンの(ss-)saddle connection diagramに現れる軌道を領域全体から抽出し、領域全体から(ss-)saddle connection diagramに現れる軌道をすべて除外して得られる連結成分に頂点を設定し、一番外側にある連結成分をルートとし、カレント成分をルートに設定し、カレント成分と互いに境界を接するような連結成分をカレント成分の子として、境界にあたる軌道に応じてラベルを割り当て,当該ラベルを所定の順序関係に従って並べ、カレント成分の子をカレント成分に設定して、子がなくなるまで繰り返すことにしてもよい。また、かかる流れパターンは、パターン語の1又は複数から開始し、流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を付与する操作を、穴の数がN個となるまで繰り返して作成された流れパターン図であることにしてもよい。
 また、正規表現作成部102dは、上述の正規表現作成工程(ステップS2)を実現するための正規表現作成手段であり、例えば、図14-Bに示す、ツリーの正規表現への変換処理等を実行することにより、グラフ表現作成部102cで作成されたグラフ表現から正規表現を作成する。
 また、語表現変換部102eは、上述の語表現変換工程(ステップS3)を実現するための語表現作成手段であり、例えば、図22-A~図22-Cに示す正規表現の語表現への変換処理等を実行して、正規表現作成部102dで作成された正規表現を語表現に変換する。語表現は、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語に対して、上記流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を、追加された穴の数だけ付与することにより形成された記号語であることにしてもよい。
 以上が、本実施形態における正規表現作成装置100の構成の一例である。なお、正規表現作成装置100は、ネットワーク130を介して外部システム120に接続されてもよい。この場合、通信制御インターフェース部104は、正規表現作成装置100とネットワーク300(またはルータ等の通信装置)との間における通信制御を行う。すなわち、通信制御インターフェース部104は、他の端末と通信回線を介してデータを通信する機能を有する。また、ネットワーク130は、正規表現作成装置100と外部システム200とを相互に接続する機能を有し、例えば、インターネット等である。
 また、外部システム120は、ネットワーク130を介して、正規表現作成装置100と相互に接続され、シミュレーション結果データや流線図データ等の各種データに関する外部データベースや、接続された情報処理装置に語表現方法を実行させるためのプログラム等を提供する機能を有する。
 ここで、外部システム120は、WEBサーバやASPサーバ等として構成していてもよい。また、外部システム120のハードウェア構成は、一般に市販されるワークステーション、パーソナルコンピュータ等の情報処理装置およびその付属装置により構成していてもよい。また、外部システム120の各機能は、外部システム120のハードウェア構成中のCPU、ディスク装置、メモリ装置、入力装置、出力装置、通信制御装置等およびそれらを制御するプログラム等により実現される。
[語表現と正規表現を使った流体中の物体の設計方法]
 例えば、橋や橋脚など流体中におかれた物体(以下設計対象と呼ぶ)の形状や配置、またその周囲の流れの制御などを通じて設計の対象にとって最適なパラメータを定めるような設計手法について,語表現と正規表現を使ってどのように設計を行うかを説明する。
 前提1:設計対象には変更可能な設計パラメータがある(形状・配置・流れの制御装置など)。
 前提2:設計対象を設計する上で、問題に応じた「最適な状態」が設定されており、その最適状態が流線構造の特徴として記述できているものとする。例えば、剥離渦を閉じこめる状態を最適状態とすると、翼であれば揚力の最大化,橋脚であれば抗力の最小化などが期待できる。
 前提3:これらの設計において、設計パラメータを変えて実験および数値計算(そしてそれにおいて得られる流線パターンの語表現・正規表現)ができるものとする。以下、二次元の流れに限定する。三次元の場合は断面をとって考えたりして二次元化することなどで設計が可能になる場合もある(三次元の場合はすべての場合でできるとは限らないことに注意する)。
 従来の最適化設計においては,初期状態から始めて、設計パラメータを経験や試行錯誤、あるいは既知の最適化手法などによる最適化が行われている。ある種の最適化が実現できても、それは局所的な最適化であったり物理的に不安定であったりして実際の設計に耐えうるかどうか分からないところがある。
 本実施の形態による設計パラメータの探索では、かかる状況に対して、最初から理想的な状況を流線パターンの語表現や正規表現として設定しておいて、それが実現されるようパラメータ領域を探索することができるようになる。こうして最適な設定を実現する設計パラメータの「候補」をできるだけ多くかつ速やかに探すことができるようになる。それらの候補から始めて、既知の最適化手法を行うことで多くのかつ実現性の高い設計パラメータが得られる可能性が高くなる。
 図24は、語表現と正規表現を使用した流体中の物体の設計方法、すなわち、上述の設計パラメータ選択工程を説明するためのフローチャートである。図24に示すフローチャートは、シミュレーション部102aにより実行される。
 図24において、まず、とりうる設計パラメータの上限と下限を定める(ステップS311)。次に、上記ステップS311で定めたパラメータ領域(設計パラメータの上限と下限で規定される領域を「パラメータ領域」と称する)を分割して、パラメータ領域からN個のパラメータの組み合わせサンプル(設計パラメータP(i=1,・・・,N))を選択する(ステップS312)。なお、サンプルの選択は、これまでの経験や計算結果などを踏まえて適当に選んでも良いし、先験的な情報がないときは、パラメータ領域を等分してサンプリングしてもよい。
 以下では、渦閉じこめ状態を表す語表現として「C」が語表現に含まれることを設計の目的とする例を説明する。各設計パラメータP(i=1,・・・,N)に対してそれぞれ流れの実験や数値計算を行う(ステップS313)。得られた各結果に語表現と正規表現の割り当てを行う(ステップS314)。その語表現が最適状態を表すICの部分語を持つ、あるいはそれに対応する正規表現を持つ構造表現か否かを判断する(ステップS315)。その語表現が最適状態を表すICの部分語を持つ,あるいはそれに対応する正規表現を持つ構造表現の場合は(ステップS315の「Yes」)、設計パラメータの候補として採用する(ステップS316)。他方、その語表現が最適状態を表すICの部分語を持つ、あるいはそれに対応する正規表現を持つ構造表現でない場合は(ステップS315の「No」)、設計パラメータの候補には採用しない。
 なお、ICなる部分語表現を持たなくても,正規表現と語表現を用いて流体遷移経路(遷移パターン)を取得する方法を使って、一回の遷移でICやそれに対応する正規表現を持つ構造に移るパターンであれば、これも設計パラメータの候補として採用することにしてもよい。
 ステップS317では、設計パラメータの候補があるか否かを判断し(ステップS317)、設計候補のパラメータがある場合には(ステップS317の「Yes」)、設計パラメータの候補に対して既存の最適化手法を適用する(ステップS319)。なお、既存の最適化手法を適用するだけでなく,最適化過程の各段階において、常に語表現や正規表現を計算してモニタリングすることで、定量的な流れの最適化に加えて、理想的な流線構造の最適化も同時に行うような最適化設計も可能である。
 他方、設計パラメータの候補がない場合には(ステップS317の「No」)、パラメータ領域を更に細かく分割して(ステップS318)、ステップS313~S317の処理を実行する。
 図25及び図26は、上記の語表現と正規表現を使用した流体中の物体の設計方法を説明するための具体例の一例を説明するための図である。図25において、一様流中におかれた物体200の背後に渦構造が二つ(渦1,渦2)存在するような状況があったとする。一様流の速さ(あるいは物体の進行速度)をUとし、物体200は角速度Gで回転できるものとする。渦の強さと位置は流れにおいて既に与えられているとして、このうち渦1を「閉じこめて」、物体200にかかる揚力を最大化し、渦2は流れに沿って物体から遠ざかることができるような装置を、「UやGのパラメータ」を探索して最適化することを考える。渦閉じこめ状態の語表現は「C」によって遠ざかる渦は「A」によって表現できるので、ターゲットとなる状況の語表現としてはIACを達成するものを本設計手法の適用によりパラメータ領域を絞り込む。
 本設計パラメータについて,Uは0から1.1まで変化でき,Gは-1.6×2πから1.6×2πまで変化できるものとする(上記ステップS311で設計パラメータの上限、下限を定める)。これらのパラメータ区間を5等分して設計パラメータ領域を分割し(上記ステップS312)する。そして、上記ステップS313~S319を実行した結果が、図26の(a)~(f)である。このパラメータ領域のうち、この設計手法によって最適化設計は「0.5≦U≦0.9,G=-0.8×2π」(図26の(a)~(c)に対応)および「0.3≦U≦0.7,G=0」(図26の(d)~(f)に対応)の間で達成されている。これらのパラメータ候補の絞り込みにより、さらにこれらのパラメータ領域の分割を繰り返して、揚力などの定量的な量の最大化をするような最良パラメータ領域を探すことで、渦閉じ込めをして揚力を最大化するための物体200の速度と各回転速度を知ることができる。
[語表現と正規表現を使った流体中の物体の設計方法2]
 語表現と正規表現を使用して、流体中に置かれた物体(以下、「設計対象」と呼ぶ)のまわりの流れを最適化する設計手法について説明する。以下では、語表現と正規表現を使用して、一様流中に置かれた物体のまわりにできる流線パターン(現在のパターン)が、所望の流線パターン(目的のパターン)となるように、物体の設計パラメータ(物理パラメータ)を設計する手法(上述の第2の設計パラメータ選択工程)の一例について説明する。以下の例では、一様流中に平板翼を置いた場合の流線パターンについて説明する。なお、第2の設計パラメータ選択工程は、正規表現のみを使用しても原理上実行可能であるが、正規表現と語表現を使用することで、より好適な設計パラメータを設計することが可能となる。
 図27及び図28、語表現と正規表現を使用した流体中の物体まわりの流れの最適化手法を説明するための説明図である。図29は、語表現と正規表現を使用した流体中の物体のまわりの流れの最適化手法、すなわち、上述の第2の設計パラメータ選択工程を説明するためのフローチャートである。図29に示すフローチャートは、シミュレーション部102aにより実行される。
 図27は、一様流Uに対して、一定の角度で平板翼300を置いた場合を示している。図28は、図27において、平板翼300のまわりに形成される流線パターンとその流線パターンの語表現及び正規表現を示している。図28(a)は一様流U中におかれた平板翼300のまわりにできた流れ(「流れ」は、流体実験又は数値計算で算出することができる)から流線パターンを抽出したもの(「パターンA」と称する)である。図28(c)は、回転流れ構造を平板翼300の上に閉じ込めるパターン(「パターンB」と称する)である。本手法では、平板翼300の設計パラメータを好適に変更して、パターンA(現在のパターン)からパターンB(目的のパターン)を実現するような設計パラメータを見つけ出し、平板翼300にかかる揚力を増加させる手法について説明する。
 図29において、まず、ステップA(ステップS400~S402)において、パターンAとパターンBに対して、語表現(特願2012-203601号(PCT/JP2013/070939号))と本実施の形態のアルゴリズムを用いて正規表現を与える。具体的には、パターンAには、IACCBという語表現と、oφ(o(o,-(-,-),+,+))という正規表現Nを付与し、パターンBには、語表現IACCBとoφ(o(o,+,-(-,-,+)))という正規表現N Bを付与する(ステップS400、図28(a)、(c)参照)。なお、パターンAとパターンBの語表現は全く同じなので、従来技術では区別ができないが、正規表現は異なるためこれを用いることにより2つの流れは区別が可能となり,これらを用いることによって、流れの最適設計の目的形状を特定することが可能となる。
 次に、パターンAの正規表現NとパターンBの正規表現Nが同じであるか否かを判断する(ステップS401)。パターンAの正規表現NとパターンBの正規表現Nが同じである場合には(ステップS401の「Yes」)、パスP=φを出力する(ステップS402)。パターンAの正規表現NとパターンBの正規表現N Bが異なる場合には(ステップS401の「No」)、ステップS403に移行する。
 ステップBでは、パターンAとパターンBの遷移の経路の決定する(ステップS403)。具体的には、流れパターンの遷移同定プログラム(特願2013-230678号(PCT/JP2014/079512))を用いて,図28(a)と図28(c)の間の同じ正規表現を高々1つしか持たない経路を決定する(ステップS403)。より具体的には、同じ頂点を高々一度しか通らない、パターンAの正規表現NからパターンBの正規表現Nへのパスを全てリストアップし、N個の経路R(i=1,2,…,N)を得る。ただし,R,…,Rは,経路の長さが短い順に並べてあるものとする(length(R)≦length(Ri+1))。各経路R(i=1,2,…,N)には中間流線パターンPij(j=1,2,…,M)があるものとする。但し、便宜上、Pi0=N,PiMi=Nであるとする。
 ステップCでは、経路Rに沿って遷移するように設計パラメータを変更する(ステップS404~S415)。以下、N個の経路R(i=1,2,…,N)に対して、i=1,…,Nの順に(経路の長さが短い順に)、経路に沿って遷移するように設計パラメータを変更して、最適な設計パラメータを求める。すなわち、パターンAからパターンBに至るN個の経路R(i=1,2,…,N)のうち、設計パラメータを変更して、パターンBの正規表現と一致することになるいずれかの経路を決定し、当該パターンBの正規表現と一致することになる変更した設計パラメータを最適な設計パラメータとして出力する。
 経路R上にあるパターンPij(但し、j=1,2,…,M)の設計パラメータを中心にして、変更可能な設計パラメータの範囲を指定し、それをK分割して、それぞれのパラメータ値をS(k=1,….K)とする。さらに、パターンPijの設計パラメータのパラメータ値をSに変更して、実験又は数値計算により得られるパターンの正規表現をPij(S)と表す。
 最終的に、Pij(S)がパターンBの語表現と正規表現N(=PiMi)に一致した場合に、このパラメータ値を最適設計を達成する設計パラメータとして出力して、当該フローを終了する。
 Sを変更して、Pij(S)が経路R上の正規表現Pi(j+1)に一致した場合に、この設計パラメータを選択して、さらに、jを「1」インクリメントして、当該選択した設計パラメータを中心にして、変更可能な設計パラメータの範囲を指定し、それをK分割して、それぞれのパラメータ値をS(k=1,….K)として、Pij(S)がパターンBの語表現と正規表現N(=PiMi)に一致するまで処理を繰り返す。
 Sを変更しても、Pij(S)が経路R上の正規表現Pi(j+1)に一致しない場合は、この経路上での遷移は不可能として、次の経路についてSを変更する。
 まず、i=1、j=1、k=1とする(ステップS404~S406)。次に、Pij(S)=Pi(j+1)であるか否かを判断する(ステップS407)。Pij(S)=Pi(j+1)である場合には(ステップS407の「Yes」)、j←j+1とし(ステップS413)、j=Mであるか否かを判断する(ステップS414)。j=Mである場合(ステップS414の「Yes」)、すなわち、Pij(S)=PiMi(=N)となる場合は、この設計パラメータが最適な設計パラメータであるので、パスP=R及びこの設計パラメータを出力して(ステップS415)、当該フローを終了する。j=Mでない場合は(ステップS414の「No」)、この設計パラメータを選択して、ステップS406に戻り、選択した設計パラメータを中心にして、変更可能な設計パラメータの範囲を指定し、それをK分割して、それぞれのパラメータ値をS(k=1,….K)として、Pij(S)がパターンBの語表現と正規表現N(=PiMi)に一致するまで処理を繰り返す。
 Pij(S)=Pi(j+1)でない場合には(ステップS407の「No」)、k←k+1とし(ステップS408)、k≦Kであるか否かを判断する(ステップS409)。k≦Kである場合には(ステップS409の「Yes」)、ステップS407に戻る。他方、k≦Kでない場合には(ステップS409の「No」)、i←i+1として(ステップS410)、i≦Nであるか否かを判断する(ステップS411)。i≦Nである場合には(ステップS411の「Yes」)、ステップS405に戻り、次の経路に対して探索を行う一方、i≦Nでない場合には(ステップS411の「No」)、全ての経路について最適化が不可能であるので、「不可能」を出力する(ステップS412)。
 一例として、ステップBで得られた経路の1つとして、最短経路である図28(b)のパターン(語表現ICCCB,正規表現oφ(o(-(-,-,+,+))))を介した場合(経路の長さM=1の場合))について説明する。ステップCにおいて,図28(a)のパターンAから設計パラメータを変更して、図28(b)の語表現と正規表現になるものを探索する。変更可能な設計パラメータを全て探索しても図28(b)となるパターンが見つからない場合は,この経路に沿った最適設計が不可能と判断して,ステップBで取得した、図28(b)の最短経路以外の経路で同じことを繰り返す。
 もし、いきなり目的とする図28(c)のパターンBを達成する設計パラメータが見つかった場合は、これが最適な設計パラメータであるので、この設計パラメータを出力してフローを終了する(上記ステップS415)。一方、図28(b)のパターンを達成する設計パラメータが見つかった場合は、この設計パラメータを選択して、さらに、この設計パラメータまわりで変更可能な設計パラメータの範囲を指定し、それを分割して、それぞれに対して実験又は数値計算を行って図28(c)のパターンBになるものがあるかどうかを探索する。図28(c)のパターンBになる設計パラメータが見つかった場合は、これが最適な設計パラメータであるので、この設計パラメータを出力してフローを終了する(上記ステップS415)。
 変更可能な設計パラメータの範囲を全て探索しても見つからない場合は、この経路での最適化は不可能と判断し、ステップBで得られた別経路を探索する。ステップBで取得した全ての経路を探索しても、図28(c)のパターンBになる設計パラメータが見つからない場合は、最適化不可能と判断する(上記ステップS412)。
[他の実施形態]
 さて、これまで本発明の実施形態について説明したが、本発明は、上述した実施形態以外にも、請求の範囲に記載した技術的思想の範囲内において種々の異なる実施形態にて実施されてよいものである。
 例えば、正規表現作成装置100がスタンドアローンの形態で処理を行う場合を一例に説明したが、正規表現作成装置100は、クライアント端末からの要求に応じて処理を行い、その処理結果を当該クライアント端末に返却するようにしてもよい。
 また、実施形態において説明した各処理のうち、自動的に行われるものとして説明した処理の全部または一部を手動的に行うこともでき、あるいは、手動的に行われるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。
 このほか、上記文献中や図面中で示した処理手順、制御手順、具体的名称、各処理の登録データや検索条件等のパラメータを含む情報、画面例、データベース構成については、特記する場合を除いて任意に変更することができる。
 また、正規表現作成装置100に関して、図示の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。
 例えば、正規表現作成装置100の各装置が備える処理機能、特に制御部102にて行われる各処理機能については、その全部または任意の一部を、CPU(Central Processing Unit)および当該CPUにて解釈実行されるプログラムにて実現してもよく、また、ワイヤードロジックによるハードウェアとして実現してもよい。尚、プログラムは、後述する記録媒体に記録されており、必要に応じて正規表現作成装置100に機械的に読み取られる。すなわち、ROMまたはHDなどの記憶部106などには、OS(Operating System)として協働してCPUに命令を与え、各種処理を行うためのコンピュータプログラムが記録されている。このコンピュータプログラムは、RAMにロードされることによって実行され、CPUと協働して制御部102を構成する。
 また、このコンピュータプログラムは、正規表現作成装置100に対して任意のネットワーク300を介して接続されたアプリケーションプログラムサーバに記憶されていてもよく、必要に応じてその全部または一部をダウンロードすることも可能である。
 また、本発明に係るプログラムを、コンピュータ読み取り可能な記録媒体に格納してもよく、また、プログラム製品として構成することもできる。ここで、この「記録媒体」とは、メモリーカード、USBメモリ、SDカード、フレキシブルディスク、光磁気ディスク、ROM、EPROM、EEPROM、CD-ROM、MO、DVD、および、Blu-ray(登録商標)Disc等の任意の「可搬用の物理媒体」を含むものとする。
 また、「プログラム」とは、任意の言語や記述方法にて記述されたデータ処理方法であり、ソースコードやバイナリコード等の形式を問わない。なお、「プログラム」は必ずしも単一的に構成されるものに限られず、複数のモジュールやライブラリとして分散構成されるものや、OS(Operating System)に代表される別個のプログラムと協働してその機能を達成するものをも含む。なお、実施形態に示した各装置において記録媒体を読み取るための具体的な構成、読み取り手順、あるいは、読み取り後のインストール手順等については、周知の構成や手順を用いることができる。
 記憶部106に格納される各種のデータベース等(シミュレーション結果ファイル106a、流線図ファイル106b、語表現ファイル106c等)は、RAM、ROM等のメモリ装置、ハードディスク等の固定ディスク装置、フレキシブルディスク、および、光ディスク等のストレージ手段であり、各種処理やウェブサイト提供に用いる各種のプログラム、テーブル、データベース、および、ウェブページ用ファイル等を格納する。
 また、正規表現作成装置100は、既知のパーソナルコンピュータ、ワークステーション等の情報処理装置として構成してもよく、また、該情報処理装置に任意の周辺装置を接続して構成してもよい。また、語表現作成装置100は、該情報処理装置に本発明の方法を実現させるソフトウェア(プログラム、データ等を含む)を実装することにより実現してもよい。
 更に、装置の分散・統合の具体的形態は図示するものに限られず、その全部または一部を、各種の付加等に応じて、または、機能負荷に応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。すなわち、上述した実施形態を任意に組み合わせて実施してもよく、実施形態を選択的に実施してもよい。
 以上説明したように、本実施の形態によれば、位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成方法であって、前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、を含むこととしたので、流れパターンと1対1に対応させることが可能な新たな表現方法を提供することが可能となるという効果を奏する。
 また、グラフ表現は、流れパターンで規定される構造安定なハミルトンベクトル場Hに対して、固有のルート付き、ラベル付き、及び向き有りのツリーT=(V,E)を割り当て(但し、Vは頂点と呼ばれる点の集合、Eは、頂点の間を結ぶエッジの集合である)、平面グラフとして可視化したものであるので、構造安定なハミルトンベクトル場に対して、固有のルート付き、ラベル付き、及び向き有りのツリーを平面グラフとして可視化したものをグラフ表現として提供することが可能となるという効果を奏する。
 また、本実施の形態によれば、グラフ表現は、親の頂点をv、その子の頂点をw,親の頂点vに割り当てられたラベルをl(v)、子の頂点wに割り当てられたラベルをl(w),vの子頂点集合Γ(v)とした場合、vの子頂点集合Γ(v)を所定の順序関係のルールに従って並び替え、w∈Γ(v)について、l(v)からl(w)への矢印を左から右に並べて描画することとしたので、親子の連結関係が視認しやすいグラフ表現を提供することが可能となるという効果を奏する。
 また、本実施の形態によれば、前記流れパターンは、(1)一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、(2)前記一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、(3)二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンOの1又は複数であることとしたので、基本となる全ての流れパターンに対して正規表現を付与することができ、具体的な流れパターンの正規表現を扱えるようになるという効果を奏する。
 また、本実施の形態によれば、正規表現作成工程で作成された正規表現を語表現に変換する語表現変換工程を備えているので、正規表現を語表現に変換することが可能となるという効果を奏する。
 また、本実施の形態の形態によれば、語表現は、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語に対して、位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を、追加された穴の数だけ付与することにより形成された記号語であることとしたので、正規表現を、基本となる全ての流れパターンと位相幾何学的に採り得る5種類の操作を規定した語表現に変換することが可能となるという効果を奏する。
 本実施の形態によれば、語表現と同様に、流れ場における構造物設計を行うにあたって、構造物に対して採り得る流れパターンを経験や直感に頼ることなく容易に扱うことができる、流れパターンの正規表現作成方法、正規表現作成装置、および、プログラムを提供することができる。また、語表現と正規表現の両方を使用することで、語表現だけでは特定の流れを限定できないものを限定することができ、流れの最適制御の理論をより発展させることができる。
 なお、語表現と正規表現は、例えば、橋脚の設計、防波堤の配置、港湾の汚染物の除去、風力発電のブレードの設計、列車のパンタグラフの構造、オイルフェンスの最適配置などのように、構造物設計や配置を伴う様々な分野において極めて有用である。また、スポーツ用品の構造設計などのようにスポーツ力学等の分野に応用することも可能である。
 100 正規表現作成装置
 102 制御部
 102a シミュレーション部
 102b 流線解析部
 102c グラフ表現作成部
 102d 正規表現作成部
 102e 語表現変換部
 104 通信制御インターフェース部
 106 記憶部
 106a シミュレーション結果ファイル
 106b 流線図ファイル
 106c 正規表現ファイル
 108 入出力制御インターフェース部
 112 入力装置
 114 出力装置
 120 外部システム
 130 ネットワーク
 200 物体
 300 平板翼

Claims (16)

  1.  位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成方法であって、
     前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、
     前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、
     を含む流れパターンの正規表現作成方法。
  2.  前記グラフ表現は、前記流れパターンで規定される構造安定なハミルトンベクトル場Hに対して、固有のルート付き、ラベル付き、及び向き有りのツリーT=(V,E)を割り当て(但し、Vは頂点と呼ばれる点の集合、Eは、頂点の間を結ぶエッジの集合である)、平面グラフとして可視化したものであることを特徴とする請求項1に記載の流れパターンの正規表現作成方法。
  3.  前記グラフ表現は、親の頂点をv、その子の頂点をw,親の頂点vに割り当てられたラベルをl(v)、子の頂点wに割り当てられたラベルをl(w),vの子頂点集合Γ(v)とした場合、vの子頂点集合Γ(v)を所定の順序関係のルールに従って並び替え、w∈Γ(v)について、l(v)からl(w)への矢印を左から右に並べて描画したものを含むことを特徴とする請求項1又は請求項2に記載の流れパターンの正規表現作成方法。
  4.  前記流れパターンは、
    (1)一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、二つのss-∂-saddle connectionをもつパターンI、
    (2)前記一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンのうちの、吸い込み湧き出し対をもち、一つのsaddle point、それを結ぶhomoclinic saddle connectionと二つのss-saddle connectionをもつパターンII、
    (3)二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンO、
     からなるパターン語の1又は複数で規定されることを特徴とする請求項1~請求項3のいずれか1つに記載の流れパターンの正規表現作成方法。
  5.  前記グラフ作成工程は、
     前記流れパターンが吸い込み湧き出し対を有する場合は、前記流れパターンを、前記吸い込み湧き出し対に最も近い反時計回りのss-orbitを含む領域が一番外側領域になるように変換し、
     前記変換した流れパターンの(ss-)saddle connection diagramに現れる軌道を領域全体から抽出し、
    領域全体から(ss-)saddle connection diagramに現れる軌道をすべて除外して得られる連結成分に頂点を設定し、一番外側にある連結成分をルートとし、
     カレント成分をルートに設定し、
    前記カレント成分と互いに境界を接するような連結成分をカレント成分の子として、境界にあたる軌道に応じてラベルを割り当て,当該ラベルを所定の順序関係に従って並べ、
     前記カレント成分の子をカレント成分に設定して、子がなくなるまで繰り返すこと、
    を含むことを特徴とする請求項1~請求項4のいずれか1つに記載の流れパターンの正規表現作成方法。
  6.  前記流れパターンは、前記パターン語の1又は複数から開始し、流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を付与する操作を、穴の数がN個となるまで繰り返して作成された流れパターン図であることを特徴とする請求項5に記載の流れパターンの正規表現作成方法。
  7.  さらに、前記正規表現作成工程で作成された正規表現を語表現に変換する語表現変換工程と、
     を含むことを特徴とする請求項1~請求項6のいずれか1つに記載の流れパターンの正規表現作成方法。
  8.  前記語表現変換工程では、ルートつき,ラベル付き,及び方向つきのツリーに対応する正規表現のうちの基本タイプの正規表現に前方置換を施すことが可能な許容正規表現を前記語表現に変換することを特徴とする請求項7に記載の流れパターンの正規表現作成方法。
  9.  前記語表現は、一つの穴を有する単連結外部領域において位相幾何学的に採り得る2種類の流れパターンに加えて、二つの穴を有する二重連結外部領域において吸い込み湧き出し対を持たないパターンを追加した、合計3種類の流れパターンを規定するパターン語に対して、前記流れパターンに一つの穴を加える場合に位相幾何学的に採り得る5種類の操作を規定した操作語のうちのいずれか一語を、追加された穴の数だけ付与することにより形成された記号語であることを特徴とする請求項7又は請求項8に記載の流れパターンの正規表現作成方法。
  10.  さらに、流体中の物体に対して設計パラメータの候補を選択する場合に、前記設計パラメータの上限及び下限を設定し、当該設計パラメータの上限と下限で規定されるパラメータ領域から複数のパラメータを選択し、選択した複数のパラメータに対して、それぞれ流れの実験及び/又は数値計算を行い、実験及び/又は数値計算の結果に対して、前記語表現及び/又は前記正規表現を割り当て、割り当てた前記語表現及び/又は前記正規表現のうち、最適状態を示す前記語表現及び/又は前記正規表現を有する設計パラメータを、前記設計パラメータの候補として選択する設計パラメータ選択工程を含むことを特徴とする請求項1~請求項9のいずれか1つに記載の流れパターンの正規表現作成方法。
  11.  前記設計パラメータ選択工程は、さらに、前記選択した設計パラメータを使用して最適化設計を行う工程を含むことを特徴とする請求項10に記載の流れパターンの正規表現作成方法。
  12.  さらに、流体中の物体に対して最適な設計パラメータを選択する場合に、現在の流れパターンの前記正規表現から目的とする流れパターンの前記正規表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの前記正規表現となる設計パラメータを選択する第2の設計パラメータ選択工程を含むことを特徴とする請求項1~請求項9のいずれか1つに記載の流れパターンの正規表現作成方法。
  13.  前記第2の設計パラメータ選択工程では、前記現在の流れパターンの前記正規表現及び語表現から目的とする流れパターンの前記正規表現及び語表現へ遷移するように設計パラメータを変更し、当該目的とする流れパターンの前記正規表現及び語表現となる設計パラメータを選択することを特徴とする請求項12に記載の流れパターンの正規表現作成方法。
  14.  前記第2の設計パラメータ選択工程は、
    前記現在の流れパターンの前記正規表現及び語表現から前記目的とする流れパターンの正規表現及び語表現へ遷移するための1又は複数の経路を取得し、
     前記取得した1又は複数の経路のうち、設計パラメータを変更して前記現在の流れパターンから遷移可能な経路を選択し、
     前記選択した経路に対して、前記変更した設計パラメータを中心にして設計パラメータを変更し、前記目的とする流れパターンの前記正規表現及び語表現となる設計パラメータを選択することを特徴とする請求項13に記載の流れパターンの正規表現作成方法。
  15.  位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成装置であって、
     前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成手段と、
     前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成手段と、
     を備えたことを特徴とする流れパターンの正規表現作成装置。
  16.  位相幾何学的にN(但し、Nは1以上の整数)個の穴を有する多重連結外部領域における流れパターンの正規表現を作成する正規表現作成装置に搭載されるプログラムであって、
     前記流れパターンに1対1に対応するグラフ表現を作成するグラフ表現作成工程と、
     前記グラフ表現作成工程で作成されたグラフ表現から正規表現を作成する正規表現作成工程と、
     をコンピュータに実行させることを特徴とするコンピュータが実行可能なプログラム。
PCT/JP2015/081402 2014-11-06 2015-11-06 流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム WO2016072515A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2016557834A JP6401288B2 (ja) 2014-11-06 2015-11-06 流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム
US15/524,897 US10540460B2 (en) 2014-11-06 2015-11-06 Flow pattern regular expression creating method, regular expression creating apparatus, and computer-executable program
EP15858028.2A EP3217303A4 (en) 2014-11-06 2015-11-06 Regular expression creation method and regular expression creation device of flow pattern, and computer-executable program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014226532 2014-11-06
JP2014-226532 2014-11-06

Publications (1)

Publication Number Publication Date
WO2016072515A1 true WO2016072515A1 (ja) 2016-05-12

Family

ID=55909238

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/081402 WO2016072515A1 (ja) 2014-11-06 2015-11-06 流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム

Country Status (4)

Country Link
US (1) US10540460B2 (ja)
EP (1) EP3217303A4 (ja)
JP (1) JP6401288B2 (ja)
WO (1) WO2016072515A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021019883A1 (ja) 2019-07-30 2021-02-04 国立研究開発法人科学技術振興機構 有限型の流れパターンの語表現装置、語表現方法、プログラム、構造物形状の学習方法および構造物設計方法
WO2021153721A1 (ja) 2020-01-30 2021-08-05 国立大学法人京都大学 流れパターンの語表現装置、語表現方法およびプログラム

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108629065B (zh) * 2017-11-09 2021-09-17 中车株洲电力机车有限公司 一种小转角受电弓弓头设计方法
CN112182922B (zh) * 2020-09-07 2022-10-21 三峡大学 方形桥墩冲刷问题的绕流流场计算方法
CN112308051B (zh) * 2020-12-29 2021-10-29 北京易真学思教育科技有限公司 文本框检测方法、装置、电子设备和计算机存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008532176A (ja) * 2005-02-28 2008-08-14 ザイ デクマ アクチボラゲット 認識グラフ
WO2014041917A1 (ja) * 2012-09-14 2014-03-20 独立行政法人科学技術振興機構 流れパターンの語表現方法、語表現装置、および、プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008532176A (ja) * 2005-02-28 2008-08-14 ザイ デクマ アクチボラゲット 認識グラフ
WO2014041917A1 (ja) * 2012-09-14 2014-03-20 独立行政法人科学技術振興機構 流れパターンの語表現方法、語表現装置、および、プログラム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIN ARAI: "Researches on fluid mechanics from perspective of mathematical physics Towards topological understanding of the structure of fluid flows", NAGARE, vol. 33, no. 1, 25 February 2014 (2014-02-25), pages 23 - 28, XP055440865, ISSN: 0286-3154 *
See also references of EP3217303A4 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021019883A1 (ja) 2019-07-30 2021-02-04 国立研究開発法人科学技術振興機構 有限型の流れパターンの語表現装置、語表現方法、プログラム、構造物形状の学習方法および構造物設計方法
WO2021153721A1 (ja) 2020-01-30 2021-08-05 国立大学法人京都大学 流れパターンの語表現装置、語表現方法およびプログラム
JP7418719B2 (ja) 2020-01-30 2024-01-22 株式会社Cardio Flow Design 流れパターンの語表現装置、語表現方法およびプログラム

Also Published As

Publication number Publication date
US10540460B2 (en) 2020-01-21
US20170323033A1 (en) 2017-11-09
EP3217303A4 (en) 2018-06-20
JPWO2016072515A1 (ja) 2017-08-17
JP6401288B2 (ja) 2018-10-10
EP3217303A1 (en) 2017-09-13

Similar Documents

Publication Publication Date Title
JP6401288B2 (ja) 流れパターンの正規表現作成方法、正規表現作成装置、および、コンピュータが実行可能なプログラム
Mackaness et al. Generalisation of geographic information: cartographic modelling and applications
Št'ava et al. Inverse procedural modeling by automatic generation of L‐systems
Yang et al. Sketch-based modeling of parameterized objects.
Scheibel et al. A Taxonomy of Treemap Visualization Techniques.
CN107248142B (zh) 一种文物碎片自动拼接方法
CN113487024A (zh) 交替序列生成模型训练方法、从文本中抽取图的方法
Schneider et al. Interactive comparison of multifield scalar data based on largest contours
Van Ham et al. Visualization of state transition graphs
De Carli et al. Procedural generation of 3D canyons
JP6440629B2 (ja) 流体遷移経路取得装置、流体遷移経路取得方法、および、プログラム
KR101635498B1 (ko) 흐름 패턴의 워드 표현 방법, 워드 표현 장치 및 프로그램
Overholser Descendent tropical mirror symmetry for $\mathbb {P}^ 2$
Yousif et al. Shape clustering using k-medoids in architectural form finding
WO2021019883A1 (ja) 有限型の流れパターンの語表現装置、語表現方法、プログラム、構造物形状の学習方法および構造物設計方法
Naitzat et al. M-Boost: Profiling and refining deep neural networks with topological data analysis
Akutsu et al. A simple linear-time algorithm for computing the centroid and canonical form of a plane graph and its applications
Chen et al. Automatic parameterisation of semantic 3D city models for urban design optimisation
Miranda Carranza Convex maps, some basic concepts and a new method to generate them
Burch et al. The aesthetics of rapidly-exploring random trees
Borsa et al. The wreath process: A totally generative model of geometric shape based on nested symmetries
Thapa High Clearance Collision-Free Paths
Liao et al. Hierarchical Automatic Power Plane Generation with Genetic Optimization and Multilayer Perceptron
Su et al. From a link-node-based network representation model to a lane-based network representation model: two-dimensional arrangements approach
Raj Depth-based visualizations for ensemble data and graphs

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15858028

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016557834

Country of ref document: JP

Kind code of ref document: A

REEP Request for entry into the european phase

Ref document number: 2015858028

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 15524897

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE