[ysjeong]

Hi everyone,

I'd like to ask all you guys in OpenFOAM4.1 new solver development, combining cht solver (chtMultiRegionFoam) and VOF solver (interPhaseChangeFoam).

For alpha 1 and alpha2 declared in createFields.H of interPhaseChangeFoam,
autoPtr<phaseChangeTwoPhaseMixture> mixture =
phaseChangeTwoPhaseMixture::New(U, phi);
volScalarField& alpha1(mixture->alpha1());
volScalarField& alpha2(mixture->alpha2());

I want to declare the PtrList for chtMultiRegionFoam as;
PtrList<autoPtr<phaseChangeTwoPhaseMixture>> mixtureFluid(fluidRegions.size());
PtrList<volScalarField> alpha1Fluid(fluidRegions.size());
PtrList<volScalarField> alpha2Fluid(fluidRegions.size());

then,
autoPtr<phaseChangeTwoPhaseMixture> mixtureFluid[i]= myphaseChangeTwoPhaseMixture::New(UFluid[i],phiFluid[i]);

alpha1Fluid.set
(
i,
new volScalarField
(
IOobject
(
"alpha1",
runTime.timeName(),
fluidRegions[i],
IOobject::READ_IF_PRESENT
),
mixtureFluid[i]->alpha1()
)
);

At this moment, compiling is perfectly working but when I want to run the case, I got the error message at the early stage of running as;

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::Boundary(Foam:imension edField<double, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary const&) at ??:?
#4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 ? at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? at ??:?



Does anyone have some ideas on solving this problem?
And, is there any differences in declaring the volScalarField and volScalarField& for pointer list form?

Any comments are welcome!!