Detailed Description

Global assembler for the monolithic scheme of the non-isothermal Richards flow.

Governing equations without vapor diffusion

The energy balance equation is given by

\[ (\rho c_p)^{eff}\dot T - \nabla (\mathbf{k}_T^{eff} \nabla T)+\rho^l c_p^l \nabla T \cdot \mathbf{v}^l = Q_T \]

with \(T\) the temperature, \((\rho c_p)^{eff}\) the effective volumetric heat capacity, \(\mathbf{k}_T^{eff} \) the effective thermal conductivity, \(\rho^l\) the density of liquid, \(c_p^l\) the specific heat capacity of liquid, \(\mathbf{v}^l\) the liquid velocity, and \(Q_T\) the point heat source. The effective volumetric heat can be considered as a composite of the contributions of solid phase and the liquid phase as

\[ (\rho c_p)^{eff} = (1-\phi) \rho^s c_p^s + S^l \phi \rho^l c_p^l \]

with \(\phi\) the porosity, \(S^l\) the liquid saturation, \(\rho^s \) the solid density, and \(c_p^s\) the specific heat capacity of solid. Similarly, the effective thermal conductivity is given by

\[ \mathbf{k}_T^{eff} = (1-\phi) \mathbf{k}_T^s + S^l \phi k_T^l \mathbf I \]

where \(\mathbf{k}_T^s\) is the thermal conductivity tensor of solid, \( k_T^l\) is the thermal conductivity of liquid, and \(\mathbf I\) is the identity tensor.

The mass balance equation is given by

\begin{eqnarray*} \left(S^l\beta - \phi\frac{\partial S}{\partial p_c}\right) \rho^l\dot p - S \left( \frac{\partial \rho^l}{\partial T} +\rho^l(\alpha_B -S) \alpha_T^s \right)\dot T\\ +\nabla (\rho^l \mathbf{v}^l) + S \alpha_B \rho^l \nabla \cdot \dot {\mathbf u}= Q_H \end{eqnarray*}

where \(p\) is the pore pressure, \(p_c\) is the capillary pressure, which is \(-p\) under the single phase assumption, \(\beta\) is a composite coefficient by the liquid compressibility and solid compressibility, \(\alpha_B\) is the Biot's constant, \(\alpha_T^s\) is the linear thermal expansivity of solid, \(Q_H\) is the point source or sink term, \(H(S-1)\) is the Heaviside function, and \( \mathbf u\) is the displacement. While this process does not contain a fully mechanical coupling, simplfied expressions can be given to approximate the latter term under certain stress conditions. The liquid velocity \(\mathbf{v}^l\) is described by the Darcy's law as

\[ \mathbf{v}^l=-\frac{{\mathbf k} k_{ref}}{\mu} (\nabla p - \rho^l \mathbf g) \]

with \({\mathbf k}\) the intrinsic permeability, \(k_{ref}\) the relative permeability, \(\mathbf g\) the gravitational force.

Definition at line 82 of file ThermoRichardsFlowProcess.h.

#include <ThermoRichardsFlowProcess.h>

Inheritance diagram for ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess:

[legend]

Collaboration diagram for ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess:

[legend]

Public Member Functions
	ThermoRichardsFlowProcess (std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const &parameters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> &&process_variables, ThermoRichardsFlowProcessData &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)

ODESystem interface
bool	isLinear () const override

Public Member Functions inherited from ProcessLib::Process
	Process (std::string name_, MeshLib::Mesh &mesh, std::unique_ptr< AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const &parameters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> &&process_variables, SecondaryVariableCollection &&secondary_variables, const bool use_monolithic_scheme=true)

void	preTimestep (std::vector< GlobalVector * > const &x, const double t, const double delta_t, const int process_id)
	Preprocessing before starting assembly for new timestep. More...

void	postTimestep (std::vector< GlobalVector * > const &x, const double t, const double delta_t, int const process_id)
	Postprocessing after a complete timestep. More...

void	postNonLinearSolver (GlobalVector const &x, GlobalVector const &xdot, const double t, double const dt, int const process_id)

void	preIteration (const unsigned iter, GlobalVector const &x) final

void	computeSecondaryVariable (double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id)
	compute secondary variables for the coupled equations or for output. More...

NumLib::IterationResult	postIteration (GlobalVector const &x) final

