(DRAFT) feat: introduce 3-velocity grid variable to support relativistic transport #109

cellio00 · 2025-09-25T21:23:59Z

implementation of 3-velocity to be stored in grid state multifab array. currently being read in through initial conditions file.

…erature units

Currently compute_dt runs every iteration. We may need to run compute_dt just once initially during production, since (rho,T,Ye) is static.

… function.

…density for mesh quantities

… antineutrino number densities

…tities to compute the correct time step

…rallel. This function will be used to compute the minimum trace of mesh quantities, and based on that, it will compute the optimum time step

…dary conditions and black holes.

… the minimum potential

…erse mean free paths in the grid center

…des a better representation of the time derivative of N than the previous 1/Opacity*c.

…and black holes in the function that computes the time step.

…e time step.

…put parameter

… it can be used when needed)

…ata generated by EMU

different to avoid repete code

… test

…dat file that can be used as an initial condition for the Emu code.

…to make the process of writing data more efficient. Check points will not be permited for the output file.

…r cell optional. It run single boxes only if the input parameters are defined otherwise it splits particles in all cells

…mber bigger than 5 digits

… different from Fermi–Dirac. This will be useful for local simulations of CFI from NSM global cells, where the equilibrium is not necessarily Fermi–Dirac.

…al_instability_test

one that had the neutrino dir inside the plt file

…ir and just take the plt directories with data

…riables. MSW input test changed in order to initialize vup components as a parameter from the parameter class.

shankar-1729 and others added 30 commits November 4, 2024 14:07

Make IMFP calculation energy dependent

7a53b35

Ensure that rho, temp and Ye are uniformally spaced.

de0bac5

Implement function to calculate maximum of IMFP over the grid (UNTESTED)

7daa907

Move multifab for IMFB inside compute_max_IMFP function, correct temp…

2c67c2e

…erature units

do not input ba, dm and ngrow inside compute_max_IMFP function

795615c

calculate max_IMFP_abs inside function compute_dt

a9d22c4

Currently compute_dt runs every iteration. We may need to run compute_dt just once initially during production, since (rho,T,Ye) is static.

Remove dt calculation from compute_dt function

66eeb3c

Adding basics comments and documentation to new functions

3480a61

Change the position of the calculation of the maximum IMFP

8458eae

Compute maximum of IMFP when IMFP method is 1 too

792c8fd

Compute the maximum IMFP and pass it as a parameter to the compute_dt…

4b43fff

… function.

Cleaning up the compute_dt function and adding documentation.

db60725

Include a generated file to compute the trace of the neutrino number …

d02f865

…density for mesh quantities

Solving issue in generated file to compute trace of mesh neutrino and…

0959355

… antineutrino number densities

Compute maximum trace of the neutrino and antineutrino from mesh quan…

b4c008b

…tities to compute the correct time step

Create a function that computes the minimum value in a MultiFab in pa…

0399808

…rallel. This function will be used to compute the minimum trace of mesh quantities, and based on that, it will compute the optimum time step

Set some if statements to handle simulations with empty periodic boun…

be783b9

…dary conditions and black holes.

Generate a new file to compute the trace on mesh quantities

d8e9f1e

Solve syntaxis mistakes and compute the minumin time step rather than…

c5c70ed

… the minimum potential

Make function compute_max_IMFP return a multifab with the maximum inv…

9703bac

…erse mean free paths in the grid center

Compute the minimum combination of TrN/Opacity*c. This quantity provi…

aada239

…des a better representation of the time derivative of N than the previous 1/Opacity*c.

Delete commented lines and handle periodic empty boundary conditions …

8a737ea

…and black holes in the function that computes the time step.

Cleaning up and addiing comments to the compute_dt function

3420c23

Add documentation for header file of Evolve.cpp file

4594186

Delete unnecessary generation of two files that compute the trace.

acbf75d

Setting the min(TrN, TrNbar) to one in case it is zero

2296774

Cleaning space

97e981e

Cleaning up compute_max_IMFP function

c1be327

Creating documentation for compute_max_IMFP function

af7a2c8

Set the new time step to agree with the new function that computes th…

08fc11f

…e time step.

erickurquilla1999 and others added 30 commits January 4, 2025 17:28

Update write particles per cell script to recive the grid index as in…

bfbba28

…put parameter

Delete file

2ab5ab3

Fix issue with number of amrex grid subdivision

cdd589e

Add line to write pressure moment in the H5 file (it is commented but…

eaa2414

… it can be used when needed)

Solve mistake in help message

44b9535

Add script compute n and f from particles data and compare with the d…

1d06781

…ata generated by EMU

Add the particle interpolator functions in a single python script

2febf93

Moving the particle interpolator class from this script to a

24be1e8

different to avoid repete code

Add scrit to plot angular distributions in cells

2cf0a6e

Change name of python script to generate constant backgorund data for…

3d28566

… test

Delete script that do not need to be in the EMU repository

52fd30c

Create script that convert the cell hdf5 files into a particle_input.…

74ab3c7

…dat file that can be used as an initial condition for the Emu code.

Dividing by volume quantites that should be written as densities

78c2085

Solve issue with amrex GpuTuple with no definition of variable types

05356df

Trying saving particle data with WritePlotFile instead of Checkpoint …

68a5b1b

…to make the process of writing data more efficient. Check points will not be permited for the output file.

Print just time step number and physical time

87aedd9

Making the input parameters of the script that split particle data pe…

d60c5ec

…r cell optional. It run single boxes only if the input parameters are defined otherwise it splits particles in all cells

Make the file that convert average data to hdf5 file read plt file nu…

f334bb5

…mber bigger than 5 digits

Sorting directories when splitting particles data per cell

a680ead

Going back to checkpoints to save particle data

d9d9e04

Generate initial conditions script for local CFI simulations

d2983d0

Add the possibility of defining an arbitrary equilibrium distribution…

4886d59

… different from Fermi–Dirac. This will be useful for local simulations of CFI from NSM global cells, where the equilibrium is not necessarily Fermi–Dirac.

Add paramter to control arbitrary equilibrium in CFI inputs_collision…

8c0b8b9

…al_instability_test

Convert data to hdf5 for all plt file. Previously is only analize the

f0bf719

one that had the neutrino dir inside the plt file

Fix issue when trying to read data from plt files, now it order the d…

0ba3dde

…ir and just take the plt directories with data

Adding attenuation paramters for every type of hamiltonian

c49fe24

Update Jenkinsfile

641fa06

Solve issue when saving pdf in collision to equilibrium test

0807ece

Modifications necessary to include a 3-velocity in the stored grid va…

fb2c16d

…riables. MSW input test changed in order to initialize vup components as a parameter from the parameter class.

added units to velocity parameters for consistency.

0c97875

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Uh oh!

(DRAFT) feat: introduce 3-velocity grid variable to support relativistic transport #109

(DRAFT) feat: introduce 3-velocity grid variable to support relativistic transport #109

cellio00 commented Sep 25, 2025

Uh oh!

Uh oh!

(DRAFT) feat: introduce 3-velocity grid variable to support relativistic transport #109

Are you sure you want to change the base?

(DRAFT) feat: introduce 3-velocity grid variable to support relativistic transport #109

Conversation

cellio00 commented Sep 25, 2025

Uh oh!

Uh oh!