Compare commits

...

69 Commits

Author SHA1 Message Date
cf7ee78e83 USER_PROBLEM: Fix plane selection for magnetic flux integration.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-10-09 17:00:39 -03:00
e615a1a341 USER_PROBLEM: Improve user refinement criterion.
Ensure that large blocks spanning beyond the `yref` parameter are
refined appropriately.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-10-09 16:28:48 -03:00
9dffcba230 Merge branch 'master' into flux-tubes 2024-10-07 10:42:09 -03:00
fd87018ebd USER_PROBLEM: Add user refinement criterion.
The criterion sets only those blocks for refinement that lie at
a distance specified by the parameter 'yref' from the current sheet.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-09-15 11:14:02 -03:00
4429843c95 USER_PROBLEM: Implement a few types of the initial perturbation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-08-17 17:55:22 -03:00
b24ab38389 USER_PROBLEM: Normalize the magnetic flux by the X area.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-08-16 16:37:11 -03:00
ab481646b6 USER_PROBLEM: Adapt the sincos profile for the periodic box.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-08-16 16:30:31 -03:00
d62d29d0ed Merge branch 'master' into flux-tubes 2024-07-30 22:25:40 -03:00
7648240902 Merge branch 'master' into flux-tubes 2024-07-30 16:54:55 -03:00
5083bce795 Merge branch 'master' into flux-tubes 2024-07-19 10:07:20 -03:00
31a5b4975e Merge branch 'master' into flux-tubes 2024-07-10 20:56:32 -03:00
35826e398d Merge branch 'master' into flux-tubes 2024-07-03 22:36:49 -03:00
ed4c83c5f4 Merge branch 'master' into flux-tubes 2024-07-01 20:11:58 -03:00
2271b69255 Merge branch 'master' into flux-tubes 2024-06-03 22:26:20 -03:00
e6d4cc51cd Merge branch 'master' into flux-tubes 2024-06-01 23:41:46 -03:00
6a178661f5 Merge branch 'master' into flux-tubes 2024-05-25 13:21:46 -03:00
58937c4ad7 Merge branch 'master' into flux-tubes 2024-05-24 19:01:41 -03:00
5372208ba5 Merge branch 'master' into flux-tubes 2024-04-20 18:01:44 -03:00
a5684f4e2f PROBLEMS: Update the parameter file for the flux-tubes problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2024-03-07 22:10:44 -03:00
db0a979057 Merge branch 'master' into flux-tubes 2024-03-07 18:21:03 -03:00
0d6f7a47ea Merge branch 'master' into flux-tubes 2023-12-31 17:46:56 -03:00
9e82eb3f6d Merge branch 'master' into flux-tubes 2023-12-28 20:43:52 -03:00
ce859bfdc1 Merge branch 'master' into flux-tubes 2023-12-26 22:22:31 -03:00
04bc74bead Merge branch 'master' into flux-tubes 2023-12-19 16:04:24 -03:00
8253aee8b8 Merge branch 'master' into flux-tubes 2023-12-18 09:35:38 -03:00
712933eecc Merge branch 'master' into flux-tubes 2023-12-14 18:21:11 -03:00
3112b82258 Merge branch 'master' into flux-tubes 2023-09-18 09:25:38 -03:00
a896b7188a USER_PROBLEM: Do not limit B at boundaries in profile 2.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2023-08-29 18:43:44 -03:00
488eb72173 Merge branch 'master' into flux-tubes 2023-08-29 17:45:00 -03:00
f50ccca3a2 Merge branch 'master' into flux-tubes 2023-08-29 17:23:32 -03:00
b624595682 Merge branch 'master' into flux-tubes 2023-08-29 10:53:47 -03:00
bb2d821df6 Merge branch 'master' into flux-tubes 2023-08-14 17:32:01 -03:00
65c41105c6 Merge branch 'master' into flux-tubes 2023-08-05 19:02:40 -03:00
6f0ecb0e42 Merge branch 'master' into flux-tubes 2023-07-27 18:44:45 -03:00
80f1fdde06 USER_PROBLEM: Implement alternative setup for flux tubes.
This setup assume two flux tubes formed by eliptic magnetic field lines.
The two flux tubes enter in contact in the midplane. Moreover, the
magnetic field lines change the polarization across the contact plane.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2023-07-13 14:56:27 -03:00
f7daa3f60d Merge branch 'master' into flux-tubes 2023-06-02 16:26:03 -03:00
b733365066 Merge branch 'master' into flux-tubes 2023-04-28 17:29:40 -03:00
2a112214b3 Merge branch 'master' into flux-tubes 2023-04-03 18:17:16 -03:00
e8e903aadf Merge branch 'master' into flux-tubes 2023-02-28 10:16:03 -03:00
2b79fb73c8 Merge branch 'master' into flux-tubes 2023-02-27 10:07:12 -03:00
22ae707b99 Merge branch 'master' into flux-tubes 2023-02-22 16:10:36 -03:00
0c37c62585 Merge branch 'master' into flux-tubes 2023-02-15 10:32:25 -03:00
1d7e17cd07 Merge branch 'master' into flux-tubes 2023-02-08 12:33:40 -03:00
614f51099b Merge branch 'master' into flux-tubes 2023-02-07 10:18:48 -03:00
f4ef9f3595 Merge branch 'master' into flux-tubes 2022-12-08 13:09:47 -03:00
021a3de384 Merge branch 'master' into flux-tubes 2022-12-06 20:24:29 -03:00
2f4aade054 Merge branch 'master' into flux-tubes 2022-12-04 20:20:31 -03:00
2561c11567 USER_PROBLEM: Estimate the thickness from the resistivity, by default.
The thikness of the current sheet is set to the square root of the
resistivity, if it is not set explictely in the parameter file.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-27 18:23:56 -03:00
5793fe7810 USER_PROBLEM: Calculate Alfvén speed using bamp only.
This is the Alfvén speed along the X direction defining the velocity
scale, time scale, and the Lundquist number.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-27 18:16:24 -03:00
20492048aa USER_PROBLEM: Correct the way the zeta parameter works.
Parameter 'zeta' controls if the thermal pressure (zeta=0.0) or
guide field (zeta=1.0) should contribute to make the total
pressure uniform.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-27 17:06:13 -03:00
ef68c12b44 Merge branch 'master' into flux-tubes 2022-11-27 16:00:15 -03:00
c438dca1fb Merge branch 'master' into flux-tubes 2022-11-25 18:21:57 -03:00
116a49821b Merge branch 'master' into flux-tubes 2022-11-25 16:08:40 -03:00
dd46ee7ee7 Merge branch 'master' into flux-tubes 2022-11-16 15:52:38 -03:00
671b1ef6da Merge branch 'master' into flux-tubes 2022-11-11 17:50:04 -03:00
83d74d627b USER_PROBLEM: Rename magnetic flux file.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-11 17:04:13 -03:00
40213a9699 Merge branch 'master' into flux-tubes 2022-11-07 19:57:29 -03:00
fbb40a4a32 USER_PROBLEM: Fix calculation of boundary resistive terms.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-05 18:51:35 -03:00
b49deb3b71 USER_PROBLEM: Fix conditions for boundary fluxes.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-05 18:48:36 -03:00
0f8aed46a1 USER_PROBLEM: Fix slice selection for the magnetic flux integration.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-04 19:30:13 -03:00
8a66bc5769 USER_PROBLEM: Integrate the magnetic flux at x=0 only.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-04 14:07:35 -03:00
711fc23f03 USER_PROBLEM: Implement reconnection rate measure.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-11-04 12:22:45 -03:00
52d6933efb Merge branch 'master' into flux-tubes 2022-10-24 10:07:36 -03:00
43ba610aef Merge branch 'master' into flux-tubes 2022-10-24 08:56:56 -03:00
9b3fcb3fe1 Merge branch 'master' into flux-tubes 2022-10-15 16:29:37 -03:00
46c9114c8f Merge branch 'master' into flux-tubes 2022-10-06 13:10:40 -03:00
a1fd7e9bc6 PROBLEMS: Add parameters file for magnetic flux tubes problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-10-04 18:20:00 -03:00
c163687a74 USER_PROBLEM: Change default parameters.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-10-04 18:18:32 -03:00
4cb96ca2b3 USER_PROBLEM: Implement flux tube reconnection problem.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
2022-10-04 13:10:38 -03:00
2 changed files with 1036 additions and 2 deletions

