fixed writing out entries in advective bc
[OpenFOAM-1.6-ext.git] / src / lagrangian / dieselSpray / spray / sprayI.H
blob46fd38f989ad060c4d3e9239eaf482586c98b46d
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4    \\    /   O peration     |
5     \\  /    A nd           | Copyright held by original author
6      \\/     M anipulation  |
7 -------------------------------------------------------------------------------
8 License
9     This file is part of OpenFOAM.
11     OpenFOAM is free software; you can redistribute it and/or modify it
12     under the terms of the GNU General Public License as published by the
13     Free Software Foundation; either version 2 of the License, or (at your
14     option) any later version.
16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
19     for more details.
21     You should have received a copy of the GNU General Public License
22     along with OpenFOAM; if not, write to the Free Software Foundation,
23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 \*---------------------------------------------------------------------------*/
27 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
29 namespace Foam
32 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
34 inline const Time& spray::runTime() const
36     return runTime_;
40 inline const fvMesh& spray::mesh() const
42     return mesh_;
46 inline const volVectorField& spray::U() const
48     return U_;
52 inline const volScalarField& spray::rho() const
54     return rho_;
58 inline const volScalarField& spray::p() const
60     return p_;
64 inline const volScalarField& spray::T() const
66     return T_;
70 inline PtrList<injector>& spray::injectors()
72     return injectors_;
76 inline const PtrList<injector>& spray::injectors() const
78     return injectors_;
82 inline const atomizationModel& spray::atomization() const
84     return atomization_;
88 inline const breakupModel& spray::breakup() const
90     return breakupModel_;
94 inline const collisionModel& spray::collisions() const
96     return collisionModel_;
100 inline const dispersionModel& spray::dispersion() const
102     return dispersionModel_;
106 inline const dragModel& spray::drag() const
108     return drag_;
112 inline const evaporationModel& spray::evaporation() const
114     return evaporation_;
118 inline const heatTransferModel& spray::heatTransfer() const
120     return heatTransfer_;
124 inline const injectorModel& spray::injection() const
126     return injectorModel_;
130 inline const wallModel& spray::wall() const
132     return wall_;
136 inline tmp<volVectorField> spray::momentumSource() const
138     tmp<volVectorField> tsource
139     (
140         new volVectorField
141         (
142             IOobject
143             (
144                 "sms",
145                 runTime_.timeName(),
146                 mesh_,
147                 IOobject::NO_READ,
148                 IOobject::NO_WRITE
149             ),
150             mesh_,
151             dimensionedVector
152             (
153                 "zero",
154                 dimensionSet(1, -2, -2, 0, 0),
155                 vector::zero
156             )
157         )
158     );
160     tsource().internalField() = sms_/runTime_.deltaT().value()/mesh_.V();
162     return tsource;
166 inline tmp<volScalarField> spray::evaporationSource(const label si) const
168     tmp<volScalarField> tsource
169     (
170         new volScalarField
171         (
172             IOobject
173             (
174                 "srhos",
175                 runTime_.timeName(),
176                 mesh_,
177                 IOobject::NO_READ,
178                 IOobject::NO_WRITE
179             ),
180             mesh_,
181             dimensionedScalar("zero", dimensionSet(1, -3, -1, 0, 0), 0.0)
182         )
183     );
185     if (isLiquidFuel_[si])
186     {
187         label fi = gasToLiquidIndex_[si];
188         tsource().internalField() =
189             srhos_[fi]/runTime_.deltaT().value()/mesh_.V();
190     }
191     else
192     {
193         scalarField s(mesh_.nCells(), 0.0);
194         tsource().internalField() = s;
195     }
197     return tsource;
201 inline tmp<volScalarField> spray::heatTransferSource() const
203     tmp<volScalarField> tsource
204     (
205         new volScalarField
206         (
207             IOobject
208             (
209                 "shs",
210                 runTime_.timeName(),
211                 mesh_,
212                 IOobject::NO_READ,
213                 IOobject::NO_WRITE
214             ),
215             mesh_,
216             dimensionedScalar("zero", dimensionSet(1, -1, -3, 0, 0), 0.0)
217         )
218     );
220     tsource().internalField() = shs_/runTime_.deltaT().value()/mesh_.V();
222     return tsource;
226 inline Random& spray::rndGen()
228     return rndGen_;
232 inline label spray::subCycles() const
234     return subCycles_;
238 inline label spray::sprayIterate() const
240     return sprayIterate_;
244 inline scalar spray::sprayRelaxFactor() const
246     return sprayRelaxFactor_;
250 inline scalar spray::minimumParcelMass() const
252     return minimumParcelMass_;
256 inline const vector& spray::g() const
258     return g_;
262 inline const liquidMixture& spray::fuels() const
264     return fuels_;
268 inline const PtrList<gasThermoPhysics>& spray::gasProperties() const
270     return gasProperties_;
274 inline const basicMultiComponentMixture& spray::composition() const
276     return composition_;
280 inline const List<label>& spray::liquidToGasIndex() const
282     return liquidToGasIndex_;
286 inline const List<label>& spray::gasToLiquidIndex() const
288     return gasToLiquidIndex_;
292 inline const List<bool>& spray::isLiquidFuel() const
294     return isLiquidFuel_;
298 inline const bool& spray::twoD() const
300     return twoD_;
304 inline const vector& spray::axisOfSymmetry() const
306     return axisOfSymmetry_;
310 inline const vector& spray::axisOfWedge() const
312     return axisOfWedge_;
316 inline const vector& spray::axisOfWedgeNormal() const
318     return axisOfWedgeNormal_;
322 inline const scalar& spray::angleOfWedge() const
324     return angleOfWedge_;
328 inline const interpolation<vector>& spray::UInterpolator() const
330     return UInterpolator_;
334 inline const interpolation<scalar>& spray::rhoInterpolator() const
336     return rhoInterpolator_;
340 inline const interpolation<scalar>& spray::pInterpolator() const
342     return pInterpolator_;
346 inline const interpolation<scalar>& spray::TInterpolator() const
348     return TInterpolator_;
352 inline vectorField& spray::sms()
354     return sms_;
358 inline const vectorField& spray::sms() const
360     return sms_;
364 inline scalarField& spray::shs()
366     return shs_;
370 inline const scalarField& spray::shs() const
372     return shs_;
376 inline PtrList<scalarField>& spray::srhos()
378     return srhos_;
382 inline const PtrList<scalarField>& spray::srhos() const
384     return srhos_;
388 inline const scalar& spray::ambientPressure() const
390     return ambientPressure_;
394 inline const scalar& spray::ambientTemperature() const
396     return ambientTemperature_;
400 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
402 } // End namespace Foam
404 // ************************************************************************* //