JP6282958B2 - Evaluation method, evaluation apparatus, and program - Google Patents

Evaluation method, evaluation apparatus, and program Download PDF

Info

Publication number
JP6282958B2
JP6282958B2 JP2014163084A JP2014163084A JP6282958B2 JP 6282958 B2 JP6282958 B2 JP 6282958B2 JP 2014163084 A JP2014163084 A JP 2014163084A JP 2014163084 A JP2014163084 A JP 2014163084A JP 6282958 B2 JP6282958 B2 JP 6282958B2
Authority
JP
Japan
Prior art keywords
particle
data set
movement
particles
time
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
JP2014163084A
Other languages
Japanese (ja)
Other versions
JP2015072679A (en
Inventor
加納 純也
純也 加納
真吾 石原
真吾 石原
通典 山本
通典 山本
尚吾 中野
尚吾 中野
松岡 佑樹
佑樹 松岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tohoku University NUC
Sumitomo Bakelite Co Ltd
Original Assignee
Tohoku University NUC
Sumitomo Bakelite Co Ltd
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 Tohoku University NUC, Sumitomo Bakelite Co Ltd filed Critical Tohoku University NUC
Priority to JP2014163084A priority Critical patent/JP6282958B2/en
Publication of JP2015072679A publication Critical patent/JP2015072679A/en
Application granted granted Critical
Publication of JP6282958B2 publication Critical patent/JP6282958B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01FMIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F23/00Mixing according to the phases to be mixed, e.g. dispersing or emulsifying
    • B01F23/60Mixing solids with solids
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01FMIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F27/00Mixers with rotary stirring devices in fixed receptacles; Kneaders
    • B01F27/80Mixers with rotary stirring devices in fixed receptacles; Kneaders with stirrers rotating about a substantially vertical axis
    • B01F27/90Mixers with rotary stirring devices in fixed receptacles; Kneaders with stirrers rotating about a substantially vertical axis with paddles or arms 
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01FMIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F2215/00Auxiliary or complementary information in relation with mixing
    • B01F2215/04Technical information in relation with mixing
    • B01F2215/0404Technical information in relation with mixing theories or general explanations of phenomena associated with mixing or generalizations of a concept by comparison of equivalent methods

Landscapes

  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本発明は、評価方法、評価装置およびプログラムに関する。   The present invention relates to an evaluation method, an evaluation apparatus, and a program.

複数の粒子を攪拌するために攪拌容器が用いられることがある。この場合、攪拌容器の内部での粒子の攪拌状態が評価されることがある。特許文献1には、離散要素法(DEM:Discrete Element Method)に基づいて粒子の攪拌状態がシミュレーションされることが記載されている。特許文献1では、シミュレーションの結果に基づいて、任意の位置での粒子の混合度が算出されている。この混合度は、サンプリングされた粒子の数と、当該サンプリングされた粒子の中の特定の種類の粒子の数と、の比によって定められる。非特許文献1にも、DEMに基づいて粒子の攪拌状態がシミュレーションされることが記載されている。   A stirring vessel may be used to stir a plurality of particles. In this case, the stirring state of the particles inside the stirring container may be evaluated. Patent Document 1 describes that the agitation state of particles is simulated based on the discrete element method (DEM: Discrete Element Method). In Patent Literature 1, the degree of mixing of particles at an arbitrary position is calculated based on the result of simulation. This degree of mixing is determined by the ratio of the number of sampled particles to the number of specific types of particles in the sampled particles. Non-Patent Document 1 also describes that the stirring state of particles is simulated based on DEM.

特開2000−214134号公報JP 2000-214134 A

Brenda Remy, Johannes G. Khinast, Benjamin J. Glasser, Chem. Eng. Sci., 66 (2011) 1811-1824.Brenda Remy, Johannes G. Khinast, Benjamin J. Glasser, Chem. Eng. Sci., 66 (2011) 1811-1824.

本発明者らは、攪拌容器で攪拌される粒子の混合状態を判断するための指標を検討した。   The present inventors examined an index for judging the mixing state of particles stirred in a stirring vessel.

本発明によれば、
コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示部に表示させる工程と、
実行する評価方法が提供される。
本発明によれば、
コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示部に表示させる工程と、
を実行する評価方法が提供される。
本発明によれば、
コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを統計的に処理された形式で表示部に表示させる工程と、
を実行し、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離である評価方法が提供される。
According to the present invention,
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
The data set, wherein each particle, and the moving amount corresponding to the particles, a step of causing display on the display unit in the associated format,
An evaluation method for performing is provided.
According to the present invention,
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
Displaying the data set on a display unit in a format in which the movement amount of each particle and the frequency corresponding to the movement amount are associated with each other;
An evaluation method for performing is provided.
According to the present invention,
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
Displaying the data set on a display in a statistically processed format;
Run
The amount of movement of each particle is the distance that the particle has moved in a projected locus obtained by orthogonally projecting the locus of the particle from the first time to the second time on a predetermined plane or a predetermined straight line. A method is provided.

本発明によれば、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示する表示部と、
を含む評価装置が提供される。
本発明によれば、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示する表示部と、
を含む評価装置が提供される。
本発明によれば、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを統計的に処理された形式で表示する表示部と、
を含み、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離である評価装置が提供される。
According to the present invention,
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display unit for displaying the data set in a format in which each particle and the movement amount corresponding to the particle are associated with each other;
An evaluation device is provided.
According to the present invention,
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display unit for displaying the data set in a format in which the amount of movement of each particle and the frequency corresponding to the amount of movement are associated with each other;
An evaluation device is provided.
According to the present invention,
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display for displaying the data set in a statistically processed format;
Including
The amount of movement of each particle is the distance that the particle has moved in a projected locus obtained by orthogonally projecting the locus of the particle from the first time to the second time on a predetermined plane or a predetermined straight line. An apparatus is provided.

本発明によれば、
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示する機能と、
を持たせるプログラムが提供される。
本発明によれば、
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示する機能と、
を持たせるプログラムが提供される。
本発明によれば、
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを統計的に処理された形式で表示する機能と、
を持たせ、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離であるプログラムが提供される。
According to the present invention,
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a format in which each particle and the movement amount corresponding to the particle are associated with each other;
A program for providing
According to the present invention,
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a format in which the movement amount of each particle and the frequency corresponding to the movement amount are associated with each other;
A program for providing
According to the present invention,
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a statistically processed format;
Hold
The moving amount of each particle is a distance that the particle has moved in a projected locus obtained by orthogonally projecting a locus traced from the first time to the second time on a predetermined plane or a straight line. Is provided.

本発明によれば、攪拌容器内における粒子の混合状態を判断するための指標が提供される。   According to the present invention, an index for determining the mixing state of particles in a stirring vessel is provided.

第1の実施形態に係る評価方法を示すフローチャートである。It is a flowchart which shows the evaluation method which concerns on 1st Embodiment. 図1の評価方法に用いられる評価装置の構成を示すブロック図である。It is a block diagram which shows the structure of the evaluation apparatus used for the evaluation method of FIG. 二粒子の接触力モデルの概念図である。It is a conceptual diagram of the contact force model of two particles. 第1の実施形態に係るステップS40で表示されるデータセットの第1例を示す図である。It is a figure which shows the 1st example of the data set displayed by step S40 which concerns on 1st Embodiment. 第1の実施形態に係るステップS40で表示されるデータセットの第2例を示す図である。It is a figure which shows the 2nd example of the data set displayed by step S40 which concerns on 1st Embodiment. 第1の実施形態に係るステップS40で表示されるデータセットの第3例を示す図である。It is a figure which shows the 3rd example of the data set displayed by step S40 which concerns on 1st Embodiment. 第1の実施形態に係るステップS40で表示されるデータセットの第4例を示す図である。It is a figure showing the 4th example of a data set displayed at Step S40 concerning a 1st embodiment. 射影軌跡を説明するための図である。It is a figure for demonstrating a projection locus | trajectory. (a)は、攪拌容器の一例を示す平面図であり、(b)は、(a)のA−A´における断面図である。(A) is a top view which shows an example of a stirring container, (b) is sectional drawing in AA 'of (a). 第3の実施形態に係る評価装置の構成を示すブロック図である。It is a block diagram which shows the structure of the evaluation apparatus which concerns on 3rd Embodiment. 図10に示した評価部の構成を示すブロック図である。It is a block diagram which shows the structure of the evaluation part shown in FIG. 粒子の移動量の分布の類似度及び粒子の移動量の分布の散布度の評価方法の一例を説明するための図である。It is a figure for demonstrating an example of the evaluation method of the similarity of the distribution of particle | grain movement amount, and the dispersion degree of distribution of particle | grain movement amount.

以下、本発明の実施の形態について、図面を用いて説明する。尚、すべての図面において、同様な構成要素には同様の符号を付し、適宜説明を省略する。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. In all the drawings, the same reference numerals are given to the same components, and the description will be omitted as appropriate.

なお、以下に示す説明において、各装置の各構成要素は、ハードウエア単位の構成ではなく、機能単位のブロックを示している。各装置の各構成要素は、任意のコンピュータのCPU、メモリ、メモリにロードされた本図の構成要素を実現するプログラム、そのプログラムを格納するハードディスクなどの記憶メディア、ネットワーク接続用インタフェースを中心にハードウエアとソフトウエアの任意の組合せによって実現される。そして、その実現方法、装置には様々な変形例がある。   In the following description, each component of each device is not a hardware unit configuration but a functional unit block. Each component of each device includes a CPU, memory, a program that realizes the components shown in the figure loaded in the memory, a storage medium such as a hard disk for storing the program, and a network connection interface. It is realized by any combination of software and software. There are various modifications of the implementation method and apparatus.

(第1の実施形態)
図1は、第1の実施形態に係る評価方法を示すフローチャートである。この評価方法は、以下の工程を含んでいる。まず、複数の種類の粒子が攪拌容器の内部で攪拌される状態がシミュレートされる(ステップS10)。次に、シミュレートされた各粒子の第1時刻から第2時刻までの間の移動量が算出される(ステップS20)。次に、所定のデータセットが生成される(ステップS30)。データセットは、各粒子と、当該粒子に対応する移動量と、を対応付けて含んでいる。次に、データセットが統計的に処理された形式で表示される(ステップS40)。
(First embodiment)
FIG. 1 is a flowchart showing an evaluation method according to the first embodiment. This evaluation method includes the following steps. First, a state in which a plurality of types of particles are stirred inside the stirring vessel is simulated (step S10). Next, the amount of movement of each simulated particle from the first time to the second time is calculated (step S20). Next, a predetermined data set is generated (step S30). The data set includes each particle and the movement amount corresponding to the particle in association with each other. Next, the data set is displayed in a statistically processed format (step S40).

ステップS40において統計的に処理された形式で表示されたデータセットは、混合状態情報として取り扱うことができる。すなわち、混合状態情報とは、複数の種類の粒子の混合状態を示す情報をいい、たとえば以下のような情報をいう。
(i)複数の種類の粒子がどの程度均一に混合されているか(複数の種類の粒子の混合の均一度)を示す情報
(ii)混合後、系内で複数の粒子が不均一に分布する場合において、各粒子がどのような状態で分布しているか(各粒子の分布状況)を示す情報
(iii)複数の種類の粒子が種類ごとにどのように分布しているか(各種類の粒子の分布状況)を示す情報
The data set displayed in the statistically processed format in step S40 can be handled as mixed state information. That is, the mixed state information refers to information indicating a mixed state of a plurality of types of particles, for example, the following information.
(I) Information indicating how uniformly a plurality of types of particles are mixed (uniformity of mixing of a plurality of types of particles) (ii) After mixing, a plurality of particles are unevenly distributed in the system In the case, information indicating how each particle is distributed (distribution state of each particle) (iii) how a plurality of types of particles are distributed for each type (for each type of particles) Information)

本実施形態では、以上の混合状態情報に基づいて、複数の種類の粒子の混合状態を評価することができる。このようにして、本実施形態では、攪拌容器内における粒子の混合状態を判断するための指標が提供される。   In the present embodiment, the mixed state of a plurality of types of particles can be evaluated based on the above mixed state information. Thus, in this embodiment, an index for determining the mixing state of particles in the stirring vessel is provided.

次に、本実施形態の評価装置100について、図2を用いて説明する。評価装置100は、図1に示される評価方法を実行する。図2は、評価装置100の構成を示すブロック図である。評価装置100は、入力部102と、シミュレーション部104と、移動量算出部106と、データセット生成部108と、表示部110と、記憶部112と、を含んでいる。入力部102は、評価装置100のユーザの指示を示す信号を、シミュレーション部104、移動量算出部106、データセット生成部108および表示部110に送出する。シミュレーション部104は、粒子の攪拌状態をシミュレートする(ステップS10)。移動量算出部106は、移動量を算出する(ステップS20)。データセット生成部108は、データセットを生成する(ステップS30)。表示部110は、データセットを統計的に処理された形式で表示する(ステップS40)。記憶部112は、シミュレーション部104、移動量算出部106、データセット生成部108および表示部110からの情報を記憶する。評価装置100の詳細は以下のとおりである。   Next, the evaluation apparatus 100 of this embodiment is demonstrated using FIG. The evaluation apparatus 100 executes the evaluation method shown in FIG. FIG. 2 is a block diagram showing the configuration of the evaluation apparatus 100. The evaluation apparatus 100 includes an input unit 102, a simulation unit 104, a movement amount calculation unit 106, a data set generation unit 108, a display unit 110, and a storage unit 112. The input unit 102 sends a signal indicating an instruction of the user of the evaluation apparatus 100 to the simulation unit 104, the movement amount calculation unit 106, the data set generation unit 108, and the display unit 110. The simulation unit 104 simulates the stirring state of the particles (Step S10). The movement amount calculation unit 106 calculates the movement amount (step S20). The data set generation unit 108 generates a data set (step S30). The display unit 110 displays the data set in a statistically processed format (step S40). The storage unit 112 stores information from the simulation unit 104, the movement amount calculation unit 106, the data set generation unit 108, and the display unit 110. Details of the evaluation apparatus 100 are as follows.

シミュレーション部104について説明する。シミュレーション部104は、複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする。本実施形態では、シミュレーション部104は、離散要素法(DEM:Discrete Element Method)に基づいて、複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする。他の例においては、シミュレーション部104は、DEMとは異なる方法(例えば、埋め込み境界法)を用いて、複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートしてもよい。   The simulation unit 104 will be described. The simulation unit 104 simulates a state in which a plurality of types of particles are stirred inside the stirring vessel. In the present embodiment, the simulation unit 104 simulates a state in which a plurality of types of particles are agitated inside the agitation vessel based on a discrete element method (DEM: Discrete Element Method). In another example, the simulation unit 104 may simulate a state in which a plurality of types of particles are agitated inside the agitation vessel using a method different from the DEM (for example, an embedded boundary method).

シミュレーション部104は、DEMにおいて、攪拌容器および粒子をモデル化する。モデル化される攪拌容器の構造および形状は、特定の構造および形状に限られない。モデル化される粒子の形状も、特定の形状に限られない。なお、本実施形態においてモデル化される粒子は、固体粒子である。   The simulation unit 104 models a stirring container and particles in the DEM. The structure and shape of the stirring vessel to be modeled is not limited to a specific structure and shape. The shape of the modeled particle is not limited to a specific shape. In addition, the particle modeled in this embodiment is a solid particle.

モデル化された粒子について、図3を用いて説明する。図3は、二粒子iおよびjの接触力モデルの概念図である。図3(a)は、粒子iと粒子jとの接触点の法線方向の接触力モデルを示している。図3(b)は、粒子iと粒子jとの接触点の接線方向の接触力モデルを示している。本実施形態では、粒子は、球状の形状を有する。DEMにおいて、互いに接触する粒子iと粒子jとの接触点の法線方向にはたらく力は、図3(a)に示されるように、バネ定数kのバネおよび減衰係数ηのダッシュポットによりモデル化される。一方、互いに接触する粒子iと粒子jとの接触点の接線方向にはたらく力は、図3(b)に示されるように、バネ定数kのバネ、減衰係数ηのダッシュポットおよび動摩擦係数μの摩擦スライダによりモデル化される。ユーザは、入力部102により、バネ定数kおよびk、減衰係数ηおよびηならびに動摩擦係数μとして適当な値を、シミュレーション部104に入力してもよい。この場合、シミュレーション部104は、入力されたパラメータに基づいて、粒子をモデル化する。 The modeled particles will be described with reference to FIG. FIG. 3 is a conceptual diagram of a contact force model of two particles i and j. FIG. 3A shows a contact force model in the normal direction of the contact point between the particle i and the particle j. FIG. 3B shows a contact force model in the tangential direction of the contact point between the particle i and the particle j. In the present embodiment, the particles have a spherical shape. In DEM, the force acting in the normal direction of the contact points between the particles i and particle j in contact with each other, as shown in FIG. 3 (a), the dashpot spring and damping coefficient eta n of the spring constant k n Modeled. On the other hand, as shown in FIG. 3B, the force acting in the tangential direction of the contact point between the particle i and the particle j that are in contact with each other includes a spring having a spring constant k t , a dashpot having a damping coefficient η t , and a dynamic friction coefficient. It is modeled by friction slider mu f. The user may input appropriate values as the spring constants k n and k t , the damping coefficients η n and η t and the dynamic friction coefficient μ f to the simulation unit 104 through the input unit 102. In this case, the simulation unit 104 models particles based on the input parameters.

