I encountered a problem with the NIntegrate function while simulating a physical situation. The expression is quite messy, but it is just an integration over a polynomial so it should be possible. I want to integrate the function over kP and Energy, and then plot the results with T varying from 0 to 10.

    simplifiednumThermalConductivityIntegrationContentPosSpin=
{-((4.24285*10^44 (1.00032*10^-30 - 
    6.371*10^-31 Energy)^2 (1.3836*10^6 + 1. Energy^6 + 
    259200. kP^4 + 
    Energy^5 (0.000999182 Sqrt[-9 + 1000000 Energy^2] - 
       1.41159 kP^2) - 
    598856. Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
     2.3104 + 1.47148 Energy - kP^2] + 
    kP^2 (-1.19771*10^6 + 
       259200. Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2]) + 
    Energy^2 (-2.306*10^11 - 
       1.26035*10^-8 Sqrt[-9 + 1000000 Energy^2] - 
       4.32*10^10 kP^4 + 
       3.32698*10^10 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] - 
       3.72529*10^-9 Sqrt[-9 + 1000000 Energy^2] Sqrt[
        2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] + 
       kP^2 (1.99619*10^11 - 
          1.44*10^10 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
           2.3104 + 1.47148 Energy - kP^2])) + 
    Energy^4 (9.35396*10^10 - 1.5 kP^4 - 
       1.68212 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] - 
       0.00146484 Sqrt[-9 + 1000000 Energy^2] Sqrt[
        2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] + 
       kP^2 (2. - 0.000403313 Sqrt[-9 + 1000000 Energy^2] - 
          1.88631 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
           2.3104 + 1.47148 Energy - kP^2])) + 
    Energy (2.91038*10^-11 - 
       7.68665*10^7 Sqrt[-9 + 1000000 Energy^2] - 
       1.44*10^7 Sqrt[-9 + 1000000 Energy^2] kP^4 + 
       1.16415*10^-10 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] - 
       3.32698*10^7 Sqrt[-9 + 1000000 Energy^2] Sqrt[
        2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] + 
       1.44*10^7 Sqrt[-9 + 1000000 Energy^2]
         kP^2 (4.6208 + 
          1. Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
           2.3104 + 1.47148 Energy - kP^2])) + 
    Energy^3 (-9.37702*10^-6 + 
       3.11799*10^7 Sqrt[-9 + 1000000 Energy^2] - 
       0.000476971 Sqrt[-9 + 1000000 Energy^2] kP^4 - 
       7.62939*10^-6 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2] + 
       kP^2 (0.0000145193 + 
          0.000976563 Sqrt[-9 + 1000000 Energy^2] + 
          
          0.00146484 Sqrt[-9 + 1000000 Energy^2] Sqrt[
           2.3104 - 1.47148 Energy - kP^2] Sqrt[
           2.3104 + 1.47148 Energy - kP^2]))) Sech[(
   5800.46 (-1.57011 + Energy))/
   T])/((1. Energy - 
    0.001 Sqrt[-9 + 1000000 Energy^2]) (1. Energy + 
    0.001 Sqrt[-9 + 1000000 Energy^2]) (Sqrt[-9 + 
      1000000 Energy^2] (-0.0023104 + 0.001 kP^2 + 
       0.003 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2]) + 
    Energy (2.3104 - 8.88178*10^-20 Sqrt[-9 + 1000000 Energy^2] - 
       1. kP^2 + 
       1. Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2])) (Sqrt[-9 + 
      1000000 Energy^2] (-0.0023104 + 0.001 kP^2 + 
       0.003 Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2]) + 
    Energy (2.3104 + 1.77636*10^-19 Sqrt[-9 + 1000000 Energy^2] - 
       1. kP^2 + 
       1. Sqrt[2.3104 - 1.47148 Energy - kP^2] Sqrt[
        2.3104 + 1.47148 Energy - kP^2])) T^2))}

Now, trying the following line, I get the NIntegrate::inumr error in the headline:

answer = NIntegrate[simplifiednumThermalConductivityIntegrationContentPosSpin,
{kP, 0, 1}, {Energy, 0, Infinity}]

I have also tried to integrate the energy only from 0 to 0.55, but that won't solve the problem.

Do any of you have any tips as to what I should try to be able to plot this? Any help will be greatly appreciated!