Module containing the physics routines for advection-reaction-diffusion equations. More...
Functions/Subroutines | |
| subroutine, public | ard_phys_init () |
Variables | |
| integer, public, protected | u_ = 1 |
| Indices of the unknowns. | |
| integer, public, protected | v_ = 2 |
| For 2 or more equations. | |
| integer, public, protected | w_ = 3 |
| For 3 or more equations. | |
| logical, public, protected | ard_particles = .false. |
| Whether particles module is added. | |
| double precision, public, protected | dtreacpar = 0.5d0 |
| Parameter with which to multiply the reaction timestep restriction. | |
| character(len=20), public, protected | equation_name = "gray-scott" |
| Name of the system to be solved. | |
| double precision, public, protected | d1 = 0.05d0 |
| Diffusion coefficient for first species (u) | |
| double precision, public, protected | d2 = 1.0d0 |
| Diffusion coefficient for second species (v) (if applicable) | |
| double precision, public, protected | d3 = 1.0d0 |
| Diffusion coefficient for third species (w) (if applicable) | |
| integer, public, protected | adv_pow = 1 |
| Power of the unknown in the advection term (1 for linear) | |
| double precision, dimension(^nd), public, protected | a1 = 0.0d0 |
| Advection coefficients for first species (u) | |
| double precision, dimension(^nd), public, protected | a2 = 0.0d0 |
| Advection coefficients for second species (v) (if applicable) | |
| double precision, dimension(^nd), public, protected | a3 = 0.0d0 |
| Advection coefficients for third species (w) (if applicable) | |
| double precision, public, protected | sb_alpha = 0.1305d0 |
| Parameters for Schnakenberg model. | |
| double precision, public, protected | sb_beta = 0.7695d0 |
| double precision, public, protected | sb_kappa = 100.0d0 |
| double precision, public, protected | gs_f = 0.046d0 |
| Feed rate for Gray-Scott model. | |
| double precision, public, protected | gs_k = 0.063d0 |
| Kill rate for Gray-Scott model. | |
| double precision, public, protected | br_a = 4.5d0 |
| Parameters for Brusselator model. | |
| double precision, public, protected | br_b = 8.0d0 |
| double precision, public, protected | br_c = 1.0d0 |
| double precision, public, protected | br_d = 1.0d0 |
| double precision, public, protected | lg_lambda = 1.0d0 |
| Parameter for logistic model (Fisher / KPP equation) | |
| double precision, public, protected | bzfn_epsilon = 1.0d0 |
| Parameters for the Field-Noyes model of the Belousov-Zhabotinsky reaction. | |
| double precision, public, protected | bzfn_delta = 1.0d0 |
| double precision, public, protected | bzfn_lambda = 1.0d0 |
| double precision, public, protected | bzfn_mu = 1.0d0 |
| double precision, public, protected | lor_r = 28.0d0 |
| Parameter for Lorenz system (Rayleigh number) | |
| double precision, public, protected | lor_sigma = 10.0d0 |
| Parameter for Lorenz system (Prandtl number) | |
| double precision, public, protected | lor_b = 8.0d0 / 3.0d0 |
| Parameter for Lorenz system (aspect ratio of the convection rolls) | |
| type(mg_bc_t), dimension(3, mg_num_neighbors), public | ard_mg_bc |
| Boundary condition information for the multigrid method. | |
Detailed Description
Module containing the physics routines for advection-reaction-diffusion equations.
This module can be seen as an extension of the reaction-diffusion (rd) module and includes the same reaction systems and more: the Gray-Scott model, the Schnakenberg model, the Brusselator model, the diffusive logistic equation, an analytical testcase from "Numerical solution of time-dependent advection- diffusion-reaction equations" by Hundsdorfer & Verwer, the Oregonator model, the extended Brusselator model, the diffusive Lorenz system and the advection- diffusion equation. See the documentation of the advection-reaction-diffusion module for more information.
An advection term can be aplied to these systems of the form: nabla( (A1/adv_pow) * u^(adv_pow) ) (for the first unknown) nabla( (A2/adv_pow) * v^(adv_pow) ) (for the second unknown, if applicable) nabla( (A3/adv_pow) * w^(adv_pow) ) (for the third unknown, if applicable)
IMEX methods are also supported. The implicit system is solved by a multigrid solver coupled into MPI-AMRVAC.
Function/Subroutine Documentation
◆ ard_phys_init()
| subroutine, public mod_ard_phys::ard_phys_init |
Variable Documentation
◆ a1
| double precision, dimension(^nd), public, protected mod_ard_phys::a1 = 0.0d0 |
Advection coefficients for first species (u)
Definition at line 63 of file mod_ard_phys.t.
◆ a2
| double precision, dimension(^nd), public, protected mod_ard_phys::a2 = 0.0d0 |
Advection coefficients for second species (v) (if applicable)
Definition at line 65 of file mod_ard_phys.t.
◆ a3
| double precision, dimension(^nd), public, protected mod_ard_phys::a3 = 0.0d0 |
Advection coefficients for third species (w) (if applicable)
Definition at line 67 of file mod_ard_phys.t.
◆ adv_pow
| integer, public, protected mod_ard_phys::adv_pow = 1 |
Power of the unknown in the advection term (1 for linear)
Definition at line 60 of file mod_ard_phys.t.
◆ ard_mg_bc
| type(mg_bc_t), dimension(3, mg_num_neighbors), public mod_ard_phys::ard_mg_bc |
Boundary condition information for the multigrid method.
Definition at line 105 of file mod_ard_phys.t.
◆ ard_particles
| logical, public, protected mod_ard_phys::ard_particles = .false. |
Whether particles module is added.
Definition at line 33 of file mod_ard_phys.t.
◆ br_a
| double precision, public, protected mod_ard_phys::br_a = 4.5d0 |
Parameters for Brusselator model.
Definition at line 80 of file mod_ard_phys.t.
◆ br_b
| double precision, public, protected mod_ard_phys::br_b = 8.0d0 |
Definition at line 81 of file mod_ard_phys.t.
◆ br_c
| double precision, public, protected mod_ard_phys::br_c = 1.0d0 |
Definition at line 82 of file mod_ard_phys.t.
◆ br_d
| double precision, public, protected mod_ard_phys::br_d = 1.0d0 |
Definition at line 83 of file mod_ard_phys.t.
◆ bzfn_delta
| double precision, public, protected mod_ard_phys::bzfn_delta = 1.0d0 |
Definition at line 90 of file mod_ard_phys.t.
◆ bzfn_epsilon
| double precision, public, protected mod_ard_phys::bzfn_epsilon = 1.0d0 |
Parameters for the Field-Noyes model of the Belousov-Zhabotinsky reaction.
Definition at line 89 of file mod_ard_phys.t.
◆ bzfn_lambda
| double precision, public, protected mod_ard_phys::bzfn_lambda = 1.0d0 |
Definition at line 91 of file mod_ard_phys.t.
◆ bzfn_mu
| double precision, public, protected mod_ard_phys::bzfn_mu = 1.0d0 |
Definition at line 92 of file mod_ard_phys.t.
◆ d1
| double precision, public, protected mod_ard_phys::d1 = 0.05d0 |
Diffusion coefficient for first species (u)
Definition at line 53 of file mod_ard_phys.t.
◆ d2
| double precision, public, protected mod_ard_phys::d2 = 1.0d0 |
Diffusion coefficient for second species (v) (if applicable)
Definition at line 55 of file mod_ard_phys.t.
◆ d3
| double precision, public, protected mod_ard_phys::d3 = 1.0d0 |
Diffusion coefficient for third species (w) (if applicable)
Definition at line 57 of file mod_ard_phys.t.
◆ dtreacpar
| double precision, public, protected mod_ard_phys::dtreacpar = 0.5d0 |
Parameter with which to multiply the reaction timestep restriction.
Definition at line 36 of file mod_ard_phys.t.
◆ equation_name
| character(len=20), public, protected mod_ard_phys::equation_name = "gray-scott" |
Name of the system to be solved.
Definition at line 39 of file mod_ard_phys.t.
◆ gs_f
| double precision, public, protected mod_ard_phys::gs_f = 0.046d0 |
Feed rate for Gray-Scott model.
Definition at line 75 of file mod_ard_phys.t.
◆ gs_k
| double precision, public, protected mod_ard_phys::gs_k = 0.063d0 |
Kill rate for Gray-Scott model.
Definition at line 77 of file mod_ard_phys.t.
◆ lg_lambda
| double precision, public, protected mod_ard_phys::lg_lambda = 1.0d0 |
Parameter for logistic model (Fisher / KPP equation)
Definition at line 86 of file mod_ard_phys.t.
◆ lor_b
| double precision, public, protected mod_ard_phys::lor_b = 8.0d0 / 3.0d0 |
Parameter for Lorenz system (aspect ratio of the convection rolls)
Definition at line 99 of file mod_ard_phys.t.
◆ lor_r
| double precision, public, protected mod_ard_phys::lor_r = 28.0d0 |
Parameter for Lorenz system (Rayleigh number)
Definition at line 95 of file mod_ard_phys.t.
◆ lor_sigma
| double precision, public, protected mod_ard_phys::lor_sigma = 10.0d0 |
Parameter for Lorenz system (Prandtl number)
Definition at line 97 of file mod_ard_phys.t.
◆ sb_alpha
| double precision, public, protected mod_ard_phys::sb_alpha = 0.1305d0 |
Parameters for Schnakenberg model.
Definition at line 70 of file mod_ard_phys.t.
◆ sb_beta
| double precision, public, protected mod_ard_phys::sb_beta = 0.7695d0 |
Definition at line 71 of file mod_ard_phys.t.
◆ sb_kappa
| double precision, public, protected mod_ard_phys::sb_kappa = 100.0d0 |
Definition at line 72 of file mod_ard_phys.t.
◆ u_
| integer, public, protected mod_ard_phys::u_ = 1 |
Indices of the unknowns.
Definition at line 28 of file mod_ard_phys.t.
◆ v_
| integer, public, protected mod_ard_phys::v_ = 2 |
For 2 or more equations.
Definition at line 29 of file mod_ard_phys.t.
◆ w_
| integer, public, protected mod_ard_phys::w_ = 3 |
For 3 or more equations.
Definition at line 30 of file mod_ard_phys.t.