DEMでは、各粒子の並進運動の方程式および回転運動の方程式が、図3(a)および図3(b)に示されるモデルに基づいて立てられる。粒子iにおける並進運動の方程式および回転運動の方程式は、それぞれ、以下の式(1)および(2)に示されるようになる。シミュレーション部104は、図3(a)および図3(b)に示されるモデルに基づいて、式(1)および(2)に示される運動方程式を立てる。

Figure 0006282958
ただし、m、I、v、ω、Rは、それぞれ、粒子iの質量、慣性モーメント、速度、角速度、半径である。Fij およびFij は、それぞれ、粒子iと粒子jとの接触によって生じる、法線方向および接線方向の力である。gは、重力加速度である。τrijは、粒子iと粒子jとの接触によって生じるトルクである。式(1)では、シミュレーション部104は、重力加速度を考慮して並進運動の方程式を立てている。重力加速度が考慮されない場合、mgの項が式(1)の右辺から除去されてもよい。式(1)および(2)と同様の式が、粒子i以外の各粒子にも立てられる。本実施形態では、シミュレーション部104がこれらの運動方程式を解く。シミュレーション部104は、差分法を用いた数値解析により運動方程式を解いてもよい。このようにしてシミュレーション部104は、粒子の攪拌状態をシミュレートする(ステップS10)。 In DEM, equations of translational motion and rotational motion of each particle are established based on the models shown in FIGS. 3 (a) and 3 (b). The equations for translational motion and rotational motion in the particle i are as shown in the following formulas (1) and (2), respectively. The simulation unit 104 establishes the equations of motion shown in the equations (1) and (2) based on the models shown in FIGS. 3 (a) and 3 (b).
Figure 0006282958
However, m i , I i , v i , ω i , and R i are the mass, moment of inertia, velocity, angular velocity, and radius of the particle i, respectively. F ij N and F ij T are normal and tangential forces generated by the contact between the particle i and the particle j, respectively. g is a gravitational acceleration. τ rij is a torque generated by the contact between the particle i and the particle j. In Expression (1), the simulation unit 104 sets up an equation for translational motion in consideration of gravitational acceleration. If the gravitational acceleration is not taken into account, term of m i g may be removed from the right hand side of equation (1). Expressions similar to Expressions (1) and (2) are established for each particle other than the particle i. In the present embodiment, the simulation unit 104 solves these equations of motion. The simulation unit 104 may solve the equation of motion by numerical analysis using a difference method. Thus, the simulation part 104 simulates the stirring state of particle | grains (step S10).

シミュレーション部104は、シミュレーションの結果を記憶部112に出力する。記憶部112は、シミュレーションの結果を記憶する。当該シミュレーションの結果は、所定時刻における各粒子の座標の情報を示している。シミュレーション部104は、シミュレーションが終了した後、シミュレーションが終了した旨を示すシミュレーション終了信号を移動量算出部106に送出する。   The simulation unit 104 outputs the simulation result to the storage unit 112. The storage unit 112 stores simulation results. The result of the simulation shows information on the coordinates of each particle at a predetermined time. After the simulation is finished, the simulation unit 104 sends a simulation end signal indicating that the simulation is finished to the movement amount calculation unit 106.

次に、移動量算出部106について説明する。移動量算出部106は、シミュレーション終了信号をシミュレーション部104から受け取る。次に、移動量算出部106は、記憶部112にアクセスする。そして移動量算出部106は、記憶部112に記憶されているシミュレーションの結果を読み取る。次に、移動量算出部106は、シミュレートされた各粒子の第1時刻から第2時刻までの間の移動量を算出する(ステップS20)。本実施形態では、各粒子の移動量は、当該粒子が第1時刻から第2時刻までの間に攪拌容器の内部で移動した距離である。第1時刻および第2時刻は、任意の時刻である。第1時刻は、例えば、攪拌が開始する時刻であってもよい。第2時刻は、例えば、攪拌が終了する時刻であってもよい。   Next, the movement amount calculation unit 106 will be described. The movement amount calculation unit 106 receives a simulation end signal from the simulation unit 104. Next, the movement amount calculation unit 106 accesses the storage unit 112. Then, the movement amount calculation unit 106 reads the simulation result stored in the storage unit 112. Next, the movement amount calculation unit 106 calculates the movement amount of each simulated particle from the first time to the second time (step S20). In the present embodiment, the movement amount of each particle is the distance that the particle has moved inside the stirring vessel from the first time to the second time. The first time and the second time are arbitrary times. The first time may be, for example, the time when stirring starts. For example, the second time may be a time at which stirring is completed.

本実施形態では、移動量算出部106は、各粒子が第1時刻から第2時刻までの間に攪拌容器の内部で移動した距離を算出する。この場合、移動量算出部106は、式(1)の解に基づいて、当該距離を算出してもよい。式(1)によって所定時刻における各粒子の座標が与えられる。移動量算出部106は、所定時刻における各粒子の座標をトレースすることで、各粒子が第1時刻から第2時刻までの間に攪拌容器の内部で移動した距離を算出することができる。   In the present embodiment, the movement amount calculation unit 106 calculates the distance that each particle has moved within the stirring vessel from the first time to the second time. In this case, the movement amount calculation unit 106 may calculate the distance based on the solution of Expression (1). Equation (1) gives the coordinates of each particle at a given time. The movement amount calculation unit 106 can calculate the distance that each particle has moved within the stirring vessel from the first time to the second time by tracing the coordinates of each particle at a predetermined time.

移動量算出部106は、各粒子の移動距離を記憶部112に出力する。記憶部112は、各粒子の移動距離を記憶する。移動量算出部106は、移動量の算出が終了した後、移動量の算出が終了した旨を示す移動量算出終了信号をデータセット生成部108に送出する。   The movement amount calculation unit 106 outputs the movement distance of each particle to the storage unit 112. The storage unit 112 stores the movement distance of each particle. After completing the calculation of the movement amount, the movement amount calculation unit 106 sends a movement amount calculation end signal indicating that the calculation of the movement amount has ended to the data set generation unit 108.

次に、データセット生成部108について説明する。データセット生成部108は、移動量算出終了信号を移動量算出部106から受け取る。次に、データセット生成部108は、記憶部112にアクセスする。そしてデータセット生成部108は、記憶部112に記憶されている各粒子の移動量を読み取る。次に、データセット生成部108は、データセットを生成する(ステップS30)。データセットは、各粒子と、当該粒子に対応する移動量と、を対応付けて含んでいる。   Next, the data set generation unit 108 will be described. The data set generation unit 108 receives a movement amount calculation end signal from the movement amount calculation unit 106. Next, the data set generation unit 108 accesses the storage unit 112. Then, the data set generation unit 108 reads the movement amount of each particle stored in the storage unit 112. Next, the data set generation unit 108 generates a data set (step S30). The data set includes each particle and the movement amount corresponding to the particle in association with each other.

データセット生成部108は、データセットを記憶部112に出力する。記憶部112は、データセットを記憶する。データセット生成部108は、データセットの生成が終了した後、データセットの生成が終了した旨を示すデータセット終了信号を表示部110に送出する。   The data set generation unit 108 outputs the data set to the storage unit 112. The storage unit 112 stores a data set. After the generation of the data set is completed, the data set generation unit 108 sends a data set end signal indicating that the generation of the data set has ended to the display unit 110.

次に、表示部110について説明する。表示部110は、データセット生成終了信号をデータセット生成部108から受け取る。次に、表示部110は、記憶部112にアクセスする。そして表示部110は、記憶部112に記憶されているデータセットを読み取る。次に、表示部110は、データセットを統計的に処理する。次に、表示部110は、データセットを統計的に処理された形式で、評価装置100のユーザに対して表示する(ステップS40)。本実施形態では、表示部110は、ディスプレイに当該データセットを表示する。統計的に処理された形式で表示されたデータセットは、上述した混合状態情報として取り扱うことができる。   Next, the display unit 110 will be described. The display unit 110 receives a data set generation end signal from the data set generation unit 108. Next, the display unit 110 accesses the storage unit 112. The display unit 110 reads the data set stored in the storage unit 112. Next, the display unit 110 statistically processes the data set. Next, the display unit 110 displays the data set in a statistically processed format to the user of the evaluation device 100 (step S40). In the present embodiment, the display unit 110 displays the data set on the display. A data set displayed in a statistically processed format can be handled as the above-described mixed state information.

次に、ステップS40の詳細について、図4から図7を用いて説明する。   Next, details of step S40 will be described with reference to FIGS.

