Skip to content

Commit

Permalink
Add new process in UnsteadyDiffusion to compute diffusion from a Magn…
Browse files Browse the repository at this point in the history 
…eticField function or from a file B.pts

Eventually add a RHS to the system to handle source

Two testcases :
  1. anisotropic source diffusion in circular B
  2. MASTU shot #45470, with list of points to compare to E2D
  • Loading branch information
Romain Futtersack committed Jun 22, 2023
1 parent d9d3584 commit 20e4d3e
Show file tree
Hide file tree
Showing 12 changed files with 16,978 additions and 0 deletions.
5,497 changes: 5,497 additions & 0 deletions example/MASTU/B.pts
View file

Large diffs are not rendered by default.

5,497 changes: 5,497 additions & 0 deletions example/MASTU/BFI.pts
View file

Large diffs are not rendered by default.

94 changes: 94 additions & 0 deletions example/MASTU/FirstOpen.pts
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,94 @@
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<POINTS DIM="2" FIELDS="">
7691.19097835 -17180.65772181
7693.49517877 -17153.76156737
7695.85075229 -17112.18905152
7698.82105168 -17056.11881256
7700.33363929 -16985.65126154
7698.92788438 -16900.87513383
7693.20678885 -16802.09166263
7681.67936513 -16689.60208112
7663.78868209 -16563.68193594
7637.29008541 -16424.86480683
7601.26469709 -16273.57685702
7555.70259278 -16110.01657350
7500.77190265 -15934.27152407
7437.28934225 -15746.25647606
7366.69043955 -15545.73908295
7290.94546980 -15332.37841415
7211.91262147 -15105.94796097
7131.91244096 -14866.10679283
7052.27771003 -14612.81170962
6973.23026707 -14346.34643097
6894.19625117 -14067.12077413
6813.66599365 -13775.72319465
6730.26759915 -13472.63589914
6642.14696569 -13158.45785129
6547.72261902 -12833.75182009
6446.14865427 -12498.99884458
6339.12271900 -12154.15922737
6232.50017934 -11799.53069894
6148.00835505 -11443.02821983
6130.56777045 -11162.12826687
6192.12792598 -10994.72083849
6418.14505920 -10806.03794015
6860.36531164 -10550.26994258
7441.75588327 -10271.51131395
8119.76583415 -9943.48004191
8864.52929423 -9535.50871925
9641.30870548 -9015.09043139
10412.67409223 -8360.68992395
11141.65608541 -7565.23405470
11784.90687820 -6628.00943808
12301.12483112 -5560.13134541
12692.61939864 -4395.28214287
13017.84032844 -3175.75502920
13317.44592866 -1926.28644529
13526.59166790 -647.70091485
13526.59166790 647.70215539
13317.44476109 1926.28775881
13017.83916087 3175.75619677
12692.61939864 4395.28360233
12301.12483112 5560.13251298
11784.90687820 6628.01060565
11141.65608541 7565.23522227
10412.67292466 8360.69109152
9641.30753791 9015.09159897
8864.52812666 9535.50988682
8119.76408280 9943.48120948
7441.75413191 10271.51248152
6860.36414406 10550.27111015
6418.14330784 10806.03910772
6192.12734220 10994.72200606
6130.56718666 11162.12826687
6148.00777127 11443.02821983
6232.50017934 11799.53186651
6339.12330278 12154.16039494
6446.14865427 12499.00001215
6547.72320280 12833.75415523
6642.14754948 13158.46018643
6730.26818293 13472.63940185
6813.66774500 13775.72786493
6894.19741874 14067.12544441
6973.23201842 14346.35343639
7052.28004517 14612.82105018
7131.91594368 14866.11730096
7211.91729175 15105.96080425
7290.95130765 15332.39476013
7366.69744498 15545.75893165
7437.29693146 15746.27865990
7500.78007564 15934.29721062
7555.71018199 16110.04459520
7601.27170251 16273.60721386
7637.29650705 16424.89749881
7663.79451994 16563.71462792
7681.68403541 16689.63710823
7693.20912399 16802.12668974
7698.92905195 16900.91016095
7700.33305550 16985.68628866
7698.81871654 17056.15500724
7695.84900093 17112.22524620
7693.49693012 17153.79776205
</POINTS>
</NEKTAR>
43 changes: 43 additions & 0 deletions example/MASTU/LowerOuter.pts
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,43 @@
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<POINTS DIM="2" FIELDS="">
7244.38976298 16700.59826025
7297.61463260 16751.36889663
7346.85808674 16796.90181168
7391.34077211 16837.08958889
7430.89280717 16873.85638449
7466.60120010 16907.52911821
7499.11979036 16937.64426457
7528.49294566 16963.71494666
7554.58289266 16986.58298281
7577.42057198 17007.10420223
7597.15134616 17025.09062605
7614.14942154 17040.73256828
7628.79542611 17054.21800775
7641.25865760 17065.57379861
7651.73818678 17074.99142240
7660.44242499 17082.65885799
7667.60722147 17088.78626807
7673.46258761 17093.51259352
7678.22218883 17096.96276440
7682.11253379 17098.94296400
7686.36424183 17102.33242121
7695.85075229 17112.18905152
7714.09462541 17130.13461037
7740.32701601 17155.78613464
7773.78666812 17188.97782932
7814.04158063 17229.60811575
7861.50565731 17276.34246086
7915.86307033 17328.46514390
7975.28015018 17386.92656729
8040.77618593 17450.85572340
8112.07563094 17519.68517298
8188.14868801 17594.26257222
8270.43438780 17674.36258000
8359.32912051 17760.94844201
8455.07574080 17854.69968429
8559.38882597 17956.37989795
8672.45635489 18067.18350763
8796.98358844 18189.71302802
</POINTS>
</NEKTAR>
94 changes: 94 additions & 0 deletions example/MASTU/ambipol.xml
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,94 @@
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<COLLECTIONS DEFAULT="auto" />
<EXPANSIONS>
<E COMPOSITE="C[10]" NUMMODES="4" TYPE="MODIFIED" FIELDS="T" />
</EXPANSIONS>

