JP2016538536A5 - - Google Patents
Download PDFInfo
- Publication number
- JP2016538536A5 JP2016538536A5 JP2016526276A JP2016526276A JP2016538536A5 JP 2016538536 A5 JP2016538536 A5 JP 2016538536A5 JP 2016526276 A JP2016526276 A JP 2016526276A JP 2016526276 A JP2016526276 A JP 2016526276A JP 2016538536 A5 JP2016538536 A5 JP 2016538536A5
- Authority
- JP
- Japan
- Prior art keywords
- polymerization
- polymerization reaction
- control parameter
- mixture
- probe
- 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.)
- Granted
Links
- 238000006116 polymerization reaction Methods 0.000 claims description 28
- 239000000523 sample Substances 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 14
- 239000000203 mixture Substances 0.000 claims description 12
- 239000002245 particle Substances 0.000 claims description 9
- 238000000691 measurement method Methods 0.000 claims description 8
- 238000004164 analytical calibration Methods 0.000 claims description 6
- 238000011088 calibration curve Methods 0.000 claims description 6
- 229920000642 polymer Polymers 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 125000000217 alkyl group Chemical group 0.000 claims description 4
- 239000011541 reaction mixture Substances 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 239000004793 Polystyrene Substances 0.000 claims description 2
- 238000010928 TGA analysis Methods 0.000 claims description 2
- 238000007563 acoustic spectroscopy Methods 0.000 claims description 2
- 125000003545 alkoxy group Chemical group 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 claims description 2
- 239000007900 aqueous suspension Substances 0.000 claims description 2
- 125000003118 aryl group Chemical group 0.000 claims description 2
- 230000001276 controlling effect Effects 0.000 claims description 2
- 230000000875 corresponding Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims description 2
- 229910052736 halogen Inorganic materials 0.000 claims description 2
- 150000002367 halogens Chemical class 0.000 claims description 2
- 229910052739 hydrogen Inorganic materials 0.000 claims description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 2
- 239000001257 hydrogen Substances 0.000 claims description 2
- 230000001678 irradiating Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 claims description 2
- 125000002496 methyl group Chemical group [H]C([H])([H])* 0.000 claims description 2
- 239000000178 monomer Substances 0.000 claims description 2
- 229920002223 polystyrene Polymers 0.000 claims description 2
- 238000002798 spectrophotometry method Methods 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 238000002411 thermogravimetry Methods 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 239000004711 α-olefin Substances 0.000 claims description 2
- 238000010191 image analysis Methods 0.000 claims 1
- 238000000985 reflectance spectrum Methods 0.000 claims 1
- 239000011324 bead Substances 0.000 description 2
- 230000012010 growth Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 230000003595 spectral Effects 0.000 description 1
Description
上記結果は、本発明の方法が上記平均粒度における変動に対して感度が高く、またそれ故に上記重合反応中のこれらの成長を、効果的にしかも期待通りに監視するのに適していることを示している。
その上、上記バッチサンプルの重合サイクルは、上記プローブの如何なるクリーニングまたはメンテナンスの介入をも必要とすることなしに、同一のオートクレーブ内で実施された。
上記キャリブレーションの検証後、本発明の方法を、上記EPSビーズの成長を、その重合中に連続的に監視するのに使用した。図5は、該重合反応中にインラインで測定されたNIRスペクトルの幾つかを示しており、本図において、約300μm(下方のスペクトルA)乃至約1.2mm(上方のスペクトルB)に広がっている、該EPSビーズの平均寸法における増加と関連付けられるスペクトルベースラインの高まりを明確に見ることができる。
本発明の好ましい態様は、下記の通りである。
〔1〕不均質相における重合反応混合物の制御パラメータを監視する方法であって、以下の工程、
(a) 該混合物の少なくとも一つのNIR反射スペクトルを得る工程、
(b) 検量線を用いて該制御パラメータの値を計算する工程であって、該検量線は、該NIR反射スペクトルと、基準測定法で測定された該制御パラメータの値とを相互に関連付けている、工程、
を含む、方法。
〔2〕前記重合反応が、一般式:R 1 R 2 C=CH 2 を持つα-オレフィンから選択される1種以上のモノマーの重合反応であり、
R 1 は水素またはメチルであり、
R 2 は、ハロゲン、C 1 −C 4 アルキル及びC 1 −C 4 アルコキシルから選択される1種以上の基で置換されていてもよい、C 1 −C 10 アルキル、C 1 −C 6 アリールから選択される基である、
前記〔1〕に記載の方法。
〔3〕前記重合反応が、水性懸濁液中でのポリスチレンの合成反応である、前記〔2〕に記載の方法。
〔4〕前記制御パラメータが、前記重合反応中に形成されるポリマー粒子の平均サイズである、前記〔1〕〜〔3〕の1つ以上に記載の方法。
〔5〕前記平均粒子サイズが、100μm〜3mm、好ましくは300μm〜1.5mmの範囲内にある、前記〔4〕に記載の方法。
〔6〕前記基準測定法が、音響分光法、分光光度法、画像解析法及び前記ポリマー粒子と、既知寸法を持つ同一のポリマーの粒子に係る標準的なサンプルとの視覚的比較から選択される、前記〔4〕又は〔5〕に記載の方法。
〔7〕前記制御パラメータが、前記重合混合物の少なくとも一つのモノマーの転化度である、前記〔1〕〜〔6〕の1つ以上に記載の方法。
〔8〕前記基準測定法が、屈折率分析又は熱重量分析から選択される、前記〔7〕に記載の方法。
〔9〕前記検量線が、多変量線形回帰法によって、複数のNIR反射スペクトルと対応する複数の前記制御パラメータの値とを相互に関連付けることにより得られ、該スペクトルが、前記重合反応の様々な進行度において、一連のサンプル重合混合物について取得され、該制御パラメータの値は、基準測定法を用いて該サンプル重合混合物について測定される、前記〔1〕〜〔8〕の1つ以上に記載の方法。
〔10〕前記多変量線形回帰法が、多重最小二乗法、部分最小二乗法、主成分による方法およびこれらの組合せから選択される、前記〔1〕〜〔9〕の1つ以上に記載の方法。
〔11〕前記NIR反射スペクトルが、波数4,000〜15,000cm -1 の範囲内、好ましくは6,000〜10,000cm -1 の範囲内で取得される、前記〔1〕〜〔10〕の1つ以上に記載の方法。
〔12〕不均質相における重合反応の制御方法であって、該反応に係る少なくとも一つの制御パラメータを、前記〔1〕〜〔11〕の1つ以上に従って監視する段階を含む、方法。
〔13〕前記〔1〕に記載の方法を実施するための装置であって、NIR反射スペクトルの取得システムを備える、少なくとも一つの重合反応器9を含み、該取得システムが、
該反応器9内に含まれる重合反応混合物8上に光放射線を照射し、かつ実質上、該照射に応答して該混合物8から反射される放射線のみを集めるための、少なくとも一つのプローブ1、
該反射された放射線を検出するための、該プローブ1と光学的に結合されている少なくとも一つの検出システム5、
を含む、前記装置。
〔14〕前記プローブ1が、少なくとも一つの測定ヘッド6及び少なくとも一つのサンプリングウインドウ7を含み、該プローブ1が、前記反応器9の壁11内に含まれ、該サンプリングウインドウ7の内側表面が、実質的に該壁11の内側表面と位置合わせされている、前記〔13〕に記載の装置。
その上、上記バッチサンプルの重合サイクルは、上記プローブの如何なるクリーニングまたはメンテナンスの介入をも必要とすることなしに、同一のオートクレーブ内で実施された。
上記キャリブレーションの検証後、本発明の方法を、上記EPSビーズの成長を、その重合中に連続的に監視するのに使用した。図5は、該重合反応中にインラインで測定されたNIRスペクトルの幾つかを示しており、本図において、約300μm(下方のスペクトルA)乃至約1.2mm(上方のスペクトルB)に広がっている、該EPSビーズの平均寸法における増加と関連付けられるスペクトルベースラインの高まりを明確に見ることができる。
本発明の好ましい態様は、下記の通りである。
〔1〕不均質相における重合反応混合物の制御パラメータを監視する方法であって、以下の工程、
(a) 該混合物の少なくとも一つのNIR反射スペクトルを得る工程、
(b) 検量線を用いて該制御パラメータの値を計算する工程であって、該検量線は、該NIR反射スペクトルと、基準測定法で測定された該制御パラメータの値とを相互に関連付けている、工程、
を含む、方法。
〔2〕前記重合反応が、一般式:R 1 R 2 C=CH 2 を持つα-オレフィンから選択される1種以上のモノマーの重合反応であり、
R 1 は水素またはメチルであり、
R 2 は、ハロゲン、C 1 −C 4 アルキル及びC 1 −C 4 アルコキシルから選択される1種以上の基で置換されていてもよい、C 1 −C 10 アルキル、C 1 −C 6 アリールから選択される基である、
前記〔1〕に記載の方法。
〔3〕前記重合反応が、水性懸濁液中でのポリスチレンの合成反応である、前記〔2〕に記載の方法。
〔4〕前記制御パラメータが、前記重合反応中に形成されるポリマー粒子の平均サイズである、前記〔1〕〜〔3〕の1つ以上に記載の方法。
〔5〕前記平均粒子サイズが、100μm〜3mm、好ましくは300μm〜1.5mmの範囲内にある、前記〔4〕に記載の方法。
〔6〕前記基準測定法が、音響分光法、分光光度法、画像解析法及び前記ポリマー粒子と、既知寸法を持つ同一のポリマーの粒子に係る標準的なサンプルとの視覚的比較から選択される、前記〔4〕又は〔5〕に記載の方法。
〔7〕前記制御パラメータが、前記重合混合物の少なくとも一つのモノマーの転化度である、前記〔1〕〜〔6〕の1つ以上に記載の方法。
〔8〕前記基準測定法が、屈折率分析又は熱重量分析から選択される、前記〔7〕に記載の方法。
〔9〕前記検量線が、多変量線形回帰法によって、複数のNIR反射スペクトルと対応する複数の前記制御パラメータの値とを相互に関連付けることにより得られ、該スペクトルが、前記重合反応の様々な進行度において、一連のサンプル重合混合物について取得され、該制御パラメータの値は、基準測定法を用いて該サンプル重合混合物について測定される、前記〔1〕〜〔8〕の1つ以上に記載の方法。
〔10〕前記多変量線形回帰法が、多重最小二乗法、部分最小二乗法、主成分による方法およびこれらの組合せから選択される、前記〔1〕〜〔9〕の1つ以上に記載の方法。
〔11〕前記NIR反射スペクトルが、波数4,000〜15,000cm -1 の範囲内、好ましくは6,000〜10,000cm -1 の範囲内で取得される、前記〔1〕〜〔10〕の1つ以上に記載の方法。
〔12〕不均質相における重合反応の制御方法であって、該反応に係る少なくとも一つの制御パラメータを、前記〔1〕〜〔11〕の1つ以上に従って監視する段階を含む、方法。
〔13〕前記〔1〕に記載の方法を実施するための装置であって、NIR反射スペクトルの取得システムを備える、少なくとも一つの重合反応器9を含み、該取得システムが、
該反応器9内に含まれる重合反応混合物8上に光放射線を照射し、かつ実質上、該照射に応答して該混合物8から反射される放射線のみを集めるための、少なくとも一つのプローブ1、
該反射された放射線を検出するための、該プローブ1と光学的に結合されている少なくとも一つの検出システム5、
を含む、前記装置。
〔14〕前記プローブ1が、少なくとも一つの測定ヘッド6及び少なくとも一つのサンプリングウインドウ7を含み、該プローブ1が、前記反応器9の壁11内に含まれ、該サンプリングウインドウ7の内側表面が、実質的に該壁11の内側表面と位置合わせされている、前記〔13〕に記載の装置。
Claims (14)
- 不均質相における重合反応混合物の制御パラメータを監視する方法であって、以下の工程、
(a) 該混合物の少なくとも一つのNIR反射スペクトルを得る工程、
(b) 検量線を用いて該制御パラメータの値を計算する工程であって、該検量線は、該NIR反射スペクトルと、基準測定法で測定された該制御パラメータの値とを相互に関連付けている、工程、
を含む、方法。 - 前記重合反応が、一般式:R1R2C=CH2を持つα-オレフィンから選択される1種以上のモノマーの重合反応であり、
R1は水素またはメチルであり、
R2は、ハロゲン、C1−C4アルキル及びC1−C4アルコキシルから選択される1種以上の基で置換されていてもよい、C1−C10アルキル、C1−C6アリールから選択される基である、
請求項1に記載の方法。 - 前記重合反応が、水性懸濁液中でのポリスチレンの合成反応である、請求項2に記載の方法。
- 前記制御パラメータが、前記重合反応中に形成されるポリマー粒子の平均サイズである、請求項1〜3のいずれか1項に記載の方法。
- 前記平均粒子サイズが、100μm〜3mmの範囲内にある、請求項4に記載の方法。
- 前記基準測定法が、音響分光法、分光光度法、画像解析法、及び前記ポリマー粒子と既知寸法を持つ同一のポリマーの粒子に係る標準的なサンプルとの視覚的比較、から選択される、請求項4又は5に記載の方法。
- 前記制御パラメータが、前記重合混合物の少なくとも一つのモノマーの転化度である、請求項1〜6のいずれか1項に記載の方法。
- 前記基準測定法が、屈折率分析又は熱重量分析から選択される、請求項7に記載の方法。
- 前記検量線が、多変量線形回帰法によって、複数のNIR反射スペクトルと対応する複数の前記制御パラメータの値とを相互に関連付けることにより得られ、該スペクトルが、前記重合反応の様々な進行度において、一連のサンプル重合混合物について取得され、該制御パラメータの値は、基準測定法を用いて該サンプル重合混合物について測定される、請求項1〜8のいずれか1項に記載の方法。
- 前記多変量線形回帰法が、多重最小二乗法、部分最小二乗法、主成分による方法およびこれらの組合せから選択される、請求項1〜9のいずれか1項に記載の方法。
- 前記NIR反射スペクトルが、波数4,000〜15,000cm-1の範囲内で取得される、請求項1〜10のいずれか1項に記載の方法。
- 不均質相における重合反応の制御方法であって、該反応に係る少なくとも一つの制御パラメータを、請求項1〜11のいずれか1項に従って監視する段階を含む、方法。
- 請求項1に記載の方法を実施するための装置であって、NIR反射スペクトルの取得システムを備える、少なくとも一つの重合反応器9を含み、該取得システムが、
該反応器9内に含まれる重合反応混合物8上に光放射線を照射し、かつ実質上、該照射に応答して該混合物8から反射される放射線のみを集めるための、少なくとも一つのプローブ1、
該反射された放射線を検出するための、該プローブ1と光学的に結合されている少なくとも一つの検出システム5、
を含む、前記装置。 - 前記プローブ1が、少なくとも一つの測定ヘッド6及び少なくとも一つのサンプリングウインドウ7を含み、該プローブ1が、前記反応器9の壁11内に含まれ、該サンプリングウインドウ7の内側表面が、実質的に該壁11の内側表面と位置合わせされている、請求項13に記載の装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ITMI2013A001916 | 2013-11-19 | ||
IT001916A ITMI20131916A1 (it) | 2013-11-19 | 2013-11-19 | Metodo per monitorare un parametro di controllo di una reazione di polimerizzazione e relativo apparato per attuare detto metodo |
PCT/IB2014/066129 WO2015075629A1 (en) | 2013-11-19 | 2014-11-18 | Method for monitoring a control parameter of a polymerization reaction and relative apparatus for implementing said method |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2016538536A JP2016538536A (ja) | 2016-12-08 |
JP2016538536A5 true JP2016538536A5 (ja) | 2017-12-28 |
JP6619736B2 JP6619736B2 (ja) | 2019-12-11 |
Family
ID=49841733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016526276A Expired - Fee Related JP6619736B2 (ja) | 2013-11-19 | 2014-11-18 | 重合反応に係る制御パラメータを監視するための方法および該方法を実施するための関連する装置 |
Country Status (17)
Country | Link |
---|---|
US (1) | US9714899B2 (ja) |
EP (1) | EP3071954B1 (ja) |
JP (1) | JP6619736B2 (ja) |
KR (1) | KR102209231B1 (ja) |
CN (1) | CN105659072B (ja) |
BR (1) | BR112016007745B1 (ja) |
ES (1) | ES2804087T3 (ja) |
HU (1) | HUE051229T2 (ja) |
IT (1) | ITMI20131916A1 (ja) |
MX (1) | MX368418B (ja) |
MY (1) | MY185715A (ja) |
PL (1) | PL3071954T3 (ja) |
PT (1) | PT3071954T (ja) |
RU (1) | RU2662039C1 (ja) |
SA (1) | SA516371040B1 (ja) |
SI (1) | SI3071954T1 (ja) |
WO (1) | WO2015075629A1 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201704447D0 (en) * | 2017-03-21 | 2017-05-03 | Asterope Ltd | Wire coating determination |
WO2019235542A1 (ja) * | 2018-06-07 | 2019-12-12 | 横河電機株式会社 | 光学分析システム及び光学分析方法 |
CN108872144B (zh) * | 2018-09-13 | 2021-02-12 | 中国农业大学 | 一种用于厌氧消化过程的在线监测装置 |
EP3711852A1 (de) * | 2019-03-21 | 2020-09-23 | Covestro Deutschland AG | Verfahren und vorrichtung zur quantitativen überwachung der zusammensetzung einer oligomer-/monomermischung |
CN117309793B (zh) * | 2023-10-25 | 2024-05-07 | 连云港石化有限公司 | 一种苯乙烯装置聚合物测定方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5604132A (en) * | 1995-01-23 | 1997-02-18 | Olin Corporation | Process flow injection analyzer and method |
BE1009667A3 (fr) * | 1995-09-25 | 1997-06-03 | Solvay | Procede et dispositif de controle de la qualite de latex de polymeres vinyliques halogenes. |
DE19805612A1 (de) * | 1998-02-12 | 1999-08-19 | Bayer Ag | Verfahren zur kontrollierten Herstellung oder Modifizierung von polymeren Produkten mittels IR-ATR-Spektroskopie |
JP2000159897A (ja) * | 1998-11-27 | 2000-06-13 | Mitsui Chemicals Inc | 高分子化合物の不良品分離方法 |
JP2002194079A (ja) * | 2000-12-22 | 2002-07-10 | Mitsubishi Gas Chem Co Inc | ポリアミドの製造方法 |
US6864331B1 (en) * | 2003-12-09 | 2005-03-08 | Fina Technology, Inc. | Process for the production of polymers |
JP5255833B2 (ja) * | 2004-03-24 | 2013-08-07 | エクソンモービル・ケミカル・パテンツ・インク | エチレン・インターポリマーの製造方法、該方法により製造されたインターポリマー、組成物、および該インターポリマーを含む電気デバイス |
CA2470887C (en) * | 2004-06-14 | 2014-10-21 | Nova Chemicals Corporation | Detection of gels in a solution polymerization |
EP2325216B1 (en) | 2006-03-13 | 2014-05-21 | LG Chem, Ltd. | Method for preparation of 1-alkene-acrylate based copolymer |
US8843324B2 (en) * | 2011-02-03 | 2014-09-23 | Nova Chemicals (International) S.A. | Double derivative NIR process control |
CN102288575B (zh) * | 2011-07-25 | 2013-03-06 | 哈尔滨工业大学 | 在线监测甲氨基阿维菌素反应的装置 |
-
2013
- 2013-11-19 IT IT001916A patent/ITMI20131916A1/it unknown
-
2014
- 2014-11-18 EP EP14806452.0A patent/EP3071954B1/en active Active
- 2014-11-18 PT PT148064520T patent/PT3071954T/pt unknown
- 2014-11-18 CN CN201480055493.7A patent/CN105659072B/zh active Active
- 2014-11-18 RU RU2016116157A patent/RU2662039C1/ru active
- 2014-11-18 PL PL14806452T patent/PL3071954T3/pl unknown
- 2014-11-18 WO PCT/IB2014/066129 patent/WO2015075629A1/en active Application Filing
- 2014-11-18 ES ES14806452T patent/ES2804087T3/es active Active
- 2014-11-18 MX MX2016005767A patent/MX368418B/es active IP Right Grant
- 2014-11-18 HU HUE14806452A patent/HUE051229T2/hu unknown
- 2014-11-18 US US15/027,179 patent/US9714899B2/en active Active
- 2014-11-18 KR KR1020167012835A patent/KR102209231B1/ko active IP Right Grant
- 2014-11-18 JP JP2016526276A patent/JP6619736B2/ja not_active Expired - Fee Related
- 2014-11-18 MY MYPI2016701687A patent/MY185715A/en unknown
- 2014-11-18 BR BR112016007745-8A patent/BR112016007745B1/pt active IP Right Grant
- 2014-11-18 SI SI201431592T patent/SI3071954T1/sl unknown
-
2016
- 2016-05-01 SA SA516371040A patent/SA516371040B1/ar unknown
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2016538536A5 (ja) | ||
Gendre et al. | Development of a process analytical technology (PAT) for in-line monitoring of film thickness and mass of coating materials during a pan coating operation | |
JP2014518390A5 (ja) | ||
RU2016116157A (ru) | Способ мониторинга контролируемого параметра реакции полимеризации и соответствующее устройство для осуществления указанного способа | |
CN103983604B (zh) | 一种基于太赫兹时域光谱技术的检测系统 | |
CN105548026A (zh) | 一种快速的郁金药材质控检测方法 | |
JP2010002327A5 (ja) | ||
CN104833654A (zh) | 利用近红外光谱法快速检测野菊花药材的方法及应用 | |
JP2011082286A5 (ja) | ||
CN104977271A (zh) | 一种红花醇沉过程有效成分近红外在线检测方法 | |
Nishii et al. | Selection of the NIR region for a regression model of the ethanol concentration in fermentation process by an online NIR and mid-IR dual-region spectrometer and 2D heterospectral correlation spectroscopy | |
CN106092945A (zh) | 食品接触材料中双酚a迁移量的在线检测方法及检测系统 | |
Szymaszek et al. | Review of quantitative and qualitative methods for monitoring photopolymerization reactions | |
CN104316490A (zh) | 利用近红外光谱法快速检测鸡血藤药材的方法及应用 | |
JP6488659B2 (ja) | コンクリート劣化診断装置 | |
CN104390930A (zh) | 一种基于漫反射红外光谱技术的木材含水率测量方法 | |
JP2012154718A (ja) | 分光分析方法および分光分析用サンプリングユニット | |
CN104819955A (zh) | 基于粒子群最小二乘支持向量机算法检测川芎方法及应用 | |
Sánchez-Paternina et al. | When “homogeneity” is expected—Theory of Sampling in pharmaceutical manufacturing | |
JP2017058217A5 (ja) | ||
RU2013128966A (ru) | Устройство и способ определения активности болезни | |
CN108051398A (zh) | 用于全透射光在农产品中探测距离的测量系统 | |
Gua et al. | Evaluation of wheat seeds by terahertz imaging | |
CN106841097B (zh) | 一种判别光学异构体旋光性的近红外光谱分析方法 | |
Lin et al. | Investigation of pharmaceutical film coating process with terahertz sensing, optical coherence tomography and numerical modelling |