1. 26 Nov, 2021 2 commits
  2. 07 Nov, 2021 1 commit
    • Hang Zhou's avatar
      Modifying the method to set char model for coal combustion. CHar model is... · 9ca8504b
      Hang Zhou authored
      Modifying the method to set char model for coal combustion. CHar model is given by "Particles: CharModel: Name:". Three char models are available: GasifOxid, CCK and OFF. When "GasifOxid" is chosen, "GasificationModel" and "OxidationModel" need to be given by "Particles: CharModel: GasificationModel: " and "Particles: CharModel: OxidationModel: ". If CCK model is chosen, the sensitivity of the model can be turn on/off by "Particles: CharModel: SensitivityOfCCK:" (default is OFF).
      9ca8504b
  3. 08 Oct, 2021 1 commit
  4. 23 Sep, 2021 1 commit
    • Hang Zhou's avatar
      1. Correct the mistake in the sensitivity of... · 5613aa2f
      Hang Zhou authored
      1. Correct the mistake in the sensitivity of `src/particles/coal/CharCombustion/CCK/CharSpeciesRHS.h`. The missing nagative sign is added.
      2. Modifying the sensitivity of CCK model. Our results have shown that using the real sensitivity or setting the sensitivity to be zero does not affect the results. So, we can set the sensitivity to be zero for most of our cases
      because the simulation will be faster. We can turn on/off this in the input file by “Particles: CharModel: SensitivityOfCCK: ON/OFF”. The default value of “SensitivityOFCCK” is OFF.
      3. The normalized number fraction of particles per gas mass (n_p,norm), instead of the real number fraction of particles per gas mass (n_p) is solved now. n_p,norm=n_p/n_p,stand. The standard value (n_p,stand) can be given in the input file by “Particles: StandardNumFrac:“. Its default value is 1e9. The reason for this is that very big sensitivity is obtained when using CCK model, because n_p is very big (In equation dn_p/dt=-n_p*S_mg/rho_g/u. The sensitivity of the numerator is dn_pdv*S_mg+n_p*dS_mgdv. The second term is big because n_p is very big. Thus, big sensitivity is obtained from this).
      5613aa2f
  5. 04 May, 2021 1 commit
  6. 30 Sep, 2020 1 commit
  7. 29 Sep, 2020 2 commits
  8. 26 Sep, 2020 1 commit
  9. 25 Sep, 2020 1 commit
  10. 23 Aug, 2020 4 commits
    • Hang Zhou's avatar
      When `Restrat` to get S-Curve, that is, using steady state as initial... · 1cfdbfc6
      Hang Zhou authored
      When `Restrat` to get S-Curve, that is, using steady state as initial conditions to generate the upper part of S-Curve for combustion, initial value of dual time step `ds` should be set based on `dsInitial` or `GESAT_MAX` givein in the input file. When restarting from one specific dual time step for steady-state running or one specific physical time step for transient running,
      1cfdbfc6
    • Hang Zhou's avatar
      1. When particle size is `UserDefinedList`, `NumbersFraction`, instead of... · 1655fb30
      Hang Zhou authored
      1. When particle size is `UserDefinedList`, `NumbersFraction`, instead of `Numbers`, is given in the input file to give the fractions of number of particles with each size. So, the value type should be double, not int.
      2. When paricle size distribution is `Uniform', 'Normal', 'LogNormal' or 'Weibull', the number of particle sizes should be given in the input file by `particle:{Size: {NumberOfSizes:}}`. The fraction of numbers of each size is `1/NumberOfSizes`.
      1655fb30
    • Hang Zhou's avatar
      Solving the problem about `stof: out of range` when setting Fields based on... · 88fd3685
      Hang Zhou authored
      Solving the problem about `stof: out of range` when setting Fields based on database for restarting. Using `std::stod`, instead of `std::stof` will solve the problem.
      88fd3685
    • Hang Zhou's avatar
      Modifying the method for giving particle loadings for coal particle. Two... · 39fd8ac8
      Hang Zhou authored
      Modifying the method for giving particle loadings for coal particle. Two methods are provided. One is giving fuel-air equivalence ratio in the input file by `Particles: FuelAirEquivalenceRatio:'. The other is giving particle loading directly in the input file by `Particles: {ParticleLoading: ParticleLoadingInflow:}'. For the second method, equivalence ratio is not required to be given in the input file.
      39fd8ac8
  11. 14 Aug, 2020 1 commit
    • Hang Zhou's avatar
      Making dual time step `ds` and physical simulation time `simulation_time` to... · 2508d0cf
      Hang Zhou authored
      Making dual time step `ds` and physical simulation time `simulation_time` to be saved by default because `ds` and `simulation_time` from the restart step are requireded to setup the initial conditions for the next time step. If they are not read from the restart step, `simulation_time` will be set to be zero, which is not what we want when restarting from one specific physical time step for transient process, and `ds` will be set to be the `dsInitial` given in the inout file, which will slow down the converging process if `dsInitial` is too small or blow up the running if `dsInitial` is too big.
      2508d0cf
  12. 11 Aug, 2020 1 commit
  13. 10 Aug, 2020 2 commits
    • Hang Zhou's avatar
      Do not need to save mole fractions of all species by default (mass fractions... · 7d750a52
      Hang Zhou authored
      Do not need to save mole fractions of all species by default (mass fractions have been saved by default).
      7d750a52
    • Hang Zhou's avatar
      Modifying the method to get mole fraction of the last species for... · b6ea6c96
      Hang Zhou authored
      Modifying the method to get mole fraction of the last species for initial/inflow condiitons when species fractions are set based on database. When species fractions are given in the input file, `pokitt::ERRORSPECN` is used (error is given when getting fraction of the last species to be negative or bigger than one.) When species fractions are given based on the database, `pokitt::CLIPSPECN' is used. (For the fraction of the last species, it is set to be `1` when value bigger than one is obtained, and is set to be `0` when negative value is obtained.). In this way, we need to be careful when writing the database to make sure the sum of the first (n-1) species fractions is equal or less than 1.
      b6ea6c96
  14. 04 Aug, 2020 1 commit
    • Hang Zhou's avatar
      Adding the function which allows setting initial or inflow conditions based on... · 80acac46
      Hang Zhou authored
      Adding the function which allows setting initial or inflow conditions based on database coresponding to the grid value. For example, when we have various surrounding temperatures for heat transfer, and want to have various inflow temperatures for each surrouding temperature, we can read the values for temperature from the database, like 'Data.txt'. We can give it by `InitialConditions: {-value: FromDatabase -variable:T}'and 'InitialInflowData: {FileName:Data.txt}' in the input file. Also, if we want to set species mole fractions based on database as well, we can give it by `InitialInflowData: {InitialSpeciesSetting: FromDataBase}' or 'InitialInflowData: {InflowSpeciesSetting: FromDataBase}', instead of specifying species in `InitialConditions` or 'InflowConditions'. (The default for `InitialSpeciesSetting` and `InflowSpeciesSetting` is `FromInputFile`. For this, species fractions need to be given in `InitialConditiobs' or 'InflowConditions' in the input file. Fractions of species not specified in the input file or the database are set to be zero.)
      80acac46
  15. 30 Jul, 2020 1 commit
    • Hang Zhou's avatar
      When checking if a speices name read from input file (for example, 'X_O2' and... · 10acb120
      Hang Zhou authored
      When checking if a speices name read from input file (for example, 'X_O2' and 'Y_O2', or 'X_O2_inflow' and 'Y_O2_inflow') is inflow condition, we need to check if the last seven terms of the name are `_inflow`. (if name.substr(name.size()-7)=="_inflow"). However, the mininum size of name is three, liek `X_H`. Checking the last seven terms will give error when the size of name is less than seven. So, I just check if the last three terms are 'low'.
      10acb120
  16. 23 Jul, 2020 2 commits
  17. 17 Jun, 2020 1 commit
    • Hang Zhou's avatar
      Deleting the weight functions used to calculate the source terms in CPD... · 1cabe72b
      Hang Zhou authored
      Deleting the weight functions used to calculate the source terms in CPD models. Because the source term is releated to the mass of labile bridge and delta, it will become zero when the mass is zero. And the source terms are set to be zero if mass is negative. So, it will not oscillate around zero like what we have for vaporization model, that is, weight function is not required.
      1cabe72b
  18. 27 May, 2020 1 commit
  19. 26 May, 2020 1 commit
    • Hang Zhou's avatar
      Modifying the method to set initial/inflow values of coal particle loading,... · 39d42a6a
      Hang Zhou authored
      Modifying the method to set initial/inflow values of coal particle loading, CPD parameters (l, delta_i and g_i), and tar and soot mass. They can be given in the input file by `Particles: {ParticleLoading: ParticleLoadingInflow: }`, `Particles:{CPDInflow:{l: deltai: gi:} CPDInit:{l: deltai: gi:}}`, and `Particles:{TarSootInflowFrac: TarSootInitFrac:}`, respectively. If they are not given in the input file, they will be set based on the equivalence ratio and coal composition.
      39d42a6a
  20. 19 May, 2020 3 commits
  21. 15 May, 2020 1 commit
    • Hang Zhou's avatar
      Add a coal type named `Given_From_InputFile`, which ultimate and proximate... · 1b31ef35
      Hang Zhou authored
      Add a coal type named `Given_From_InputFile`, which ultimate and proximate analysis are given in the input file by "CoalComposiiton: UltimateAnalysis: {C: H: O: N: S:} ProximateAnalysis:{Moisture: Ash: Volatile: FixedC:}. Also, its CPD data is interpolated based on the CPD data of six exist coal types.
      1b31ef35
  22. 11 May, 2020 2 commits
  23. 05 May, 2020 2 commits
  24. 28 Apr, 2020 4 commits
  25. 15 Apr, 2020 1 commit
    • Hang Zhou's avatar
      1. Modifying the tag name for heat tranfer coefficient and surroundings... · be042414
      Hang Zhou authored
      1. Modifying the tag name for heat tranfer coefficient and surroundings temperature for calculating heat transfer between gas and surrounding: Using lower case with words separated by underscores, instead of upper case.
      2. Modifying the input files for test in `rtest` directory based on the new method giving heat transfer parameters between gas and surroundings.
      be042414
  26. 10 Apr, 2020 1 commit