<CONDITIONS>
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyDiffusion" />
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder3" />
</SOLVERINFO>

<PARAMETERS>
<P> NumSteps = 200 </P>
<P> TimeStep = 5e4 </P>
<P> IO_InfoSteps = 10 </P>
<P> IO_CheckSteps = 10 </P>
<P> epsilon = 1.0 </P>
<P> k_par = 1e4 </P>
<P> k_perp = 1e0 </P>
<P> TCore = 300.0 </P>
</PARAMETERS>

<VARIABLES>
<V ID="0"> T </V>
</VARIABLES>

<BOUNDARYREGIONS>
<B ID="0"> C[1] </B>
<B ID="1"> C[2] </B>
<B ID="2"> C[3] </B>
<B ID="3"> C[4] </B>
<B ID="4"> C[5] </B>
<B ID="5"> C[6] </B>
<B ID="6"> C[7] </B>
<B ID="7"> C[8] </B>
<B ID="8"> C[9] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<!--Core -->
<REGION REF="0">
<D VAR="T" VALUE="TCore" />
</REGION>
<!--Lower outer -->
<REGION REF="1">
<N VAR="T" VALUE="0.0" />
<!--<R VAR="T" VALUE="0.0" PRIMCOEFF="1.0/k_par"/>-->
</REGION>
<!--Upper outer -->
<REGION REF="2">
<N VAR="T" VALUE="0.0" />
<!--<R VAR="T" VALUE="0.0" PRIMCOEFF="1.0/k_par"/>-->
</REGION>
<!--Upper inner -->
<REGION REF="3">
<N VAR="T" VALUE="0.0" />
<!--<R VAR="T" VALUE="0.0" PRIMCOEFF="1.0/k_par"/>-->
</REGION>
<!--Lower inner -->
<REGION REF="4">
<N VAR="T" VALUE="0.0" />
<!--<R VAR="T" VALUE="0.0" PRIMCOEFF="1.0/k_par"/>-->
</REGION>
<!--LFS SOL -->
<REGION REF="5">
<D VAR="T" VALUE="0.0" />
</REGION>
<!--HFS SOL -->
<REGION REF="6">
<D VAR="T" VALUE="0.0" />
</REGION>
<!--Bottom PFR -->
<REGION REF="7">
<D VAR="T" VALUE="0.0" />
</REGION>
<!--Top PFR -->
<REGION REF="8">
<D VAR="T" VALUE="0.0" />
</REGION>
</BOUNDARYCONDITIONS>-->

<FUNCTION NAME="InitialConditions">
<E VAR="T" VALUE="1.0" />
</FUNCTION>

<FUNCTION NAME="MagneticField">
<F VAR="Bx,By,Bz" FILE="B.pts" />
</FUNCTION>

</CONDITIONS>
</NEKTAR>
Loading

0 comments on commit 20e4d3e

Please sign in to comment.