Numerical Method: Heat Conductions of White Dwarfs, Neutron Stars, Black Holes and Derivation of Perfect Black Body Radiation of Black Holes

Abstract

In this paper, the authors plan to derive the famous perfect black body radiation formula in order to suit the statement for black holes that M. Hawking had ever desired. Based on the calculation resulted by numerical method, a black hole’s thermal conductivity is k = 1.909 W·m−1·K−1, in a condition of a black hole undergoes on thermal conduction and information saved or escaped such complies with AdS/CFT correspondence. This can be generalized to the study of electronic engineering considering that the material where computer information can conserve maximum capability (Bakenstein bound) can tolerate and contribute to this exploration returned as to reference study in cosmology and astrophysics. In another respect, neutron stars possess interestingly the critical point of phases: 4He (I/II) and 3He superflow or superglass in our calculations by the same physical principle worked out and this quite supports the other studies.

Share and Cite:

Su, H.-T.H. and Lee, P.-H. (2025) Numerical Method: Heat Conductions of White Dwarfs, Neutron Stars, Black Holes and Derivation of Perfect Black Body Radiation of Black Holes. Journal of High Energy Physics, Gravitation and Cosmology, 11, 203-208. doi: 10.4236/jhepgc.2025.112016.

1. Introduction

Heat transfer is an important subject in science and is one of the cores of classical physics. This subject is not only used in chemical and mechanical engineering, but also widely used in astrophysics [1]-[9]. Heat transfer is a discipline of thermal engineering that concerns the generation, use, conversion, and exchange of thermal energy between physical systems. It is classified into various mechanisms, such as thermal conduction, thermal convection, thermal radiation, and transfer of energy by changes of phases. Engineers also consider the transfer of mass of differing chemical species (mass transfer in the form of advection), either cold or hot, to achieve heat transfer. While these mechanisms have distinct characteristics, they often occur simultaneously in the same system. The famous physician M. Hawking once submitted the idea of black hole evaporation, which later became known as the Hawking radiation formula: T H = c 3 8πGM k B , but this formula only shows the relationship between the black hole temperature and its gravity parameters. Upon investigation, it was found that there is no statement of any Boltzmann-Planck black body radiation formula in references [10]-[16]. This paper uses deductive methods to improve the discussion of many previous studies that used numerical analysis methods to find the key magic number “4”, thereby connecting it with Bekenstein Bound and finding conductive photometric transfer to be used in a descriptive general equation for heat conduction in compact stars (e.g. the equation for thermal convection applies to neutron stars; the equation to heat conduction applies to black holes). This article is a numerical method and not a theoretical/phenomenological method; astronomy and physics bias a few chemical engineering.

2. Discussions

2.1. Bekenstein Bound

The discussion sections begin with the famous Bekenstein bound, as to be addressed to black-body radiation, eventually. The advantage is quickly forced Hawking radiation connecting with thermal conduction handling object of black holes. Of course, where Newton method is used to unfold the frustration of numerical method when undergoing approximately acquired the routine of Triple-Alpha Process. The black hole’s entropy is denoted as well-known:

S BH = k B A 4 l p 2 (1)

Via entropy derivations, we obtain that the parameter 4T represents in Equation (1). This origins from the thermal conductivity possesses 4-directions a r 1 , a r 2 , a r 3 , a r 4 :

S= | ΔQ | T = 1 4ΔT | i=1 4 Q i ( r 0 ) Q 0 ( R ) |= E B.E. 4( T T 0 ) 1 4 (2)

where E B.E. =M c 2 and T = up 4ΔT . Followings resulted in

S BH = M c 2 4ΔT (3)

Apply M= r c 2 2G substituted into Equation (3). Therefore

S BH = 4π r 2 = A BH 4ΔT c/ ( 8πr ) G/ c 3 = 4π r 2 4 c/ ( 8πrΔT ) =[ k B ] G/ c 3 2 , S BH = k B A 4 l p 2 (4)

This reproduction is very important since it helps one to think the physical sense of form of black hole that required main sequence stars ( 25 M Θ ) raised-up the temperature T4T in the evolution process of collapse.

2.2. The Heat Transfer of White Dwarf

