Data for "A fast 1D model of active magnetic regeneration with a compressible working fluid"

The is the data set for the paper [Diamantopoulos, T., Moreno-Ramírez, L. M., Franco, V. and Bjørk, R., A fast 1D model of active magnetic regeneration with a compressible working fluid, International Journal of Refrigeration, 178, 359-366, 2025]. The DOI for the publication is 10.1016/j.ijrefrig.2025.07.005
These data are located on the repository at data.dtu.dk, with DOI:10.11583/DTU.29604440 and URL: https://doi.org/10.11583/DTU.29604440
The file Field_and_flow_100e-3_P_1.mat contains the information to plot the field and flow profiles, Fig. 1 in the paper. The file is in MATLAB file format. The field is plotted as plot(z,H5T(z)), i.e. the distance metric is z and the field is saved as an interpolation function in the variable H5T in tesla. The flow profile is plotted as plot(t_node,mdot), i.e. the time variable is t_node and the flow amplitude is mdot in kg/s.
The file DyCoNi.mat contains the specific heat capacity (Cp_tot) in J kg^-1 K^-1, the specific magnetization (mag) in A m^2 kg^-1 and entropy (S) in J kg^-1 K^-1 as tables as function of the magnetic field (Barr) in tesla and the temperature (Tarr) in kelvin. This information is plotted in Fig. 2 in the paper.
The results are structured in six folders, named after the two methods (IMEX and SST), the pressure (1bar and 5bar) and whether the mass flow rate (mdot) or CFL number (CFL) has been varied. Each folder contains a set of MATLAB files. In general these files are named as follows: dyconiCFL_A_freq100mHz_d_particle5e4_mdotBg_PCDT20Nx120 where A, B and C are the CFL number, the mass flow rate in g/s and the pressure in bar, respectively. Thus the file named dyconiCFL_0.1_freq100mHz_d_particle5e4_mdot10g_P1DT20Nx120 has a CFL number of 0.1, a mass flow rate of 10 g/s and a pressure of 1 bar. The remaining numbers, which are not varied in this study, are the frequency of operation, the particle size and the numerical discretization.
Each file is a MATLAB file that contains the following variables
CFL_limit- The CFL limit for the simulation
COP-The computed COP during the different iterations until convergence. Is stated as zero after the model has converged.
D-The diameter of the AMR device.
dE_cycle-The change in energy in a cycle until convergence. Is stated as zero after the model has converged.
err-The numerical error in a cycle until convergence. Is stated as zero after the model has converged.
frac-Always equal to 1.
freq-The frequency of operation in Hz.
L-The length of the AMR device.
mdot-The flow rate profile over a cycle.
mdot_amp-The amplitude of the flow rate.
Pcool-The cooling power during the different iterations until convergence. Is stated as zero after the model has converged.
pf-The pressure within the AMR as function of position (120) and time in a cycle (85001).
Pheat-The released heat during the different iterations until convergence. Is stated as zero after the model has converged.
porosity-The porosity of the AMR.
rhof-The density of the fluid within the AMR as function of position (120) and time in a cycle (85001).
Tf-The temperature of the fluid within the AMR as function of position (120) and time in a cycle (85001).
time-Time during a single cycle.
time_wc-Computation time of a single cycle until convergence. Is stated as zero after the model has converged.
Ts-The temperature of the solid within the AMR as function of position (120) and time in a cycle (85001).
Uavg-The average flow velocity.
v-The flow velocity of the fluid within the AMR as function of position (120) and time in a cycle (85001).
x-Spatial position along the AMR.
To replicate Fig. 4, plot the last value in the Pcool array that is not zero as function of mdot_amp for the files in the folder IMEX-1bar-mdot_variation and SST-1bar-mdot_variation.
The file 0D_Results_DyCoNi_P1.mat contains the results of the 0D model in MATLAB format. For this file, the relevant variables to plot to replicate Fig. 4 are mdot_amp and Pcool_vec.
To replicate Fig. 5, the CFL_limit and Pcool array are the variables to plot from the IMEX-1bar-CFL_variation and SST-1bar-CFL_variation for 1 bar (or alternatively for 5 bar).
To replicate Fig. 6, the mdot_amp and sum(time_wc) (the total computation time) are the variables to plot from the IMEX-1bar-CFL_variation and SST-1bar-CFL_variation for 1 bar (or alternatively for 5 bar).