/// <summary> /// default plotting /// </summary> protected override void PlotCurrentState(double phystime, TimestepNumber timestepNo, int superSampling = 0) { string caseStr = ""; DGField[] Fields = new DGField[0]; SinglePhaseField phiDist = ComputeDistanceField(); Fields = Fields.Cat(this.phi, this.mu, this.Velocity, this.gradPhi0, this.Curvature, this.DCurvature, phiDist); BoSSS.Solution.Tecplot.Tecplot.PlotFields(Fields, "Phasefield-" + timestepNo + caseStr, phystime, superSampling); }
//protected override IEnumerable<DGField> InstantiateParameterFields() { // DGField[] ParameterFields = Velocity.ToArray().Cat(phi0); // return ParameterFields; //} public override IEnumerable <DGField> InstantiateResidualFields() { DGField[] ResidualFields = new DGField[] { phi_Resi }; switch (this.Control.ModTyp) { case PhasefieldControl.ModelType.modelB: ResidualFields = ResidualFields.Cat(mu_Resi); break; case PhasefieldControl.ModelType.modelA: case PhasefieldControl.ModelType.modelC: default: break; } if (this.Control.CurvatureCorrectionType != PhasefieldControl.CurvatureCorrection.None) { ResidualFields.Cat(curvature_Resi); } return(ResidualFields); }
protected override IEnumerable <DGField> InstantiateSolutionFields() { DGField[] SolutionFields = new DGField[] { phi }; switch (this.Control.ModTyp) { case PhasefieldControl.ModelType.modelB: SolutionFields = SolutionFields.Cat(mu); break; case PhasefieldControl.ModelType.modelA: case PhasefieldControl.ModelType.modelC: default: break; } if (this.Control.CurvatureCorrectionType != PhasefieldControl.CurvatureCorrection.None) { SolutionFields.Cat(Curvature); } return(SolutionFields); }