|
MFC
Exascale flow solver
|
|
MFC
Exascale flow solver
|
Tracks Lagrangian bubbles and couples their dynamics to the Eulerian flow via volume averaging. More...
Functions/Subroutines | |
| impure subroutine | s_initialize_bubbles_el_module (q_cons_vf, bc_type) |
| Initializes the lagrangian subgrid bubble solver. | |
| impure subroutine | s_read_input_bubbles (q_cons_vf, bc_type) |
| Read initial bubble data from input files. | |
| impure subroutine | s_add_bubbles (inputbubble, q_cons_vf, bub_id) |
| Add a new bubble from input data for a fresh start. | |
| impure subroutine | s_restart_bubbles (bub_id, save_count) |
| Restore bubble data from a restart file. | |
| subroutine | s_compute_bubble_el_dynamics (q_prim_vf, bc_type, stage) |
| Contains the bubble dynamics subroutines. | |
| subroutine | s_compute_bubbles_el_source (q_cons_vf, q_prim_vf, rhs_vf) |
| Compute the Lagrangian bubble source terms and add them to the RHS. | |
| subroutine | s_compute_cson_from_pinf (q_prim_vf, pinf, cell, rhol, gamma, pi_inf, cson) |
| Compute the speed of sound from a given driving pressure. | |
| subroutine | s_smear_voidfraction (bc_type) |
| Smear the bubble effects onto the Eulerian grid. | |
| subroutine | s_get_pinf (bub_id, q_prim_vf, ptype, f_pinfl, cell, preterm1, term2, romega) |
| Compute the bubble driving pressure p_inf. | |
| impure subroutine | s_update_lagrange_tdv_rk (q_prim_vf, bc_type, stage) |
| Update Lagrangian bubble variables using TVD Runge-Kutta time stepping. | |
| impure subroutine | s_enforce_el_bubbles_boundary_conditions (q_prim_vf) |
| subroutine | s_locate_cell (pos, cell, scoord) |
| Locate the cell index for a given physical position. | |
| impure subroutine | s_transfer_data_to_tmp () |
| Transfer data into the temporal variables. | |
| logical function | particle_in_domain (pos_part) |
| Determine if a bubble position lies within the current MPI subdomain including ghost cells. | |
| logical function | particle_in_domain_physical (pos_part) |
| Determine if a Lagrangian bubble is within the physical domain excluding ghost cells. | |
| subroutine | s_gradient_dir (q, dq, dir) |
| Compute the gradient of a scalar field using second-order central differences on a non-uniform grid. | |
| impure subroutine | s_write_lag_particles (qtime) |
| Write Lagrangian bubble state data at each time step. | |
| impure subroutine | s_open_lag_bubble_evol () |
| Open the file to write the evolution of the lagrangian bubbles on each time step. | |
| impure subroutine | s_write_lag_bubble_evol (qtime) |
| Write on each time step the changes of the lagrangian bubbles. | |
| impure subroutine | s_close_lag_bubble_evol |
| subroutine | s_open_void_evol |
| impure subroutine | s_write_void_evol (qtime) |
| Write some useful statistics related to the volume fraction of the particles (void fraction). | |
| subroutine | s_close_void_evol |
| impure subroutine | s_write_restart_lag_bubbles (t_step) |
| Write restart files for the Lagrangian bubble solver. | |
| subroutine | s_calculate_lag_bubble_stats () |
| Compute the maximum and minimum radius of each bubble. | |
| impure subroutine | s_open_lag_bubble_stats () |
| impure subroutine | s_write_lag_bubble_stats () |
| Write the maximum and minimum radius of each bubble. | |
| subroutine | s_close_lag_bubble_stats |
| impure subroutine | s_copy_lag_bubble (dest, src) |
| Copy a lag bubble from a src index to a destination index. | |
| impure subroutine | s_finalize_lagrangian_solver () |
Variables | |
| integer, dimension(:,:), allocatable | lag_id |
| Global and local IDs. | |
| real(wp), dimension(:), allocatable | bub_r0 |
| Initial bubble radius. | |
| real(wp), dimension(:), allocatable | rmax_stats |
| Maximum radius. | |
| real(wp), dimension(:), allocatable | rmin_stats |
| Minimum radius. | |
| real(wp), dimension(:), allocatable | gas_mg |
| Bubble's gas mass. | |
| real(wp), dimension(:), allocatable | gas_betat |
| heatflux model (Preston et al., 2007) | |
| real(wp), dimension(:), allocatable | gas_betac |
| massflux model (Preston et al., 2007) | |
| real(wp), dimension(:), allocatable | bub_dphidt |
| subgrid velocity potential (Maeda & Colonius, 2018) | |
| real(wp), dimension(:,:), allocatable | gas_p |
| Pressure in the bubble. | |
| real(wp), dimension(:,:), allocatable | gas_mv |
| Vapor mass in the bubble. | |
| real(wp), dimension(:,:), allocatable | intfc_rad |
| Bubble radius. | |
| real(wp), dimension(:,:), allocatable | intfc_vel |
| Velocity of the bubble interface. | |
| real(wp), dimension(:,:,:), allocatable | mtn_pos |
| Bubble's position. | |
| real(wp), dimension(:,:,:), allocatable | mtn_posprev |
| Bubble's previous position. | |
| real(wp), dimension(:,:,:), allocatable | mtn_vel |
| Bubble's velocity. | |
| real(wp), dimension(:,:,:), allocatable | mtn_s |
| Bubble's computational cell position in real format. | |
| real(wp), dimension(:,:), allocatable | intfc_draddt |
| Time derivative of bubble's radius. | |
| real(wp), dimension(:,:), allocatable | intfc_dveldt |
| Time derivative of bubble's interface velocity. | |
| real(wp), dimension(:,:), allocatable | gas_dpdt |
| Time derivative of gas pressure. | |
| real(wp), dimension(:,:), allocatable | gas_dmvdt |
| Time derivative of the vapor mass in the bubble. | |
| real(wp), dimension(:,:,:), allocatable | mtn_dposdt |
| Time derivative of the bubble's position. | |
| real(wp), dimension(:,:,:), allocatable | mtn_dveldt |
| Time derivative of the bubble's velocity. | |
| integer, private | lag_num_ts |
| Number of time stages in the time-stepping scheme. | |
| real(wp) | rmax_glb |
| real(wp) | rmin_glb |
| Maximum and minimum bubbe size in the local domain. | |
| type(scalar_field), dimension(:), allocatable | q_beta |
| Projection of the lagrangian particles in the Eulerian framework. | |
| type(scalar_field), dimension(:), allocatable | kahan_comp |
| Kahan compensation for q_beta accumulation. | |
| integer | q_beta_idx |
| Size of the q_beta vector field. | |
| integer, parameter | lag_evol_id = 11 |
| integer, parameter | lag_stats_id = 12 |
| integer, parameter | lag_void_id = 13 |
| integer, dimension(:), allocatable | keep_bubble |
| integer, dimension(:,:), allocatable | wrap_bubble_loc |
| integer, dimension(:,:), allocatable | wrap_bubble_dir |
Tracks Lagrangian bubbles and couples their dynamics to the Eulerian flow via volume averaging.
| logical function m_bubbles_el::particle_in_domain | ( | real(wp), dimension(3), intent(in) | pos_part | ) |
Determine if a bubble position lies within the current MPI subdomain including ghost cells.
Definition at line 4144 of file m_bubbles_EL.fpp.f90.
| logical function m_bubbles_el::particle_in_domain_physical | ( | real(wp), dimension(3), intent(in) | pos_part | ) |
Determine if a Lagrangian bubble is within the physical domain excluding ghost cells.
Definition at line 4197 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_add_bubbles | ( | real(wp), dimension(8), intent(in) | inputbubble, |
| type(scalar_field), dimension(sys_size), intent(in) | q_cons_vf, | ||
| integer, intent(in) | bub_id ) |
Add a new bubble from input data for a fresh start.
Definition at line 2133 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_calculate_lag_bubble_stats |
Compute the maximum and minimum radius of each bubble.
Definition at line 4686 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_close_lag_bubble_evol |
| subroutine m_bubbles_el::s_close_lag_bubble_stats |
| subroutine m_bubbles_el::s_close_void_evol |
| subroutine m_bubbles_el::s_compute_bubble_el_dynamics | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_prim_vf, |
| type(integer_field), dimension(1:num_dims,1:2), intent(in) | bc_type, | ||
| integer, intent(in) | stage ) |
Contains the bubble dynamics subroutines.
Definition at line 2392 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_compute_bubbles_el_source | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_cons_vf, |
| type(scalar_field), dimension(sys_size), intent(inout) | q_prim_vf, | ||
| type(scalar_field), dimension(sys_size), intent(inout) | rhs_vf ) |
Compute the Lagrangian bubble source terms and add them to the RHS.
Definition at line 2625 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_compute_cson_from_pinf | ( | type(scalar_field), dimension(sys_size), intent(in) | q_prim_vf, |
| real(wp), intent(in) | pinf, | ||
| integer, dimension(3), intent(in) | cell, | ||
| real(wp), intent(in) | rhol, | ||
| real(wp), intent(in) | gamma, | ||
| real(wp), intent(in) | pi_inf, | ||
| real(wp), intent(out) | cson ) |
Compute the speed of sound from a given driving pressure.
Definition at line 2865 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_copy_lag_bubble | ( | integer, intent(in) | dest, |
| integer, intent(in) | src ) |
Copy a lag bubble from a src index to a destination index.
Definition at line 4798 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_enforce_el_bubbles_boundary_conditions | ( | type(scalar_field), dimension(sys_size), intent(in) | q_prim_vf | ) |
Definition at line 3739 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_finalize_lagrangian_solver |
| subroutine m_bubbles_el::s_get_pinf | ( | integer, intent(in) | bub_id, |
| type(scalar_field), dimension(sys_size), intent(in) | q_prim_vf, | ||
| integer, intent(in) | ptype, | ||
| real(wp), intent(out) | f_pinfl, | ||
| integer, dimension(3), intent(out) | cell, | ||
| real(wp), intent(out), optional | preterm1, | ||
| real(wp), intent(out), optional | term2, | ||
| real(wp), intent(out), optional | romega ) |
Compute the bubble driving pressure p_inf.
Getting p_cell in terms of only the current cell by interpolation
Obtain bilinear interpolation coefficients, based on the current location of the bubble.
Perform bilinear interpolation
Obtain biquadratic interpolation coefficients, based on the current location of the bubble.
Perform biquadratic interpolation
Obtaining the cell volume
Update values
Definition at line 3044 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_gradient_dir | ( | real(stp), dimension(idwbuff(1)%beg:,idwbuff(2)%beg:,idwbuff(3)%beg:), intent(inout) | q, |
| real(stp), dimension(idwbuff(1)%beg:,idwbuff(2)%beg:,idwbuff(3)%beg:), intent(inout) | dq, | ||
| integer, intent(in) | dir ) |
Compute the gradient of a scalar field using second-order central differences on a non-uniform grid.
Definition at line 4213 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_initialize_bubbles_el_module | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_cons_vf, |
| type(integer_field), dimension(1:num_dims,1:2), intent(in) | bc_type ) |
Initializes the lagrangian subgrid bubble solver.
Definition at line 488 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_locate_cell | ( | real(wp), dimension(3), intent(in) | pos, |
| integer, dimension(3), intent(inout) | cell, | ||
| real(wp), dimension(3), intent(out) | scoord ) |
Locate the cell index for a given physical position.
Definition at line 4011 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_open_lag_bubble_evol |
Open the file to write the evolution of the lagrangian bubbles on each time step.
Definition at line 4382 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_open_lag_bubble_stats |
| subroutine m_bubbles_el::s_open_void_evol |
| impure subroutine m_bubbles_el::s_read_input_bubbles | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_cons_vf, |
| type(integer_field), dimension(1:num_dims,1:2), intent(in) | bc_type ) |
Read initial bubble data from input files.
Definition at line 1993 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_restart_bubbles | ( | integer, intent(inout) | bub_id, |
| integer, intent(inout) | save_count ) |
Restore bubble data from a restart file.
Definition at line 2238 of file m_bubbles_EL.fpp.f90.
| subroutine m_bubbles_el::s_smear_voidfraction | ( | type(integer_field), dimension(1:num_dims,1:2), intent(in) | bc_type | ) |
Smear the bubble effects onto the Eulerian grid.
Definition at line 2935 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_transfer_data_to_tmp |
Transfer data into the temporal variables.
Definition at line 4094 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_update_lagrange_tdv_rk | ( | type(scalar_field), dimension(sys_size), intent(in) | q_prim_vf, |
| type(integer_field), dimension(1:num_dims,1:2), intent(in) | bc_type, | ||
| integer, intent(in) | stage ) |
Update Lagrangian bubble variables using TVD Runge-Kutta time stepping.
Definition at line 3369 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_write_lag_bubble_evol | ( | real(wp), intent(in) | qtime | ) |
Write on each time step the changes of the lagrangian bubbles.
Definition at line 4409 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_write_lag_bubble_stats |
Write the maximum and minimum radius of each bubble.
Definition at line 4760 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_write_lag_particles | ( | real(wp), intent(in) | qtime | ) |
Write Lagrangian bubble state data at each time step.
Definition at line 4353 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_write_restart_lag_bubbles | ( | integer, intent(in) | t_step | ) |
Write restart files for the Lagrangian bubble solver.
Definition at line 4548 of file m_bubbles_EL.fpp.f90.
| impure subroutine m_bubbles_el::s_write_void_evol | ( | real(wp), intent(in) | qtime | ) |
Write some useful statistics related to the volume fraction of the particles (void fraction).
Definition at line 4459 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::bub_dphidt |
subgrid velocity potential (Maeda & Colonius, 2018)
Definition at line 362 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::bub_r0 |
Initial bubble radius.
Definition at line 344 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::gas_betac |
massflux model (Preston et al., 2007)
Definition at line 361 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::gas_betat |
heatflux model (Preston et al., 2007)
Definition at line 360 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::gas_dmvdt |
Time derivative of the vapor mass in the bubble.
Definition at line 411 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::gas_dpdt |
Time derivative of gas pressure.
Definition at line 410 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::gas_mg |
Bubble's gas mass.
Definition at line 359 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::gas_mv |
Vapor mass in the bubble.
Definition at line 377 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::gas_p |
Pressure in the bubble.
Definition at line 376 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::intfc_draddt |
Time derivative of bubble's radius.
Definition at line 408 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::intfc_dveldt |
Time derivative of bubble's interface velocity.
Definition at line 409 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::intfc_rad |
Bubble radius.
Definition at line 378 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:), allocatable m_bubbles_el::intfc_vel |
Velocity of the bubble interface.
Definition at line 379 of file m_bubbles_EL.fpp.f90.
| type(scalar_field), dimension(:), allocatable m_bubbles_el::kahan_comp |
Kahan compensation for q_beta accumulation.
Definition at line 442 of file m_bubbles_EL.fpp.f90.
| integer, dimension(:), allocatable m_bubbles_el::keep_bubble |
Definition at line 460 of file m_bubbles_EL.fpp.f90.
| integer, parameter m_bubbles_el::lag_evol_id = 11 |
Definition at line 457 of file m_bubbles_EL.fpp.f90.
| integer, dimension(:,:), allocatable m_bubbles_el::lag_id |
Global and local IDs.
Definition at line 343 of file m_bubbles_EL.fpp.f90.
|
private |
Number of time stages in the time-stepping scheme.
Definition at line 426 of file m_bubbles_EL.fpp.f90.
| integer, parameter m_bubbles_el::lag_stats_id = 12 |
Definition at line 458 of file m_bubbles_EL.fpp.f90.
| integer, parameter m_bubbles_el::lag_void_id = 13 |
Definition at line 459 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_dposdt |
Time derivative of the bubble's position.
Definition at line 412 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_dveldt |
Time derivative of the bubble's velocity.
Definition at line 413 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_pos |
Bubble's position.
Definition at line 392 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_posprev |
Bubble's previous position.
Definition at line 393 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_s |
Bubble's computational cell position in real format.
Definition at line 395 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:,:,:), allocatable m_bubbles_el::mtn_vel |
Bubble's velocity.
Definition at line 394 of file m_bubbles_EL.fpp.f90.
| type(scalar_field), dimension(:), allocatable m_bubbles_el::q_beta |
Projection of the lagrangian particles in the Eulerian framework.
Definition at line 441 of file m_bubbles_EL.fpp.f90.
| integer m_bubbles_el::q_beta_idx |
Size of the q_beta vector field.
Definition at line 443 of file m_bubbles_EL.fpp.f90.
| real(wp) m_bubbles_el::rmax_glb |
Definition at line 439 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::rmax_stats |
Maximum radius.
Definition at line 345 of file m_bubbles_EL.fpp.f90.
| real(wp) m_bubbles_el::rmin_glb |
Maximum and minimum bubbe size in the local domain.
Definition at line 439 of file m_bubbles_EL.fpp.f90.
| real(wp), dimension(:), allocatable m_bubbles_el::rmin_stats |
Minimum radius.
Definition at line 346 of file m_bubbles_EL.fpp.f90.
| integer, dimension(:,:), allocatable m_bubbles_el::wrap_bubble_dir |
Definition at line 461 of file m_bubbles_EL.fpp.f90.
| integer, dimension(:,:), allocatable m_bubbles_el::wrap_bubble_loc |
Definition at line 461 of file m_bubbles_EL.fpp.f90.