The material in a white dwarf no longer undergoes fusion reactions, so the star has no source of energy. As a result, it cannot support itself by the heat generated by fusion against gravity-collapse, but is supported only by electron degeneracy pressure, causing it to be extremely dense. The physics of degeneracy yields a maximum mass for a non-rotating white dwarf, the Chandrasekhar limit, approximately 1.44 M , beyond which it cannot be supported by electron degeneracy pressure.

White dwarf masses < 10 M is located on in general case. The main mode of heat transfer belongs to thermal radiation, but to a certain extent, the problem of white dwarfs can be regarded as through conductive photometric transfer.

dT dr = 1 k Q r = 1 k l 4π r 2 (5)

Therefore

dT dr = 1 k l 4π r 2 , T 0 T 4dT = l kπ r=R r= r 0 1 r 2 dr (6)

Obviously, we get “four times original temperature expanded in radial directions”

4( T T 0 )= l kπ ( 1 r 0 1 R ) (7)

This equation contained magic number 4 such actually appeals one that black hole Hawking radiation could be derived the corresponding equation of black hole’s black body spectrum (i.e., Stephan-Boltzmann function) M. Hawking desired. As such, this equation is flavored and helpful. Of course above deduction is based on the conclusion that resulted in Bekenstein bound.

2.3. The Heat Transfer of Neutron Stars

The majority mode of heat transfer on the neutron star ( 1.5 M Θ <M3.2 M Θ ) surface belongs to the thermal convection. Indicated local heat flux is heat convention and is expressed as

Q ˙ 2 Q r =h( T s T )=h( T T 0 ),T> T 0 (8)

Therefore

h( T T 0 )=l 2 1 4π r 2 = l 2π r 3 ,r= r 0 0( core ),r=R( surface ) (9)

And [2]2,

lim T T 0 h ¯ ( T T 0 ) = ( L.H. ) Q ˙ 6π r 2 (10)

Followings,

h ¯ lim T T 0 Q ˙ 4( 3π r 2 /2 )( T T 0 ) = 1 4( T T 0 ) lim T T 0 Q ˙ 3π r 2 /2 , 2 3 ( T T 0 ) Q ˙ h ¯ ( 4π r 2 ) ,0<rR 1 3 ( T T 0 ) Q ˙ h( 4π r 2 ) ,0<rR (11)

Due to 4π 0.5+43!/2 , 4πh=h( 0.5+43!/2 )=1/3 , thus we have Heat Transfer Coefficients h= 1 3( 0.5+43!/2 ) 0.027 [W/m2K]1. By this thread, see Triple-Alpha Process. And due to Equation (10), such both lead that the surface R to the inner-core r 0 of neutron stars the heat convention corresponds to SUPERFLOW coefficients (i.e., SUPERFLUIDITY). Supposed that its physical behaviors are like 4He (I/II) [or superglass/supersolid] in regions of strong magnetic fields, otherwise 3He in the outer core through 9 km of neutron stars (i.e., regions of weak magnetic fields) [1] [2]. This mathematical deduction by factorial ( n! ) for cosmology-physics is well to use, because it can reduce the complexities of deduction process.

2.4. The Heat Transfer of Black Holes

Based on the previous deductions and the logics, we can think that the one of the three mode heat transfers such can be worked out the famous black hole’s perfect black body radiation formula that Hawking desired. Here the thermal conduction is the definite option as deduction and does not violate any fundamental physical principle. In detailed:

4( T T 0 )= l kπ ( 1 r 0 1 R ), r 0 0( singularity ),R= r s (12.1)

4( T T 0 )= lim r 0 0 M l kπ ( 1 r 0 c 2 2GM ), (12.2)

4πk( T T 0 )= lim r 0 0 M l( 1 r 0 c 2 2GM ), (12.3)

4! T 0 T T 0 T T 0 T T 0 T ( T T 0 ) d 4 T = M M M M lim r 0 0 M l( 1 r 0 c 2 2GM ) d 4 M , K( T 4 T 0 4 )= M M M M lim r 0 0 M l( 1 r 0 c 2 2GM ) d 4 M ,K=const= A BH εσ A BH σε( T 4 T 0 4 )= M M M M lim r 0 0 M l( 1 r 0 c 2 2GM ) d 4 M } (12.4)

For black holes l=0 and via observation, the terms inside this bracket resulted in . As such, we obtained M lim r 0 0 M l( 1 r 0 c 2 2GM )dM =finite (e.g. power of back-body radiation). So far, the famous black hole’s perfect black body radiation formula Hawking desired is now derived:

P BH = A BH σε( T 4 T 0 4 ) (13)

Based on the above calculation, via investigation of tables of material thermal conductivities, we found that k = 1.909 W·m−1·K−1 (thermal conductivity) corresponds to material made of china/porcelains. As a result, we can think that black holes if were in touch, would be felt like ordinary porcelains, and only if gravity were not taken into consideration (actually this is impossible). As such, it would be such interesting. This thermal property is exactly corresponding to Ads/CFT correspondence since as a material of porcelains, it saves information or lets information escape is not strange.

3. Conclusion

Although this article does not discuss heat transfer in chemical engineering, the paper uses a lot of knowledge in this field as a research subject. We eventually derived the perfect black body radiation formula of black holes that Hawking desired. This formula is very important because it not only reproduces Planck’s blackbody radiation formula but also marks the correctness of the beginning of modern physics, and combines this correctness with new physics, adding further impetus to the creation of new physics. Furthermore, this paper deduced that the material of the surface thermal convection of the neutron star is 4He, and found that it is a superglass, which is consistent with the research results of other studies.

4. Justification

Indicated Equation (13), PBH where B is denoted as Boltzmann/Bekenstein/Black, and H is denoted as Hawking/Henry/Hole, in attitudes of solemnity. Honorifics orders of names are stepped by the year-times and then as to list matters.

Acknowledgements

We would like to thank Dr. Ben Jiu for his useful suggestion to this paper's development.

NOTES

*The authors contribution equaled.

#Corresponding author.

1i.e., Newton method is used to solve h approximately.

2Interestingly, Equation (10) can be rewritten as D= lim η k B T 6πη r 2 r= lim η k B T 6πηr (the famous Einstein relation) where Q= k B T . And the indicated diverging viscosity η simply gives the superglass phase.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Geankoplis, C.J. (2003) Transport Processes and Separation Principles. 4th Edition, Prentice Hall.
[2] Abedin, E.Z. (2023) What Are the Basic Concepts of Engineering Thermodynamics?
[3] New Jersey Institute of Technology, Chemical Engineering Department (2010) B.S. Chemical Engineering.
[4] Lienhard, J.H. and Lienhard, J.H.V. (2019) A Heat Transfer Textbook. 5th Edition, Dover Pub., 3.
[5] Welty, J.R., Wicks, C.E. and Wilson, R.E. (1976) Fundamentals of Momentum, Heat, and Mass Transfer. 2nd Edition, Wiley.
[6] Faghri, A., Zhang, Y.W. and Howell, J. (2010) Advanced Heat and Mass Transfer. Global Digital Press.
[7] Taylor, R.A., Phelan, P.E., Otanicar, T., Prasher, R.S. and Phelan, B.E. (2012) Socioeconomic Impacts of Heat Transfer Research. International Communications in Heat and Mass Transfer, 39, 1467-1473.[CrossRef]
[8] Abbott, J.M., Smith, H.C. and Van Ness, M.M. (2005) Introduction to Chemical Engineering Thermodynamics. 7th Edition, McGraw-Hill.
[9] Çengel, Y. (2003) Heat Transfer: A Practical Approach. 2nd Edition, McGraw-Hill.
[10] Rose, C. (2013) A Conversation with Dr. Stephen Hawking & Lucy Hawking.
[11] Hawking, S.W. (1974) Black Hole Explosions? Nature, 248, 30-31.[CrossRef]
[12] Hamilton, A. (2021) A Black Hole Is a Waterfall of Space.
[13] Levin, J. (2020) Black Hole Survival Guide. Alfred A. Knopf, Penguin Random House.
[14] Susskind, L. and Lindesay, J. (2004) An Introduction to Black Holes, Information and the String Theory Revolution—The Holographic Universe. World Scientific Publishing Co. Pte. Ltd.[CrossRef]
[15] Bekenstein, J.D. (1972) Black Holes and the Second Law. Lettere Al Nuovo Cimento Series 2, 4, 737-740.[CrossRef]
[16] Hawking, S. (1988) A Brief History of Time. Bantam Books.

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.