64
problems/flux-tubes.in Normal file
View File

@ -0,0 +1,64 @@
# problem name and parameters
#
problem = "flux-tubes"
profile = "sincos"
beta = 2.00d+00
bgui = 0.00d+00
delta = 3.16d-03
# diffusion coefficients
#
viscosity = 1.00d-05
resistivity = 1.00d-05
glm_source_terms = "none"
# random number generator
#
gentype = "same"
# physics
#
equation_system = "mhd"
equation_of_state = "iso"
# methods
#
time_advance = "ssprk3(2)4"
riemann_solver = "roe"
reconstruction = "ocmp5"
# mesh parameters
#
xmin = -5.00d-01
xmax = 5.00d-01
ymin = -5.00d-01
ymax = 5.00d-01
# refinement control
#
ncells = 32
minlev = 1
maxlev = 9
refinement_variables = "vort jabs"
currmin = 5.0d-02
currmax = 2.0d-01
vortmin = 5.0d-02
vortmax = 2.0d-01
# boundary conditions
#
xlbndry = "user"
xubndry = "user"
ylbndry = "user"
yubndry = "user"
# runtime control parameters
#
tmax = 2.0d+00
cfl = 5.0d-01
# data output control
#
snapshot_interval = 1.0d-01
compression_format = "zstd"
compression_level = 19

File diff suppressed because it is too large Load Diff