39 lines
1.2 KiB
Matlab
39 lines
1.2 KiB
Matlab
function [Physic, Material, Port1Solver, Port2Solver] = configPhysicMaterialFem4()
|
||
%CONFIGPHYSICMATERIALFEM4 Fem4 / MeshData2x 基准物理与材料参数
|
||
% 对应 COMSOL Fem4.mph：脊形波导，传播方向 z，长度 4 um（MeshData2x 实测 z in [0,4um]）
|
||
% 材料与 MainFem2.m 一致：SiO2 + Si + 掺杂 Si(σ=5000)
|
||
|
||
Physic.lam0 = 1.55e-6;
|
||
Physic.c_const = 299792458;
|
||
Physic.k0 = 2*pi/Physic.lam0;
|
||
Physic.omega = 2*pi*Physic.c_const/Physic.lam0;
|
||
|
||
Material.epsilonr = [4; 11.9; 11.9; 11.9];
|
||
Material.mur = [1; 1; 1; 1];
|
||
Material.sigma = [0; 0; 5000; 0];
|
||
Material.epsilon0 = 8.854187817e-12;
|
||
Material.temp = 1/Physic.k0*120*pi;
|
||
|
||
Material.epsilonrF = [0 0 4 4 0 0 0 11.9 11.9 0 ...
|
||
0 0 11.9 11.9 0 0 0 0 0 11.9 ...
|
||
11.9 0 0 0]';
|
||
Material.murF = ones(24, 1);
|
||
Material.sigmaF = [0 0 0 0 0 0 0 0 0 0 0 0 5000 5000 ...
|
||
0 0 0 0 0 0 0 0 0 0]';
|
||
|
||
Physic.Port1 = [3 8 13 20];
|
||
Physic.Port2 = [4 9 14 21];
|
||
Physic.PEC = [1 2 5 24];
|
||
|
||
Port1Solver.targetNeff = 2.6;
|
||
Port1Solver.modeNum = 1;
|
||
Port1Solver.lam0 = Physic.lam0;
|
||
Port1Solver.type = 1;
|
||
|
||
Port2Solver.targetNeff = 2.6;
|
||
Port2Solver.modeNum = 1;
|
||
Port2Solver.lam0 = Physic.lam0;
|
||
Port2Solver.type = 2;
|
||
|
||
end
|