図4は、ステップS40で表示されるデータセットの第1例を示す図である。ステップS40では、図4に示されるように、データセットが統計的に処理されてテーブル形式で表示される。図4のテーブルでは、各粒子と、当該粒子に対応する移動量と、が対応付けられている。さらに図4では、各粒子と、当該粒子に対応する移動量と、を示すセル(例えば、図4の「粒子A2」および「DA2」のセル)が、粒子の種類(粒子群)ごとに分類されている。表示部110は、図4に示されるテーブルを評価装置100のユーザに対して表示する。評価装置100のユーザは、図4のテーブルによって、各粒子の移動距離を俯瞰的に把握することができる。このようにして評価装置100のユーザは、図4のテーブルに基づいて、粒子の混合状態を評価する。 FIG. 4 is a diagram illustrating a first example of the data set displayed in step S40. In step S40, as shown in FIG. 4, the data set is statistically processed and displayed in a table format. In the table of FIG. 4, each particle is associated with the amount of movement corresponding to the particle. Further, in FIG. 4, cells (for example, “particle A2” and “D A2 ” cells in FIG. 4) indicating each particle and the movement amount corresponding to the particle are provided for each type of particle (particle group). It is classified. The display unit 110 displays the table shown in FIG. 4 to the user of the evaluation apparatus 100. The user of the evaluation apparatus 100 can grasp the movement distance of each particle from a bird's-eye view using the table of FIG. In this way, the user of the evaluation apparatus 100 evaluates the particle mixing state based on the table of FIG.

図5は、ステップS40で表示されるデータセットの第2例を示す図である。ステップS40では、図5に示されるように、データセットが統計的に処理されてテーブル形式で表示される。図5のテーブルでは、各粒子が移動量に基づいて分類されている度数分布表が表示されている。図5の度数分布表は、例えば、移動距離がd1[cm]からd2[cm]までの範囲に入る粒子の頻度がfA1[%]であること(図5の粒子群Aの一行目)を意味している。図5では、粒子の種類(粒子群)ごとに度数分布表が生成されている。他の例では、複数の種類の粒子群ごとに度数分布表が生成されてもよいし、すべての種類の粒子を含む度数分布表が生成されてもよい。表示部110は、図5に示されるテーブルを評価装置100のユーザに対して表示する。評価装置100のユーザは、図5のテーブルによって、各粒子の移動距離を俯瞰的に把握することができる。このようにして評価装置100のユーザは、図5のテーブルに基づいて、粒子の混合状態を評価する。 FIG. 5 is a diagram illustrating a second example of the data set displayed in step S40. In step S40, as shown in FIG. 5, the data set is statistically processed and displayed in a table format. In the table of FIG. 5, a frequency distribution table in which each particle is classified based on the movement amount is displayed. In the frequency distribution table of FIG. 5, for example, the frequency of particles that fall within the range of the moving distance from d1 [cm] to d2 [cm] is f A1 [%] (first row of particle group A in FIG. 5). Means. In FIG. 5, a frequency distribution table is generated for each type (particle group) of particles. In another example, a frequency distribution table may be generated for each of a plurality of types of particle groups, or a frequency distribution table including all types of particles may be generated. The display unit 110 displays the table shown in FIG. 5 to the user of the evaluation apparatus 100. The user of the evaluation apparatus 100 can grasp the movement distance of each particle from a bird's-eye view using the table of FIG. In this way, the user of the evaluation apparatus 100 evaluates the particle mixing state based on the table of FIG.

図6は、ステップS40で表示されるデータセットの第3例を示す図である。ステップS40では、図6に示されるように、データセットが統計的に処理されてグラフィカル形式で表示される。図6のグラフでは、各粒子と、当該粒子に対応する移動量と、が対応付けられている。さらに図6では、各粒子を示すプロット(例えば、図6の粒子A2を示す白四角プロット)が、粒子の種類(粒子群)ごとに分類されている。表示部110は、図6に示されるグラフを評価装置100のユーザに対して表示する。評価装置100のユーザは、図6のグラフによって、各粒子の移動距離を俯瞰的に把握することができる。このようにして評価装置100のユーザは、図6のグラフに基づいて、粒子の混合状態を評価する。   FIG. 6 is a diagram illustrating a third example of the data set displayed in step S40. In step S40, the data set is statistically processed and displayed in a graphical format, as shown in FIG. In the graph of FIG. 6, each particle is associated with a movement amount corresponding to the particle. Further, in FIG. 6, plots indicating each particle (for example, a white square plot indicating the particle A <b> 2 in FIG. 6) are classified for each particle type (particle group). The display unit 110 displays the graph shown in FIG. 6 to the user of the evaluation device 100. The user of the evaluation apparatus 100 can grasp the movement distance of each particle from a bird's-eye view with the graph of FIG. In this way, the user of the evaluation apparatus 100 evaluates the particle mixing state based on the graph of FIG.

図7は、ステップS40で表示されるデータセットの第4例を示す図である。ステップS40では、図7に示されるように、データセットが統計的に処理されてグラフィカル形式で表示される。図7のグラフでは、各粒子が移動量に基づいて分類されているヒストグラムが表示されている。図7では、粒子の種類(粒子群)ごとにヒストグラムが生成されている。他の例では、複数の種類の粒子群ごとにヒストグラムが生成されてもよいし、すべての種類の粒子を含むヒストグラムが生成されてもよい。表示部110は、図7に示されるグラフを評価装置100のユーザに対して表示する。評価装置100のユーザは、図7のグラフによって、各粒子の移動距離を俯瞰的に把握することができる。このようにして評価装置100のユーザは、図7のグラフに基づいて、粒子の混合状態を評価する。   FIG. 7 is a diagram illustrating a fourth example of the data set displayed in step S40. In step S40, as shown in FIG. 7, the data set is statistically processed and displayed in a graphical format. In the graph of FIG. 7, a histogram in which each particle is classified based on the amount of movement is displayed. In FIG. 7, a histogram is generated for each type of particle (particle group). In another example, a histogram may be generated for each of a plurality of types of particle groups, or a histogram including all types of particles may be generated. The display unit 110 displays the graph shown in FIG. 7 to the user of the evaluation apparatus 100. The user of the evaluation apparatus 100 can grasp the moving distance of each particle from a bird's-eye view with the graph of FIG. In this way, the user of the evaluation apparatus 100 evaluates the particle mixing state based on the graph of FIG.

本実施形態では、以上のように統計的に処理された形式で表示されたデータセットは、混合状態情報として取り扱うことができる。このようにして、本実施形態では、攪拌容器内における粒子の混合状態を判断するための指標が提供される。   In the present embodiment, the data set displayed in the statistically processed format as described above can be handled as mixed state information. Thus, in this embodiment, an index for determining the mixing state of particles in the stirring vessel is provided.

以上、本発明の実施形態について説明したが、上記の実施形態に限らず種々の変形が可能である。すなわち、本実施形態では、流体中(例えば、気体雰囲気および液体中)における粒子の混合を取り扱うことができる。たとえば、本実施形態では、気相(例えば、空気)中の粒子の混合を例に挙げて説明したが、混合の場は気相に限らず液相であってもよい。たとえば、溶媒中の粒子の混合にも適用できる。この場合、溶媒となる液体の物性として、粘性等を考慮してシミュレーションする方法を採用してもよい。   As mentioned above, although embodiment of this invention was described, not only said embodiment but various deformation | transformation are possible. That is, in the present embodiment, mixing of particles in a fluid (for example, in a gas atmosphere and a liquid) can be handled. For example, in the present embodiment, the mixing of particles in a gas phase (for example, air) has been described as an example, but the mixing field is not limited to the gas phase, and may be a liquid phase. For example, it can be applied to mixing particles in a solvent. In this case, a simulation method may be employed in consideration of viscosity and the like as the physical properties of the liquid serving as the solvent.

(第2の実施形態)
第1の実施形態では、粒子の移動量を、3次元空間内の実際の軌跡に沿った移動量としたが、本実施形態では、射影軌跡において当該粒子が移動した距離を移動量とする。具体的には、各粒子が第1時刻から第2時刻までに辿った軌跡を所定の平面または所定の直線に正射影させた軌跡に基づいて移動量を算出する。
(Second Embodiment)
In the first embodiment, the amount of movement of the particles is the amount of movement along the actual trajectory in the three-dimensional space. However, in this embodiment, the distance traveled by the particles in the projected locus is the amount of movement. Specifically, the movement amount is calculated based on a trajectory obtained by orthogonally projecting a trajectory that each particle traces from the first time to the second time on a predetermined plane or a predetermined straight line.