void	initialize ()

void	setInitialConditions (std::vector< GlobalVector * > &process_solutions, std::vector< GlobalVector * > const &process_solutions_prev, double const t, int const process_id)

MathLib::MatrixSpecifications	getMatrixSpecifications (const int process_id) const override

void	setCoupledSolutionsForStaggeredScheme (CoupledSolutionsForStaggeredScheme *const coupled_solutions)

void	updateDeactivatedSubdomains (double const time, const int process_id)

bool	isMonolithicSchemeUsed () const

virtual void	setCoupledTermForTheStaggeredSchemeToLocalAssemblers (int const)

virtual void	extrapolateIntegrationPointValuesToNodes (const double, std::vector< GlobalVector * > const &, std::vector< GlobalVector * > &)

void	preAssemble (const double t, double const dt, GlobalVector const &x) final

void	assemble (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) final

void	assembleWithJacobian (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) final

std::vector< NumLib::IndexValueVector< GlobalIndexType > > const *	getKnownSolutions (double const t, GlobalVector const &x, int const process_id) const final

virtual NumLib::LocalToGlobalIndexMap const &	getDOFTable (const int) const

MeshLib::Mesh &	getMesh () const

std::vector< std::reference_wrapper< ProcessVariable > > const &	getProcessVariables (const int process_id) const

SecondaryVariableCollection const &	getSecondaryVariables () const

std::vector< std::unique_ptr< IntegrationPointWriter > > const *	getIntegrationPointWriter (MeshLib::Mesh const &mesh) const

virtual Eigen::Vector3d	getFlux (std::size_t, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &) const

virtual void	solveReactionEquation (std::vector< GlobalVector * > &, std::vector< GlobalVector * > const &, double const, double const, NumLib::EquationSystem &, int const)

Private Types
using	LocalAssemblerIF = LocalAssemblerInterface

