NavierStokesDictionary.C

Go to the documentation of this file.

#include "NavierStokesDictionary.H"

namespace navierstokes {
  
  void CreateDictionaryConfiguration(int numDim,int gridType,int numScalars,bool withDV,bool withAuxData,
                                     bool withLegacyData,bool withGridMetrics,const std::string &configPrefix,
                                     fixtures::ConfigurationType &inConfig)
  {
    std::ostringstream advanceFieldsStream;
    advanceFieldsStream << "rho rhoV rhoE";
    if(numScalars > 0)
      advanceFieldsStream << " scalarVars";
    std::string advancerFields(advanceFieldsStream.str());
    
    
    std::ostringstream allFields;
    std::string myConfigPrefix;
    if(!configPrefix.empty()){
      myConfigPrefix = configPrefix+std::string(":");
    }
    inConfig.SetParameter(myConfigPrefix+std::string("Advancer:Fields"),advancerFields);
    allFields << paramNames << " ";
    inConfig.SetParameter(myConfigPrefix+std::string("Param:Fields"),paramNames);
    std::ostringstream fieldListStream;
    fieldListStream << cvNames << " ";
    if(withLegacyData){
      for(int iDim = 0;iDim < numDim+2;iDim++){
        fieldListStream << "cv0" << iDim+1 << " ";
      }
    }

    
    
    std::list<std::string> legacyAuxVars;
    if(numScalars > 0){
      fieldListStream << scalarNames << " ";
      if(withLegacyData){
        for(int iScalar = 0;iScalar < numScalars;iScalar++){
          std::ostringstream scalarStream;
          scalarStream << "aux" << (iScalar < 9 ? "0" : "")
                       << iScalar+1;
          fieldListStream << scalarStream.str() << " ";
          legacyAuxVars.push_back(scalarStream.str());
        }
      }
    }
    fieldListStream << timeFieldNames << " ";
    if(withDV)
      fieldListStream << dvNames << " ";
    if(withAuxData) {
      fieldListStream << auxDataNames << " ";
    }
    if(withGridMetrics)
      fieldListStream << "gridMetrics ";
    std::string stateFields(fieldListStream.str());
    allFields << stateFields << " ";
    pcpp::io::RenewStream(fieldListStream);
    inConfig.SetParameter(myConfigPrefix+std::string("State:Fields"),stateFields);


    inConfig.SetParameter(myConfigPrefix+std::string("Field:simTime:Meta"),"1 d 8 time");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:inputCFL:Meta"),"1 d 8 stabilityno");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:simCFL:Meta"),"1 d 8 stabilityno");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:inputDT:Meta"),"1 d 8 time");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:simDT:Meta"),"1 d 8 time");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:gridTime:Meta"),"1 m 8 time");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:gridDT:Meta"),"1 m 8 time");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:gridCFL:Meta"),"1 m 8 stabilityno");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:gamma:Meta"),"1 s 8 shratio");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:power:Meta"),"1 s 8 powerLawTransport");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:beta:Meta"),"1 s 8 powerLawTransport");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:bulkViscFac:Meta"),"1 s 8 powerLawTransport");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:sigmaDissipation:Meta"),"1 s 8 artDissParam");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:sigmaDilatation:Meta"),"1 s 8 artDissParam");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:dilatationCutoff:Meta"),"1 s 8 artDissParam");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refRe:Meta"),"1 s 8 ReynoldsNo");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refPr:Meta"),"1 s 8 PrandtlNo");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:Flag:Meta"),"2 s 4 fluidOptions");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:nonDimensional:Meta"),"1 s 4 non-dimensional");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refRho:Meta"),"1 m 8 density");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refLength:Meta"),"1 m 8 length");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refPressure:Meta"),"1 m 8 pressure");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:refTemperature:Meta"),"1 m 8 temperature");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:imask:Meta"),"1 n 4 imask");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:rho:Meta"),"1 n 8 density");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:rhoE:Meta"),"1 n 8 energy");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:pressure:Meta"),"1 n 8 pressure");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:temperature:Meta"),"1 n 8 temperature");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:rhom1:Meta"),"1 n 8 specificVolume");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:cv01:Meta"),"1 n 8 density");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:cv02:Meta"),"1 n 8 x-momentum");
    inConfig.SetParameter(myConfigPrefix+std::string("Field:cv03:Meta"),"1 n 8 y-momentum");
    if(numDim > 2){
      inConfig.SetParameter(myConfigPrefix+std::string("Field:cv04:Meta"),"1 n 8 z-momentum");
      inConfig.SetParameter(myConfigPrefix+std::string("Field:cv05:Meta"),"1 n 8 energy");
    } else {
      inConfig.SetParameter(myConfigPrefix+std::string("Field:cv04:Meta"),"1 n 8 energy");
    }
    if(!legacyAuxVars.empty()){
      std::list<std::string>::iterator auxVarIt = legacyAuxVars.begin();
      while(auxVarIt != legacyAuxVars.end()){
        inConfig.SetParameter(myConfigPrefix+std::string("Field:"+*auxVarIt+":Meta"),"1 n 8 massfrac");
        auxVarIt++;
      }
    }
    if(withGridMetrics){
      std::ostringstream metricMetaStream;
      if(gridType == 0){
        metricMetaStream << numDim << " g 8 gridMetric";
      } else if(gridType == 1){
        metricMetaStream << numDim << " n 8 gridMetric";
      } else {
        metricMetaStream << numDim*numDim << " n 8 gridMetric";
      }
      inConfig.SetParameter(myConfigPrefix+std::string("Field:gridMetrics:Meta"),metricMetaStream.str());
    }


    std::ostringstream Ostr;
    Ostr << numDim << " n 4 sconn";
    inConfig.SetParameter(myConfigPrefix+std::string("Field:sconn:Meta"),Ostr.str());
    pcpp::io::RenewStream(Ostr);
    Ostr << numDim << " n 8 momentum";
    inConfig.SetParameter(myConfigPrefix+std::string("Field:rhoV:Meta"),Ostr.str());
    pcpp::io::RenewStream(Ostr);
    Ostr << numDim << " n 8 velocity";
    inConfig.SetParameter(myConfigPrefix+std::string("Field:velocity:Meta"),Ostr.str());
    pcpp::io::RenewStream(Ostr);
  }
}