InterfaceCompositionModels.C 12.6 KB
Newer Older
1 2 3 4
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
OpenFOAM bot's avatar
OpenFOAM bot committed
5
    \\  /    A nd           | www.openfoam.com
6
     \\/     M anipulation  |
OpenFOAM bot's avatar
OpenFOAM bot committed
7 8
-------------------------------------------------------------------------------
    Copyright (C) 2017 OpenCFD Ltd.
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

\*---------------------------------------------------------------------------*/

#include "addToRunTimeSelectionTable.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "thermoPhysicsTypes.H"

#include "rhoConst.H"
#include "perfectFluid.H"
#include "Boussinesq.H"

#include "pureMixture.H"
#include "multiComponentMixture.H"
#include "reactingMixture.H"
#include "SpecieMixture.H"

#include "rhoThermo.H"
#include "rhoReactionThermo.H"
#include "heRhoThermo.H"

#include "solidThermo.H"
#include "heSolidThermo.H"
#include "solidThermoPhysicsTypes.H"

#include "kineticGasEvaporation.H"
52
#include "Lee.H"
53
#include "interfaceHeatResistance.H"
54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{
    typedef
        constTransport
        <
            species::thermo
            <
                hConstThermo
                <
                    rhoConst<specie>
                >,
                sensibleEnthalpy
            >
70
        > constRhoHThermoPhysics;
71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93


    typedef
        constTransport
        <
            species::thermo
            <
                hConstThermo
                <
                    Boussinesq<specie>
                >,
                sensibleEnthalpy
            >
        > BoussinesqFluidEThermoPhysics;
};


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{
    using namespace meltingEvaporationModels;

94
    //NOTE: First thermo (from) and second otherThermo (to)
95 96

    // kineticGasEvaporation model definitions
97 98 99

        // From pure liquid (rhoConst) to a multi-component gas incomp phase
        makeInterfaceContSpecieMixtureType
100 101 102 103 104
        (
            kineticGasEvaporation,
            heRhoThermo,
            rhoThermo,
            pureMixture,
105 106 107 108 109
            constRhoHThermoPhysics,
            heRhoThermo,
            rhoReactionThermo,
            multiComponentMixture,
            constIncompressibleGasHThermoPhysics
110 111
        );

Andrew Heather's avatar
Andrew Heather committed
112 113
        // From pure liquid (BoussinesqFluid) to a multi-component gas incomp
        // phase
114 115 116 117 118 119
        makeInterfaceContSpecieMixtureType
        (
            kineticGasEvaporation,
            heRhoThermo,
            rhoThermo,
            pureMixture,
120
            BoussinesqFluidEThermoPhysics,
121 122 123 124 125 126 127
            heRhoThermo,
            rhoReactionThermo,
            multiComponentMixture,
            constIncompressibleGasHThermoPhysics
        );


128
        // From pure liquid (rhoConst) to pure gas (incompressible ideal gas)
129 130 131 132 133 134
        makeInterfacePureType
        (
            kineticGasEvaporation,
            heRhoThermo,
            rhoThermo,
            pureMixture,
135
            constRhoHThermoPhysics,
136 137 138
            heRhoThermo,
            rhoThermo,
            pureMixture,
139
            constIncompressibleGasHThermoPhysics
140 141
        );

142
        // From pure liquid (const rho) to pure gas (rhoConst) gas
143 144 145 146 147 148
        makeInterfacePureType
        (
            kineticGasEvaporation,
            heRhoThermo,
            rhoThermo,
            pureMixture,
149
            constRhoHThermoPhysics,
150 151 152
            heRhoThermo,
            rhoThermo,
            pureMixture,
153
            constRhoHThermoPhysics
154 155 156
        );


157
        // From pure liquid (Boussinesq) to pure gas (incompressible ideal gas)
158 159
        makeInterfacePureType
        (
160
            kineticGasEvaporation,
161 162 163
            heRhoThermo,
            rhoThermo,
            pureMixture,
164
            BoussinesqFluidEThermoPhysics,
165 166 167
            heRhoThermo,
            rhoThermo,
            pureMixture,
168
            constIncompressibleGasHThermoPhysics
169 170
        );

171
        // From pure liquid (Boussinesq) to pure gas (rho const)
172 173
        makeInterfacePureType
        (
174
            kineticGasEvaporation,
175 176 177
            heRhoThermo,
            rhoThermo,
            pureMixture,
178
            BoussinesqFluidEThermoPhysics,
179 180 181
            heRhoThermo,
            rhoThermo,
            pureMixture,
182
            constRhoHThermoPhysics
183 184 185
        );


186
    // Lee model definitions
187

188
        // From pure phase (rho const) to phase (rho const)
189 190
        makeInterfacePureType
        (
191
            Lee,
192 193 194
            heRhoThermo,
            rhoThermo,
            pureMixture,
195
            constRhoHThermoPhysics,
196 197 198
            heRhoThermo,
            rhoThermo,
            pureMixture,
199
            constRhoHThermoPhysics
200 201
        );

202 203
        // From pure phase (rho const) to phase (Boussinesq)
        makeInterfacePureType
204
        (
205
            Lee,
206 207 208
            heRhoThermo,
            rhoThermo,
            pureMixture,
209
            constRhoHThermoPhysics,
210 211 212
            heRhoThermo,
            rhoThermo,
            pureMixture,
213
            BoussinesqFluidEThermoPhysics
214 215
        );

216 217

        // From pure phase (solidThermo) to phase (Boussinesq)
218 219
        makeInterfacePureType
        (
220 221 222
            Lee,
            heSolidThermo,
            solidThermo,
223
            pureMixture,
224
            hConstSolidThermoPhysics,
225 226 227
            heRhoThermo,
            rhoThermo,
            pureMixture,
228
            BoussinesqFluidEThermoPhysics
229 230
        );

231
        // From pure phase (solidThermo) to phase (rho const)
232 233
        makeInterfacePureType
        (
234
            Lee,
235 236 237
            heSolidThermo,
            solidThermo,
            pureMixture,
238
            hConstSolidThermoPhysics,
239 240 241
            heRhoThermo,
            rhoThermo,
            pureMixture,
242
            constRhoHThermoPhysics
243 244
        );

245
        // From pure phase (all-poly solidThermo) to phase (ico-rho)
246 247 248 249 250 251 252 253 254 255 256 257 258
        makeInterfacePureType
        (
            Lee,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hPolyTranspPolyIcoSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            icoPoly8HThermoPhysics
        );

259 260 261 262 263 264 265 266 267 268 269 270 271 272
        // From pure phase (exp-Transp, hPower solidThermo) to phase (ico-rho)
        makeInterfacePureType
        (
            Lee,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hPowerSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            icoPoly8HThermoPhysics
        );

273 274

        // From pure phase (const rho) to multi phase (incomp ideal gas)
275 276
        makeInterfaceContSpecieMixtureType
        (
277
            Lee,
278 279 280
            heRhoThermo,
            rhoThermo,
            pureMixture,
281
            constRhoHThermoPhysics,
282 283 284 285 286 287
            heRhoThermo,
            rhoReactionThermo,
            multiComponentMixture,
            constIncompressibleGasHThermoPhysics
        );

288 289 290

        // From pure phase (Boussinesq) to phase (solidThermo)
        makeInterfacePureType
291
        (
292
            Lee,
293
            heRhoThermo,
294 295 296 297 298 299 300
            rhoThermo,
            pureMixture,
            BoussinesqFluidEThermoPhysics,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hConstSolidThermoPhysics
301 302
        );

303
        // From pure phase (rho const) to phase (solidThermo)
304 305
        makeInterfacePureType
        (
306
            Lee,
307 308 309
            heRhoThermo,
            rhoThermo,
            pureMixture,
310 311 312
            constRhoHThermoPhysics,
            heSolidThermo,
            solidThermo,
313
            pureMixture,
314
            hConstSolidThermoPhysics
315
        );
316

317 318 319 320 321 322 323
        //From pure liquid phase (ico-rho) to pure phase (exp-Transp, hPower solidThermo)
        makeInterfacePureType
        (
            Lee,
            heRhoThermo,
            rhoThermo,
            pureMixture,
324 325 326 327
            icoPoly8HThermoPhysics,
            heSolidThermo,
            solidThermo,
            pureMixture,
328 329 330
            hPowerSolidThermoPhysics
        );

331

332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477

        // interfaceHeatResistance model definitions

        // From pure phase (rho const) to phase (rho const)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics
        );

        // From pure phase (rho const) to phase (Boussinesq)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            BoussinesqFluidEThermoPhysics
        );


        // From pure phase (solidThermo) to phase (Boussinesq)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hConstSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            BoussinesqFluidEThermoPhysics
        );

        // From pure phase (solidThermo) to phase (rho const)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hConstSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics
        );

        // From pure phase (all-poly solidThermo) to phase (ico-rho)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hPolyTranspPolyIcoSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            icoPoly8HThermoPhysics
        );

        // From pure phase (exp-Transp, hPower solidThermo) to phase (ico-rho)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hPowerSolidThermoPhysics,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            icoPoly8HThermoPhysics
        );


        // From pure phase (const rho) to multi phase (incomp ideal gas)
        makeInterfaceContSpecieMixtureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics,
            heRhoThermo,
            rhoReactionThermo,
            multiComponentMixture,
            constIncompressibleGasHThermoPhysics
        );


        // From pure phase (Boussinesq) to phase (solidThermo)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            BoussinesqFluidEThermoPhysics,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hConstSolidThermoPhysics
        );

        // From pure phase (rho const) to phase (solidThermo)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            constRhoHThermoPhysics,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hConstSolidThermoPhysics
        );

        //From pure liquid phase (ico-rho) to pure phase (exp-Transp, hPower solidThermo)
        makeInterfacePureType
        (
            interfaceHeatResistance,
            heRhoThermo,
            rhoThermo,
            pureMixture,
            icoPoly8HThermoPhysics,
            heSolidThermo,
            solidThermo,
            pureMixture,
            hPowerSolidThermoPhysics
        );

478 479
}

Andrew Heather's avatar
Andrew Heather committed
480

481
// ************************************************************************* //