Dear Lorenzo,
this looks reasonable at first glance, however, the TwoPImmiscible fluid
system can only handle two components in total - your code will probably
not compile due to a static_assert.
Do you consider a water/polymer mixture at a fixed concentration? If
yes, then you could create some sort of pseudo-component (such as "Air",
which is actually a mixture of multiple components) which you could then
just plug into OnePLiquid.
If you actually want to balance water and the polymer explicitly, you
will probably have to write your own fluid system
where you can then also specify your phase and component indices.
Best regards
Kilian
On 01/27/2020 10:54 AM, lc wrote:
Hello,
I refine a bit my question.
Is it corrent, respect to the definition of fluid system that I need,
to write:
// In the following we define our fluid properties.
template<class TypeTag>
struct FluidSystem<TypeTag, TTag::FiveSpotProblemTypeTag>
{
// We define a convenient shortcut to the property Scalar:
using Scalar = GetPropType<TypeTag, Properties::Scalar>;
// First, we create a fluid system that consists of one liquid
water phase. We use the simple
// description of water, which means we do not use tabulated
values but more general equations of state.
using WettingPhase = FluidSystems::LiquidPhaseTwoC<Scalar,
myH2O<Scalar>, myPolymer<Scalar> >;
// Second, we create another fluid system consisting of a liquid
phase as well, the Trichlorethene (DNAPL) phase:
using NonwettingPhase = FluidSystems::OnePLiquid<Scalar,
myOil<Scalar> >;
// Third, we combine both fluid systems in our final fluid system
which consist of two immiscible liquid phases:
using type = FluidSystems::TwoPImmiscible<Scalar, WettingPhase,
NonwettingPhase>;
};
If, yes, then, I'm a bit confused on how to define the indexes for
phases and components.
Best regards,
Lorenzo
On 27.01.2020 12:35, lc wrote:
Dear Dumux-ers,
I'd like to ask, how to create a fluid system with 2 phases: water
and oil and 2 components, water and "polymer" only in the water phase?
The definietion of the "polymer" component is not important for the
moment, it will be obtained overloading some other component and
changing the property values.
Thank you,
Kind regards,
Lorenzo Campoli
--
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Kilian Weishaupt M.Sc.
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
