Dear Etienne,
my guess is that ElementVolumeVariables is not default-constructible.
You therefore need to initialize it in the constructor of
fvlocalassemblerbase.hh as done for the other elementVolVars.
If that is, for any reason, not possible consider using
std::optional<ElementVolumeVariables>.
Best wishes
Kilian
On 10.03.21 10:05, Etienne Ahusborde wrote:
Dear Dumux Users,
With my phD student we are interested in the numerical simulation of
convective dissolution of CO2.
We would like to implement a second order scheme for the time derivative :
$$\frac{\partial Storage}{\partial t}\approx \frac{3
Storage^{n+1}-4Storage^n + Storage^{n-1}}{2\Delta t}$$
We identified that we have to modify (at least) the files
fvlocalassemblerbase.hh and fvlocalresidual.hh.
We have to create a variable "prevprevElemVolVars_" to store
$Storage^{n-1}$ in fvlocalassemblerbase.hh
When we just try ro declare ElementVolumeVariables prevElemVolVars_
in fvlocalassemblerbase.hh as a private variable (similarly to
ElementVolumeVariables curElemVolVars_ and ElementVolumeVariables
prevElemVolVars_) without any use of this variable, we obtain the
following error during the compilation :
/home/eahusbor/DUMUX_PYRENE/DUMUX_3.1/dumux/dumux/assembly/fvlocalassemblerbase.hh:107:37:
error: no matching function for call to
‘Dumux::CCTpfaElementVolumeVariables<Dumux::CCGridVolumeVariables<Dumux::CCTpfaDefaultGridVolumeVariablesTraits<Dumux::HeterogeneousProblem<Dumux::Properties::TTag::HeterogeneousCCTpfa>,
Dumux::TwoPTwoCVolumeVariablesImplementation<Dumux::TwoPNCVolumeVariablesTraits<Dumux::SwitchablePrimaryVariables<Dune::FieldVector<double,
5>, int>, Dumux::FluidSystems::BrineCO2<double,
Dumux::HeterogeneousCO2Tables::CO2Tables,
Dumux::Components::TabulatedComponent<Dumux::Components::H2O<double>,
true>, Dumux::FluidSystems::BrineCO2DefaultPolicy<true, true> >,
Dumux::CompositionalFluidState<double,
Dumux::FluidSystems::BrineCO2<double,
Dumux::HeterogeneousCO2Tables::CO2Tables,
Dumux::Components::TabulatedComponent<Dumux::Components::H2O<double>,
true>, Dumux::FluidSystems::BrineCO2DefaultPolicy<true, true> > >,
Dumux::SolidSystems::InertSolidPhase<double,
Dumux::Components::Constant<1, double> >,
Dumux::InertSolidState<double,
Dumux::SolidSystems::InertSolidPhase<double,
Dumux::Components::Constant<1, double> > >, double,
Dumux::TwoPNCModelTraits<4, true, true, Dumux::TwoPFormulation::p0s1,
4> >, false, true> >, false>, false>::CCTpfaElementVolumeVariables()’
107 | , elementIsGhost_(elementIsGhost)
We don't understand this issue.
Can somebody help us?
Thanks in advance
Regards
Etienne
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Kilian Weishaupt
Institut für Wasser- und Umweltsystemmodellierung (IWS)
Lehrstuhl für Hydromechanik und Hydrosystemmodellierung
Universität Stuttgart, Pfaffenwaldring 61, 70569 Stuttgart
Email: [email protected]
Telefon: 0049 711 685-60461 ** fax: 0049-711-685-60430
http://www.hydrosys.uni-stuttgart.de
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