射影軌跡について、図8を用いて説明する。図8は、射影軌跡を説明するための図である。図8では、xyz直交座標系において粒子が第1時刻から第2時刻までに所定の軌跡を辿っている。第1例において、射影軌跡は、xy平面に正射影された軌跡である。第2例においては、射影軌跡は、z軸に正射影された軌跡である。なお、図8の粒子を示す座標は、座標Mにおいてzの最小値をとる。このようにして、粒子は、z軸に正射影された軌跡において、zが減少する方向に移動した後、zが増加する方向に移動している。各粒子の移動量が射影軌跡において当該粒子が移動した距離になる場合、移動量算出部106は、シミュレーションの結果に対して、図8に示される座標系を適用する。そして、移動量算出部106は、各粒子が射影軌跡において移動した距離を算出する。   The projection locus will be described with reference to FIG. FIG. 8 is a diagram for explaining the projected trajectory. In FIG. 8, the particles follow a predetermined locus from the first time to the second time in the xyz rectangular coordinate system. In the first example, the projected trajectory is a trajectory that is orthogonally projected on the xy plane. In the second example, the projected trajectory is a trajectory orthogonally projected on the z axis. In addition, the coordinate which shows the particle | grains of FIG. In this way, the particles move in the direction in which z increases after moving in the direction in which z decreases in the locus orthogonally projected on the z axis. When the movement amount of each particle is the distance traveled by the particle in the projected locus, the movement amount calculation unit 106 applies the coordinate system shown in FIG. 8 to the simulation result. Then, the movement amount calculation unit 106 calculates the distance that each particle has moved in the projected locus.

図8に示されるxyz座標系の方向は、攪拌容器との関係で決められてもよい。具体例について、図9を用いて説明する。図9(a)は、攪拌容器の一例(攪拌容器200)を示す平面図である。図9(b)は、図9(a)のA−A´における断面図である。   The direction of the xyz coordinate system shown in FIG. 8 may be determined in relation to the stirring vessel. A specific example will be described with reference to FIG. FIG. 9A is a plan view showing an example of the stirring container (stirring container 200). FIG. 9B is a cross-sectional view taken along the line AA ′ of FIG.

攪拌容器200は、羽根202を含んでいる。攪拌容器200では、羽根202は、重力加速度gと直交する方向に回転している。シミュレーション部104は、攪拌容器200のモデル化において、x軸方向、y軸方向およびz軸方向を、図9(a)および図9(b)に示されるように規定する。図9(a)および(b)では、z軸方向は、羽根202の回転軸方向に沿っている。一方、x軸およびy軸方向は、羽根202の回転面方向に沿っている。移動量算出部106は、当該座標系において、各粒子が射影軌跡において移動した距離を算出する。   The stirring container 200 includes a blade 202. In the stirring vessel 200, the blade 202 rotates in a direction orthogonal to the gravitational acceleration g. The simulation unit 104 defines the x-axis direction, the y-axis direction, and the z-axis direction as shown in FIGS. 9A and 9B in modeling the stirring vessel 200. 9A and 9B, the z-axis direction is along the rotation axis direction of the blade 202. On the other hand, the x-axis and y-axis directions are along the rotation surface direction of the blade 202. The movement amount calculation unit 106 calculates the distance that each particle has moved in the projected locus in the coordinate system.

本実施形態のデータセット生成部108および表示部110は、第1の実施形態のデータセット生成部108および表示部110と同様に動作する。本実施形態のデータセット生成部108は、第1の実施形態のデータセット生成部108と同様に、データセットを生成する。データセットは、各粒子と、当該粒子に対応する移動量と、を対応付けて含んでいる。本実施形態の表示部110は、第1の実施形態の表示部110と同様に、データセットを統計的に処理することで得られる混合状態情報を、評価装置100のユーザに対して表示する。ただし、本実施形態は、各粒子の移動量が射影軌跡において当該粒子が移動した距離である点について、第1の実施形態と異なる。   The data set generation unit 108 and the display unit 110 of the present embodiment operate in the same manner as the data set generation unit 108 and the display unit 110 of the first embodiment. The data set generation unit 108 of the present embodiment generates a data set in the same manner as the data set generation unit 108 of the first embodiment. The data set includes each particle and the movement amount corresponding to the particle in association with each other. Similar to the display unit 110 of the first embodiment, the display unit 110 of the present embodiment displays the mixed state information obtained by statistically processing the data set to the user of the evaluation apparatus 100. However, the present embodiment is different from the first embodiment in that the amount of movement of each particle is the distance traveled by the particle in the projected locus.

本実施形態では、xyz直交座標系を採用し、射影軌跡を求めたが、これに限らず極座標系を採用し、射影軌跡を求めてもよい   In the present embodiment, the xyz orthogonal coordinate system is adopted and the projection locus is obtained. However, the present invention is not limited thereto, and the polar locus system may be adopted to obtain the projection locus.

(第3の実施形態)
図10は、第3の実施形態に係る評価装置100の構成を示すブロック図であり、第1の実施形態の図2に対応する。本実施形態に係る評価装置100は、表示部110に代わって評価部114を有している点を除いて、第1の実施形態に係る評価装置100と同様の構成である。詳細を後述するように、評価部114は、粒子の混合状態を評価するために用いられる。
(Third embodiment)
FIG. 10 is a block diagram illustrating a configuration of the evaluation apparatus 100 according to the third embodiment, and corresponds to FIG. 2 of the first embodiment. The evaluation apparatus 100 according to the present embodiment has the same configuration as that of the evaluation apparatus 100 according to the first embodiment, except that an evaluation unit 114 is provided instead of the display unit 110. As will be described in detail later, the evaluation unit 114 is used to evaluate the mixed state of particles.

図11は、図10に示した評価部114の構成を示すブロック図である。本図に示すように、評価部114は、類似度評価部116及び散布度評価部118を有している。   FIG. 11 is a block diagram illustrating a configuration of the evaluation unit 114 illustrated in FIG. As shown in the figure, the evaluation unit 114 includes a similarity evaluation unit 116 and a distribution degree evaluation unit 118.

類似度評価部116及び散布度評価部118は、データセット生成部108又は記憶部112から、上記したデータセットを受け取る。この場合、データセットは、一の種類の粒子の移動量の分布を示す情報及び他の種類の粒子の移動量の分布を示す情報を含んでいる(例えば、図7)。   The similarity evaluation unit 116 and the dispersion degree evaluation unit 118 receive the above-described data set from the data set generation unit 108 or the storage unit 112. In this case, the data set includes information indicating the distribution of movement amounts of one type of particles and information indicating the distribution of movement amounts of other types of particles (for example, FIG. 7).

類似度評価部116は、一の種類の粒子の移動量の分布と他の種類の粒子の移動量の分布の類似度を評価する。散布度評価部118は、一の種類の粒子の移動量の分布の散布度及び他の種類の粒子の移動量の分布の散布度を評価する。   The similarity evaluation unit 116 evaluates the similarity between the movement amount distribution of one type of particles and the movement amount distribution of other types of particles. The distribution degree evaluation unit 118 evaluates the distribution degree of the distribution of the movement amount of one type of particles and the distribution degree of the distribution of the movement amount of other types of particles.

評価部114は、上記した類似度に基づいて、一の種類の粒子と他の種類の粒子の混合状態を評価する。さらに、評価部114は、上記した散布度を考慮して、上記した混合状態を評価する。具体的には、評価部114は、一の種類の粒子の上記した散布度及び他の種類の粒子の上記した散布度のそれぞれが基準値以下である場合、上記した混合状態を評価する。   The evaluation unit 114 evaluates the mixed state of one type of particles and other types of particles based on the above-described similarity. Furthermore, the evaluation unit 114 evaluates the above-described mixed state in consideration of the above-described spreading degree. Specifically, the evaluation unit 114 evaluates the above-described mixed state when each of the above-described spreading degree of one type of particles and the above-described spreading degree of other types of particles is equal to or less than a reference value.

本発明者らが検討したところ、一の種類の粒子と他の種類の粒子の混合の均一性は、上記した類似度及び上記した散布度と以下の関係を有する可能性が高いことが明らかとなった。   As a result of examination by the present inventors, it is clear that the uniformity of mixing of one type of particles and other types of particles is highly likely to have the following relationship with the above-mentioned similarity and the above-mentioned degree of dispersion. became.

第1に、一の種類の粒子の上記した散布度及び他の種類の粒子の上記した散布度のいずれもが低い場合、上記した類似度と上記した均一性が強い相関関係を有することが明らかとなった。言い換えると、この場合、上記した類似度が高いほど上記した均一性も高くなる傾向がある。   First, it is clear that the above-mentioned similarity and the above-mentioned uniformity have a strong correlation when both the above-mentioned degree of dispersion of one kind of particles and the above-mentioned degree of dispersion of other kinds of particles are low. It became. In other words, in this case, the higher the similarity is, the higher the uniformity is.

