Flow123d  JS_before_hm-1008-g3dab983
heat_model.cc
Go to the documentation of this file.
1 /*!
2  *
3  * Copyright (C) 2015 Technical University of Liberec. All rights reserved.
4  *
5  * This program is free software; you can redistribute it and/or modify it under
6  * the terms of the GNU General Public License version 3 as published by the
7  * Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
8  *
9  * This program is distributed in the hope that it will be useful, but WITHOUT
10  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11  * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
12  *
13  *
14  * @file heat_model.cc
15  * @brief Discontinuous Galerkin method for equation of transport with dispersion.
16  * @author Jan Stebel
17  */
18 
19 #include "input/input_type.hh"
20 #include "mesh/mesh.h"
21 #include "mesh/accessors.hh"
22 //#include "transport/transport_operator_splitting.hh"
23 #include "heat_model.hh"
24 #include "tools/unit_si.hh"
25 #include "coupling/balance.hh"
26 
27 
28 
29 using namespace std;
30 using namespace Input::Type;
31 
32 
33 
34 
35 
36 
37 
38 
40  return Selection("Heat_BC_Type", "Types of boundary conditions for heat transfer model.")
41  .add_value(bc_inflow, "inflow",
42  "Default heat transfer boundary condition.\n"
43  "On water inflow (($(q_w \\le 0)$)), total energy flux is given by the reference temperature 'bc_temperature'. "
44  "On water outflow we prescribe zero diffusive flux, "
45  "i.e. the energy flows out only due to advection.")
46  .add_value(bc_dirichlet, "dirichlet",
47  "Dirichlet boundary condition (($T = T_D $)).\n"
48  "The prescribed temperature (($T_D$)) is specified by the field 'bc_temperature'.")
49  .add_value(bc_total_flux, "total_flux",
50  "Total energy flux boundary condition.\n"
51  "The prescribed incoming total flux can have the general form (($\\delta(f_N+\\sigma_R(T_R-T) )$)), "
52  "where the absolute flux (($f_N$)) is specified by the field 'bc_flux', "
53  "the transition parameter (($\\sigma_R$)) by 'bc_robin_sigma', "
54  "and the reference temperature (($T_R$)) by 'bc_temperature'.")
55  .add_value(bc_diffusive_flux, "diffusive_flux",
56  "Diffusive flux boundary condition.\n"
57  "The prescribed incoming energy flux due to diffusion can have the general form (($\\delta(f_N+\\sigma_R(T_R-T) )$)), "
58  "where the absolute flux (($f_N$)) is specified by the field 'bc_flux', "
59  "the transition parameter (($\\sigma_R$)) by 'bc_robin_sigma', "
60  "and the reference temperature (($T_R$)) by 'bc_temperature'.")
61  .close();
62 }
63 
64 
66 {
67  *this+=bc_type
68  .name("bc_type")
69  .description(
70  "Type of boundary condition.")
71  .units( UnitSI::dimensionless() )
72  .input_default("\"inflow\"")
73  .input_selection( get_bc_type_selection() )
75  *this+=bc_dirichlet_value
76  .name("bc_temperature")
77  .description("Boundary value of temperature.")
78  .units( UnitSI().K() )
79  .input_default("0.0")
80  .flags_add(in_rhs);
81  *this+=bc_flux
82  .disable_where(bc_type, { bc_dirichlet, bc_inflow })
83  .name("bc_flux")
84  .description("Flux in Neumann boundary condition.")
85  .units( UnitSI().kg().m().s(-1).md() )
86  .input_default("0.0")
87  .flags_add(FieldFlag::in_rhs);
88  *this+=bc_robin_sigma
89  .disable_where(bc_type, { bc_dirichlet, bc_inflow })
90  .name("bc_robin_sigma")
91  .description("Conductivity coefficient in Robin boundary condition.")
92  .units( UnitSI().m(4).s(-1).md() )
93  .input_default("0.0")
94  .flags_add(FieldFlag::in_rhs & FieldFlag::in_main_matrix);
95 
96  *this+=init_temperature
97  .name("init_temperature")
98  .description("Initial temperature.")
99  .units( UnitSI().K() )
100  .input_default("0.0");
101 
102  *this+=porosity
103  .name("porosity")
104  .description("Porosity.")
105  .units( UnitSI::dimensionless() )
106  .input_default("1.0")
107  .flags_add(in_main_matrix & in_time_term)
108  .set_limits(0.0);
109 
110  *this+=water_content
111  .name("water_content")
112  .units( UnitSI::dimensionless() )
113  .input_default("1.0")
114  .flags_add(input_copy & in_main_matrix & in_time_term);
115 
116  *this += flow_flux.name("flow_flux")
117  .flags( FieldFlag::input_copy )
118  .flags_add(in_time_term & in_main_matrix & in_rhs);
119 
120  *this+=fluid_density
121  .name("fluid_density")
122  .description("Density of fluid.")
123  .units( UnitSI().kg().m(-3) )
124  .input_default("1000")
125  .flags_add(in_main_matrix & in_time_term);
126 
127  *this+=fluid_heat_capacity
128  .name("fluid_heat_capacity")
129  .description("Heat capacity of fluid.")
130  .units( UnitSI::J() * UnitSI().kg(-1).K(-1) )
131  .flags_add(in_main_matrix & in_time_term);
132 
133  *this+=fluid_heat_conductivity
134  .name("fluid_heat_conductivity")
135  .description("Heat conductivity of fluid.")
136  .units( UnitSI::W() * UnitSI().m(-1).K(-1) )
137  .flags_add(in_main_matrix)
138  .set_limits(0.0);
139 
140 
141  *this+=solid_density
142  .name("solid_density")
143  .description("Density of solid (rock).")
144  .units( UnitSI().kg().m(-3) )
145  .flags_add(in_time_term);
146 
147  *this+=solid_heat_capacity
148  .name("solid_heat_capacity")
149  .description("Heat capacity of solid (rock).")
150  .units( UnitSI::J() * UnitSI().kg(-1).K(-1) )
151  .flags_add(in_time_term);
152 
153  *this+=solid_heat_conductivity
154  .name("solid_heat_conductivity")
155  .description("Heat conductivity of solid (rock).")
156  .units( UnitSI::W() * UnitSI().m(-1).K(-1) )
157  .flags_add(in_main_matrix)
158  .set_limits(0.0);
159 
160  *this+=disp_l
161  .name("disp_l")
162  .description("Longitudinal heat dispersivity in fluid.")
163  .units( UnitSI().m() )
164  .input_default("0.0")
165  .flags_add(in_main_matrix);
166 
167  *this+=disp_t
168  .name("disp_t")
169  .description("Transverse heat dispersivity in fluid.")
170  .units( UnitSI().m() )
171  .input_default("0.0")
172  .flags_add(in_main_matrix);
173 
174  *this+=fluid_thermal_source
175  .name("fluid_thermal_source")
176  .description("Density of thermal source in fluid.")
177  .units( UnitSI::W() * UnitSI().m(-3) )
178  .input_default("0.0")
179  .flags_add(in_rhs);
180 
181  *this+=solid_thermal_source
182  .name("solid_thermal_source")
183  .description("Density of thermal source in solid.")
184  .units( UnitSI::W() * UnitSI().m(-3) )
185  .input_default("0.0")
186  .flags_add(in_rhs);
187 
188  *this+=fluid_heat_exchange_rate
189  .name("fluid_heat_exchange_rate")
190  .description("Heat exchange rate of source in fluid.")
191  .units( UnitSI().s(-1) )
192  .input_default("0.0")
193  .flags_add(in_rhs);
194 
195  *this+=solid_heat_exchange_rate
196  .name("solid_heat_exchange_rate")
197  .description("Heat exchange rate of source in solid.")
198  .units( UnitSI().s(-1) )
199  .input_default("0.0")
200  .flags_add(in_rhs);
201 
202  *this+=fluid_ref_temperature
203  .name("fluid_ref_temperature")
204  .description("Reference temperature of source in fluid.")
205  .units( UnitSI().K() )
206  .input_default("0.0")
207  .flags_add(in_rhs);
208 
209  *this+=solid_ref_temperature
210  .name("solid_ref_temperature")
211  .description("Reference temperature in solid.")
212  .units( UnitSI().K() )
213  .input_default("0.0")
214  .flags_add(in_rhs);
215 
216  *this+=cross_section
217  .name("cross_section")
218  .units( UnitSI().m(3).md() )
219  .flags(input_copy & in_time_term & in_main_matrix);
220 
221  *this+=output_field
222  .name("temperature")
223  .description("Temperature solution.")
224  .units( UnitSI().K() )
225  .flags(equation_result);
226 }
227 
228 
229 
230 IT::Record HeatTransferModel::get_input_type(const string &implementation, const string &description)
231 {
232  return IT::Record(
233  std::string(ModelEqData::name()) + "_" + implementation,
234  description + " for heat transfer.")
237  .declare_key("balance", Balance::get_input_type(), Default("{}"),
238  "Settings for computing balance.")
239  .declare_key("output_stream", OutputTime::get_input_type(), Default("{}"),
240  "Parameters of output stream.");
241 }
242 
243 
245 {
246  // Return empty selection just to provide model specific selection name and description.
247  // The fields are added by TransportDG using an auxiliary selection.
248  return IT::Selection(
249  std::string(ModelEqData::name()) + "_DG_output_fields",
250  "Selection of output fields for Heat Transfer DG model.");
251 }
252 
253 
255  AdvectionProcessBase(mesh, in_rec)
256 {
257  time_ = new TimeGovernor(in_rec.val<Input::Record>("time"));
258  ASSERT( time_->is_default() == false ).error("Missing key 'time' in Heat_AdvectionDiffusion_DG.");
259  substances_.initialize({""});
260 
261  output_stream_ = OutputTime::create_output_stream("heat", in_rec.val<Input::Record>("output_stream"), time().get_unit_string());
262  //output_stream_->add_admissible_field_names(in_rec.val<Input::Array>("output_fields"));
263 
264  balance_ = std::make_shared<Balance>("energy", mesh_);
265  balance_->init_from_input(in_rec.val<Input::Record>("balance"), *time_);
266  // initialization of balance object
267  subst_idx = {balance_->add_quantity("energy")};
268  balance_->units(UnitSI().m(2).kg().s(-2));
269 }
270 
271 
273 {
274  output_stream_->write_time_frame();
275 }
276 
277 
279  const ElementAccessor<3> &ele_acc,
280  std::vector<double> &mm_coef)
281 {
282  vector<double> elem_csec(point_list.size()),
283  por(point_list.size()),
284  f_rho(point_list.size()),
285  s_rho(point_list.size()),
286  f_c(point_list.size()),
287  s_c(point_list.size());
288 
289  data().cross_section.value_list(point_list, ele_acc, elem_csec);
290  data().porosity.value_list(point_list, ele_acc, por);
291  data().fluid_density.value_list(point_list, ele_acc, f_rho);
292  data().fluid_heat_capacity.value_list(point_list, ele_acc, f_c);
293  data().solid_density.value_list(point_list, ele_acc, s_rho);
294  data().solid_heat_capacity.value_list(point_list, ele_acc, s_c);
295 
296  for (unsigned int i=0; i<point_list.size(); i++)
297  mm_coef[i] = elem_csec[i]*(por[i]*f_rho[i]*f_c[i] + (1.-por[i])*s_rho[i]*s_c[i]);
298 }
299 
300 
302  const std::vector<arma::vec3> &velocity,
303  const ElementAccessor<3> &ele_acc,
306 {
307  const unsigned int qsize = point_list.size();
308  std::vector<double> f_rho(qsize), f_cap(qsize), f_cond(qsize),
309  s_cond(qsize), por(qsize), csection(qsize), disp_l(qsize), disp_t(qsize);
310 
311  data().fluid_density.value_list(point_list, ele_acc, f_rho);
312  data().fluid_heat_capacity.value_list(point_list, ele_acc, f_cap);
313  data().fluid_heat_conductivity.value_list(point_list, ele_acc, f_cond);
314  data().solid_heat_conductivity.value_list(point_list, ele_acc, s_cond);
315  data().disp_l.value_list(point_list, ele_acc, disp_l);
316  data().disp_t.value_list(point_list, ele_acc, disp_t);
317  data().porosity.value_list(point_list, ele_acc, por);
318  data().cross_section.value_list(point_list, ele_acc, csection);
319 
320  for (unsigned int k=0; k<qsize; k++) {
321  ad_coef[0][k] = velocity[k]*f_rho[k]*f_cap[k];
322 
323  // dispersive part of thermal diffusion
324  // Note that the velocity vector is in fact the Darcian flux,
325  // so to obtain |v| we have to divide vnorm by porosity and cross_section.
326  double vnorm = arma::norm(velocity[k], 2);
327  if (fabs(vnorm) > 0)
328  for (int i=0; i<3; i++)
329  for (int j=0; j<3; j++)
330  dif_coef[0][k](i,j) = ((velocity[k][i]/vnorm)*(velocity[k][j]/vnorm)*(disp_l[k]-disp_t[k]) + disp_t[k]*(i==j?1:0))
331  *vnorm*f_rho[k]*f_cond[k];
332  else
333  dif_coef[0][k].zeros();
334 
335  // conductive part of thermal diffusion
336  dif_coef[0][k] += csection[k]*(por[k]*f_cond[k] + (1.-por[k])*s_cond[k])*arma::eye(3,3);
337  }
338 }
339 
340 
342  const ElementAccessor<3> &ele_acc,
343  std::vector<std::vector<double> > &init_values)
344 {
345  data().init_temperature.value_list(point_list, ele_acc, init_values[0]);
346 }
347 
348 
350  arma::uvec &bc_types)
351 {
352  // Currently the bc types for HeatTransfer are numbered in the same way as in TransportDG.
353  // In general we should use some map here.
354  bc_types = { data().bc_type.value(ele_acc.centre(), ele_acc) };
355 }
356 
357 
358 void HeatTransferModel::get_flux_bc_data(unsigned int index,
359  const Armor::array &point_list,
360  const ElementAccessor<3> &ele_acc,
361  std::vector< double > &bc_flux,
362  std::vector< double > &bc_sigma,
363  std::vector< double > &bc_ref_value)
364 {
365  data().bc_flux.value_list(point_list, ele_acc, bc_flux);
366  data().bc_robin_sigma.value_list(point_list, ele_acc, bc_sigma);
367  data().bc_dirichlet_value[index].value_list(point_list, ele_acc, bc_ref_value);
368 
369  // Change sign in bc_flux since internally we work with outgoing fluxes.
370  for (auto f : bc_flux) f = -f;
371 }
372 
374  const Armor::array &point_list,
375  const ElementAccessor<3> &ele_acc,
376  std::vector< double > &bc_sigma)
377 {
378  data().bc_robin_sigma.value_list(point_list, ele_acc, bc_sigma);
379 }
380 
381 
383  const ElementAccessor<3> &ele_acc,
384  std::vector<std::vector<double> > &sources_value,
385  std::vector<std::vector<double> > &sources_density,
386  std::vector<std::vector<double> > &sources_sigma)
387 {
388  const unsigned int qsize = point_list.size();
389  std::vector<double> por(qsize), csection(qsize), f_rho(qsize), s_rho(qsize), f_cap(qsize), s_cap(qsize),
390  f_source(qsize), s_source(qsize), f_sigma(qsize), s_sigma(qsize), f_temp(qsize), s_temp(qsize);
391  data().porosity.value_list(point_list, ele_acc, por);
392  data().cross_section.value_list(point_list, ele_acc, csection);
393  data().fluid_density.value_list(point_list, ele_acc, f_rho);
394  data().solid_density.value_list(point_list, ele_acc, s_rho);
395  data().fluid_heat_capacity.value_list(point_list, ele_acc, f_cap);
396  data().solid_heat_capacity.value_list(point_list, ele_acc, s_cap);
397  data().fluid_thermal_source.value_list(point_list, ele_acc, f_source);
398  data().solid_thermal_source.value_list(point_list, ele_acc, s_source);
399  data().fluid_heat_exchange_rate.value_list(point_list, ele_acc, f_sigma);
400  data().solid_heat_exchange_rate.value_list(point_list, ele_acc, s_sigma);
401  data().fluid_ref_temperature.value_list(point_list, ele_acc, f_temp);
402  data().solid_ref_temperature.value_list(point_list, ele_acc, s_temp);
403 
404  sources_density[0].resize(point_list.size());
405  sources_sigma[0].resize(point_list.size());
406  sources_value[0].resize(point_list.size());
407  for (unsigned int k=0; k<point_list.size(); k++)
408  {
409  sources_density[0][k] = csection[k]*(por[k]*f_source[k] + (1.-por[k])*s_source[k]);
410  sources_sigma[0][k] = csection[k]*(por[k]*f_rho[k]*f_cap[k]*f_sigma[k] + (1.-por[k])*s_rho[k]*s_cap[k]*s_sigma[k]);
411  if (fabs(sources_sigma[0][k]) > numeric_limits<double>::epsilon())
412  sources_value[0][k] = csection[k]*(por[k]*f_rho[k]*f_cap[k]*f_sigma[k]*f_temp[k]
413  + (1.-por[k])*s_rho[k]*s_cap[k]*s_sigma[k]*s_temp[k])/sources_sigma[0][k];
414  else
415  sources_value[0][k] = 0;
416  }
417 }
418 
419 
421  const ElementAccessor<3> &ele_acc,
422  std::vector<std::vector<double> > &sources_sigma)
423 {
424  const unsigned int qsize = point_list.size();
425  std::vector<double> por(qsize), csection(qsize), f_rho(qsize), s_rho(qsize), f_cap(qsize), s_cap(qsize),
426  f_source(qsize), s_source(qsize), f_sigma(qsize), s_sigma(qsize), f_temp(qsize), s_temp(qsize);
427  data().porosity.value_list(point_list, ele_acc, por);
428  data().cross_section.value_list(point_list, ele_acc, csection);
429  data().fluid_density.value_list(point_list, ele_acc, f_rho);
430  data().solid_density.value_list(point_list, ele_acc, s_rho);
431  data().fluid_heat_capacity.value_list(point_list, ele_acc, f_cap);
432  data().solid_heat_capacity.value_list(point_list, ele_acc, s_cap);
433  data().fluid_heat_exchange_rate.value_list(point_list, ele_acc, f_sigma);
434  data().solid_heat_exchange_rate.value_list(point_list, ele_acc, s_sigma);
435  sources_sigma[0].resize(point_list.size());
436  for (unsigned int k=0; k<point_list.size(); k++)
437  {
438  sources_sigma[0][k] = csection[k]*(por[k]*f_rho[k]*f_cap[k]*f_sigma[k] + (1.-por[k])*s_rho[k]*s_cap[k]*s_sigma[k]);
439  }
440 }
441 
442 
444 {}
445 
446 
447 
448 
TimeGovernor & time()
Definition: equation.hh:149
vector< unsigned int > subst_idx
List of indices used to call balance methods for a set of quantities.
Definition: heat_model.hh:261
Field< 3, FieldValue< 3 >::Scalar > solid_heat_capacity
Heat capacity of solid.
Definition: heat_model.hh:134
static const Input::Type::Record & get_input_type()
Main balance input record type.
Definition: balance.cc:50
Field< 3, FieldValue< 3 >::Scalar > fluid_density
Density of fluid.
Definition: heat_model.hh:126
unsigned int size() const
Definition: armor.hh:718
void compute_advection_diffusion_coefficients(const Armor::array &point_list, const std::vector< arma::vec3 > &velocity, const ElementAccessor< 3 > &ele_acc, std::vector< std::vector< arma::vec3 > > &ad_coef, std::vector< std::vector< arma::mat33 > > &dif_coef) override
Definition: heat_model.cc:301
static constexpr Mask in_main_matrix
A field is part of main "stiffness matrix" of the equation.
Definition: field_flag.hh:49
Field< 3, FieldValue< 3 >::Scalar > disp_t
Transversal heat dispersivity.
Definition: heat_model.hh:140
Field< 3, FieldValue< 3 >::Scalar > cross_section
Pointer to DarcyFlow field cross_section.
Definition: heat_model.hh:155
static const Input::Type::Record & get_input_type()
The specification of output stream.
Definition: output_time.cc:37
Class Input::Type::Default specifies default value of keys of a Input::Type::Record.
Definition: type_record.hh:61
static IT::Selection get_output_selection()
Definition: heat_model.cc:244
void initialize(const Input::Array &in_array)
Read from input array.
Definition: substance.cc:58
void get_flux_bc_sigma(unsigned int index, const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< double > &bc_sigma) override
Return transition coefficient for flux b.c.
Definition: heat_model.cc:373
Field< 3, FieldValue< 3 >::Scalar > fluid_heat_exchange_rate
Heat exchange rate in fluid.
Definition: heat_model.hh:146
static std::shared_ptr< OutputTime > create_output_stream(const std::string &equation_name, const Input::Record &in_rec, std::string unit_str)
This method delete all object instances of class OutputTime stored in output_streams vector...
Definition: output_time.cc:184
Field< 3, FieldValue< 3 >::Scalar > solid_heat_exchange_rate
Heat exchange rate in solid.
Definition: heat_model.hh:148
Definition: mesh.h:78
BCMultiField< 3, FieldValue< 3 >::Scalar > bc_dirichlet_value
Dirichlet boundary condition for temperature.
Definition: heat_model.hh:112
#define ASSERT(expr)
Allow use shorter versions of macro names if these names is not used with external library...
Definition: asserts.hh:347
BCField< 3, FieldValue< 3 >::Enum > bc_type
Type of boundary condition (see also BC_Type)
Definition: heat_model.hh:110
void compute_init_cond(const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< std::vector< double > > &init_values) override
Definition: heat_model.cc:341
static const Input::Type::Selection & get_bc_type_selection()
Definition: heat_model.cc:39
Field< 3, FieldValue< 3 >::Scalar > fluid_heat_capacity
Heat capacity of fluid.
Definition: heat_model.hh:128
Basic time management functionality for unsteady (and steady) solvers (class Equation).
BCField< 3, FieldValue< 3 >::Scalar > bc_flux
Flux value in total/diffusive flux b.c.
Definition: heat_model.hh:114
void get_bc_type(const ElementAccessor< 3 > &ele_acc, arma::uvec &bc_types) override
Definition: heat_model.cc:349
Field< 3, FieldValue< 3 >::Scalar > porosity
Porosity of solid.
Definition: heat_model.hh:120
Field< 3, FieldValue< 3 >::Scalar > fluid_thermal_source
Thermal source in fluid.
Definition: heat_model.hh:142
static Input::Type::Record & record_template()
Template Record with common keys for derived equations.
Definition: equation.cc:35
void compute_source_coefficients(const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< std::vector< double > > &sources_conc, std::vector< std::vector< double > > &sources_density, std::vector< std::vector< double > > &sources_sigma) override
Definition: heat_model.cc:382
arma::vec::fixed< spacedim > centre() const
Computes the barycenter.
virtual Record & derive_from(Abstract &parent)
Method to derive new Record from an AbstractRecord parent.
Definition: type_record.cc:196
Field< 3, FieldValue< 3 >::Scalar > fluid_heat_conductivity
Heat conductivity of fluid.
Definition: heat_model.hh:130
static Input::Type::Abstract & get_input_type()
Common specification of the input record for secondary equations.
void compute_mass_matrix_coefficient(const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< double > &mm_coef) override
Definition: heat_model.cc:278
#define FMT_UNUSED
Definition: posix.h:75
Field< 3, FieldValue< 3 >::Scalar > fluid_ref_temperature
Reference temperature in fluid.
Definition: heat_model.hh:150
void get_flux_bc_data(unsigned int index, const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< double > &bc_flux, std::vector< double > &bc_sigma, std::vector< double > &bc_ref_value) override
Return data for diffusive or total flux b.c.
Definition: heat_model.cc:358
static constexpr Mask input_copy
Definition: field_flag.hh:44
Field< 3, FieldValue< 3 >::Scalar > init_temperature
Initial temperature.
Definition: heat_model.hh:118
Field< 3, FieldValue< 3 >::Scalar > solid_thermal_source
Thermal source in solid.
Definition: heat_model.hh:144
static UnitSI & W()
Returns Watt.
Definition: unit_si.cc:45
Accessor to the data with type Type::Record.
Definition: accessors.hh:291
const Ret val(const string &key) const
static UnitSI & J()
Returns Joule.
Definition: unit_si.cc:40
Mesh * mesh_
Definition: equation.hh:218
Selection & add_value(const int value, const std::string &key, const std::string &description="", TypeBase::attribute_map attributes=TypeBase::attribute_map())
Adds one new value with name given by key to the Selection.
virtual Value::return_type const & value(const Point &p, const ElementAccessor< spacedim > &elm) const
Definition: field.hh:434
std::shared_ptr< Balance > balance_
object for calculation and writing the mass balance to file.
Definition: equation.hh:230
Record & declare_key(const string &key, std::shared_ptr< TypeBase > type, const Default &default_value, const string &description, TypeBase::attribute_map key_attributes=TypeBase::attribute_map())
Declares a new key of the Record.
Definition: type_record.cc:503
Field< 3, FieldValue< 3 >::Scalar > solid_heat_conductivity
Heat conductivity of solid.
Definition: heat_model.hh:136
SubstanceList substances_
Transported substances.
Definition: heat_model.hh:258
static constexpr Mask in_rhs
A field is part of the right hand side of the equation.
Definition: field_flag.hh:51
virtual void value_list(const Armor::array &point_list, const ElementAccessor< spacedim > &elm, std::vector< typename Value::return_type > &value_list) const
Definition: field.hh:448
Record & copy_keys(const Record &other)
Copy keys from other record.
Definition: type_record.cc:216
void compute_sources_sigma(const Armor::array &point_list, const ElementAccessor< 3 > &ele_acc, std::vector< std::vector< double > > &sources_sigma) override
Definition: heat_model.cc:420
Field< 3, FieldValue< 3 >::Scalar > solid_ref_temperature
Reference temperature in solid.
Definition: heat_model.hh:152
const Selection & close() const
Close the Selection, no more values can be added.
virtual ModelEqData & data()=0
Derived class should implement getter for ModelEqData instance.
~HeatTransferModel() override
Definition: heat_model.cc:443
Field< 3, FieldValue< 3 >::Scalar > solid_density
Density of solid.
Definition: heat_model.hh:132
Discontinuous Galerkin method for equation of transport with dispersion.
HeatTransferModel(Mesh &mesh, const Input::Record in_rec)
Definition: heat_model.cc:254
Record type proxy class.
Definition: type_record.hh:182
BCField< 3, FieldValue< 3 >::Scalar > bc_robin_sigma
Transition coefficient in total/diffusive flux b.c.
Definition: heat_model.hh:116
std::shared_ptr< OutputTime > output_stream_
Definition: heat_model.hh:263
Class for representation SI units of Fields.
Definition: unit_si.hh:40
static UnitSI & dimensionless()
Returns dimensionless unit.
Definition: unit_si.cc:55
Template for classes storing finite set of named values.
void output_data() override
Write computed fields.
Definition: heat_model.cc:272
Field< 3, FieldValue< 3 >::Scalar > disp_l
Longitudal heat dispersivity.
Definition: heat_model.hh:138
TimeGovernor * time_
Definition: equation.hh:219