Discontinuous Galerkin method for equation of transport with dispersion. More...
#include "system/sys_profiler.hh"#include "transport/transport_dg.hh"#include "io/output_time.hh"#include "quadrature/quadrature_lib.hh"#include "fem/mapping_p1.hh"#include "fem/fe_values.hh"#include "fem/fe_p.hh"#include "fem/fe_rt.hh"#include "fem/dh_cell_accessor.hh"#include "fields/field_fe.hh"#include "fields/fe_value_handler.hh"#include "la/linsys_PETSC.hh"#include "flow/mh_dofhandler.hh"#include "transport/advection_diffusion_model.hh"#include "transport/concentration_model.hh"#include "transport/heat_model.hh"#include "coupling/balance.hh"#include "fields/multi_field.hh"#include "fields/generic_field.hh"#include "input/factory.hh"#include "fields/equation_output.hh"#include "mesh/long_idx.hh"#include "mesh/accessors.hh"Go to the source code of this file.
Macros | |
| #define | AVERAGE(i, k, side_id) (fe_values[sd[side_id]]->shape_value(i,k)*0.5) |
| #define | WAVERAGE(i, k, side_id) (arma::dot(dif_coef_edg[sd[side_id]][sbi][k]*fe_values[sd[side_id]]->shape_grad(i,k),nv)*omega[side_id]) |
| #define | JUMP(i, k, side_id) ((side_id==0?1:-1)*fe_values[sd[side_id]]->shape_value(i,k)) |
Functions | |
| double | elem_anisotropy (ElementAccessor< 3 > e) |
Detailed Description
Discontinuous Galerkin method for equation of transport with dispersion.
* Copyright (C) 2015 Technical University of Liberec. All rights reserved.
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License version 3 as published by the Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
Definition in file transport_dg.cc.
Macro Definition Documentation
| #define AVERAGE | ( | i, | |
| k, | |||
| side_id | |||
| ) | (fe_values[sd[side_id]]->shape_value(i,k)*0.5) |
| #define JUMP | ( | i, | |
| k, | |||
| side_id | |||
| ) | ((side_id==0?1:-1)*fe_values[sd[side_id]]->shape_value(i,k)) |
| #define WAVERAGE | ( | i, | |
| k, | |||
| side_id | |||
| ) | (arma::dot(dif_coef_edg[sd[side_id]][sbi][k]*fe_values[sd[side_id]]->shape_grad(i,k),nv)*omega[side_id]) |
Function Documentation
| double elem_anisotropy | ( | ElementAccessor< 3 > | e | ) |
