diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel/basicChemistryModel.H b/src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel/basicChemistryModel.H
index 553fc1a01..eb927e71b 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel/basicChemistryModel.H
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel/basicChemistryModel.H
@@ -175,27 +175,6 @@ public:
 
                 //- Return the heat release, i.e. enthalpy/sec [m2/s3]
                 virtual tmp<volScalarField> dQ() const = 0;
-
-                //- Return the mass fraction of the *specieName
-                //- 28.Nov.2017 Jinwoo Park
-                virtual tmp<volScalarField> Yspecie(const word *specieName) {}
-
-                //- Return the mass fraction of the *specieName
-                //- 28.Nov.2017 Jinwoo Park
-                virtual tmp<volScalarField> MoleFracSpecie(const word *specieName) {}
-
-                //- Return the epsilon
-                virtual tmp<volScalarField> epsilon() {}
-
-                //- Return the k
-                virtual tmp<volScalarField> k() {}
-
-                //- Return the U field
-                virtual tmp<volVectorField> U() {}
-
-                //- Return the turbulent kinematic diffusivity of chemical species (Le=1 is assumed)
-                virtual tmp<volScalarField> nut() {}
-
 };
 
 
diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C
index 6807df96e..9747e979c 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C
@@ -585,249 +585,6 @@ Foam::chemistryModel<CompType, ThermoType>::Sh() const
     return tSh;
 }
 
-// * * * * * * * quoted by JW 20170320   (13) * * * * * * * * * * * * //
-
-
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volScalarField>
-Foam::chemistryModel<CompType, ThermoType>::Yspecie(const word *specieName) const
-{
-    tmp<volScalarField> tYspecie
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "Yspecie",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedScalar("zero", dimless, 0.0)
-        )
-    );
-
-    scalarField& Yspecie = tYspecie.ref();
-
-    scalar initializationchecker=0.0;
-
-    forAll(Yspecie, celli)
-    {
-        if (Yspecie[celli] !=0.0){
-            initializationchecker = 1.0;
-        }
-    }
-
-    if (initializationchecker != 0.0){
-        Info<<"mass fraction field failed to be initialized!"<<endl;
-    }
-
-    char SpecieIndex = dynamic_cast<const reactingMixture<ThermoType>&>(this->thermo()).species()[*specieName];
-
-
-    forAll(Yspecie, celli)
-    {
-        Yspecie[celli] = Y_[SpecieIndex][celli];
-    }
-
-    return tYspecie;
-}
-
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volScalarField>
-Foam::chemistryModel<CompType, ThermoType>::MoleFracSpecie(const word *specieName) const
-{
-    tmp<volScalarField> tMoleFracSpecie
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "MoleFracSpecie",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedScalar("zero", dimless, 0.0)
-        )
-    );
-
-    scalarField& MoleFracSpecie = tMoleFracSpecie.ref();
-
-    scalar initializationchecker=0.0;
-
-    forAll(MoleFracSpecie, celli)
-    {
-        if (MoleFracSpecie[celli] !=0.0){
-            initializationchecker = 1.0;
-        }
-    }
-
-    if (initializationchecker != 0.0){
-        Info<<"mole fraction field failed to be initialized!"<<endl;
-    }
-
-    char SpecieIndex = dynamic_cast<const reactingMixture<ThermoType>&>(this->thermo()).species()[*specieName];
-
-    forAll(MoleFracSpecie, celli)
-    {
-
-        scalar totMol(0.0);
-
-        for (label i=0; i<nSpecie_; i++)
-        {
-            totMol += Y_[i][celli]/specieThermo_[i].W();
-        }
-
-        MoleFracSpecie[celli] = (Y_[SpecieIndex][celli]/specieThermo_[SpecieIndex].W())/totMol;
-    }
-
-    return tMoleFracSpecie;
-}
-
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volScalarField>
-Foam::chemistryModel<CompType, ThermoType>::epsilon() const
-{
-    tmp<volScalarField> tepsilon
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "epsilon_temp",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedScalar("zero", dimless, 0.0)
-        )
-    );
-
-    scalarField& epsilon = tepsilon.ref();
-
-    const volScalarField& epsilon_ref = this->db().objectRegistry::lookupObject<volScalarField>("epsilon");
-
-    forAll(epsilon, celli)
-    {
-        epsilon[celli] = epsilon_ref[celli];
-    }
-
-    return tepsilon;
-}
-
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volScalarField>
-Foam::chemistryModel<CompType, ThermoType>::k() const
-{
-    tmp<volScalarField> tk
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "k_temp",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedScalar("zero", dimless, 0.0)
-        )
-    );
-
-    scalarField& k = tk.ref();
-
-    const volScalarField& k_ref = this->db().objectRegistry::lookupObject<volScalarField>("k");
-
-    forAll(k, celli)
-    {
-        k[celli] = k_ref[celli];
-    }
-
-    return tk;
-}
-
-// * * * * * quoted by JW 20170331(15) * * * * * * //
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volVectorField>
-Foam::chemistryModel<CompType, ThermoType>::U() const
-{
-    tmp<volVectorField> tU
-    (
-        new volVectorField
-        (
-            IOobject
-            (
-                "U_temp",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedVector("U_temp", dimLength/dimTime, vector(0.0, 0.0, 0.0))
-        )
-    );
-
-    vectorField& U = tU.ref();
-
-    const volVectorField& U_ref = this->db().objectRegistry::lookupObject<volVectorField>("U");
-
-    forAll(U, celli)
-    {
-        U[celli] = U_ref[celli];
-    }
-
-    return tU;
-}
-
-
-template<class CompType, class ThermoType>
-Foam::tmp<Foam::volScalarField>
-Foam::chemistryModel<CompType, ThermoType>::nut() const
-{
-    tmp<volScalarField> tnut
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "nut_temp",
-                this->mesh_.time().timeName(),
-                this->mesh_,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
-            this->mesh_,
-            dimensionedScalar("zero", pow(dimLength,2)/dimTime, 0.0)
-        )
-    );
-
-    scalarField& nut = tnut.ref();
-
-    const volScalarField& nut_ref = this->db().objectRegistry::lookupObject<volScalarField>("nut");
-
-    forAll(nut, celli)
-    {
-        nut[celli] = nut_ref[celli];
-    }
-
-    return tnut;
-}
 
 template<class CompType, class ThermoType>
 Foam::tmp<Foam::volScalarField>
diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.H b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.H
index d5343adea..720ad96b9 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.H
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.H
@@ -231,29 +231,6 @@ public:
             //- Return the heat release, i.e. enthalpy/sec [kg/m2/s3]
             virtual tmp<volScalarField> dQ() const;
 
-            //- Return mass fraction of *specieName
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volScalarField> Yspecie(const word *specieName) const;
-
-            //- Return mole fraction of *specieName
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volScalarField> MoleFracSpecie(const word *specieName) const;
-
-            //- Return the epsilon
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volScalarField> epsilon() const;
-
-            //- Return the turbulent kinetic energy
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volScalarField> k() const;
-
-            //- Return the velocity field
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volVectorField> U() const;
-
-            //- Return the turbulent kinematic diffusivity of chemical species (Le=1)
-            //- 28.Nov.2017 Jinwoo Park
-            virtual tmp<volScalarField> nut() const;
 
         // ODE functions (overriding abstract functions in ODE.H)
 
diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModels.C b/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModels.C
index b45ebb318..c331cd1c7 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModels.C
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModels.C
@@ -111,23 +111,6 @@ namespace Foam
         rhoChemistryModel,
         icoPoly8EThermoPhysics
     );
-
-    // Chemistry moldels based on absoluteEnthalpy
-    makeChemistryModel
-    (
-        chemistryModel,
-        rhoChemistryModel,
-        gasHaThermoPhysics
-    );
-
-    // Chemistry moldels based on absoluteInternalEnergy
-    makeChemistryModel
-    (
-        chemistryModel,
-        rhoChemistryModel,
-        gasEaThermoPhysics
-    );
-
 }
 
 // ************************************************************************* //
diff --git a/src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/makeChemistrySolvers.C b/src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/makeChemistrySolvers.C
index 03ee392a3..069a78282 100644
--- a/src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/makeChemistrySolvers.C
+++ b/src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/makeChemistrySolvers.C
@@ -88,14 +88,6 @@ namespace Foam
         incompressibleGasEThermoPhysics
     );
     makeChemistrySolverTypes(rhoChemistryModel, icoPoly8EThermoPhysics);
-
-    // Chemistry solvers based on absoluteEnthalpy
-    makeChemistrySolverTypes(rhoChemistryModel, gasHaThermoPhysics);
-    // makeChemistrySolverTypes(rhoChemistryModel, gasHaThermoPhysics);
-
-    // Chemistry solvers based on absoluteInternalEnergy
-    makeChemistrySolverTypes(rhoChemistryModel, gasEaThermoPhysics);
-    // makeChemistrySolverTypes(rhoChemistryModel, gasEaThermoPhysics);
 }