JP7358863B2 - シミュレーション方法、及びシミュレーション装置 - Google Patents
シミュレーション方法、及びシミュレーション装置 Download PDFInfo
- Publication number
- JP7358863B2 JP7358863B2 JP2019164978A JP2019164978A JP7358863B2 JP 7358863 B2 JP7358863 B2 JP 7358863B2 JP 2019164978 A JP2019164978 A JP 2019164978A JP 2019164978 A JP2019164978 A JP 2019164978A JP 7358863 B2 JP7358863 B2 JP 7358863B2
- Authority
- JP
- Japan
- Prior art keywords
- reaction
- gas phase
- solid phase
- individual
- fixed carbon
- 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
Links
Images
Landscapes
- Manufacture And Refinement Of Metals (AREA)
Description
反応炉の一端側から供給した原料鉱石を他端側に向かって移動させながら、移動の途中で揮発分及び固定炭素の少なくとも一方を含む燃焼用材料を投入し、前記原料鉱石を前記他端側から供給された燃焼ガスと接触させて、乾燥させると共に還元を行うシミュレーション方法であって、
前記燃焼用材料を前記揮発分と前記固定炭素とに分配する分配工程と、
前記燃焼ガスを含む第1ガス相と分配された前記揮発分を含む第2ガス相のうち、反応速度が異なる第2ガス相ごとに、前記原料鉱石と分配された前記固定炭素とを含む固相との平衡反応に寄与する、第2ガス相の個別反応量をそれぞれ計算して、前記第2ガス相の前記平衡反応に寄与する第2ガス相個別反応分をそれぞれ求める第2ガス相個別反応量の計算工程と、
それぞれの前記第2ガス相個別反応分のうち、平衡計算で同一の化合物として取り扱う揮発分を合算して、揮発分合算反応分を求める揮発分の合算工程と、
前記原料鉱石を含む第1固相と分配された前記固定炭素を含む第2固相のうち、反応速度が異なる第2固相ごとに、前記平衡反応に寄与する第2固相の個別反応量をそれぞれ計算して、前記第2固相の前記平衡反応に寄与する第2固相個別反応分をそれぞれ求める第2固相個別反応量の計算工程と、
それぞれの前記第2固相個別反応分のうち、平衡計算で同一の化合物として取り扱う固定炭素を合算して、固定炭素合算反応分を求める固定炭素の合算工程と、
前記揮発分合算反応分を含むガス相と前記固定炭素合算反応分を含む固相とが平衡状態に達した時のそれぞれの熱量の変化及び流量を少なくとも計算する平衡反応計算工程と、
を含む。
図1は、一実施形態に係るシミュレーション方法が適用されるロータリーキルンの概略構成を示す。図1に示すように、ロータリーキルン1は、回転自在で略円筒形状のキルン本体11と、キルン本体11の途中に設けられる燃焼用材料供給管12とを有する。
次に、一実施形態に係るシミュレーション方法が適用される、一実施形態に係るシミュレーション装置について説明する。図2は、一実施形態に係るシミュレーション装置の機能を示すブロック図である。なお、図2では、シミュレーション装置が、ロータリーキルン1内を複数の領域に分割した時の、燃焼用材料供給管12から燃焼用材料が供給される領域Aにおける単位操作モデルとして説明する。また、図2では、領域Aのうち、燃焼ガスが吹き込まれる領域を領域A+1とし、原料鉱石が装入される領域を領域A-1とする。
以下の説明において、
「第1ガス相」とは、一方の隣接する領域(領域A+1)から領域Aに流入する燃焼ガスである。燃焼ガスは、揮発分の他に、酸素や二酸化炭素、一酸化炭素、水素など平衡反応に寄与するガス、さらに反応に寄与しない不活性な物質(例えば、窒素等)等を含んでもよい。
「第2ガス相」とは、燃焼用材料から分配された揮発分と第1ガス相とが合算されたガス相である。
「第3ガス相」とは、後述する、第2ガス相の混合計算モデルM6-1で生じるガス相である。
「第1固相」とは、一方の隣接する領域(領域A-1)から領域Aに流入する原料鉱石である。原料鉱石は、酸化ニッケルや酸化鉄等の鉱石中化合物や、固定炭素の他に、後述する平衡反応計算モデルM3において反応に寄与しない不活性な物質等を含んでもよい。
「第2固相」とは、燃焼用材料から分配された固定炭素と第1固相とが合算された固相である。
「第3固相」とは、後述する、第2固相の混合計算モデルM6-2で生じる固相である。
次に、一実施形態に係るシミュレーション装置を用いて、一実施形態に係るシミュレーション方法について説明する。一実施形態に係るシミュレーション方法は、図1に示すような構成を有するロータリーキルン1において、ロータリーキルン1の装入端14A側から供給した原料鉱石を排出端14B側に向かって移動させながら、移動の途中から燃焼用材料を投入し、原料鉱石を排出端14B側に設けられるバーナー16から供給された燃焼ガスと接触させて、乾燥させると共に還元を行う場合のシミュレーション方法である。
次に、シミュレーション装置のハードウェア構成の一例について説明する。図6は、シミュレーション装置のハードウェア構成図である。図6に示すように、シミュレーション装置20は、例えば、情報処理装置(コンピュータ)で構成され、物理的には、演算処理部であるCPU(Central Processing Unit:プロセッサ)21と、主記憶装置であるRAM(Random Access Memory)22及びROM(Read Only Memory)23と、補助記憶装置24と、入出力インタフェース25と、出力装置である表示装置26等を含むコンピュータシステムとして構成することができる。これらは、バス27で相互に接続されている。なお、補助記憶装置24及び表示装置26は、外部に設けられていてもよい。
すなわち、シミュレーションプログラムは、
反応炉の一端側から供給した原料鉱石を他端側に向かって移動させながら、移動の途中で揮発分及び固定炭素の少なくとも一方を含む燃焼用材料を投入し、前記原料鉱石を前記他端側から供給された燃焼ガスと接触させて、乾燥させると共に還元を行うシミュレーションをコンピュータに実行させるプログラムであって、
前記燃焼用材料を前記揮発分と前記固定炭素とに分配する分配モデルと、
前記燃焼ガスを含む第1ガス相と分配された前記揮発分を含む第2ガス相のうち、反応速度が異なる第2ガス相ごとに、前記原料鉱石と分配された前記固定炭素とを含む固相との平衡反応に寄与する、第2ガス相の個別反応量をそれぞれ計算して、前記第2ガス相の前記平衡反応に寄与する第2ガス相個別反応分をそれぞれ求める第2ガス相個別反応量の計算モデルと、
それぞれの前記第2ガス相個別反応分のうち、平衡計算で同一の化合物として取り扱う揮発分を合算して、揮発分合算反応分を求める揮発分の合算モデルと、
前記原料鉱石を含む第1固相と分配された前記固定炭素を含む第2固相のうち、反応速度が異なる第2固相ごとに、前記平衡反応に寄与する第2固相の個別反応量をそれぞれ計算して、前記第2固相の前記平衡反応に寄与する第2固相個別反応分をそれぞれ求める第2固相個別反応量の計算モデルと、
それぞれの前記第2固相個別反応分のうち、平衡計算で同一の化合物として取り扱う固定炭素を合算して、固定炭素合算反応分を求める固定炭素の合算モデルと、
前記揮発分合算反応分を含むガス相と前記固定炭素合算反応分を含む固相とが平衡状態に達した時のそれぞれの熱量の変化及び流量を少なくとも計算する平衡反応計算モデルと、
を少なくともコンピュータに実行させるプログラムを用いることができる。
11 キルン本体
14 バーナー
20 シミュレーション装置
M1 燃焼用材料の分配モデル(分配モデル)
M2-1 第2ガス相個別反応量の計算モデル
M2-2 第2固相個別反応量の計算モデル
M3-1 揮発分の合算モデル
M3-2 固定炭素の合算モデル
M4 平衡反応計算モデル
M5-1 揮発分合算生成分の分解モデル
M5-2 固定炭素合算生成分の分解モデル
M6-1 第2ガス相の混合計算モデル
M6-2 第2固相の混合計算モデル
G1 第1ガス相
G2 第2ガス相
G3 第3ガス相
S1 第1固相
S2 第2固相
S3 第3固相
VM11・・・VM2M 揮発分個別反応分
VM1Mtotal 揮発分合算反応分
VM21・・・VM2M 揮発分個別生成分
VM2Mtotal 揮発分合算生成分
FC21・・・FC2M 固定炭素個別反応分
FC1Mtotal 固定炭素合算反応分
FC2Mtotal 固定炭素合算生成分
Claims (5)
- ロータリーキルンのキルン本体の一端側である開口端部側から供給した原料鉱石を他端側である排出端部側に向かって、前記キルン本体の回転に連れて、前記キルン本体内を移動させながら、移動の途中に揮発分及び固定炭素の少なくとも一方を含む燃焼用材料を前記キルン本体内に投入し、前記原料鉱石を、前記他端側である前記排出端部に設置したバーナで発生させた燃焼ガスと向流接触させる際の、物質の挙動を解析するシミュレーション方法であって、
前記燃焼用材料を前記揮発分と前記固定炭素とに分配する分配工程と、
前記燃焼ガスである第1ガス相と分配された前記揮発分とからなる第2ガス相のうち、反応速度が異なる第2ガス相ごとに、前記原料鉱石と分配された前記固定炭素とを含む固相との平衡反応に寄与する、第2ガス相の個別反応量をそれぞれ計算して、前記第2ガス相の前記平衡反応に寄与する第2ガス相個別反応分をそれぞれ求める第2ガス相個別反応量の計算工程と、
それぞれの前記第2ガス相個別反応分のうち、平衡計算で同一の化合物として取り扱う揮発分を合算して、揮発分合算反応分を求める揮発分の合算工程と、
前記原料鉱石である第1固相と分配された前記固定炭素とからなる第2固相のうち、反応速度が異なる第2固相ごとに、前記平衡反応に寄与する前記第2固相の個別反応量をそれぞれ計算して、前記第2固相の前記平衡反応に寄与する第2固相個別反応分をそれぞれ求める第2固相個別反応量の計算工程と、
それぞれの前記第2固相個別反応分のうち、平衡計算で同一の化合物として取り扱う前記固定炭素を合算して、固定炭素合算反応分を求める固定炭素の合算工程と、
前記揮発分合算反応分を含むガス相と前記固定炭素合算反応分を含む固相とが平衡状態に達した時のそれぞれの熱量の変化及び流量を少なくとも計算する平衡反応計算工程と、
を含むシミュレーション方法。 - 前記平衡反応計算工程で前記第2ガス相個別反応分と前記第2固相個別反応分とが反応することで生じる第2ガス相個別生成分に含まれる揮発分合算生成分を、前記反応速度が異なる前記第2ガス相の個別反応量の割合に分解し、第2ガス相個別生成分を得る揮発分合算生成分の分解工程と、
前記平衡反応計算工程で前記第2ガス相個別反応分と前記第2固相個別反応分とが反応することで生じる第2固相個別生成分に含まれる固定炭素合算生成分を、前記反応速度の異なる第2固相の個別反応量の割合に分解し、第2固相個別生成分を得る固定炭素合算生成分の分解工程と、
を含む請求項1に記載のシミュレーション方法。 - 前記第2ガス相個別反応量の計算工程で、前記第2ガス相個別反応分の和以外の未反応分を第2ガス相個別未反応分とし、
前記第2固相個別反応量の計算工程で、前記第2固相個別反応分の和以外の未反応分を、第2固相個別未反応分とし、
前記第2ガス相個別未反応分と、前記平衡反応計算工程で前記第2ガス相個別反応分と前記第2固相個別反応分とが反応することで生じる第2ガス相個別生成分とを混合した第3ガス相の流量を少なくとも計算する第1の混合計算工程と、
前記第2固相個別未反応分と、前記平衡反応計算工程で前記第2ガス相個別反応分と前記第2固相個別反応分とが反応することで生じる第2固相個別生成分とを混合した第3固相の流量を少なくとも計算する第2の混合計算工程と、
を含む請求項1又は2に記載のシミュレーション方法。 - 前記キルン本体内をその長軸方向に沿って複数の領域に分割すると仮定した時、
シミュレーション方法は、前記複数の領域に、前記原料鉱石の流れ又は前記燃焼ガスの流れに沿って繰り返し行い、所定の領域における計算値とその領域における前回の計算値との差が所定の範囲内に収まるまで繰り返し行う請求項1~3の何れか一項に記載のシミュレーション方法。 - ロータリーキルンのキルン本体の一端側である開口端部側から供給した原料鉱石を他端側である排出端部側に向かって、前記キルン本体の回転に連れて、前記キルン本体内を移動させながら、移動の途中に揮発分及び固定炭素の少なくとも一方を含む燃焼用材料を前記キルン本体内に投入し、前記原料鉱石を、前記他端側である前記排出端部に設置したバーナで発生させた燃焼ガスと向流接触させる際の、物質の挙動を解析するシミュレーション装置であって、
前記燃焼用材料を前記揮発分と前記固定炭素とに分配する分配モデルと、
前記燃焼ガスである第1ガス相と分配された前記揮発分とからなる第2ガス相のうち、反応速度が異なる第2ガス相ごとに、前記原料鉱石と分配された前記固定炭素とを含む固相との平衡反応に寄与する、第2ガス相の個別反応量をそれぞれ計算して、前記第2ガス相の前記平衡反応に寄与する第2ガス相個別反応分をそれぞれ求める第2ガス相個別反応量の計算モデルと、
それぞれの前記第2ガス相個別反応分のうち、平衡計算で同一の化合物として取り扱う揮発分を合算して、揮発分合算反応分を求める揮発分の合算モデルと、
前記原料鉱石である第1固相と分配された前記固定炭素とからなる第2固相のうち、反応速度が異なる第2固相ごとに、前記平衡反応に寄与する前記第2固相の個別反応量をそれぞれ計算して、前記第2固相の前記平衡反応に寄与する第2固相個別反応分をそれぞれ求める第2固相個別反応量の計算モデルと、
それぞれの前記第2固相個別反応分のうち、平衡計算で同一の化合物として取り扱う前記固定炭素を合算して、固定炭素合算反応分を求める固定炭素の合算モデルと、
前記揮発分合算反応分を含むガス相と前記固定炭素合算反応分を含む固相とが平衡状態に達した時のそれぞれの熱量の変化及び流量を少なくとも計算する平衡反応計算モデルと、
を備えるシミュレーション装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019164978A JP7358863B2 (ja) | 2019-09-10 | 2019-09-10 | シミュレーション方法、及びシミュレーション装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019164978A JP7358863B2 (ja) | 2019-09-10 | 2019-09-10 | シミュレーション方法、及びシミュレーション装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021042428A JP2021042428A (ja) | 2021-03-18 |
JP7358863B2 true JP7358863B2 (ja) | 2023-10-11 |
Family
ID=74862102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019164978A Active JP7358863B2 (ja) | 2019-09-10 | 2019-09-10 | シミュレーション方法、及びシミュレーション装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7358863B2 (ja) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006016216A (ja) | 2004-06-30 | 2006-01-19 | Mitsubishi Materials Corp | シミュレーション方法、シミュレーションプログラムおよび記憶媒体 |
JP2009103357A (ja) | 2007-10-23 | 2009-05-14 | Sumitomo Metal Mining Co Ltd | ロータリーキルンの操業方法。 |
JP2014141719A (ja) | 2013-01-25 | 2014-08-07 | Hyuga Seirensho:Kk | ロータリーキルンの操業方法 |
KR101930680B1 (ko) | 2018-11-26 | 2018-12-18 | 주식회사 에스엔엔씨 | 타이어 분말 칩을 이용한 니켈 제련장치 및 제련방법 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03183719A (ja) * | 1989-12-11 | 1991-08-09 | Nippon Steel Corp | 製鋼工程におけるマンガン鉱石の還元促進方法 |
-
2019
- 2019-09-10 JP JP2019164978A patent/JP7358863B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006016216A (ja) | 2004-06-30 | 2006-01-19 | Mitsubishi Materials Corp | シミュレーション方法、シミュレーションプログラムおよび記憶媒体 |
JP2009103357A (ja) | 2007-10-23 | 2009-05-14 | Sumitomo Metal Mining Co Ltd | ロータリーキルンの操業方法。 |
JP2014141719A (ja) | 2013-01-25 | 2014-08-07 | Hyuga Seirensho:Kk | ロータリーキルンの操業方法 |
KR101930680B1 (ko) | 2018-11-26 | 2018-12-18 | 주식회사 에스엔엔씨 | 타이어 분말 칩을 이용한 니켈 제련장치 및 제련방법 |
Also Published As
Publication number | Publication date |
---|---|
JP2021042428A (ja) | 2021-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mujumdar et al. | Rotary Cement Kiln Simulator (RoCKS): Integrated modeling of pre-heater, calciner, kiln and clinker cooler | |
Nishioka et al. | Sustainable aspects of CO 2 ultimate reduction in the steelmaking process (COURSE50 Project), part 1: Hydrogen reduction in the blast furnace | |
Zhuo et al. | CFD study of hydrogen injection through tuyeres into ironmaking blast furnaces | |
Marias et al. | Modelling of a rotary kiln for the pyrolysis of aluminium waste | |
Gaurav et al. | Analysis of temperature profile and% metallization in rotary kiln of sponge iron process through CFD | |
Neto et al. | Optimizing the coke oven process by adjusting the temperature of the combustion chambers | |
JP7358863B2 (ja) | シミュレーション方法、及びシミュレーション装置 | |
Liu et al. | Computational fluid dynamics study of re-blowin operation in an ironmaking blast furnace | |
JP7358862B2 (ja) | シミュレーション方法、及びシミュレーション装置 | |
JP7358865B2 (ja) | シミュレーション方法、及びシミュレーション装置 | |
JP7358864B2 (ja) | シミュレーション方法、及びシミュレーション装置 | |
Ryan et al. | Packed bed combustion of char particles: experiments and an ash model | |
Zhang et al. | CFD-DEM simulations of municipal solid waste gasification in a pilot-scale direct-melting furnace | |
Farajzadeh et al. | Numerical investigation of the impact of injecting coke oven gas on raceway evolution in an ironmaking blast furnace | |
Zhang et al. | Numerical investigation of hydrogen-rich gas and pulverized coal injection in the raceway of a blast furnace with lower carbon emissions | |
Xu et al. | A novel CFD-DEM-DPM modelling of fluid-particles-fines reacting flows | |
JP2022144573A (ja) | 炉内反応の計算装置及び炉内反応の計算方法 | |
JP2023004584A (ja) | 炉内反応の計算装置及び炉内反応の計算方法 | |
JP2022143250A (ja) | 炉内反応の計算装置及び炉内反応の計算方法 | |
JP2022143251A (ja) | 炉内反応の計算装置及び炉内反応の計算方法 | |
JP2023103125A (ja) | 炉内反応の計算装置及び炉内反応の計算方法 | |
Runkana | Model-based optimization of industrial gas-solid reactors | |
Yu et al. | Numerical study of transient thermochemical states inside an ironmaking blast furnace: Impacts of blast temperature drop and recovery | |
JP4501554B2 (ja) | シミュレーションプログラムおよび記憶媒体 | |
van Staden et al. | Damring formation during rotary kiln chromite pre-reduction: Effects of pulverized carbonaceous fuel selection and partial pellet melting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220418 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230216 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230307 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230428 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230801 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230809 |
|
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: 20230829 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230911 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7358863 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |