Commit 27b32292 authored by frey_m's avatar frey_m

add AMR test 2

parent 295258c9
"Neighboring bunch AMR test where we track two bunches by three turns. It uses the ML Poisson solver of AMReX."
"Neighboring bunch AMR test where we track two bunches by three turns. It uses the ML Poisson solver of AMReX and one AMR level."
stat "rms_x" last 1E-7 AmrRingCyclotron-bunch-0000.smb # first bunch
stat "rms_y" last 1E-7 AmrRingCyclotron-bunch-0000.smb
stat "rms_z" last 1E-7 AmrRingCyclotron-bunch-0000.smb
......
Option, ECHO=FALSE;
Option, PSDUMPFREQ=4500;
Option, PSDUMPFRAME=BUNCH_MEAN;
Option, SPTDUMPFREQ = 1000000;
Option, ENABLEHDF5=TRUE;
Option, PSDUMPEACHTURN=TRUE;
Option, REPARTFREQ=5;
Option, ECHO=FALSE;
Option, STATDUMPFREQ=1;
Option, CZERO=FALSE;
Option, VERSION=20000;
Option, SCsolveFreq=2;
Option, AMR=TRUE;
Option, AMR_REGRID_FREQ=10;
Option, AMR_YT_DUMP_FREQ=1000000;
Option, MEMORYDUMP=TRUE;
Option, TELL=TRUE;
Title,string="OPAL-cycl: the first turn acceleration in PSI 590MeV Ring";
REAL Edes=.072;
REAL gamma=(Edes+PMASS)/PMASS;
REAL beta=sqrt(1-(1/gamma^2));
REAL gambet=gamma*beta;
REAL P0 = gamma*beta*PMASS;
REAL brho = (PMASS*1.0e9*gambet) / CLIGHT;
//value,{gamma,brho,Edes,beta,gambet};
REAL phi01=139.4281;
REAL phi02=phi01+180.0;
REAL phi04=phi01;
REAL phi05=phi01+180.0;
REAL phi03=phi01+10.0;
REAL volt1st=0.9;
REAL volt3rd=0.9*4.0*0.112;
REAL turns =3;
REAL nstep=360;
REAL frequency=50.650;
REAL frequency3=3.0*frequency;
ring: Cyclotron, TYPE="RING", CYHARMON=6, PHIINIT=0.0,
PRINIT=-0.000174, RINIT=2130.0, SYMMETRY=8.0, RFFREQ=frequency,
FMAPFN="s03av.nar";
rf0: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=35.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi01;
rf1: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=125.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi02;
rf2: RFCavity, VOLT=volt3rd, FMAPFN="Cav3.dat", TYPE="SINGLEGAP",
FREQ=frequency3,RMIN = 1900.0, RMAX = 4500.0, ANGLE=170.0, PDIS = 452.0,
GAPWIDTH = 250.0, PHI0=phi03;
rf3: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=215.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi04;
rf4: RFCavity, VOLT=volt1st, FMAPFN="Cav1.dat", TYPE="SINGLEGAP",
FREQ=frequency, RMIN = 1900.0, RMAX = 4500.0, ANGLE=305.0, PDIS = 416.0,
GAPWIDTH = 220.0, PHI0=phi05;
l1: Line = (ring,rf0,rf1,rf2,rf3,rf4);
Dist1:DISTRIBUTION, TYPE=gauss,
sigmax = 2.0e-03,
sigmapx = 1.0e-7,
corrx = 0.0,
sigmay = 2.0e-03,
sigmapy = 1.0e-7,
corry = 0.0,
sigmat = 2.0e-03,
sigmapt = 3.394e-4,
corrt=0.0;
Fs1:FIELDSOLVER, FSTYPE=ML, MX=32, MY=32, MT=32,
PARFFTX=true, PARFFTY=true, PARFFTT=true,
BCFFTX=dirichlet, BCFFTY=dirichlet, BCFFTT=dirichlet,
BBOXINCR=20, AMR_MAXLEVEL=2, AMR_MAXGRIDX=32, AMR_MAXGRIDY=32,
AMR_MAXGRIDZ=32, AMR_BFX=16, AMR_BFY=16, AMR_BFZ=16,
AMR_REFX=2, AMR_REFY=2, AMR_REFZ=2, AMR_DOMAIN_RATIO={1.0, 0.75, 0.75},
AMR_TAGGING=CHARGE_DENSITY, AMR_DENSITY=1.0e-9;
beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=32*32*32, BCURRENT=2.0E-3, CHARGE=1.0, BFREQ=frequency;
Select, Line=l1;
TRACK,LINE=l1, BEAM=beam1, MAXSTEPS=nstep*turns, STEPSPERTURN=nstep,TIMEINTEGRATOR="RK-4";
run, method = "CYCLOTRON-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1,
MBMODE=FORCE, TURNS=2, MB_BINNING=GAMMA, MB_ETA=0.25;
endtrack;
Stop;
#!/bin/bash
mpirun -np 4 $OPAL_EXE_PATH/opal AmrRingCyclotron.in --info 2 --warn 0 "$@" 2>&1
21
0.0e+00 0.63013 1.08
5.0e-02 0.68626 1.08
1.0e-01 0.73857 1.00
1.5e-01 0.78677 9.20e-1
2.0e-01 0.83060 8.30e-1
2.5e-01 0.86979 7.35e-1
3.0e-01 0.90416 6.36e-1
3.5e-01 0.93349 5.34e-1
4.0e-01 0.95762 4.29e-1
4.5e-01 0.97643 3.21e-1
5.0e-01 0.98981 2.12e-1
5.5e-01 0.99766 1.01e-1
6.0e-01 0.99998 -9.31e-3
6.5e-01 0.99673 -1.20e-1
7.0e-01 0.98793 -2.30e-1
7.5e-01 0.97363 -3.40e-1
8.0e-01 0.95393 -4.47e-1
8.5e-01 0.92891 -5.52e-1
9.0e-01 0.89871 -6.53e-1
9.5e-01 0.86352 -7.51e-1
1.0e+00 0.82354 -7.51e-1
\ No newline at end of file
21
0.0e+00 0.109737 2.91
5.0e-02 0.258532 2.91
1.0e-01 0.401421 2.76
1.5e-01 0.535168 2.55
2.0e-01 0.656550 2.27
2.5e-01 0.762840 1.95
3.0e-01 0.851785 1.58
3.5e-01 0.921153 1.17
4.0e-01 0.969231 7.40e-1
4.5e-01 0.995244 2.92e-1
5.0e-01 0.998514 -1.65e-1
5.5e-01 0.978652 -6.20e-1
6.0e-01 0.936425 -1.05
6.5e-01 0.872880 -1.47
7.0e-01 0.789161 -1.85
7.5e-01 0.687305 -2.19
8.0e-01 0.569821 -2.48
8.5e-01 0.439238 -2.71
9.0e-01 0.298508 -2.88
9.5e-01 0.150978 -2.98
1.0e+00 0.222050 -2.98
\ No newline at end of file
48a73a591194bf8e10c6067a0fe7d3d6 AmrRingCyclotron-bunch-0000.smb
7c6efedc03028394d4e23f3b56c60457 AmrRingCyclotron-bunch-0001.smb
2c767d77866390de5d72bbef4ef4353d AmrRingCyclotron.out
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment