Corsica调研与工作进展 .pptVIP

Corsica调研与工作进展 报告人：邢哲 导师：肖炳甲 研究员 报告内容 1、Corsica简介 2、BASIS初步 3、Corsica构建平衡 4、Corsica应用调研 5、工作进展和计划 一、Corsica简介 Corsica是由美国劳伦斯利弗莫尔国家实验室（LLNL）开发的一个综合的托卡马克放电模拟程序。开发的初衷是通过Corsica实现对托卡马克中所有物理过程的综合模拟，分析放电试验中的物理现象，优化放电控制，也可以用于优化托卡马克设计，比如ITER Corsica代码耦合二维平衡计算和一维输运模拟，可以使用自由边界平衡解或者固定边界平衡解模拟放电平衡演化。目前的Corsica III 加入了二维的边界输运和湍流输运模型 Corsica包含的程序包可以数值计算等离子体平衡、模拟等离子体中的输运（如芯部输运、中心-边界输运、湍流输运等）、分析MHD稳定性、辅助加热等等。 Corsica的程序包还包括常用数学计算、多项式拟合、Bessel函数、快速傅里叶变换、图形绘制、BASIS开发环境等。 Corsica中还预留与Mdsplus和外部控制器的接口 二、BASIS初步 Corsica的一个优点是代码的适应性好，用户可以通过交互界面自由的定义变量，编译程序 使用体验类似Matlab ；甚至可以根据需要向Corsica加入编译好的程序包，方便国际合作。 这一优势源于Corsica基于BASIS开发环境，BASIS提供了交互式语言解释器，内部包括大量嵌入的函数、错误恢复、宏、存取变量、标准化I/O等。 BASIS语言是一种解释性语言，语法类似Fortran90语言。用户可以自定义变量和函数，之后保存为一个脚本文件.bas 类似Matlab中.m文件 举例：test.bas # caltrans test.bas read tokamak.bas # Variables used by weqdsk... shotName probname;shotTime 1 # Create a nominal direct-solve eq from scratch... ds trim probname //.inp # Converge to q 0 1.1, q a 2.9, betaN 1... probid DCON //trim probname // direct-solve nctot 3 vo [ctroy,qsrf 1 ,qsrf msrf ] vo0 [1,1.1,2.9] vi [betaj,alfa 0 ,plcm] x0 [betaj,alfa 0 ,plcm] ihy 20; run # Save the direct-solve eq... saveq trim probname //.sav # Make some plots... layout;profiles # Write direct-solve EQDSK files... weqdsk ag # Standard a- and g-files quit Corsica中读写数据 1 使用标准I/O设备读写文本文件 读入data1.dat： 写入data2.dat： integer n 100 character*1 tab char 9 real big_array 2,n io basopen data2.dat,w integer io basopen data1.dat,r io big_array do $i 1,n call basclose io io format time $i ,10,3,1 \ real time n big_array 1, tab format current $i ,16,8,2 real current n big_array 2, enddo call basclose io 使用PFB设备读写二进制文件.pfb: 创建some_data.pfb： real time 0 , current 0 time : shotTime # shotTime is a caltrans variable current : placur # code variable placur holds Itor create some_data.pfb write time,current close 之后可以使用： open some_data.pfb ls #查看some_data.pfb中的内容 三、Corsica构建平衡 Corsica构建平衡的启动模式有三种： 1、从零开始构建平衡 2、从EQDSK文件启