JP3704195B2 - 超電導体の電磁界シミュレーション方法 - Google Patents
超電導体の電磁界シミュレーション方法 Download PDFInfo
- Publication number
- JP3704195B2 JP3704195B2 JP08653496A JP8653496A JP3704195B2 JP 3704195 B2 JP3704195 B2 JP 3704195B2 JP 08653496 A JP08653496 A JP 08653496A JP 8653496 A JP8653496 A JP 8653496A JP 3704195 B2 JP3704195 B2 JP 3704195B2
- Authority
- JP
- Japan
- Prior art keywords
- current density
- minute region
- region
- minute
- superconducting material
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Measuring Magnetic Variables (AREA)
- Superconductor Devices And Manufacturing Methods Thereof (AREA)
Description
【発明の属する技術分野】
本発明は、例えばMRI、リニアモーターカー等の超電導機器において、磁気シールドおよび超電導マグネットなどに使用される超電導体内外の電磁界シミュレーションをおこなう方法に関するものである。
【0002】
【従来の技術】
近年では、MRI、磁気浮上搬送装置など超電導を利用した様々な装置が考案、実用化されている。超電導の応用分野でも特に超電導線材をコイル状に巻き、通電する超電導マグネットとしての利用がほとんどであった。しかし、最近では超電導体をシート状、またはバルク状に加工し、磁界発生源から磁気的に隔絶する磁気シールドとしての利用や、着磁して超電導マグネットとしての利用といった超電導体の磁気的な側面を利用した応用、実用化が広がりつつある。これらの応用では、そのシールド能力や磁界発生能力は、超電導特性、特に超電導体に流れうる最大の電流密度(臨界電流密度;Jc)とその形状によって様々に変化する。それゆえ、装置設計の段階から超電導の特性と形状を考慮しておく必要がある。しかし、実験でのシミュレーションは経済性及び時間の観点から頻繁に試行実験を行うことは難しい。そこで、実験にかわり計算による超電導体の電磁界シミュレーションによって、必要とする特性や形状をあらかじめ予測し、絞った開発を行うことのできる手法が求められていた。
【0003】
シールド体や着磁マグネットとして利用する場合、外部の磁場の大きさや方向の変化に対して、超電導体内部にはその磁場の大きさに応じて電流が自発的に流れ、この電流密度値が超電導体に流すことのできる最大の電流密度値(臨界電流密度)が流れる。これを臨界状態と呼ぶが、超電導シミュレーション手法での目標は、超電導体内での電流分布を臨界状態に近づけるように解くことである。つまり、臨界状態とは超電導体内の磁界変動を出来るだけ抑えるべく、自発的に電流を臨界電流密度で流している状態である。これを模擬する手法は、大きく分けて2種類ある。一つは、直接求める変数を磁気ベクトルポテンシャル(A)とする方法(1)であり、もう一つは電流ベクトルポテンシャル(T)とする方法(2)である。以下にマクスウェル方程式から導かれる基礎方程式をあげておく。
【数1】
【0004】
ここでJ0 、H0 はそれぞれ外部から与えられる電流密度(強制電流密度)、磁界の強さ(強制磁界)を表し、σ、μはそれぞれ電気伝導度、透磁率を表し、tは時間を表している。この式から分かるように、前者は電流密度を外部から条件が与えられる項(強制項)として代入することが可能であるため、単に臨界電流密度の値を代入するだけで、臨界状態を模擬できる反面、軸対称でない三次元問題を解く場合にそのままでは電流の保存則が成り立たない欠点があり、三次元問題を解く場合の障害となっている。というのも、本来求める解が磁場に関する変数、磁気ベクトルポテンシャル(A)であり、電流密度ではないので、電流密度J0 は、そのAの外部条件を与える単なるパラメータとして扱われ、その連続性は保証されないのである。
【0005】
一方、本発明もこの範疇に含められる後者の方法は、計算する行列が大きくなり計算時間の点で不利ではあるが、求める解が電流密度に関する変数、電流ベクトルポテンシャル(T)であるので、解を求める際に連続性は考慮されて解かれ、三次元問題でも電流の保存則は成立するため任意形状で超電導体を扱うことが比較的容易である。しかし、電流密度を強制項として代入することができないため、個々の微小領域に流れる電流密度を該領域の電気伝導度をもって間接的に操作する必要がある。
【0006】
この方法での従来の手法を図2に基づき以下に述べる。全体として超電導状態の理想状態(電流密度いくらでも流れる状態)を出発点とし、徐々に臨界電流密度の制限を加えてるシミュレーション方法をとる。それゆえ初期値として理想状態に近くなるよう電気伝導度(σ)を無限大としたいが、計算機上では扱うことができないので、計算機で扱える最大の値に近い値を代入してVを解くことからシミュレーションを開始する。このVから各電流密度(J)分布をもとめ、Jcを超えた微小領域の電気伝導度(σ)を減少させて、再度計算する。以後同様に電流密度(J)分布をもとめ、Jcを超えた微小領域の電気伝導度(σ)を減少させて、Jcを超えた微小領域がなくなるまで繰り返し計算する。
【0007】
この場合、問題点が2つある。1つは電気伝導度(σ)の初期値によって最終結果が異なることと、もう1つは、各ループの際にJcを超えた微小領域しか扱わないため、一度でもJc以下に低下した微小領域は、本来Jcとするべきであるにもかかわらず、取り残されてしまい、微小領域の分け方で解が異なってしまう点である。前者はσの初期値が計算機の扱える最大値と関係するため、個々の計算機の能力に左右され、収束結果も異なる可能性がある。また、後者では一般に微小領域が全体に小さく取られていれば解の精度が高まるはずが、それとは逆の現象が起こる可能性を有している。従来の方法ではこの微小領域が小さくなればなるほど、図4(B)に示すがごとく微小領域の平均電流密度として扱うので、そこに流れるJがJcより極端に大きくなり、電気伝導度(σ)の更新幅が大きくなり、更新したσを用い、計算したJはJcよりもかなり低い値になり、計算の対象から外され、本来Jcの値を取る領域が大きく外れたJのままとなり、結果的に微小領域が小さい方が計算機精度が悪くなる欠点をもっていた。
【0008】
以上のように問題点としては第1に初期値によって最終結果が異なる点、第2に微小領域の分け方で解が異なってしまう点の2つがあった。
【0009】
【発明が解決しようとする課題】
本発明は、従来方法の電気伝導度の初期値の任意性をなくし、微小領域分割の大きさによる変動が小さく、計算精度を高め、かつ計算機資源を有効に利用することができる電磁界シミュレーション方法を提供することにある。
【0010】
【課題を解決するための手段】
本発明は、印加磁場に対して超電導材料の電流密度と、超電導材料内部およびその周囲の磁界分布を予測する目的で、超電導材料を微小領域に分割し、該微小領域に流れる電流密度Jを後述の(4)式で記述するポテンシャルVで表現し、該超電導材料の電流密度が臨界電流密度を超えないよう繰り返し計算する電磁界シミュレーションにおいて、繰り返し計算の1回目に電気抵抗率を含まない支配方程式、具体的には後述の(5)式で該ポテンシャルを算出して各微小領域の電流密度を計算する第1手法と、求めた各微小領域の電流密度と該超電導材料の臨界電流密度値とを比較し、臨界電流密度値を超えた微小領域を計算機の記憶領域に記憶し、以降継続して保持する第2手法と、記憶した微小領域については、該微小領域が電気抵抗率ゼロであるとき、該微小領域の電場と電流密度から該微小領域の電気抵抗率を求めておき、臨界電流密度値Jcを超えた電流密度Jの微小領域はその電気抵抗率ρをρ=ρ・J/Jcにより増加し、臨界電流密度値Jcより低い電流密度Jの微小領域はその電気抵抗率ρをρ=ρ・J/Jcにより減少させ、電気抵抗率を含む支配方程式、具体的には後述の(4)式で該ポテンシャルを算出し、再度各微小領域の電流密度を計算する第3手法と、求めた各微小領域の電流密度と該超電導材料の臨界電流密度値とを比較し、該微小領域の電流密度が臨界電流密度を超えた微小領域がある場合は第2手法以降の処理を繰り返し、ない場合は収束とみなし計算を終了する第4手法を特徴とする電磁界シミュレーション方法である。
【0011】
本発明は前述の第1手法から第4手法までを全て満たすことを特徴とする電磁界シミュレーション方法である。
【0012】
本発明で述べる超電導材料とは、超電導体のみで構成された材料だけでなく、超電導線材や超電導多層材料といった金属材料を含むものでも良い。該金属材料には銅、銅合金材料や、銀、銀合金などが代表例としてあげられる。
【0013】
本発明で述べる電流を記述するポテンシャルは、次の(3)および(4)式が考えられるが、本発明では(4)式に限定することとする。(4)式は(3)式に電流バーJが2方向の成分しか流れない条件(例えば、超電導体のある面内で流れ、厚み方向への電流を無視できる場合を課した場合に得られる式で、上述の制約を受けるものの、(4)式の変数Vはスカラー量、つまり変数が1つとして扱えるため簡便であり、計算速度や計算機資源が有効に使える。ここでのバーnは微小領域の表面法線単位ベクトルである。金属層と超電導層を何層か重ねた材料や酸化物超電導体のように異方性の高く、ある結晶方位については電流が流れにくい性質の材料をシミュレートする場合には(4)式の変数Vが有効である。
【数2】
【0014】
本発明で述べる電気抵抗率を含む支配方程式は(5)式の有限要素回路法と呼ばれる支配方程式である。
【数3】
【0015】
(5)式の各変数は表1に示す標記に従う。但しtは時間である。
【0016】
【表1】
【0017】
本発明で述べた電気抵抗率を含まない支配方程式とは、(7)式を指す。これは電気抵抗率が零としてマクスウェル方程式から定式化した式である。
【数4】
【0018】
本発明で述べた臨界電流密度値を超えた微小領域を計算機の記憶領域に記憶し、保持することとは、臨界電流密度値を超えた微小領域がどれなのかを計算機のメモリ上で記憶させ、一度超えた微小領域は第4手法で収束と判定されるまで保持することを意味する。具体的には微小領域に分割した数(要素数;Ne)の配列を用意し、初期設定としてある値(例えば0)を全ての配列に代入し、臨界電流密度値を超えた微小領域の通し番号に対応する配列に初期設定と違う別の値(例えば1)を代入し、一度超えた微小領域は第4手法で収束と判定されるまで保持することを指す。
【0019】
本発明での微小領域の電場Escは、印加磁界に対応する磁気ベクトルポテンシャルA0 と、各微小領域の電流密度による反磁界に対応する礎気ベクトルポテンシャルAexで表される。具体的にEscは(15)式および(16)式で表される。
【数5】
ここでμ0 、π、hはそれぞれ透磁率、円周率、超電導材料の厚みを指し、j(r′)はr′の位置での電流密度を指す。また、積分は超電導材料の微小領域の表面積分を意味する。
【0020】
本発明での該微小領域の電場Escと電流密度Jから該微小領域の電気抵抗率ρを求める方法は、オームの法則に従い、(17)式より求めることを指す。
【数6】
【0021】
本発明での臨界電流密度値を超えた電流密度の微小領域はその電気抵抗率を増加し、臨界電流密度値より低い電流密度の微小領域はその電気抵抗率を減少させる方法とは、ρ=ρ・J/Jcにより、電流密度Jが臨界電流密度Jcに収束するように増減させて、繰り返し計算することを指す。
【0022】
【発明の実施の形態】
本発明を図1に基づき説明する。ここで第1手法から第4手法までは図1における(a)から(d)の計算ステップに対応する。前述のように従来の手法での問題点は初期値の任意性と微小領域の細分化に伴う解の精度悪化の2点である。これらは同形状で同特性をもつ超電導材料の計算結果が計算上の初期設定値や領域区分の方法によって異なることを意味する。本発明では、電気抵抗率をゼロとみなした時間に依存しない支配方程式で簡単に計算することで初期値を設定することなく、初期値の任意性の問題解決を図っている。これが第1手法をおこなう理由である。
【0023】
また、領域の区分の方法で異なる最大の原因は、微小領域の電流密度がJc近傍の値になるべきところ、図4(B)のようにJc近傍の値とかけ離れた値に留まってしまう点である。このため第2手法において第3手法で操作する微小領域をチェックし、それを以降の繰り返し計算の際にも保持し、一度でも臨界電流密度を超えたことのある微小領域は第3手法において電気抵抗率を増減させて、電流密度Jが臨界電流密度Jc近傍の値に収まるよう計算している。これを第4手法において、繰り返し計算することにより、微小領域の細分化に伴う解の精度悪化を避けることができ、本来の領域の区分の微細化を伴った計算精度向上が図れる。これが第2手法、第3手法および第4手法をおこなう理由である。
【0024】
また、第3手法では微小領域で初めて臨界電流密度を超えた場合は、該領域の電気抵抗率を該領域の電場と電流密度から求めることで、適当な電気抵抗値を得ることができ、計算収束を速め、計算機資源の有効利用ができる。
このように計算精度を高め、計算機資源を有効に利用することができる。さらに、この解析結果を用いて超電導材料を使用した装置設計などに対して経済的、時間的にも削減できる利点を持つ。
【0025】
【実施例】
円筒加工したCu/Nb/Nb−Tiの30層の多層材でNb−Ti層の1層の厚みが平均18μmの材料について円筒軸に垂直に磁界を印加した場合(図7) について有限要素回路法をもちいて解析を行った。解析に利用した微小領域の区分は半径方向に1分割、同径方向に10分割、長さ方向に8分割した。微小領域の形状は三角形であり、微小領域要素数は160である。多層材のサイズは内径20mm、長さ20mm、厚さ1mmである。超電導体の臨界電流密度を通電法で測定し、その結果を材料特性値として用いた。計算結果と実験値を比較した図が図8である。このように実験値と定性的に10%の相対誤差で一致する。また本解析の収束までのループ回数は平均13回であり、従来の方法では平均24回であった。約半分で計算ことができた。加えて、微小領域分割数を6倍に増やしても、5%の精度で同じ値に収束した。このように計算機資源の有効利用が図れ、超電導体を含んだ材料の、形状と特性値を盛り込んでの電磁気的な応答能力を精度よく予測することが可能となった。
【0026】
【発明の効果】
本発明を用いることにより、収束までのループ回数を低減し、かつ初期値を与える必要がないため計算解の初期値による任意性がなく一意の解が得られ、精度を高める上で必要な微小領域の細分化をおこなっても、収束結果の変わらない安定した解が得られた。これにより計算機資源の有効利用が図れ、超電導体を含んだ材料の、形状と特性値を盛り込んでの電磁気的な応答能力を精度よく予測できる効果がある。
【図面の簡単な説明】
【図1】本発明の方法の流れ図の例。
【図2】従来方法の流れ図。
【図3】従来方法で超電導板に垂直な磁界が印加された場合に初期値σにより得られる連続した電流分布図。
【図4】従来方法での微小領域区分が大きい場合の電流分布と繰り返し計算による電流分布の変化(A)。微小領域区分が大きい場合の電流分布と繰り返し計算による電流分布の変化(B)を示す図。
【図5】本発明の方法で超電導板に垂直な磁界が印加された場合に初期値を必要としない第1手法により得られる連続した電流分布図。
【図6】本発明方法での微小領域区分が小さい場合の電流分布と繰り返し計算による電流分布の変化図。
【図7】円筒加工したCu/Nb/Nb−Tiの多層材の形状、配置と磁界測定点。
【図8】円筒加工したCu/Nb/Nb−Tiの多層材の計算結果(予測値)と実験結果(実験値)との比較を示す図。
【符号の説明】
(a) 請求項の第1手法に対応する流れ図の箇所。
(b) 請求項の第2手法に対応する流れ図の箇所。
(c) 請求項の第3手法に対応する流れ図の箇所。
(d) 請求項の第4手法に対応する流れ図の箇所。
(21) 微小領域。
(22) 印加磁界の方向。
(23) 超電導材料。
(24) 初期値の電気電導度σを用いて得た電流密度分布。
(25) 電流を変化させる方向と大きさを概念的に示す矢印。
(26) 超電導材料の臨界電流密度を示す点線
(27) 電気電導度σを減少させる微小領域の範囲
(28) 極端なJc以下の微小領域が発生した箇所
(29) 電気抵抗率を含まない支配方程式から得た電流密度分布。
(30) 電気抵抗率を増加させる微小領域の範囲。
Claims (1)
- 印加磁場に対して超電導材料の電流密度と、超電導材料内部およびその周囲の磁界分布を予測する目的で、超電導材料を微小領域に分割し、該微小領域に流れる電流密度Jを下記(4)式で記述するポテンシャルVで表現し、該超電導材料の電流密度が臨界電流密度を超えないよう繰り返し計算する電磁界シミュレーションにおいて、繰り返し計算の1回目に電気抵抗率を含まない下記(7)式で表される支配方程式で該ポテンシャルを算出して各微小領域の電流密度を計算する第1手法と、求めた各微小領域の電流密度と該超電導材料の臨界電流密度値とを比較し、臨界電流密度値を超えた微小領域を計算機の記憶領域に記憶し、以降継続して保持する第2手法と、記憶した微小領域については、該微小領域が電気抵抗率ゼロであるとき、該微小領域の電場と電流密度から該微小領域の電気抵抗率を求めておき、臨界電流密度値Jcを超えた電流密度Jの微小領域はその電気抵抗率ρをρ=ρ・J/Jcにより増加し、臨界電流密度値Jcより低い電流密度Jの微小領域はその電気抵抗率ρをρ=ρ・J/Jcにより減少させ、電気抵抗率を含む下記(5)式で表される支配方程式で該ポテンシャルを算出し、再度各微小領域の電流密度を計算する第3手法と、求めた各微小領域の電流密度と該超電導材料の臨界電流密度値とを比較し、該微小領域の電流密度が臨界電流密度を超えた微小領域がある場合は第2手法以降の処理を繰り返し、ない場合は収束とみなし計算を終了する第4手法を特徴とする超電導体の電磁界シミュレーション方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP08653496A JP3704195B2 (ja) | 1996-04-09 | 1996-04-09 | 超電導体の電磁界シミュレーション方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP08653496A JP3704195B2 (ja) | 1996-04-09 | 1996-04-09 | 超電導体の電磁界シミュレーション方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09283319A JPH09283319A (ja) | 1997-10-31 |
JP3704195B2 true JP3704195B2 (ja) | 2005-10-05 |
Family
ID=13889676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP08653496A Expired - Lifetime JP3704195B2 (ja) | 1996-04-09 | 1996-04-09 | 超電導体の電磁界シミュレーション方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3704195B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101768211B1 (ko) * | 2016-11-01 | 2017-08-17 | 한국기초과학지원연구원 | 초전도코일 상변화 시 온도예측 시뮬레이션 장치 및 그 방법 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11048846B2 (en) * | 2018-06-12 | 2021-06-29 | International Business Machines Corporation | Surface participation analysis of superconducting qubits with the boundary element method |
-
1996
- 1996-04-09 JP JP08653496A patent/JP3704195B2/ja not_active Expired - Lifetime
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101768211B1 (ko) * | 2016-11-01 | 2017-08-17 | 한국기초과학지원연구원 | 초전도코일 상변화 시 온도예측 시뮬레이션 장치 및 그 방법 |
Also Published As
Publication number | Publication date |
---|---|
JPH09283319A (ja) | 1997-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lee et al. | Modeling eddy currents with boundary conditions by using Coulomb's law and the method of images | |
Lee et al. | Distributed multipole model for design of permanent-magnet-based actuators | |
Adeyeye et al. | Size dependence of the magnetoresistance in submicron FeNi wires | |
White et al. | Evidence for interaction effects in the low-temperature resistance rise in ultrathin metallic wires | |
KR910001860B1 (ko) | 핵자기 공명 영상화용 횡 그레디언트 자계코일구성 및 그 제조방법 | |
Mykola et al. | A–V formulation for numerical modelling of superconductor magnetization in true 3D geometry | |
Adeyeye et al. | Magnetostatic interactions and magnetization reversal in ferromagnetic wires | |
DeGeorge et al. | Multiphase resistivity model for magnetic nanocomposites developed for high frequency, high power transformation | |
JP3704195B2 (ja) | 超電導体の電磁界シミュレーション方法 | |
Gilchrist et al. | Screening effect of Ohmic and superconducting planar thin films | |
Webb et al. | The low-frequency performance of H-phi and T-Omega methods using edge elements for 3D eddy current problems | |
Ying et al. | Theory of the Diamagnetic Interactions among the Conduction Electrons in Metals | |
Webb et al. | Adaptive improvement of magnetic fields using hierarchal tetrahedral finite elements | |
Long et al. | A new modeling technique for simulating 3-D arbitrary conductor-magnet structures for RFIC applications | |
Akhmeteli et al. | Superconducting and resistive tilted coil magnets for generation of high and uniform transverse magnetic field | |
Li et al. | Numerical modelling of soldered superconducting REBCO stacks of tapes suggests strong reduction in cross-field demagnetization | |
Machado | Axisymmetric eddy current problems using a special hybrid FE/BE method | |
JP2941516B2 (ja) | 直流重畳の鉄損計算装置 | |
Chen et al. | Analytical method for two dimensional current crowding effect in magnetic tunnel junctions | |
Kurzawa et al. | Analysis of Eddy current system using equivalent multi-branch Foster Circuit and Edge Element Method | |
JP2984108B2 (ja) | 直流重畳のインダクタンス計算装置 | |
Rodger et al. | Calculation of transient 3D eddy currents in nonlinear media verification using a rotational test rig | |
Longhitano | Electromagnetic devices modeling with material constitutive laws adapted to high temperature/frequency | |
Lavis et al. | The effect of sub-lattice order in binary alloys with one magnetic component. II | |
Tang et al. | Computation of 3-D open boundary eddy-current fields in large transformers using T,/spl psi/-/spl phi//sub m/method coupled with GTM |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050208 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050215 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050413 |
|
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: 20050719 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050722 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080729 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090729 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090729 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100729 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110729 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120729 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130729 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130729 Year of fee payment: 8 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130729 Year of fee payment: 8 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130729 Year of fee payment: 8 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130729 Year of fee payment: 8 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |