---

Hallo miteinander,

ich bin dabei die Krümmung in interfaceProperties.C (v6) anders zu berechnen, als Alternative zum Standard CSF-Modell (man soll aber mit einem Dict-File zwischen den Modellen switchen können).

Mein Problem: Die "Sharp Surface Force" Methode nach Raeini et al. ermittelt am Ende eine Krümmung an den Seitenflächen (also vom Typus surfaceScalarField statt volScalarField), nennen wir sie Kf_, aber folgende Verrechnung mit sigma erfordert regulär K_ als volScalarField:

Code:

---

Foam::tmp<Foam::volScalarField>
Foam::interfaceProperties::sigmaK() const
{
    return sigmaPtr_->sigma()*K_;
}

---

Ich habe das hier versucht:

Code:

---

Foam::tmp<Foam::volScalarField>
Foam::interfaceProperties::sigmaK() const
{
    return sigmaPtr_->sigma();
}

---

...und die Oberflächenkraftsberechnung so umgeschrieben:

Code:

---

Foam::tmp<Foam::surfaceScalarField>
Foam::interfaceProperties::surfaceTensionForce() const
{
    return fvc::interpolate(sigmaK())*Kf_*fvc::snGrad(alpha1_);
}

---

...aber es kompiliert nicht, sigma alleine interpoliert macht hier Probleme.

Meine Frage wäre: Wie kann ich effektiv eine Krümmung an der Seitenfläche ermitteln und es mit sigma verrechnen? C++ habe ich noch nicht komplett "gemeistert", womit ich sicherlich einige entscheidende Aspekt hier nicht berücksichtigt habe.

Danke im voraus!

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hey,

das ist ganz klar wieso das nicht geht, da Kf_ ein (wahrscheinlich) privates Objekt der Klasse interfaceProperties ist. Du hast damit keinen Zugriff in der Klasse surfaceTensionForce.

Ich möchte jetzt keine Aussage zu Deinem Wissensstand zu c++ nennen aber das sei mal dahingestellt.

Da Du auch keine Fehlermeldung vom Kompiler lieferst, gehe ich mal (ohne den Code jetzt zu prüfen - dafür hab ich grad keine Zeit) davon aus, dass Du das von mir beschrieben Problem hast.

Demnach musst Du eine Public Funktion in interfaceProperties erstellen, die den Wert von Kf_ zurückgibt.

Quasi so in der Art:

interfaceProperties.H

Code:

---

    public:
  
            const volScalarField& Kf();

---

interfaceProperties.C

Code:

---

const volScalarField& Foam::interfaceProperties::Kf()
{
    return Kf_;
}

---

Vllt sind die Namespaces jetzt nicht ganz richtig aber das siehst Du ja dann.
Im Anschluss, sofern surfaceTensionForce einen Bezug zur interfaceProperties hat, kannst Du folgendes machen:

Code:

---

Foam::tmp<Foam::surfaceScalarField>
Foam::interfaceProperties::surfaceTensionForce() const
{
    return fvc::interpolate(sigmaK())*fvc::interpolate(Kf())*fvc::snGrad(alpha1_);
}

---

Anmerkung:

• Du hast einmal K_ und dann wieder Kf_, Fehler oder gewollt? Ich geh gerade  davon aus das es beidesmal Kf_ ist.
  
• Meine Analyse ist nur Oberflächlich ohne das ich irgendwas von dem verstehe was Du machen willst
  
• Codeschnippsel sind ggf. falsch
  
• Da ich nicht weiß ob Kf_ und K_ gleich sind und welchen Typ diese darstellen, musst Du natürlich sicherstellen das Du nicht Gummibärchen mit Plastikspielzeug vertauscht. Du kannst keine Multiplication von surfaceScalarField*volScalarField machen. Das ist einfach nonsense.
  

------------------
Glück Auf,
Tobi

OpenFOAM® Community - Knowledge Base

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP