changed to better initialisation of table

This commit is contained in:
andy
2009-03-26 12:08:45 +00:00
parent d504063ce1
commit ab2fe81805
2 changed files with 109 additions and 121 deletions

View File

@ -29,119 +29,110 @@ License
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
const Foam::scalar
Foam::radiation::blackBodyEmission::
emissivePowerTable[Foam::radiation::blackBodyEmission::nTableEntries][2] =
{
{1000.0, 0.00},
{1100.0, 0.00},
{1200.0, 0.00},
{1300.0, 0.00},
{1400.0, 0.01},
{1500.0, 0.01},
{1600.0, 0.02},
{1700.0, 0.03},
{1800.0, 0.04},
{1900.0, 0.05},
{2000.0, 0.07},
{2100.0, 0.08},
{2200.0, 0.10},
{2300.0, 0.12},
{2400.0, 0.14},
{2500.0, 0.16},
{2600.0, 0.18},
{2700.0, 0.21},
{2800.0, 0.23},
{2900.0, 0.25},
{3000.0, 0.27},
{3100.0, 0.30},
{3200.0, 0.32},
{3300.0, 0.34},
{3400.0, 0.36},
{3500.0, 0.38},
{3600.0, 0.40},
{3700.0, 0.42},
{3800.0, 0.44},
{3900.0, 0.46},
{4000.0, 0.48},
{4100.0, 0.50},
{4200.0, 0.52},
{4300.0, 0.53},
{4400.0, 0.55},
{4500.0, 0.56},
{4600.0, 0.58},
{4700.0, 0.59},
{4800.0, 0.61},
{4900.0, 0.62},
{5000.0, 0.63},
{5100.0, 0.65},
{5200.0, 0.66},
{5300.0, 0.67},
{5400.0, 0.68},
{5500.0, 0.69},
{5600.0, 0.70},
{5700.0, 0.71},
{5800.0, 0.72},
{5900.0, 0.73},
{6000.0, 0.74},
{6100.0, 0.75},
{6200.0, 0.75},
{6300.0, 0.76},
{6400.0, 0.77},
{6500.0, 0.78},
{6600.0, 0.78},
{6700.0, 0.79},
{6800.0, 0.80},
{6900.0, 0.80},
{7000.0, 0.81},
{7100.0, 0.81},
{7200.0, 0.82},
{7300.0, 0.82},
{7400.0, 0.83},
{7500.0, 0.83},
{7600.0, 0.84},
{7700.0, 0.84},
{7800.0, 0.85},
{7900.0, 0.85},
{8000.0, 0.86},
{8100.0, 0.86},
{8200.0, 0.86},
{8300.0, 0.87},
{8400.0, 0.87},
{8500.0, 0.87},
{8600.0, 0.88},
{8700.0, 0.88},
{8800.0, 0.88},
{8900.0, 0.89},
{9000.0, 0.89},
{9100.0, 0.89},
{9200.0, 0.90},
{9300.0, 0.90},
{9400.0, 0.90},
{9500.0, 0.90},
{9600.0, 0.91},
{9700.0, 0.91},
{9800.0, 0.91},
{9900.0, 0.91},
{10000.0, 0.92}
};
// * * * * * * * * * * * * * Private member funstions * * * * * * * * * * * //
Foam::List<Foam::Tuple2<Foam::scalar, Foam::scalar> >
Foam::radiation::blackBodyEmission::tableToList() const
{
List<Tuple2<scalar, scalar> > newList(nTableEntries);
forAll(newList, i)
{
newList[i].first() = emissivePowerTable[i][0];
newList[i].second() = emissivePowerTable[i][1];
}
return newList;
}
const Foam::List<Foam::Tuple2<Foam::scalar, Foam::scalar> >
Foam::radiation::blackBodyEmission::emissivePowerTable
(
IStringStream
(
"( 1000 0.00032)"
"( 1100 0.00091)"
"( 1200 0.00213)"
"( 1300 0.00432)"
"( 1400 0.00779)"
"( 1500 0.01280)"
"( 1600 0.01972)"
"( 1700 0.02853)"
"( 1800 0.03934)"
"( 1900 0.05210)"
"( 2000 0.06672)"
"( 2100 0.08305)"
"( 2200 0.10088)"
"( 2300 0.12002)"
"( 2400 0.14025)"
"( 2500 0.16135)"
"( 2600 0.18311)"
"( 2700 0.20535)"
"( 2800 0.22788)"
"( 2900 0.25055)"
"( 3000 0.27322)"
"( 3100 0.29576)"
"( 3200 0.31809)"
"( 3300 0.34009)"
"( 3400 0.36172)"
"( 3500 0.38290)"
"( 3600 0.40359)"
"( 3700 0.42375)"
"( 3800 0.44336)"
"( 3900 0.46240)"
"( 4000 0.48085)"
"( 4100 0.49872)"
"( 4200 0.51599)"
"( 4300 0.53267)"
"( 4400 0.54877)"
"( 4500 0.56429)"
"( 4600 0.57925)"
"( 4700 0.59366)"
"( 4800 0.60753)"
"( 4900 0.62088)"
"( 5000 0.63372)"
"( 5100 0.64606)"
"( 5200 0.65794)"
"( 5300 0.66935)"
"( 5400 0.68033)"
"( 5500 0.69087)"
"( 5600 0.70101)"
"( 5700 0.71076)"
"( 5800 0.72012)"
"( 5900 0.72913)"
"( 6000 0.73778)"
"( 6100 0.74610)"
"( 6200 0.75410)"
"( 6300 0.76180)"
"( 6400 0.76920)"
"( 6500 0.77631)"
"( 6600 0.78316)"
"( 6700 0.78975)"
"( 6800 0.79609)"
"( 6900 0.80219)"
"( 7000 0.80807)"
"( 7100 0.81373)"
"( 7200 0.81918)"
"( 7300 0.82443)"
"( 7400 0.82949)"
"( 7500 0.83436)"
"( 7600 0.83906)"
"( 7700 0.84359)"
"( 7800 0.84796)"
"( 7900 0.85218)"
"( 8000 0.85625)"
"( 8100 0.86017)"
"( 8200 0.86396)"
"( 8300 0.86762)"
"( 8400 0.87115)"
"( 8500 0.87456)"
"( 8600 0.87786)"
"( 8700 0.88105)"
"( 8800 0.88413)"
"( 8900 0.88711)"
"( 9000 0.88999)"
"( 9100 0.89277)"
"( 9200 0.89547)"
"( 9300 0.89807)"
"( 9400 0.90060)"
"( 9500 0.90304)"
"( 9600 0.90541)"
"( 9700 0.90770)"
"( 9800 0.90992)"
"( 9900 0.91207)"
"(10000 0.91415)"
"(12000 0.94505)"
"(15000 0.96893)"
"(20000 0.98555)"
"(30000 0.99529)"
"(40000 0.99792)"
"(50000 0.99890)"
)()
);
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
@ -154,7 +145,7 @@ Foam::radiation::blackBodyEmission::blackBodyEmission
:
table_
(
tableToList(),
emissivePowerTable,
interpolationTable<scalar>::CLAMP,
"blackBodyEmissivePower"
),

View File

@ -28,6 +28,9 @@ Class
Description
Class black body emission
Table of black body emissive power taken from:
Modest, "Radiative Heat Transfer", pp.775-777, 1993
SourceFiles
blackBodyEmission.C
@ -58,11 +61,8 @@ class blackBodyEmission
{
public:
//- Number of entries in the table
static const label nTableEntries = 91;
//- Static table of black body emissive power
static const scalar emissivePowerTable[nTableEntries][2];
static const List<Tuple2<scalar, scalar> > emissivePowerTable;
private:
@ -87,9 +87,6 @@ private:
// Private member functions
//- Convert the static table into a list of Tuple2<scalar, scalar>
List<Tuple2<scalar, scalar> > tableToList() const;
scalar fLambdaT(const scalar lambdaT) const;