本発明者らが検討したところ、上記した相関関係の理由として、一の種類のいずれの粒子も他の種類のいずれの粒子もほぼ同じ距離を移動していることが可能性として挙がった。上記した散布度が低い場合、一の種類の粒子においてはいずれの粒子もほぼ同じ距離を移動している。同様に、他の種類の粒子においてはいずれの粒子もほぼ同じ距離を移動している。さらに上記した類似度が高い場合、一の種類の上記した分布と他の種類の上記した分布の重なり合いが大きい。言い換えると、一の種類の粒子も他の種類の粒子もほぼ同じ距離を移動している。このような場合、上記した均一性が高くなる可能性が高い。   As a result of the study by the present inventors, as a reason for the above-described correlation, it was mentioned that any one type of particles and any other type of particles may move substantially the same distance. When the above-mentioned degree of dispersion is low, all the particles are moving the same distance in one kind of particles. Similarly, in the other types of particles, all of the particles move approximately the same distance. Further, when the above-described similarity is high, the overlap between one type of the above-described distribution and another type of the above-described distribution is large. In other words, one type of particle and another type of particle are moving approximately the same distance. In such a case, the above-described uniformity is likely to be high.

第2に、一の種類の粒子の上記した散布度及び他の種類の粒子の上記した散布度のいずれか一方又は両方が高い場合、上記した類似度と上記した均一性の相関関係は弱くなることが明らかとなった。言い換えると、この場合、上記した類似度が高くなっても、上記した均一性が高くならない場合が多い。   Second, when one or both of the above-described degree of dispersion of one kind of particles and the above-mentioned degree of dispersion of other kinds of particles are high, the correlation between the above-described similarity and the above-described uniformity is weakened. It became clear. In other words, in this case, the above-described uniformity is often not increased even when the above-described similarity is increased.

本発明者らが検討したところ、上記した相関関係の理由として、粒子の移動距離のばらつきが上記した均一性に影響を与えていることが可能性として挙がった。上記した類似度が高い場合、一の種類の上記した分布の中心と他の種類の上記した分布の中心がほぼ同じ位置になる。この場合においても一の種類の粒子の上記した散布度が高いとき、一の種類の粒子には、上記した中心から移動量のずれた粒子がある程度の割合で含まれている。そしてこのような粒子が上記した均一性を低くしている可能性が高い。   As a result of the study by the present inventors, the reason for the above-described correlation was listed as the possibility that the variation in the moving distance of the particles affects the above-described uniformity. When the above-described similarity is high, the center of one type of the above-described distribution and the center of the other type of the above-described distribution are substantially at the same position. Even in this case, when the above-described degree of dispersion of one kind of particles is high, the one kind of particles contains particles having a certain amount of movement from the center described above. And it is highly possible that such particles have lowered the above-described uniformity.

以上によれば、一の種類の粒子の上記した散布度及び他の種類の粒子の上記した散布度のそれぞれが基準値以下である場合に、上記した均一性を上記した類似度に基づいて正確に評価することができるといえる。   According to the above, when each of the above-mentioned distribution degree of one kind of particles and the above-mentioned degree of dispersion of other kinds of particles is equal to or less than a reference value, the above-described uniformity is accurately based on the above-described similarity degree. It can be said that it can be evaluated.

図12は、粒子の移動量の分布の類似度及び粒子の移動量の分布の散布度の評価方法の一例を説明するための図である。本図(a)は、一の種類の粒子(粒子群A)の移動距離の分布(確率密度)及び他の種類の粒子(粒子群B)の移動距離の分布(確率密度)を示している。本図(b)は、粒子群Aの上記した分布の累積確率及び粒子群Bの上記した分布の累積確率を示している。なお、粒子群Aの上記した累積確率は、粒子群Aの上記した確率密度の積分によって与えられる。同様に、粒子群Bの上記した累積確率は、粒子群Bの上記した確率密度の積分によって与えられる。   FIG. 12 is a diagram for explaining an example of a method for evaluating the degree of similarity of the distribution of particle movement amounts and the degree of distribution of the distribution of particle movement amounts. This figure (a) has shown the distribution (probability density) of the movement distance of one kind of particle | grains (particle group A), and the distribution (probability density) of the movement distance of another kind of particle | grains (particle group B). . This figure (b) has shown the cumulative probability of the above-mentioned distribution of the particle group A, and the cumulative probability of the above-mentioned distribution of the particle group B. FIG. The cumulative probability of the particle group A is given by the integral of the probability density of the particle group A. Similarly, the cumulative probability of the particle group B is given by the integral of the probability density of the particle group B.

本図(b)に示すように、粒子群Aの上記した分布と粒子群Bの上記した分布の類似度は、粒子群Aの上記した累積確率と粒子群Bの上記した累積確率を比較することにより、評価することができる。詳細には、粒子群Aの累積確率と粒子群Bの累積確率の差を各移動距離で算出する。そして上記した差の絶対値の最大値(コルモゴロフ−スミルノフ検定のD値)に基づいて上記した類似度を評価する。この場合、上記した最大値が小さいほど粒子群Aの上記した分布と粒子群Bの上記した分布の類似度が高いといえる。これにより、上記した最大値が所定の基準値以下であるか否かを判断することにより、粒子群Aと粒子群Bの混合状態を評価することができる。   As shown in FIG. 4B, the similarity between the above-described distribution of the particle group A and the above-described distribution of the particle group B compares the above-described cumulative probability of the particle group A with the above-described cumulative probability of the particle group B. This can be evaluated. Specifically, the difference between the cumulative probability of the particle group A and the cumulative probability of the particle group B is calculated for each moving distance. Then, the degree of similarity is evaluated based on the maximum absolute value of the difference (D value of Kolmogorov-Smirnov test). In this case, it can be said that the smaller the maximum value is, the higher the degree of similarity between the above distribution of the particle group A and the above distribution of the particle group B is. Thereby, the mixed state of the particle group A and the particle group B can be evaluated by determining whether or not the above-described maximum value is equal to or less than a predetermined reference value.

ただし、粒子群Aの上記した分布と粒子群Bの上記した分布の類似度の評価方法は、上記した例に限定されるものではない。例えば、本図(b)に示す例では、粒子群Aの累積確率と粒子群Bの累積確率の差の平均値を算出してもよい。この場合、上記した平均値が所定の基準値以下であるか否かを判断することにより、粒子群Aと粒子群Bの混合状態を評価することができる。   However, the method for evaluating the similarity between the above-described distribution of the particle group A and the above-described distribution of the particle group B is not limited to the above-described example. For example, in the example shown in FIG. 5B, an average value of the difference between the cumulative probability of the particle group A and the cumulative probability of the particle group B may be calculated. In this case, the mixed state of the particle group A and the particle group B can be evaluated by determining whether or not the above average value is equal to or less than a predetermined reference value.

本図(b)に示すように、粒子群Aの分布の散布度は、粒子群Aの上記した累積確率の傾きに基づいて評価することができる。この場合、例えば、累積確率が第1値以上第2値以下である範囲における累積確率の傾きの平均値に基づいて、散布度を評価することができる。なお、第1値及び第2値は、それぞれ、例えば、中央値(50%)未満の値(例えば、25%)及び中央値(50%)より大きい値(例えば、75%)である。他の例として、累積確率が基準確率値である位置における累積確率の傾きに基づいて、散布度を評価することができる。この場合、基準確率値は、例えば、累積確率の傾き(累積確率の一次導関数の値)が最大値をとる位置における累積確率の値であり、又は適当な特定の値(例えば、50%)である。なお、粒子群Bの分布の散布度も、上記と同様にして評価することができる。   As shown in FIG. 5B, the distribution degree of the distribution of the particle group A can be evaluated based on the slope of the cumulative probability of the particle group A described above. In this case, for example, the distribution degree can be evaluated based on the average value of the slopes of the cumulative probabilities in the range where the cumulative probability is greater than or equal to the first value and less than or equal to the second value. The first value and the second value are, for example, a value (for example, 25%) less than the median value (50%) and a value (for example, 75%) greater than the median value (50%), respectively. As another example, the distribution degree can be evaluated based on the slope of the cumulative probability at a position where the cumulative probability is the reference probability value. In this case, the reference probability value is, for example, the value of the cumulative probability at the position where the slope of the cumulative probability (the value of the first derivative of the cumulative probability) takes the maximum value, or an appropriate specific value (for example, 50%). It is. The distribution degree of the particle group B distribution can also be evaluated in the same manner as described above.

ただし、粒子群Aの上記した分布の散布度及び粒子群Bの上記した分布の散布度の評価方法は、上記した例に限定されるものではない。例えば、本図(a)に示すように、分布の標準偏差に基づいて、散布度を評価してもよい。   However, the evaluation method of the distribution degree of the particle group A and the distribution degree of the particle group B described above is not limited to the above example. For example, as shown in FIG. 5A, the distribution degree may be evaluated based on the standard deviation of the distribution.

以上、図面を参照して本発明の実施形態について述べたが、これらは本発明の例示であり、上記以外の様々な構成を採用することもできる。   As mentioned above, although embodiment of this invention was described with reference to drawings, these are the illustrations of this invention, Various structures other than the above are also employable.