Private Member Functions
void	initializeConcreteProcess (NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
	Process specific initialization called by initialize(). More...

void	setInitialConditionsConcreteProcess (std::vector< GlobalVector * > &x, double const t, int const) override

void	assembleConcreteProcess (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override

void	assembleWithJacobianConcreteProcess (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) override

void	postTimestepConcreteProcess (std::vector< GlobalVector * > const &x, double const t, double const dt, const int process_id) override

void	computeSecondaryVariableConcrete (double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id) override

std::vector< NumLib::LocalToGlobalIndexMap const * >	getDOFTables (const int number_of_processes) const

Private Attributes
std::vector< MeshLib::Node * >	_base_nodes

std::unique_ptr< MeshLib::MeshSubset const >	_mesh_subset_base_nodes

ThermoRichardsFlowProcessData	_process_data

std::vector< std::unique_ptr< LocalAssemblerIF > >	_local_assemblers

GlobalSparsityPattern	_sparsity_pattern_with_linear_element

MeshLib::PropertyVector< double > *	_heat_flux = nullptr

MeshLib::PropertyVector< double > *	_hydraulic_flow = nullptr

Additional Inherited Members
Public Types inherited from ProcessLib::Process
using	NonlinearSolver = NumLib::NonlinearSolverBase

using	TimeDiscretization = NumLib::TimeDiscretization

Public Attributes inherited from ProcessLib::Process
std::string const	name

Protected Member Functions inherited from ProcessLib::Process
NumLib::Extrapolator &	getExtrapolator () const

NumLib::LocalToGlobalIndexMap const &	getSingleComponentDOFTable () const

void	initializeProcessBoundaryConditionsAndSourceTerms (const NumLib::LocalToGlobalIndexMap &dof_table, const int process_id)

virtual void	constructDofTable ()

void	constructMonolithicProcessDofTable ()

void	constructDofTableOfSpecifiedProcessStaggeredScheme (const int specified_prosess_id)

virtual std::tuple< NumLib::LocalToGlobalIndexMap *, bool >	getDOFTableForExtrapolatorData () const

Protected Attributes inherited from ProcessLib::Process
MeshLib::Mesh &	_mesh

std::unique_ptr< MeshLib::MeshSubset const >	_mesh_subset_all_nodes

std::unique_ptr< NumLib::LocalToGlobalIndexMap >	_local_to_global_index_map

SecondaryVariableCollection	_secondary_variables

VectorMatrixAssembler	_global_assembler

const bool	_use_monolithic_scheme

CoupledSolutionsForStaggeredScheme *	_coupled_solutions

unsigned const	_integration_order

std::vector< std::unique_ptr< IntegrationPointWriter > >	_integration_point_writer

GlobalSparsityPattern	_sparsity_pattern

std::vector< std::vector< std::reference_wrapper< ProcessVariable > > >	_process_variables

std::vector< BoundaryConditionCollection >	_boundary_conditions

Member Typedef Documentation

◆ LocalAssemblerIF

using ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::LocalAssemblerIF = LocalAssemblerInterface

private

Definition at line 106 of file ThermoRichardsFlowProcess.h.

Constructor & Destructor Documentation

◆ ThermoRichardsFlowProcess()

ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::ThermoRichardsFlowProcess	(	std::string	name,
		MeshLib::Mesh &	mesh,
		std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&	jacobian_assembler,
		std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const &	parameters,
		unsigned const	integration_order,
		std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> &&	process_variables,
		ThermoRichardsFlowProcessData &&	process_data,
		SecondaryVariableCollection &&	secondary_variables,
		bool const	use_monolithic_scheme
	)

Definition at line 26 of file ThermoRichardsFlowProcess.cpp.

     : Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
               integration_order, std::move(process_variables),
               std::move(secondary_variables), use_monolithic_scheme),
       _process_data(std::move(process_data))
 {
     _heat_flux = MeshLib::getOrCreateMeshProperty<double>(
         mesh, "HeatFlux", MeshLib::MeshItemType::Node, 1);
  
     _hydraulic_flow = MeshLib::getOrCreateMeshProperty<double>(
         mesh, "HydraulicFlow", MeshLib::MeshItemType::Node, 1);
  
     // TODO (naumov) remove ip suffix. Probably needs modification of the mesh
     // properties, s.t. there is no "overlapping" with cell/point data.
     // See getOrCreateMeshProperty.
     _integration_point_writer.emplace_back(
         std::make_unique<IntegrationPointWriter>(
             "saturation_ip", 1 /*n components*/, integration_order,
             _local_assemblers, &LocalAssemblerIF::getSaturation));
  
     _integration_point_writer.emplace_back(
         std::make_unique<IntegrationPointWriter>(
             "porosity_ip", 1 /*n components*/, integration_order,
             _local_assemblers, &LocalAssemblerIF::getPorosity));
 }

References _heat_flux, _hydraulic_flow, ProcessLib::Process::_integration_point_writer, _local_assemblers, ProcessLib::ThermoRichardsFlow::LocalAssemblerInterface::getPorosity(), ProcessLib::ThermoRichardsFlow::LocalAssemblerInterface::getSaturation(), and MeshLib::Node.

Member Function Documentation

◆ assembleConcreteProcess()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::assembleConcreteProcess	(	const double	t,
		double const	dt,
		std::vector< GlobalVector * > const &	x,
		std::vector< GlobalVector * > const &	xdot,
		int const	process_id,
		GlobalMatrix &	M,
		GlobalMatrix &	K,
		GlobalVector &	b
	)

overrideprivatevirtual

Implements ProcessLib::Process.

Definition at line 181 of file ThermoRichardsFlowProcess.cpp.

 {
     DBUG("Assemble the equations for ThermoRichardsFlowProcess.");
  
     std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
         dof_table = {std::ref(*_local_to_global_index_map)};
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     // Call global assembler for each local assembly item.
     GlobalExecutor::executeSelectedMemberDereferenced(
         _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
         pv.getActiveElementIDs(), dof_table, t, dt, x, xdot, process_id, M, K,
         b);
 }

References ProcessLib::Process::_global_assembler, _local_assemblers, ProcessLib::Process::_local_to_global_index_map, ProcessLib::VectorMatrixAssembler::assemble(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::Process::getProcessVariables().

◆ assembleWithJacobianConcreteProcess()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::assembleWithJacobianConcreteProcess	(	const double	t,
		double const	dt,
		std::vector< GlobalVector * > const &	x,
		std::vector< GlobalVector * > const &	xdot,
		const double	dxdot_dx,
		const double	dx_dx,
		int const	process_id,
		GlobalMatrix &	M,
		GlobalMatrix &	K,
		GlobalVector &	b,
		GlobalMatrix &	Jac
	)

overrideprivatevirtual

Implements ProcessLib::Process.

Definition at line 199 of file ThermoRichardsFlowProcess.cpp.

 {
     std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
         dof_tables;
  
     DBUG(
         "Assemble the Jacobian of ThermoRichardsFlow for the monolithic "
         "scheme.");
     dof_tables.emplace_back(*_local_to_global_index_map);
  
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     GlobalExecutor::executeSelectedMemberDereferenced(
         _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
         _local_assemblers, pv.getActiveElementIDs(), dof_tables, t, dt, x, xdot,
         dxdot_dx, dx_dx, process_id, M, K, b, Jac);
  
     auto copyRhs = [&](int const variable_id, auto& output_vector)
     {
         transformVariableFromGlobalVector(b, variable_id, dof_tables[0],
                                           output_vector, std::negate<double>());
     };
  
     copyRhs(0, *_heat_flux);
     copyRhs(1, *_hydraulic_flow);
 }

References ProcessLib::Process::_global_assembler, _heat_flux, _hydraulic_flow, _local_assemblers, ProcessLib::Process::_local_to_global_index_map, ProcessLib::VectorMatrixAssembler::assembleWithJacobian(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), ProcessLib::Process::getProcessVariables(), and NumLib::transformVariableFromGlobalVector().

◆ computeSecondaryVariableConcrete()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::computeSecondaryVariableConcrete	(	double const	t,
		double const	dt,
		std::vector< GlobalVector * > const &	x,
		GlobalVector const &	x_dot,
		int const	process_id
	)

overrideprivatevirtual

Reimplemented from ProcessLib::Process.

Definition at line 249 of file ThermoRichardsFlowProcess.cpp.

 {
     if (process_id != 0)
     {
         return;
     }
     DBUG(
         "Compute the secondary variables for "
         "ThermoRichardsFlowProcess.");
     auto const dof_tables = getDOFTables(x.size());
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LocalAssemblerIF::computeSecondaryVariable, _local_assemblers,
         pv.getActiveElementIDs(), dof_tables, t, dt, x, x_dot, process_id);
 }

References _local_assemblers, ProcessLib::LocalAssemblerInterface::computeSecondaryVariable(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), getDOFTables(), and ProcessLib::Process::getProcessVariables().

◆ getDOFTables()

std::vector< NumLib::LocalToGlobalIndexMap const * > ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::getDOFTables ( const int number_of_processes ) const

private

Definition at line 269 of file ThermoRichardsFlowProcess.cpp.

 {
     std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
     dof_tables.reserve(number_of_processes);
     std::generate_n(std::back_inserter(dof_tables), number_of_processes,
                     [&]() { return _local_to_global_index_map.get(); });
     return dof_tables;
 }

References ProcessLib::Process::_local_to_global_index_map.

Referenced by computeSecondaryVariableConcrete(), and postTimestepConcreteProcess().

◆ initializeConcreteProcess()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::initializeConcreteProcess	(	NumLib::LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh,
		unsigned const	integration_order
	)

overrideprivatevirtual

Process specific initialization called by initialize().

Implements ProcessLib::Process.

Definition at line 62 of file ThermoRichardsFlowProcess.cpp.

 {
     using nlohmann::json;
  
     const int process_id = 0;
     const int variable_id = 0;
     ProcessLib::createLocalAssemblers<ThermoRichardsFlowLocalAssembler>(
         mesh.getDimension(), mesh.getElements(), dof_table,
         getProcessVariables(process_id)[variable_id]
             .get()
             .getShapeFunctionOrder(),
         _local_assemblers, mesh.isAxiallySymmetric(), integration_order,
         _process_data);
  
     auto add_secondary_variable = [&](std::string const& name,
                                       int const num_components,
                                       auto get_ip_values_function)
     {
         _secondary_variables.addSecondaryVariable(
             name,
             makeExtrapolator(num_components, getExtrapolator(),
                              _local_assemblers,
                              std::move(get_ip_values_function)));
     };
  
     add_secondary_variable("velocity", mesh.getDimension(),
                            &LocalAssemblerIF::getIntPtDarcyVelocity);
  
     add_secondary_variable("saturation", 1,
                            &LocalAssemblerIF::getIntPtSaturation);
  
     add_secondary_variable("porosity", 1, &LocalAssemblerIF::getIntPtPorosity);
  
     add_secondary_variable("dry_density_solid", 1,
                            &LocalAssemblerIF::getIntPtDryDensitySolid);
  
     _process_data.element_saturation = MeshLib::getOrCreateMeshProperty<double>(
         const_cast<MeshLib::Mesh&>(mesh), "saturation_avg",
         MeshLib::MeshItemType::Cell, 1);
  
     _process_data.element_porosity = MeshLib::getOrCreateMeshProperty<double>(
         const_cast<MeshLib::Mesh&>(mesh), "porosity_avg",
         MeshLib::MeshItemType::Cell, 1);
  
     // Set initial conditions for integration point data.
     for (auto const& ip_writer : _integration_point_writer)
     {
         // Find the mesh property with integration point writer's name.
         auto const& name = ip_writer->name();
         if (!mesh.getProperties().existsPropertyVector<double>(name))
         {
             continue;
         }
         auto const& mesh_property =
             *mesh.getProperties().template getPropertyVector<double>(name);
  
         // The mesh property must be defined on integration points.
         if (mesh_property.getMeshItemType() !=
             MeshLib::MeshItemType::IntegrationPoint)
         {
             continue;
         }
  
         auto const ip_meta_data = getIntegrationPointMetaData(mesh, name);
  
         // Check the number of components.
         if (ip_meta_data.n_components !=
             mesh_property.getNumberOfGlobalComponents())
         {
             OGS_FATAL(
                 "Different number of components in meta data ({:d}) than in "
                 "the integration point field data for '{:s}': {:d}.",
                 ip_meta_data.n_components, name,
                 mesh_property.getNumberOfGlobalComponents());
         }
  
         // Now we have a properly named vtk's field data array and the
         // corresponding meta data.
         std::size_t position = 0;
         for (auto& local_asm : _local_assemblers)
         {
             std::size_t const integration_points_read =
                 local_asm->setIPDataInitialConditions(
                     name, &mesh_property[position],
                     ip_meta_data.integration_order);
             if (integration_points_read == 0)
             {
                 OGS_FATAL(
                     "No integration points read in the integration point "
                     "initial conditions set function.");
             }
             position += integration_points_read * ip_meta_data.n_components;
         }
     }
  
     // Initialize local assemblers after all variables have been set.
     GlobalExecutor::executeMemberOnDereferenced(&LocalAssemblerIF::initialize,
                                                 _local_assemblers,
                                                 *_local_to_global_index_map);
 }

◆ isLinear()

bool ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::isLinear ( ) const

inlineoverride

Definition at line 102 of file ThermoRichardsFlowProcess.h.

102 { return false; }

◆ postTimestepConcreteProcess()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::postTimestepConcreteProcess	(	std::vector< GlobalVector * > const &	x,
		double const	t,
		double const	dt,
		const int	process_id
	)

overrideprivatevirtual

Reimplemented from ProcessLib::Process.

Definition at line 230 of file ThermoRichardsFlowProcess.cpp.

 {
     if (process_id != 0)
     {
         return;
     }
  
     DBUG("PostTimestep ThermoRichardsFlowProcess.");
  
     auto const dof_tables = getDOFTables(x.size());
  
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LocalAssemblerIF::postTimestep, _local_assemblers,
         pv.getActiveElementIDs(), dof_tables, x, t, dt);
 }

References _local_assemblers, DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), getDOFTables(), ProcessLib::Process::getProcessVariables(), and ProcessLib::LocalAssemblerInterface::postTimestep().

◆ setInitialConditionsConcreteProcess()

void ProcessLib::ThermoRichardsFlow::ThermoRichardsFlowProcess::setInitialConditionsConcreteProcess	(	std::vector< GlobalVector * > &	x,
		double const	t,
		int const	process_id
	)

overrideprivatevirtual

Reimplemented from ProcessLib::Process.

Definition at line 166 of file ThermoRichardsFlowProcess.cpp.

 {
     if (process_id != 0)
     {
         return;
     }
     DBUG("SetInitialConditions ThermoRichardsFlowProcess.");
  
     GlobalExecutor::executeMemberOnDereferenced(
         &LocalAssemblerIF::setInitialConditions, _local_assemblers,
         *_local_to_global_index_map, *x[process_id], t, _use_monolithic_scheme,
         process_id);
 }