|
[Sponsors] |
November 6, 2013, 05:50 |
Read out value of keyword in boundaryField
|
#1 |
Member
Join Date: May 2012
Posts: 55
Rep Power: 14 |
Hello foamers,
I've added a new boundary condition (modified constantAlphaContactAngle) for calculating a temperature dependent contact angle. This BC writes out the current contact angle in the write function via: Code:
os.writeKeyword("theta0") << calcMic(0) << token::END_STATEMENT << nl; Code:
const volScalarField::GeometricBoundaryField& abf = alpha1_.boundaryField(); const fvBoundaryMesh& boundary = mesh.boundary(); forAll(boundary, patchi) { if (isA<alphaContactAngleFvPatchScalarField>(abf[patchi])) { alphaContactAngleFvPatchScalarField& acap = const_cast<alphaContactAngleFvPatchScalarField&> ( refCast<const alphaContactAngleFvPatchScalarField> ( abf[patchi] ) ); fvsPatchVectorField& nHatp = nHatfv.boundaryField()[patchi]; const scalarField theta ( convertToRad*acap.theta(U_.boundaryField()[patchi], nHatp) ); Info << "######## acap ======" << acap << endl; } } My Question: How can I get access to the value of the keyword theta0? thank you! |
|
November 11, 2013, 09:44 |
|
#2 |
Member
Join Date: May 2012
Posts: 55
Rep Power: 14 |
Does anybody have an idea?
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Commercial meshers] fluentMeshToFoam multidomain mesh conversion problem | Attesz | OpenFOAM Meshing & Mesh Conversion | 12 | May 2, 2013 10:52 |
999999 (../../src/mpsystem.c@1123):mpt_read: failed:errno = 11 | UDS_rambler | FLUENT | 2 | November 22, 2011 09:46 |
Problem in running ICEM grid in Openfoam | Tarak | OpenFOAM | 6 | September 9, 2011 17:51 |
[Gmsh] GmshToFoam keyword patch0 is undefined | steve999 | OpenFOAM Meshing & Mesh Conversion | 5 | September 14, 2008 14:45 |
Phase locked average in run time | panara | OpenFOAM | 2 | February 20, 2008 14:37 |