以下、参考形態の例を付記する。
1.複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記複数の種類のうちの第1種類の粒子の前記移動量の分布と前記複数の種類のうちの第2種類の粒子の前記移動量の分布の類似度を評価する工程と、
前記類似度に基づいて前記第1種類の粒子と前記第2種類の粒子の混合状態を評価する工程と、
を含む評価方法。
2.1.に記載の評価方法において、
前記類似度を評価する工程では、
前記第1種類の粒子の前記移動量の前記分布の累積確率と前記第2種類の粒子の前記移動量の前記分布の累積確率を比較することにより、前記類似度を評価する評価方法。
3.1.又は2.に記載の評価方法において、
前記第1種類の粒子の前記移動量の分布の散布度及び前記第2種類の粒子の前記移動量の分布の散布度を評価する工程をさらに含み、
前記第1種類の粒子の前記散布度及び前記第2種類の粒子の前記散布度のそれぞれが基準値以下である場合に前記第1種類の粒子と前記第2種類の粒子の混合状態を評価する評価方法。
4.複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記複数の種類のうちの第1種類の粒子の前記移動量の分布と前記複数の種類のうちの第2種類の粒子の前記移動量の分布の類似度を評価する類似度評価部と、
前記類似度に基づいて前記第1種類の粒子と前記第2種類の粒子の混合状態を評価する評価部と、
を含む評価装置。
5.コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記複数の種類のうちの第1種類の粒子の前記移動量の分布と前記複数の種類のうちの第2種類の粒子の前記移動量の分布の類似度を評価する機能と、
前記類似度に基づいて前記第1種類の粒子と前記第2種類の粒子の混合状態を評価する機能と、
を持たせるプログラム。
以下、参考形態の例をさらに付記する。
1. 複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを統計的に処理された形式で表示する工程と、
を含む評価方法。
2. 1.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記データセットが統計的に処理されてグラフィカル形式で表示される評価方法。
3. 2.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられたグラフが表示される評価方法。
4. 3.に記載の評価方法であって、
前記グラフでは、前記各粒子を示すプロットが前記粒子の種類ごとに分類されている評価方法。
5. 2.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記各粒子が前記移動量に基づいて分類されているヒストグラムが表示される評価方法。
6. 1.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記データセットが統計的に処理されてテーブル形式で表示される評価方法。
7. 6.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられたテーブルが表示される評価方法。
8. 7.に記載の評価方法であって、
前記テーブルでは、前記各粒子と、当該粒子に対応する前記移動量と、を示すセルが前記粒子の種類ごとに分類されている評価方法。
9. 6.に記載の評価方法であって、
前記データセットを統計的に処理された形式で表示する前記工程では、前記各粒子が前記移動量に基づいて分類されている度数分布表が表示される評価方法。
10. 1.から9までのいずれか一つに記載の評価方法であって、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までの間に前記攪拌容器の内部で移動した距離である評価方法。
11. 1.から9までのいずれか一つに記載の評価方法であって、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離である評価方法。
12. 複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを統計的に処理された形式で表示する表示部と、
を含む評価装置。
13. コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを統計的に処理された形式で表示する機能と、
を持たせるプログラム。
Hereinafter, examples of the reference form will be added.
1. Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Evaluating the similarity between the movement amount distribution of the first type of the plurality of types of particles and the movement amount distribution of the second type of the plurality of types of particles;
Evaluating a mixed state of the first type particles and the second type particles based on the similarity;
Evaluation method including
2.1. In the evaluation method described in
In the step of evaluating the similarity,
An evaluation method for evaluating the similarity by comparing a cumulative probability of the distribution of the movement amount of the first type of particles with a cumulative probability of the distribution of the movement amount of the second type of particles.
3.1. Or 2. In the evaluation method described in
Further comprising the step of evaluating the distribution degree of the movement amount distribution of the first type particles and the distribution degree of the movement amount distribution of the second type particles,
The mixing state of the first-type particles and the second-type particles is evaluated when each of the first-type particles and the second-type particles is less than a reference value. Evaluation method.
4). A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A similarity evaluation unit that evaluates the similarity between the movement amount distribution of the first type of the plurality of types of particles and the movement amount distribution of the second type of the plurality of types of particles;
An evaluation unit that evaluates a mixed state of the first type particles and the second type particles based on the similarity;
Evaluation device including
5. A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of evaluating the similarity between the movement amount distribution of the first type of particles of the plurality of types and the movement amount distribution of the second type of particles of the plurality of types;
A function of evaluating a mixed state of the first type particles and the second type particles based on the similarity;
A program to give
Hereinafter, examples of the reference form will be additionally described.
1. Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
Displaying the data set in a statistically processed format;
Evaluation method including
2. 1. The evaluation method described in
In the step of displaying the data set in a statistically processed format, the evaluation method wherein the data set is statistically processed and displayed in a graphical format.
3. 2. The evaluation method described in
In the step of displaying the data set in a statistically processed format, an evaluation method in which a graph in which each particle is associated with the amount of movement corresponding to the particle is displayed.
4). 3. The evaluation method described in
In the graph, an evaluation method in which plots indicating the respective particles are classified for each type of the particles.
5. 2. The evaluation method described in
In the step of displaying the data set in a statistically processed format, an evaluation method in which a histogram in which each particle is classified based on the amount of movement is displayed.
6). 1. The evaluation method described in
In the step of displaying the data set in a statistically processed format, the evaluation method wherein the data set is statistically processed and displayed in a table format.
7). 6). The evaluation method described in
In the step of displaying the data set in a statistically processed format, an evaluation method in which a table in which each particle is associated with the amount of movement corresponding to the particle is displayed.
8). 7). The evaluation method described in
In the table, an evaluation method in which cells indicating the respective particles and the movement amount corresponding to the particles are classified for each type of the particles.
9. 6). The evaluation method described in
In the step of displaying the data set in a statistically processed format, an evaluation method in which a frequency distribution table in which the particles are classified based on the amount of movement is displayed.
10. 1. The evaluation method according to any one of 1 to 9,
The amount of movement of each particle is an evaluation method that is the distance that the particle has moved inside the stirring vessel from the first time to the second time.
11. 1. The evaluation method according to any one of 1 to 9,
The amount of movement of each particle is the distance that the particle has moved in a projected locus obtained by orthogonally projecting the locus of the particle from the first time to the second time on a predetermined plane or a predetermined straight line. Method.
12 A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display for displaying the data set in a statistically processed format;
Evaluation device including
13. A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a statistically processed format;
A program to give

100 評価装置
102 入力部
104 シミュレーション部
106 移動量算出部
108 データセット生成部
110 混合状態情報表示部
112 記憶部
114 評価部
116 類似度評価部
118 散布度評価部
DESCRIPTION OF SYMBOLS 100 Evaluation apparatus 102 Input part 104 Simulation part 106 Movement amount calculation part 108 Data set production | generation part 110 Mixed state information display part 112 Storage part 114 Evaluation part 116 Similarity evaluation part 118 Scatter degree evaluation part

Claims (18)

コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示部に表示させる工程と、
実行する評価方法。
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
The data set, wherein each particle, and the moving amount corresponding to the particles, a step of causing display on the display unit in the associated format,
Evaluation method to perform .
請求項1に記載の評価方法であって、
前記データセットを表示する前記工程では、前記データセットが統計的に処理されてグラフィカル形式で表示される評価方法。
The evaluation method according to claim 1,
Wherein in the data set table Shimesuru the process, evaluation process that the data set is displayed in a statistically processed by graphically.
請求項2に記載の評価方法であって、
前記データセットを表示する前記工程では、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられたグラフが表示される評価方法。
The evaluation method according to claim 2,
Wherein in the data set table Shimesuru said step, said each particle, the evaluation methods and the moving amount corresponding to the particles, the graph associated is displayed.
請求項3に記載の評価方法であって、
前記グラフでは、前記各粒子を示すプロットが前記粒子の種類ごとに分類されている評価方法。
The evaluation method according to claim 3, wherein
In the graph, an evaluation method in which plots indicating the respective particles are classified for each type of the particles.
請求項2に記載の評価方法であって、
前記データセットを表示する前記工程では、前記各粒子が前記移動量に基づいて分類されているヒストグラムが表示される評価方法。
The evaluation method according to claim 2,
Wherein in the data set table Shimesuru the process, evaluation process wherein each particle histograms are classified on the basis of the movement amount is displayed.
請求項1に記載の評価方法であって、
前記データセットを表示する前記工程では、前記データセットが統計的に処理されてテーブル形式で表示される評価方法。
The evaluation method according to claim 1,
Wherein in the data set table Shimesuru the process, evaluation process that the data set is displayed in a statistically processed by table format.
請求項6に記載の評価方法であって、
前記データセットを表示する前記工程では、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられたテーブルが表示される評価方法。
The evaluation method according to claim 6, wherein
Wherein in the data set table Shimesuru said step, said each particle, the evaluation methods and the moving amount corresponding to the particles, there is associated table is displayed.
請求項7に記載の評価方法であって、
前記テーブルでは、前記各粒子と、当該粒子に対応する前記移動量と、を示すセルが前記粒子の種類ごとに分類されている評価方法。
The evaluation method according to claim 7,
In the table, an evaluation method in which cells indicating the respective particles and the movement amount corresponding to the particles are classified for each type of the particles.
コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示部に表示させる工程と、
実行する評価方法。
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
The data set, and the moving amount of each of the particles, the frequency corresponding to the movement amount, a step of causing display on the display unit in the associated format,
Evaluation method to perform .
請求項に記載の評価方法であって、
前記データセットを表示する前記工程では、前記各粒子が前記移動量に基づいて分類されている度数分布表が表示される評価方法。
The evaluation method according to claim 9 , comprising:
Wherein in the data set table Shimesuru the process, the evaluation method of each particle frequency distribution table are classified on the basis of the movement amount is displayed.
請求項1から9までのいずれか一項に記載の評価方法であって、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までの間に前記攪拌容器の内部で移動した距離である評価方法。
The evaluation method according to any one of claims 1 to 9, wherein
The amount of movement of each particle is an evaluation method that is the distance that the particle has moved inside the stirring vessel from the first time to the second time.
コンピュータが、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする工程と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する工程と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する工程と、
前記データセットを統計的に処理された形式で表示部に表示させる工程と、
を実行し、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離である評価方法。
Computer
Simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
Calculating the amount of movement of each simulated particle from a first time to a second time;
Generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
Displaying the data set on a display in a statistically processed format;
Run
The movement amount of each particle is an evaluation of the distance traveled by the particle in a projected locus in which the locus traced from the first time to the second time is orthogonally projected onto a predetermined plane or a predetermined straight line. Method.
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示する表示部と、
を含む評価装置。
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display unit for displaying the data set in a format in which each particle and the movement amount corresponding to the particle are associated with each other;
Evaluation device including
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示する表示部と、
を含む評価装置。
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display unit for displaying the data set in a format in which the amount of movement of each particle and the frequency corresponding to the amount of movement are associated with each other;
Evaluation device including
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートするシミュレーション部と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する移動量算出部と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成するデータセット生成部と、
前記データセットを統計的に処理された形式で表示する表示部と、
を含み、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離である評価装置。
A simulation unit that simulates a state in which a plurality of types of particles are stirred inside the stirring vessel;
A movement amount calculation unit that calculates a movement amount of each simulated particle from the first time to the second time;
A data set generation unit for generating a data set including each particle and the movement amount corresponding to the particle in association with each other;
A display for displaying the data set in a statistically processed format;
Only including,
The amount of movement of each particle is the distance that the particle has moved in a projected locus obtained by orthogonally projecting the locus of the particle from the first time to the second time on a predetermined plane or a predetermined straight line. apparatus.
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを、前記各粒子と、当該粒子に対応する前記移動量と、が対応付けられた形式で表示する機能と、
を持たせるプログラム。
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a format in which each particle and the movement amount corresponding to the particle are associated with each other;
A program to give
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを、前記各粒子の前記移動量と、当該移動量に対応する頻度と、が対応付けられた形式で表示する機能と、
を持たせるプログラム。
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a format in which the movement amount of each particle and the frequency corresponding to the movement amount are associated with each other;
A program to give
コンピュータを評価装置として機能させるためのプログラムであって、
前記コンピュータに、
複数の種類の粒子が攪拌容器の内部で攪拌される状態をシミュレートする機能と、
シミュレートされた前記各粒子の第1時刻から第2時刻までの間の移動量を算出する機能と、
前記各粒子と、当該粒子に対応する前記移動量と、を対応付けて含むデータセットを生成する機能と、
前記データセットを統計的に処理された形式で表示する機能と、
を持たせ、
前記各粒子の移動量は、当該粒子が前記第1時刻から前記第2時刻までに辿った軌跡が所定の平面または所定の直線に正射影された射影軌跡において当該粒子が移動した距離であるプログラム。
A program for causing a computer to function as an evaluation device,
In the computer,
A function of simulating a state in which a plurality of types of particles are stirred inside the stirring vessel;
A function of calculating the amount of movement of each simulated particle from the first time to the second time;
A function of generating a data set including each particle and the amount of movement corresponding to the particle in association with each other;
A function of displaying the data set in a statistically processed format;
It was given,
The moving amount of each particle is a distance that the particle has moved in a projected locus obtained by orthogonally projecting a locus traced from the first time to the second time on a predetermined plane or a straight line. .
JP2014163084A 2013-09-04 2014-08-08 Evaluation method, evaluation apparatus, and program Active JP6282958B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014163084A JP6282958B2 (en) 2013-09-04 2014-08-08 Evaluation method, evaluation apparatus, and program

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2013182847 2013-09-04
JP2013182847 2013-09-04
JP2014163084A JP6282958B2 (en) 2013-09-04 2014-08-08 Evaluation method, evaluation apparatus, and program

Publications (2)

Publication Number Publication Date
JP2015072679A JP2015072679A (en) 2015-04-16
JP6282958B2 true JP6282958B2 (en) 2018-02-21

Family

ID=52628222

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014163084A Active JP6282958B2 (en) 2013-09-04 2014-08-08 Evaluation method, evaluation apparatus, and program

Country Status (2)

Country Link
JP (1) JP6282958B2 (en)
WO (1) WO2015033739A1 (en)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000214134A (en) * 1998-11-17 2000-08-04 Sharp Corp Method for simulating behavior of particle in container and behavior-simulating apparatus and its program recording medium
JP2007017380A (en) * 2005-07-11 2007-01-25 Fuji Xerox Co Ltd Particle simulation device, particle simulation method and computer program
JP5190251B2 (en) * 2007-11-19 2013-04-24 住友化学株式会社 Medium behavior simulation device of rotary rocking mill device
JP2010079493A (en) * 2008-09-25 2010-04-08 Fuji Xerox Co Ltd Particle behavior analyzer and program
JP2013029790A (en) * 2011-07-29 2013-02-07 Ricoh Co Ltd Toner production method and toner

Also Published As

Publication number Publication date
WO2015033739A1 (en) 2015-03-12
JP2015072679A (en) 2015-04-16

Similar Documents

Publication Publication Date Title
CN111859770B (en) Particle behavior simulation method and particle behavior simulation system
AU2017227323B2 (en) Particle simulation device, particle simulation method, and particle simulation program
US9477794B2 (en) Simulation method for high polymer material
US20130238302A1 (en) Simulation method for high polymer material
US7702490B1 (en) Method and system for adaptive mesh-free shell structures
Smeets et al. Modeling contact interactions between triangulated rounded bodies for the discrete element method
EP2808812A1 (en) Analysis device and simulation method
US11188692B2 (en) Turbulent boundary layer modeling via incorporation of pressure gradient directional effect
Ramos et al. Cargo dynamic stability in the container loading problem-a physics simulation tool approach
JP2011081530A (en) Discrete element method analysis simulation method for particle model, discrete element method analysis simulation program and discrete element method analysis simulation device
EP3809308A1 (en) Method for automatic calculation of axial cooling fan shroud circular opening size
Lu et al. A framework for problem standardization and algorithm comparison in multibody system
JP6282958B2 (en) Evaluation method, evaluation apparatus, and program
JP7246636B2 (en) Information processing device, particle simulator system, and particle simulator method
JP2016024177A (en) Method for creating model for specific substance analysis, computer program for creating model for specific substance analysis, method for simulating specific substance and computer program for simulating specific substance
JP6746971B2 (en) Composite material analysis method and computer program for composite material analysis
Hung et al. Automatic clustering method for real-time construction simulation
Mondesire et al. Physics engine benchmarking in three-dimensional virtual world simulation
EP3809307A1 (en) Method for automatic detection of axial cooling fan rotation direction
Araújo et al. A part complexity measurement method supporting 3D printing
Krstić et al. Using unity to study properties of a two-dimensional ideal gas
JP2002109445A (en) Behavior simulation method of particles in vessel, behavior simulation system and computer-readable recording medium
US20220414298A1 (en) Simulation method, simulation apparatus, and non-transitory computer readable medium storing program
WO2023157235A1 (en) Arithmetic device
JP5464160B2 (en) Particle presence field analysis apparatus, particle presence field analysis method, and computer program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20161221

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20161221

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20171219

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180112

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180123

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180125

R150 Certificate of patent or registration of utility model

Ref document number: 6282958

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250