Next Article in Journal
A Pareto–Pontryagin Maximum Principle for Optimal Control
Next Article in Special Issue
Lifetime and Dynamics of Natural Orbits around Titan
Previous Article in Journal
On the Exiting Patterns of Multivariate Renewal-Reward Processes with an Application to Stochastic Networks
Previous Article in Special Issue
Trajectories Derived from Periodic Orbits around the Lagrangian Point L1 and Lunar Swing-Bys: Application in Transfers to Near-Earth Asteroids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Enhanced Calculation Method of the Heat Rejection System of a Free-Piston Stirling Engine (FPSE) Operating on the Moon

1
Academy of Engineering, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Str., 117198 Moscow, Russia
2
Gas Turbine Technologies Laboratory, Joint Institute for High Temperatures, Russian Academy of Sciences, 13 2 Izhorskaya Str., 125412 Moscow, Russia
3
Nauka Power Technology LLC, 42 1 Bolshoy Boulevard, 121205 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(6), 1168; https://doi.org/10.3390/sym14061168
Submission received: 3 April 2022 / Revised: 15 May 2022 / Accepted: 1 June 2022 / Published: 6 June 2022
(This article belongs to the Special Issue Advances in Mechanics and Control)

Abstract

:
In this paper, an enhanced calculation method of a heat rejection system operating on the moon is presented. This was taken into consideration in the developed calculation method and in the propagation of heat fluxes with the radiation of the removed heat. The developed method made it possible to effectively evaluate the capabilities of various refrigerants and choose the radiator parameters and the refrigerant flow regime in a less time-consuming process and with minimal deviations (<5%) compared to the previously developed two-dimensional radiator model by the authors. A comparative analysis was carried out for two refrigerants: helium and liquid ammonia. It has been established that when using liquid ammonia, there are more possibilities for varying the geometric parameters of the radiator. The use of liquid ammonia as a refrigerant made it possible to reduce the power spent on pumping the refrigerant through the radiator. Using helium, the power for pumping the refrigerant was N R = 5.1 W during a turbulent flow Re = 4500. On the other hand, the power for pumping liquid ammonia was N R = 0.27 W. In addition, using liquid ammonia increased the heat flux radiated by the radiator pipe by 3.9 times, which made it to possible to increase the fin width and reduce the length of the radiator pipe.

1. Introduction

Mankind has always been interested in space exploration. Providing a highly reliable power system plays a critical role in the success of space missions. The space missions of this era are evolving to a whole new level, which includes the colonization of new planets and the construction of bases with a permanent presence of humans. The first mission on the list of modern space missions is the colonization of the Moon. Russia, USA, China, Japan and other countries are planning to build a base on the Moon within the next 15–20 years [1,2,3].
The discovery of water at the Moon’s poles about 10 years ago and the possibility of extracting rare earth gases such as helium 3 (a non-radioactive isotope) revived interest in the Moon. In addition, it is believed that the colonization of the Moon will provide people a useful experience in preparation for the exploration of Mars. Lunar bases will have high power requirements to support activities such as scientific experiments, mining and processing, astronomical observations and surface exploration. Therefore, providing a reliable source of energy at these bases is considered extremely important.
A key component of space nuclear power systems is the prime mover, which converts thermal energy into electrical energy. Among the well-known nuclear converters with a high potential for application in space missions is the free-piston Stirling engine (FPSE). In July 2018, NASA acknowledged the FPSE as the most reliable heat engine in the history of civilization. The FPSE unit began operation in the Glenn research center in 2002 and has been operating continuously for more than 14 years without any signs of degradation [4].
Due to the extremely low temperatures and harsh environmental conditions in space, the heat rejection process is considered to be quite complex and implies high demands on the design. Therefore, designing a highly efficient heat rejection system is crucial for the nominal performance and maximum efficiency of the power generation system. Heat rejection systems for space nuclear power systems should be coordinated with the thermodynamic cycles of the power plants in such a way as to maximize the performance of the space nuclear power system while maintaining the specific mass of the heat rejection system and the total specific gravity of the power generation system (kg/kW) to a minimum [5].
Several heat rejection systems have been implemented in space, one of which is widely used—the fluid loop system. This system transfers heat from the equipment, in our case the power generation system, to the radiators, which then rejects the heat into free space. This system can be a mechanically pumped single-phase circuit or a two-phase heat pipe circuit.
Heat pipes have the significant advantage of being completely passive with no moving parts, making them exceptionally suitable for use in space. A heat pipe is a thin hollow pipe filled with a liquid appropriate for the temperature range in which it is intended to operate. At the hot end, the liquid is in the vapor phase and tries to fill the pipe by passing through the pipe to the cold end, where it gradually condenses into the liquid phase. The walls of the pipe, or corresponding channels cut into the pipe, are filled with wick material which returns the liquid by surface tension to the hot end where it is re-evaporated and recirculated. Heat pipes are most commonly used because of their lightness and high thermal conductivity.
Currently, heat pipes are often used for the thermal control of spacecrafts. The heat pipe evaporator is mounted on the components that need to be cooled, and the condenser is mounted on the heatsink panel to dissipate heat. However, heat pipes have a number of limitations compared to a mechanically pumped loop heat rejection system. These limitations are mainly related to the amount of heat rejected and the flexibility of heat transfer control.
On the other hand, the mechanically pumped loop heat rejection system ensures efficient transfer of large amounts of heat at a controlled heat transfer rate and operating temperature. Several works discuss mechanically pumped single-phase liquid loop heat rejection systems and their potential for future space missions. In one study [6], the experimental and analytical results of a heat rejection system consisting of a single-phase mechanically pumped loop and a space radiator are presented. The prototype was developed for future crewed exploration missions to provide a large amount of heat dissipation capability. In [7], a study on refrigerant selection for a high-temperature mechanically pumped fluid loop for space applications is presented; it describes the trade study used to select the high-temperature working fluid for the system and the initial development testing of loop components. In [8], a mini mechanically pumped loop designed for the thermal control of small satellites is presented. In [9,10], a highly self-adaptive cold plate for the single-phase mechanically pumped fluid loop for spacecraft thermal management and hybrid system are studied.
In addition to mechanically pumped single-phase liquid loop systems, heat pipes are also actively being studied. In [11], a heat rejection technology based on heat pipes for space application is presented, where the feasibility of the proposed model is described. In [12], a titanium–water heat pipe radiator for the thermal management of space fission power systems, such as the Kilopower system, is studied. The testing results demonstrated that the developed radiator is able to transfer the required power at the working temperature of 400 K under space-like testing conditions. In [13], a space power system consisting of FPSEs based on potassium heat pipes is presented. In [14,15], sodium variable conductance heat pipes with a Carbon–Carbon radiator for radioisotope Stirling systems are studied.
The active research and development of heat rejection systems for space applications confirms the importance of the studies conducted in this paper. In mathematical modeling, the symmetry of the distribution of heat fluxes as a result of heat conduction along the radiator fins from the radiator pipes and radiation heat flows from two symmetrical surfaces of the radiator was taken into account. In this work, we consider a mechanically pumped single-phase heat rejection system, which is a continuation of our previous work [16]. The calculation method in our previous work took into account the influence of uneven temperature distribution over the surface of the radiator, variable speed and temperature of the refrigerant in the radiator pipes. This is considered crucial while designing a radiator for space applications. However, the main design parameters of the radiator based on the choice of the refrigerant were determined using a time-consuming iterative process, which is discussed in detail in the work itself.
The authors have continued research in this area and developed an enhanced calculation method of a heat rejection system operating on the moon, which allows the time to be decreased for determining the main design parameters of the radiator and the most effective refrigerant. The enhanced method made it possible to effectively evaluate the capabilities of various refrigerants and choose the radiator parameters and the refrigerant flow regime. This was possible by deriving formulas based on acceptable assumptions, which are explained in detail in the next section. In addition, a comparative analysis is presented for two refrigerants: helium and liquid ammonia. It is important to note that the calculation method is not limited to a specific type of refrigerant. Any refrigerant can be considered in the calculation method; however, the operating range of the refrigerant should meet the requirements of the heat rejection system. The refrigerant operating temperature and pressure ranges before encountering a phase change, thermal capacity, thermal conductivity and viscosity should be taken into consideration when choosing a refrigerant.
Helium and liquid ammonia were considered in this work due to their attractive thermophysical properties and suitable operational conditions for the considered FPSE heat rejection system. In a previously published work of the authors, a detailed comparative study of several refrigerants was presented [17]. Helium gas is known for its high heat-transfer capabilities, one of the main reasons it is used as the working fluid inside the FPSE. On the other hand, the choice of liquid ammonia as a refrigerant is due to its proven advantages in space thermal management systems. Advantages of ammonia include: high thermal capacity, wide range of operating temperatures, low weight when compared to water and its viscosity which leads to minimum pumping power through pipes. The International Space Station’s (ISS) active thermal control systems include a liquid-ammonia coolant loop along the station’s main truss which keeps the station’s electricity-generating solar panels cool [18].
This work is organized as follows: First, an overview of the lunar power plant proposed is described. Second, the main steps of the enhanced calculation method are presented. The dependences of the change in the geometrical parameters of the radiator on the thermodynamic properties of the refrigerant and the flow regime are established. Then, the calculation results are presented. An example of selecting the dimensions of a radiator with a given specific radiation flux per unit area using helium and liquid ammonia is presented. Finally, the conclusions are provided.

2. Materials and Methods

A detailed description of the heat rejection system of the FPSE operating on the Moon has already been presented in one of the published works of the authors [19]. However, it is also presented in this work to give the reader a better understanding of the calculation method and the results presented in this work. An overview diagram of the FPSE and the heat rejection system is shown in Figure 1. The engine model considered in this work was developed by Microgen Engine Corporation and is the same model implemented by Nauka-Power Technology LLC in its autonomous off-grid power solutions for hard to reach and remote facilities which require a reliable source of power [20,21]. This particular FPSE model has a high potential for application in future space missions.
The proposed power plant considers supplying heat from the nuclear reactor to the FPSE. The FPSE has two main moving parts: the displacer and the working piston. Both operate in a closed helium environment and are not mechanically coupled to each other. Helium shuttles from the expansion zone (heater) to the compression zone (cooler) due to the temperature difference and cyclical movements of the displacer and power piston. Between the expansion and compression zones there is a porous heat regenerator, which increases the efficiency of the cooling and heat processes of the helium shuttling between the zones. The mechanical energy of the reciprocating movement of the working piston is transformed into electrical energy through a single-phase synchronous linear generator, the inductor of which with permanent magnets is connected to the working piston. Several FPSEs can be powered by the nuclear reactor in the system, however, for convenience, only one FPSE was considered in this paper.
The FPSE is a combination of heat exchangers: a heater, a regenerator and a cooler, which form an internal circuit. In the cooler, heat is rejected from the helium to the outer circuit. The outer circuit (heat rejection system) consists of a pump and a radiator, in which a refrigerant circulates. The heat transferred from the internal circuit to the outer circuit is rejected into free space by the radiator.
Structurally, the radiator is a set of radiant panels operating according to a parallel scheme of inclusion in the radiator’s hydraulic path, connected together by a system of inlet and outlet collectors. The structural form of the radiating panel can be made in the form of parallel pipes, to which the fins are rigidly attached [22]. In this work, two options for the radiator were considered (Figure 2), in which the diameter of the pipes is smaller (option 1) or greater (option 2) than the thickness of the radiator fin. Each radiator pipe has two fins for radiating heat into free space: one from the left side and another from the right side. The fins are symmetric with respect to the pipe. Appendix A provides more detailed information about the design features of the radiator options.
The heat transfer processes occurring in the outer circuit (heat rejection system) are convective heat transfer ( q c o n v ) and radiation ( q r a d ) . Convective heat transfer occurs in the radiator pipes while the refrigerant is passing through them. Heat is transferred from the refrigerant to the pipes’ walls. Afterwards, heat is radiated from the radiator surface into free space [23,24]. Considering the two heat transfer processes occurring, the following condition must be satisfied:
q L = q r a d = q c o n v
q L = q s n p
To calculate the radiator, we must specify the total heat flux that must be rejected ( q s ); temperature at the inlet and outlet of the radiator pipe ( T i n ; T o u t ); and the thermodynamic parameters of the refrigerant ( μ i n ; μ o u t ;   ρ i n ; ρ o u t ; h i n   ; h o u t   ; λ ).
Next, the total flow of the refrigerant through the radiator is determined:
G s = q s Δ h L    
where the change in stagnation enthalpy through one pipe:
Δ h L   = h i n   h o u t    
Through one radiator pipe, the refrigerant mass flow is as follows:
G r e f = G s n p  
The refrigerant mass flow must be constant and is determined as follows:
G r e f = ρ i n · w i n · π · d 2 4
Heat q l transferred as a result of convection:
q c o n v = α · S t · ( T h T w )
where the surface area for heat transfer:
S t = π · d · L R
The heat transferred to the walls of the radiator pipe is further distributed along the fins of the radiator by thermal conduction and then the heat is rejected by radiation from the surface of the radiator q R :
q R = σ · ε · F R   · ( T w 4 T S 4 )  
where the radiation area:
F R = 2 · L R · B R
It is worth noting that the radiation is considered from two surfaces of the radiator. Detailed information about the calculation of the width of one section of the radiation surface B R and the total width of the radiator B C is given in Appendix A.
Based on the condition of heat balance ( q c o n v = q R ), the thermal state of the radiator fin was determined. As already mentioned, the enhanced calculation method of the heat rejection system allows the time to be decreased for determining the main design parameters of the radiator and the most effective refrigerant. The enhanced calculation method is carried out in two stages:
  • Preliminary calculation to determine the main geometric parameters of the radiator based on the refrigerant chosen using derived formula.
  • Calculation of the radiator based on the result of the preliminary calculations. At this stage, the calculation is carried out using the previously developed method which takes into consideration the uneven distribution of temperatures over the surface of the radiator, variable speed and temperature of the refrigerant inside the pipe. The method for determining the temperature state of the radiator in a two-dimensional formulation is described in detail in the work [19].
The amount of heat radiated from the surfaces of the radiator fins at a constant ambient temperature is determined by the surface area and the temperature level in the radiator fins. Moreover, in formula (9), the temperatures are to the fourth power ( T R 4     T S 4 ) ; therefore, the greater the difference between the surface temperature TR and ambient temperature T S , the more efficient the heat rejection. At a given ambient temperature T S , the specific heat flux per unit area q F depends on the fin surface temperature TR.
For the preliminary calculations, it was assumed that the temperature of the fin is constant and equal to the temperature of the inner surface of the pipe ( T R = T w ), and the temperature of the refrigerant in the pipe is equal to the average temperature:
T m = T i n + T o u t 2
Subsequent studies have shown the validity of this assumption. A temperature difference coefficient for convective heat transfer K T = T w T h is introduced for convenience. The closer the wall temperature T w to the refrigerant temperature T h , the greater the temperature difference between the radiator fin TR and the ambient temperature T S ; therefore, the greater the heat flux that will be rejected per unit area q F . Theoretically, the minimum total radiator area for heat rejection of a given heat flux q s will be when T R = T h :
F R a m i n = q s σ · ε · ( T h 4 T S 4 )  
Then, the maximum heat flux emitted from the radiator surface:
q F m a x = q s F R a m i n = σ · ε · ( T h 4 T S 4 )
According to these formulas, it is possible to determine the limiting values of the radiator area and heat flux for evaluating the results obtained. In practice, it is impossible to achieve them, since a temperature difference is necessary for the heat rejection process, therefore T h > T w and K T < 1.0 .
The total area of the radiation surface of the radiator is as follows:
F R a = 2 · L R · B R a = 2 · L R · B R · n p
F R a = q s σ · ε · ( T w 4 T S 4 )  
A formula for determining the pipe diameter d for a given mean Reynolds number R e m is derived. The Reynolds number of the refrigerant at the pipe inlet can be determined as follows:
R e i n = 2 · R e m μ i n μ o u t + 1
and
R e i n = w i n · d ν i n
Further, we perform rearrangement of the formulas:
w i n · d = R e i n · ν i n
The total flow rate G S is found from the formulas (5) and (6):
G S = ρ i n · w i n · π · d 2 4 · n p
Substituting (18) into (19):
G S = ρ i n · π 4 · n p · R e i n · ν i n · d
Considering:
μ i n = ρ i n · ν i n  
and
μ v = μ i n + μ o u t μ i n · μ o u t  
We obtain a formula showing the relationship between the pipe diameter and the number of radiator pipes:
d · n p = 2 · G S · μ v π · R e m  
Or considering formula (3), we obtain:
d · n p = 2 · q s · μ v π · Δ h L   · R e m  
We combine all the parameters of the initial data for calculation in this formula and introduce a coefficient for convenience:
A 1 = 2 · q s · μ v π · Δ h L    
We obtain:
d · n p = A 1 R e m  
For a given number of pipes, the formula for calculating the diameter of the radiator pipe is as follows:
d = A 1 n p · R e m  
It is possible to set the diameter of the radiator pipe and determine the number of pipes:
n p = A 1 d · R e m  
In this case, the number of pipes must be an integer. This can be achieved by appropriate adjustment of the diameter d. Then, we derive the formula for calculating the width of one section of the radiation surface from ( q r a d = q c o n v ) and ( α · d = N u · λ ):
B R = π 2   · σ · ε   · T h · ( 1 K T ) ( K T · T h )   4 T S 4   · λ ·   N u  
A necessary condition (check Appendix A) in this case is:
B S > d + 2 · t  
The thickness of the pipe (t) is chosen based on the strength calculation of the structure. The total width of the radiator was calculated as follows:
B R a = B R · n p  
Ceteris paribus, setting the coefficient K T or wall temperature T w = T h · K T , we can determine the width of the radiator fin B R , providing the temperature difference ( Δ T = T h T w ). When specifying the number of radiator pipes, we define:
  • Radiation surface area:
    F R = q s σ · ε · n p · ( T w 4 T S 4 )    
  • Radiator pipe length:
    L R = F R 2 · B R  
The formula for determining the length of the radiator can be obtained by substituting into formula (33): the formulas for determining the radiator’s radiation area F R (32); width of the radiation surface for one fin section B R (29); and number of pipes n p (28):
L R = d 2   · 1 T h · ( 1 K T ) · Δ h L   ·   R e m μ v · λ ·   N u  
3.
Specific heat flux of radiator surface:
q F = q s F R · · n p = q s F R a · = σ · ε · ( T w 4 T S 4 )  
The total radiation area of the radiator F R a   and the specific heat flux q F , taking into account the assumptions made for a given heat flux q s and ambient temperature T s , depends on the pipe wall temperature T w and coefficient of temperature difference K T . The dimensional formulas B R and L R are determined by the number of pipes n p , pipe diameter d, thermodynamic parameters of the refrigerant and flow regime Re.
Another important factor for evaluating the efficiency of a radiator design is the power spent on pumping the refrigerant through the radiator. By changing the number of radiator pipes, it is possible to achieve the required geometric dimensions with acceptable power losses for pumping the refrigerant through the radiator. The calculation of the hydraulic resistance inside the radiator was carried out using the following method.
The power spent on pumping the refrigerant through the radiator was calculated as follows:
N R = G s · Δ P R η · ρ h    
The Darcy friction factor for laminar flow:
f = 64 R e m  
For turbulent flow:
f = [ 0.79 l n ( R e m 8 ) ] 2  
The pressure drop was calculated as follows:
Δ P R = f L R d · ρ h · w h 2 2  
Rearranging:
N R = G s · f η   · L R d · w h 2 2  
Considering that the refrigerant velocity:
w h = R e m · ν m d  
We obtain:
N R = G s · f η   · L R d 3 · R e m 2 · ν   m 2 2  
For a laminar flow, this formula can be rearranged, taking into account formula (37):
N R = 32 · G s · ν   m 2 η   · L R d 3 · R e m  
or
N R = A N R · L R d 3 · R e m  
where the pumping power coefficient:
A N R = 32 · G s · ν   m 2 η    
An increase in the pipe length L R and a decrease in the pipe diameter d lead to an increase in the power for pumping the refrigerant. The pipe length L R depends on the diameter d (34), and the diameter d is related to the number of pipes n p (28).
The formulas include the Nusselt number Nu and Reynolds Re, which characterize the features of heat transfer and refrigerant flow [24,25].
  • Nusselt number for laminar flow when P e · d L R < 100 :
    N u = 3.66
    where Peclet number:
    P e = R e · P r
  • Transitional flow from laminar to turbulent when 2300 ≤ Re < 4000:
For the condition 0.5 < Pr < 1.5:
N u = 0.0214 R e 0.8 100 P r 0.4 [ 1 + ( d L R ) 2 3 ]
For the condition 1.5 < Pr < 500:
N u = 0.012 R e 0.87 280 P r 0.4 [ 1 + ( d L R ) 2 3 ]
3.
Nusselt number for turbulent flow:
N u = f 8 · R e · P r R 1 + R 2 f 8 P r n 1  
where Prandtl number:
P r = γ a  
R 1 = 1 , R 2 = 12.7 , n = 2 3  
The above formulas made it possible to study the influence of the flow regime on the geometrical parameters of the radiator and the power spent on pumping the refrigerant depending on the given specific heat flux emitted by the radiator for various refrigerants.
Based on the studies carried out, the main design parameters of the radiator can be preliminarily determined ( d ;   B R ;   L R ;   n p ) . After conducting the preliminary calculation using the above formulas, the obtained results are refined using the previously developed calculation method of the radiator (two-dimensional model), which takes into account the uneven distribution of temperatures over the area of the radiator fin [19].
To compare two refrigerants with the same initial data ( q s ;   T i n / T o u t ; Δ T ) , the following formulas are proposed (subscripts 1 and 2 are used to indicate the refrigerant):
  • From formula (24), we obtain:
    d 1 · n p 1 d 2 · n p 2 = μ v 1 · Δ h L 2   ·   R e m 2 μ v 2 · Δ h L 1   ·   R e m 1
For flows with the same mean value of Reynolds number R e m :
d 1 · n p 1 d 2 · n p 2 = μ v 1 · Δ h L 2   μ v 2 · Δ h L 1  
or
d 1 · n p 1 d 2 · n p 2 = K μ   K h    
where the viscosity and enthalpy ratios:
K μ = μ v 1 μ v 2 ;   K h = Δ h L 1   Δ h L 2    
2.
From formula (29):
B R 1 B R 2 = λ 1 · N u 1   λ 2 · N u 2   = B R 1 B R 2 = K λ · N u 1   N u 2    
where the radiation surface ratio:
K λ = λ 1   λ 2    
3.
From formula (34):
L R 1 L R 2 = d 1 · λ 2 · N u 2   d 2 · λ 1 · N u 1   · μ v 2 · Δ h L 1   ·   R e m 1 μ v 1 · Δ h L 2   ·   R e m 2  
4.
Total required radiation surface, an important parameter for evaluating the heat rejection efficiency of the refrigerant:
S t a = S t · n p = q s · d · 1 T h · ( 1 K T ) · 1 λ ·   N u  
d 1 · n p 1 d 2 · n p 2 = μ v 1 · Δ h L 2   ·   R e m 2 μ v 2 · Δ h L 1   ·   R e m 1  
S t a 1 S t a 2 = d 1 · λ 2 · N u 2   d 2 · λ 1 · N u 1    
For laminar flows ( P e · d L R < 100 ; Nu = 3.66), the formulas are simplified:
B R 1 B R 2 = K λ  
S t a 1 S t a 2 = d 1 · λ 2   d 2 · λ 1   = d 1   d 2   · 1   K λ  
L R 1 L R 2 = d 1 · R e m 1   d 2 · R e m 2   · K h   K λ · K μ
For flows with the same mean value of Reynolds number R e m :
L R 1 L R 2 = d 1   d 2   · K h   K λ · K μ  

3. Results and Discussion

The initial data for calculating the radiator are:
  • The FPSE parameters:
    • The amount of heat that must be rejected from the working fluid in the FPSE;
    • The minimum temperature of the working fluid at the cold side of the FPSE.
  • The refrigerant used in the heat rejection system; thermodynamic properties.
  • Ambient temperature.
The ambient temperature in this study is of great importance due to its direct influence on the efficiency of the heat rejection process. Temperatures on the Moon are extreme. Ranges of temperature changes relative to their location on the Moon and a parametric study of the energy systems of the lunar base can be found in [26,27]. The location of the lunar power plant considered in this work was at one of the poles, where the maximum ambient temperature reaches 200 K (−73 °C), which is considered the most suitable for lunar bases according to several space agencies. At lower ambient temperatures, the difference between the temperature of the fins and the ambient temperature increases, resulting in an increase in the efficiency of heat rejection.
Helium and liquid ammonia were considered as refrigerants; their main properties (gas constant, enthalpy, heat capacity Cp, dynamic viscosity, thermal conductivity) were set according to the data from these works [28,29]. The main parameters of the FPSE considered in the calculation are shown in Table 1.
To compare the efficiency of refrigerants, the following conditions were set constant:
  • The heat flow that must be rejected in the radiator from the FPSE q s = 2500 W;
  • Refrigerant temperatures at the inlet and outlet of the radiator pipe T i n / T o u t = 280/260 ( T m = T h = 270   K ) ;
  • Average pipe wall temperature T w = 269   K , consequently K T = T w T h = 0.996 ;
  • Refrigerant pressure in the heat rejection loop system—1 MPa;
  • Fin thickness h R = 0.005 m.
If the specified temperature difference is provided as Δ T = T h T w = 1   K , then the total radiation surface area of the radiator, regardless of the refrigerant used and the flow regime, will be constant (14): F R a = 14.3   m 2 .
Initially, helium was considered as a refrigerant for the calculation. During a laminar flow ( P e · d L R < 100 ; Nu = 3.66) and based on the given initial data, the fin width remains constant (29) B S = B R = 0.0052   m (Option 1 in Appendix A).
When specifying the pipe diameter, the condition B S > d + 2 · t must be taken into account. For the example under consideration, with a pipe diameter of 2.0 mm and a thickness t = 1.0 mm, the distance between two pipes b r 1 is only 1.2 mm.
Depending on the considered pipe diameter, we determine the number of pipes n p , pipe length L R , as well as the total width of the radiator B a and the power spent on pumping the refrigerant N R for different flow rates R e . The results for a pipe of a diameter d = 2.0 mm are shown in Figure 3, Figure 4 and Figure 5. With a decrease in the refrigerant flow rate, the number of pipes increases and, consequently, the total width of the radiator, while the length of the pipe and the pumping power decrease. During a transitional flow (2300 < Re < 4000), the value of the Nu number increases, the heat rejection from the refrigerant increases, which consequently affects the geometric dimensions of the radiator.
For a pipe with a diameter d = 2.0   mm at R e = 500   B R = 0.0052   m   , the following results were obtained:
  • Number of pipes: n p = 1648 ;
  • Total width of the radiator: B C = B R a = 8.53   m ;
  • Pipe length: L R = 0.84   m ;
  • Power spent on pumping the refrigerant: N R = 5.1   W .
It is worth noting that with an increase in the number of pipes, the dimensions of the distribution collectors of the radiator increase. This would increase the total hydraulic losses at the inlet and outlet of the pipes.
During turbulent flows (Re > 4000), the heat transfer from the refrigerant to the pipe walls increases, therefore providing the required specific heat radiation flux per unit area, and the width of the fin increases.
For the same pipe diameter d = 2.0   mm at R e = 4500 B R = 0.0268   m   , the following results were obtained:
  • Number of pipes: n p = 183 ;
  • Total width of the radiator: B C = B R a = 4.91   m ;
  • Pipe length: L R = 1.45   m ;
  • Power spent on pumping the refrigerant: N R = 222.6   W .
Analyzing the results, it can be seen that the number of pipes decreased. However, the power required to pump the refrigerant through the radiator increased to a high value which is unacceptable. The design can be significantly improved by increasing the fin width of the radiator B R , and increasing the pipe diameter d, which leads to minimizing the power required to pump the refrigerant N R (Option 2 in Appendix A). For instance, the increase in diameter to d = 10.8   mm leads to the following results at R e = 4500 :
  • Number of pipes: n p = 34   ;
  • Total width of the radiator: B C = 0.81   m ;
  • Total width of the radiation surface: B R a = 0.91   m ;
  • Pipe length: L R = 7.83   m ;
  • Power spent on pumping the refrigerant: N R = 7.7   W .
Similar calculations for liquid ammonia were carried out and the results were compared. During a laminar flow ( P e · d L R < 100 ; Nu = 3.66) and based on the given initial data, the fin width remains constant B R = 0.0202   m . For the given initial data, A 1 = 196.2 . The calculation results are shown in Figure 6, Figure 7 and Figure 8.
For a pipe with a diameter d = 2.0   mm at R e = 500 , the following results were obtained:
  • Number of pipes: n p = 196 ;
  • Total width of the radiator: B C = B R a = 3.97   m ;
  • Pipe length: L R = 1.8   m ;
  • Power spent on pumping the refrigerant: N R = 0.009   W .
During turbulent flows, the heat transfer from the refrigerant to the pipe walls increases. At R e = 4500 , the fin width B R increases ( B R = 0.143   m ) . The fin width B R allows the pipe diameter to increase.
The graphs of the dependences of the number of pipes n p , pipe lengths L R , as well as the total width of the radiation surface of the radiator B R a , and the power spent on pumping liquid ammonia N R for various pipes diameters d at Re = 4500 are shown in Figure 9, Figure 10 and Figure 11.
For a pipe with a diameter d = 10.9   mm at R e = 4500   B R = 0.143   m , the following results were obtained:
Number of pipes: n p = 4 ;
  • Total width of radiator: B C = 0.56   m ;
  • Total width of the radiation surface: B R a = 0.57   m ;
  • Pipe length: L R = 12.5   m ;
  • Power spent on pumping the refrigerant: N R = 0.0094   W .
For a pipe with a diameter d = 5.5   mm at R e = 4500   B R = 0.143   m , the following results were obtained:
  • Number of pipes: n p = 8 ;
  • Total width of radiator: B C = 0.138   m ;
  • Total width of the radiation surface: B R a = 1.144   m ;
  • Pipe length:   L R = 6.25   m ;
  • Power spent on pumping the refrigerant: N R = 0.038   W .
Analyzing the results, it can be seen that although the length of the pipe has decreased by 2 times, the power for pumping the refrigerant increased; however, the values remain insignificant.
The conducted studies show that when using liquid ammonia, there are more possibilities for varying the geometric parameters of the radiator and the power for pumping the refrigerant can be significantly reduced. The final decision related to the suitable radiator design can be made based on the operating conditions of the lunar power plant.
Next, the influence of the thermodynamic properties of the refrigerant on the main geometrical parameters of the radiator is studied using the previously derived formulas (53)–(66) for two refrigerants:
  • Helium;
  • Ammonia.
From the formula (55) for a flow with the same R e m and given initial conditions K μ = μ v 1 μ v 2   = 9.49; K h = Δ h L 1   Δ h L 2   = 1.12:
d 1 · n p 1 = K μ   K h   · d 2 · n p 2  
If pipes of the same diameter are used, then:
n p 1 = K μ   K h   · n p 2  
While using helium, the number of pipes n p 1 will be 8.4 greater, compared to using liquid ammonia.
The fin width of the radiator during laminar flow ( P e · d L R < 100 ; Nu = 3.66):
B R 1 = K λ · B R 2  
where the radiation surface ratio:
K λ = λ 1   λ 2     =   0.256
The total heat transfer surface of the pipes:
S t a 1 = S t a 2   K λ    
To reject the required heat from the FPSE using helium as a refrigerant, the total heat radiation surface area of the pipes will be 3.9 times larger, i.e., the specific heat flux emitted by the pipe surface will be 3.9 times smaller, compared to using liquid ammonia. Therefore, the fin width is 3.9 times smaller when using helium, since less heat is radiated from the surface. Liquid ammonia allows the heat flux emitted by the radiator pipe section to increase. Therefore, it is possible to increase the width of the fin and reduce the length of the pipe. When using helium, the number of pipes and fin sections increases.
The results of the preliminary stage of calculation can be summed as follows. Taking into account the accepted assumptions, the total area of the radiator depends only on the temperature difference ( Δ T = T h T w ). If the difference is provided, then the area and the specific heat flux will be the same using any refrigerant. The thermodynamic properties affect the design parameters of the radiator (Pipe diameter d; number of pipes n p ; pipe length L R , total radiator width B a and the pumping power of the refrigerant through the radiator N R ).
It has been established that the use of liquid ammonia as a refrigerant makes it possible to reduce the power for pumping the refrigerant through the radiator. The formulas (67)–(71) can also be used to compare other refrigerants.
After performing the first stage of the enhanced calculation method, we move onto the next stage and perform calculations using the previously developed method (two-dimensional model). A comparative analysis for laminar (Re = 500) and turbulent (Re = 4500) flows for two refrigerants—helium and liquid ammonia—was carried out. The comparison results are presented in Table 2, Table 3, Table 4, Table 5 and Table 6 below.
The results of calculations using the previously developed two-dimensional radiator model by the authors confirmed the acceptable accuracy of the results obtained at the preliminary stage of the enhanced method presented in this work.
When considering a laminar flow (Re = 500), the maximum deviation in the radiator area is 0.83%, NR—1.12%. For turbulent flows (Re = 4500), the error increases:
-
For liquid ammonia:
  • Total area of the radiator—2.3%;
  • NR—2.4%.
-
For helium:
  • Total area of the radiator—1.0%;
  • NR—4.4%
The conducted studies proved the effectiveness of the preliminary calculations using the derived formulas presented in this work. The enhanced method made it possible to effectively evaluate the capabilities of various refrigerants and choose the main design radiator parameters and the refrigerant flow regime in a less time-consuming process and with minimal deviations (<5%), compared to the previously developed two-dimensional radiator model by the authors.

4. Conclusions

Thermal management plays a major role in determining the success of space missions. This work presents an enhanced calculation method for a radiator of the heat rejection system of the FPSE designed to operate on the Moon. The enhanced calculation method allows less time for determining the main design parameters of the radiator and the most effective refrigerant compared to a previously developed calculation method by the authors. The major highlights and conclusions of this work are summarized as follows:
  • The enhanced method was validated by a comparative analysis using helium or liquid ammonia as a refrigerant in the heat rejection system of the FPSE designed to operate on the Moon within the temperature range of 260–280 K.
  • The dependences of the change in the geometrical parameters of the radiator on the thermodynamic properties of the refrigerant and the flow regime were established.
  • An example of selecting the dimensions of a radiator with a given specific radiation flux per unit area using helium and liquid ammonia was presented.
  • It has been established that when using liquid ammonia, there are more possibilities for varying the geometric parameters of the radiator. The use of liquid ammonia as a refrigerant makes it possible to:
    • Reduce the power spent pumping the refrigerant through the radiator. For helium, even with an increase in pipe diameter to d = 10.8 mm with a turbulent flow Re = 4500, the power for pumping the refrigerant was N R = 5.1 W; on the other hand, while using liquid ammonia N R = 0.27 W.
    • Increase the heat flow per unit area of the radiator pipe by 3.9 times, which allows the fin width to be increased and reduces the length of the radiator pipe.
  • The conducted studies proved the effectiveness of the preliminary calculations using the derived formulas presented in this work. The enhanced method made it possible to effectively evaluate the capabilities of various refrigerants and choose the radiator parameters and the refrigerant flow regime in a less time-consuming process and with minimal deviations (<5%) compared to the previously developed two-dimensional radiator model by the authors. The largest deviations were noted for turbulent flows (Re = 4500):
    • For liquid ammonia:
      • Total area of the radiator—2.3%;
      • NR—2.4%.
    • For helium:
      • Total area of the radiator—1.0%;
      • NR—4.4%.

Author Contributions

Conceptualization, S.S. and Y.A.; Formal analysis, S.S. and H.K.; Investigation, S.S., Y.A. and H.K.; Methodology, S.S. and Y.A.; Supervision, S.S., M.S. and I.T.; Validation, M.S., Y.A. and I.T.; Visualization, H.K.; Writing—original draft, S.S. and Y.A.; Writing—review and editing, S.K. and H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This paper has been supported by the RUDN University Strategic Academic Leadership Program.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

q r a d Amount of heat released due to radiation, J
q c o n v Amount of heat released due to convection, J
q L Amount of heat that must be rejected from the refrigerant through one pipe, J
q s Amount of heat that must be rejected from the radiator, J
n p The number of radiator pipes required for the heat rejection process
Δ h L   Change in stagnation enthalpy through one pipe, J
Δ h   Change in specific enthalpy through one pipe, J/kg
G s Refrigerant mass flow rate through radiator, kg/s
h i n   Stagnation enthalpy at the inlet, J
h o u t   Stagnation enthalpy at the outlet, J
G r e f Refrigerant mass flow rate through one pipe, kg/s
ρ i n Refrigerant density at the pipe inlet, kg m 3
w i n   Refrigerant flow speed at the pipe inlet, m/s
d The inner diameter of the radiator pipe, m
α Heat transfer coefficient, W m 2 · K
StSurface area for heat transfer, m2
T h Refrigerant temperature, K
T m Refrigerant mean temperature, K
T w Pipe wall temperature, K
LRPipe length, m
q R Heat flux as a result of radiation, W
σ Stefan–Boltzmann constant,   W m 2   · K 4
ε Emissivity coefficient
F R Radiator radiation area, m2
T S   Ambient temperature, K
BRFin width, m
T R Fin surface temperature, K
F R a m i n The minimum total area of radiator radiation for heat rejection of a given heat flux, m2
q F m a x Maximum heat flux per unit area of radiator radiation, J
F R a Total radiator radiation area, m2
B R a Total width of the radiation surface of the radiator, m
R e m Mean value of Reynolds
R e i n Reynolds number at the inlet
μ i n Refrigerant dynamic viscosity at the inlet, Pa.s
μ o u t Refrigerant dynamic viscosity at the outlet, Pa.s
ν i n Refrigerant kinematic viscosity at the inlet, m 2 s  
μ v Viscosity coefficient
A 1 Coefficient for convenience
K T Coefficient of temperature difference during convective heat exchange between the refrigerant and the pipe wall
B S Distance between two pipes’ centers, m
tPipe wall thickness, m
b r 1 Fin width (related to a specific design option), m
h R Fin thickness, m
q F Specific heat flux per unit area of radiator surface, J/m2
N R The power spent on pumping the refrigerant through the radiator, W
Δ P R Pressure drop due to friction in the radiator pipes, Pa
ρ h Refrigerant density in the radiator, kg m 3
η Efficiency of the refrigerant pump
f Friction factor
w h Refrigerant velocity, m/s
ν m Refrigerant mean kinematic viscosity, m 2 s  
A N R Pumping power coefficient
N u Nusselt number
λ Refrigerant thermal conductivity, W m   · K
R e Reynolds number
P e Peclet number
PrPrandtl Number
aThermal diffusivity, m2/s
γ Kinematic viscosity, m2/s
K μ Viscosity ratio of the two refrigerants
K h Enthalpy ratio of the two refrigerants
K λ Radiation surface ratio of the two refrigerants
S t a Total required radiation surface, m2
B C Total width of radiator, m

Appendix A. Calculation of the Width of the Radiator Radiation Surface

Each radiator pipe has two fins for radiating heat into free space: one from the left side, and another from the right side. The fins are symmetric with respect to the pipe. When calculating the width of the radiator’s radiation surface, two options were considered (Figure A1 and Figure A3).
Figure A1. Radiator design (option 1).
Figure A1. Radiator design (option 1).
Symmetry 14 01168 g0a1
The radiation surface width was calculated as follows:
B S = B R = d + 2 t + b r 1  
Figure A2. (a) Radiator design (option 2); (b) Calculation of radiator design (option 2).
Figure A2. (a) Radiator design (option 2); (b) Calculation of radiator design (option 2).
Symmetry 14 01168 g0a2aSymmetry 14 01168 g0a2b
For option 2, the radiation surface width was calculated as follows:
B R = b r 1 + 2 A D 1 + A C  
where
B S = b r 1 + 2 O D
A C π d + 2 t = β 2 π
A C = d + 2 t 2 β
α = a r c s i n h R d + 2 t
β = π 2 α
O D = d 2 + t
A D 1 = B D = O D O B = d 2 + t d 2 + t c o s α = d 2 + t 1 c o s α  
Rearranging them, the following formulas are obtained:
B R = b r 1 + d + 2 t 1 c o s α + β  
B S = b r 1 + d + 2 t
Figure A3. Total dimensions of the radiator (option 2).
Figure A3. Total dimensions of the radiator (option 2).
Symmetry 14 01168 g0a3
The total width of the radiator was calculated as follows:
B C = B S · n p  
For option 2, the total width of the radiation surface of the radiator B R a and the total width of the radiator B C ( B R a > B C ) are different.
When using option 2 with protruding pipes’ radiation surfaces, it must be taken into account that the amount of heat determined by the Stefan–Boltzmann law is reflected in different directions. The change in radiation in individual directions is determined by Lambert’s law, which is strictly defined for a completely black body. For rough bodies, this law is confirmed only for angles of deviation from the normal to the radiation surface up to 60° [30]. Lambert’s law is used to calculate radiant heat transfer between surfaces of finite dimensions. The scheme of option 2 shown in Figure A4 is simplified. From protruding surfaces, radiation can reach the surface of the fins (a–b) and adjacent radiator pipes (b–c).
Figure A4. Radiation from the surface of the radiator (option 2).
Figure A4. Radiation from the surface of the radiator (option 2).
Symmetry 14 01168 g0a4
To avoid such situations in a real design, it is necessary to perform special smooth transitions (fillets) from the fin to the pipe as shown in Figure A5.
Figure A5. Radiation from the enhanced surface of the radiator (option 2).
Figure A5. Radiation from the enhanced surface of the radiator (option 2).
Symmetry 14 01168 g0a5

References

  1. Russia and China Team Up to Build a Moon Base|Science|The Guardian. Available online: https://www.theguardian.com/science/2021/jun/25/russia-china-team-up-build-moon-base (accessed on 20 January 2022).
  2. NASA’s Artemis Base Camp on the Moon Will Need Light, Water, Elevation|NASA. Available online: https://www.nasa.gov/feature/goddard/2021/nasa-s-artemis-base-camp-on-the-moon-will-need-light-water-elevation (accessed on 20 January 2022).
  3. ESA—Future Moon Base. Available online: https://www.esa.int/ESA_Multimedia/Images/2018/11/Future_Moon_base (accessed on 20 January 2022).
  4. Stirling Convertor Sets 14-Year Continuous Operation Milestone|NASA. Available online: https://www.nasa.gov/feature/glenn/2020/stirling-convertor-sets-14-year-continuous-operation-milestone (accessed on 20 January 2022).
  5. Heat Rejection Technologies for Nuclear Systems|NASA SBIR & STTR Program Homepage. Available online: https://sbir.nasa.gov/content/heat-rejection-technologies-nuclear-systems (accessed on 2 March 2022).
  6. Iwata, N.; Nakanoya, S.; Nakamura, N.; Takeda, N.; Tsutsui, F. Thermal Performance Evaluation of Space Radiator for Single-Phase Mechanically Pumped Fluid Loop. J. Spacecr. Rockets 2022, 59, 225–235. [Google Scholar] [CrossRef]
  7. Paris, A.D.; Bhandari, P.; Birur, G.C. High Temperature Mechanically Pumped Fluid Loop for Space Applications -Working Fluid Selection. SAE Tech. Pap. 2004, 113, 892–898. [Google Scholar]
  8. Van Es, J.; van den Berg, T.; van Vliet, A.; Ganzeboom, T.; Brouwer, H.; Elvik, S. Mini Mechanically Pumped Loop Modelling and Design for Standardized CubeSat Thermal Control. In Proceedings of the 50th International Conference on Environmental Systems (ICES 2021), Lisbon, Portugal, 12–15 July 2021. [Google Scholar]
  9. Wang, J.X.; Li, Y.Z.; Zhang, Y.; Li, J.X.; Mao, Y.F.; Ning, X.W. A Hybrid Cooling System Combining Self-Adaptive Single-Phase Mechanically Pumped Fluid Loop and Gravity-Immune Two-Phase Spray Module. Energy Convers. Manag. 2018, 176, 194–208. [Google Scholar] [CrossRef]
  10. Wang, J.X.; Li, Y.Z.; Zhang, H.S.; Wang, S.N.; Liang, Y.H.; Guo, W.; Liu, Y.; Tian, S.P. A Highly Self-Adaptive Cold Plate for the Single-Phase Mechanically Pumped Fluid Loop for Spacecraft Thermal Management. Energy Convers. Manag. 2016, 111, 57–66. [Google Scholar] [CrossRef]
  11. Lee, K.-L.; Li, Y.; Guzek, B.J.; Kadambi, J.R.; Kamotani, Y. Compact Heat Rejection System Utilizing Integral Variable Conductance Planar Heat Pipe Radiator for Space Application. Gravitational Sp. Res. 2015, 3, 30–41. [Google Scholar] [CrossRef]
  12. Lee, K.L.; Tarau, C.; Anderson, W.G.; Beard, D. Titanium-Water Heat Pipe Radiators for Space Fission Power System Thermal Management. Microgravity Sci. Technol. 2020, 32, 453–464. [Google Scholar] [CrossRef]
  13. Lin, M.; Mou, J.; Chi, C.; Hong, G.; Ge, P.; Hu, G. A Space Power System of Free Piston Stirling Generator Based on Potassium Heat Pipe. Front. Energy 2020, 14, 1–10. [Google Scholar] [CrossRef]
  14. Tarau, C.; Anderson, W.G.; Miller, W.O.; Ramirez, R. Sodium VCHP with Carbon-Carbon Radiator for Radioisotope Stirling Systems. AIP Conf. Proc. 2010, 1208, 42–54. [Google Scholar]
  15. Anderson, W.G.; Tarau, C. Variable Conductance Heat Pipes for Radioisotope Stirling Systems. AIP Conf. Proc. 2008, 969, 679–688. [Google Scholar]
  16. Hertzberg, A. Thermal Management in Space. Available online: https://space.nss.org/settlement/nasa/spaceresvol2/thermalmanagement.html (accessed on 21 January 2021).
  17. Antipov, Y.A.; Smirnov, S.V.; Oshchepkov, P.P.; Khalife, H.S. Design Features of a Power Plant Based on a Stirling Engine Working on the Moon. Adv. Astronaut. Sci. 2021, 174, 833–847. [Google Scholar]
  18. Badger, B. Ammonia: The Refrigerant of the Future. Available online: https://www.energy-learning.com/index.php/opinion/112-ammonia-the-refrigerant-of-the-future (accessed on 21 January 2021).
  19. Smirnov, S.V.; Sinkevich, M.V.; Antipov, Y.A.; Khalife, H.S. A Calculation Method of a Heat Rejection System in a Lunar Power Plant Consisting of a Free-Piston Stirling Engine (FPSE). Acta Astronaut. 2021, 180, 46–57. [Google Scholar] [CrossRef]
  20. Microgen Engine Corporation—Microgen. Available online: https://www.microgen-engine.com/ (accessed on 5 February 2022).
  21. Evogress—Innovative Power Generation Technologies. Available online: http://www.evogress.com/en/ (accessed on 5 February 2022).
  22. Ivanov, V.L.; Fisher, Y.V.; Korniychuk, S.P.; Kisilev, N.A. Space Radiant Heat Exchanger for Heat Elimination. Sci. Educ. 2011, 10, 12. [Google Scholar]
  23. Hesselgreaves, R.E.; Law, D.R. Compact Heat Exchangers; Butterworth-Heinemann: Oxford, UK, 2016. [Google Scholar]
  24. Tsvetkov, F.F. Problem Book in Heat and Mass Transfer; Publishing House MEI: Moscow, Russia, 2008. [Google Scholar]
  25. Petukhova, B.S.; Shilkov, V.K. Handbook of Heat Exchangers; Energoatomizdat: Moscow, Russia, 1987. [Google Scholar]
  26. Williams, J.P.; Paige, D.A.; Greenhagen, B.T.; Sefton-Nash, E. The Global Surface Temperatures of the Moon as Measured by the Diviner Lunar Radiometer Experiment. Icarus 2017, 283, 300–325. [Google Scholar] [CrossRef] [Green Version]
  27. What Is the Temperature on the Moon?|Space. Available online: https://www.space.com/18175-moon-temperature.html (accessed on 12 June 2021).
  28. Golubev, I.F.; Kiyashova, V.P. Thermophysical Properties of Ammonia; Publishing House of Standards: Moscow, Russia, 1978. [Google Scholar]
  29. Thermophysical Properties of Fluid Systems. Available online: https://webbook.nist.gov/chemistry/fluid/ (accessed on 12 June 2021).
  30. Mikheev, M.A.; Mikheeva, I.M. Basics of Heat Transfer; Energia: Moscow, Russia, 1977. [Google Scholar]
Figure 1. Overview scheme of the lunar power plant.
Figure 1. Overview scheme of the lunar power plant.
Symmetry 14 01168 g001
Figure 2. Two options of the radiator design.
Figure 2. Two options of the radiator design.
Symmetry 14 01168 g002
Figure 3. Dependence of the number of pipes on the Reynolds number (Refrigerant—He).
Figure 3. Dependence of the number of pipes on the Reynolds number (Refrigerant—He).
Symmetry 14 01168 g003
Figure 4. Dependence of the pipe length and the total width of the radiator on the flow regime (Refrigerant—He).
Figure 4. Dependence of the pipe length and the total width of the radiator on the flow regime (Refrigerant—He).
Symmetry 14 01168 g004
Figure 5. Dependence of the pumping power of the refrigerant through the radiator on the flow regime (Refrigerant—He).
Figure 5. Dependence of the pumping power of the refrigerant through the radiator on the flow regime (Refrigerant—He).
Symmetry 14 01168 g005
Figure 6. Dependence of the number of pipes on the Reynolds number (Refrigerant—NH3).
Figure 6. Dependence of the number of pipes on the Reynolds number (Refrigerant—NH3).
Symmetry 14 01168 g006
Figure 7. Dependence of the pipe length and the total width of the radiator on the flow regime (Refrigerant—NH3).
Figure 7. Dependence of the pipe length and the total width of the radiator on the flow regime (Refrigerant—NH3).
Symmetry 14 01168 g007
Figure 8. Dependence of the pumping power of the refrigerant through the radiator on the flow regime (Refrigerant—NH3).
Figure 8. Dependence of the pumping power of the refrigerant through the radiator on the flow regime (Refrigerant—NH3).
Symmetry 14 01168 g008
Figure 9. Dependence of the number of pipes on the pipe diameter (Refrigerant—NH3).
Figure 9. Dependence of the number of pipes on the pipe diameter (Refrigerant—NH3).
Symmetry 14 01168 g009
Figure 10. Dependence of the pipe length and the total width of the radiator on the pipe diameter (Refrigerant—NH3).
Figure 10. Dependence of the pipe length and the total width of the radiator on the pipe diameter (Refrigerant—NH3).
Symmetry 14 01168 g010
Figure 11. Dependence of the pumping power of the refrigerant through the radiator on the pipe diameter (Refrigerant—NH3).
Figure 11. Dependence of the pumping power of the refrigerant through the radiator on the pipe diameter (Refrigerant—NH3).
Symmetry 14 01168 g011
Table 1. The main parameters of the FPSE.
Table 1. The main parameters of the FPSE.
ParameterValueUnit
Output Power1100W
Total thermal power3600W
Amount of heat to be rejected2500W
Efficiency30%
Maximum cycle temperature600°C
Minimum cycle temperature50°C
Table 2. Comparative analysis for liquid ammonia (Re = 500; d = 0.002 m).
Table 2. Comparative analysis for liquid ammonia (Re = 500; d = 0.002 m).
ParameterUnit2D ModelEnhanced MethodDeviation, %
Pipe lengthm1.8211.81.15
Total area of the radiatorm27.217.150.83
Power for pumping refrigerant through the radiatorW0.00890.0091.12
Table 3. Comparative analysis for liquid ammonia (Re = 4500; d = 0.0055 m).
Table 3. Comparative analysis for liquid ammonia (Re = 4500; d = 0.0055 m).
ParameterUnit2D ModelEnhanced MethodDeviation, %
Pipe lengthm6.396.252.1
Total area of the radiatorm27.287.112.3
Power for pumping refrigerant through the radiatorW0.03710.0382.4
Table 4. Comparative analysis for helium (Re = 500; d = 0.002 m).
Table 4. Comparative analysis for helium (Re = 500; d = 0.002 m).
ParameterUnit2D ModelEnhanced MethodDeviation, %
Pipe lengthm0.8440.840.7
Total area of the radiatorm27.20367.150.74
Power for pumping refrigerant through the radiatorW5.095.080.27
Table 5. Comparative analysis for helium (Re = 4500; d = 0.002 m).
Table 5. Comparative analysis for helium (Re = 4500; d = 0.002 m).
ParameterUnit2D ModelEnhanced MethodDeviation, %
Pipe lengthm1.4731.460.88
Total area of the radiatorm27.22427.151.03
Power for pumping refrigerant through the radiatorW217.3795222.62.4
Table 6. Comparative analysis for helium (Re = 4500; d = 0.0102 m).
Table 6. Comparative analysis for helium (Re = 4500; d = 0.0102 m).
ParameterUnit2D ModelEnhanced MethodDeviation, %
Pipe lengthm7.97.830.9
Total area of the radiatorm26.376.31
Power for pumping refrigerant through the radiatorW7.377.74.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Smirnov, S.; Sinkevich, M.; Antipov, Y.; Tsarkov, I.; Kupreev, S.; Khalife, H. An Enhanced Calculation Method of the Heat Rejection System of a Free-Piston Stirling Engine (FPSE) Operating on the Moon. Symmetry 2022, 14, 1168. https://doi.org/10.3390/sym14061168

AMA Style

Smirnov S, Sinkevich M, Antipov Y, Tsarkov I, Kupreev S, Khalife H. An Enhanced Calculation Method of the Heat Rejection System of a Free-Piston Stirling Engine (FPSE) Operating on the Moon. Symmetry. 2022; 14(6):1168. https://doi.org/10.3390/sym14061168

Chicago/Turabian Style

Smirnov, Sergey, Mikhail Sinkevich, Yuri Antipov, Igor Tsarkov, Sergei Kupreev, and Hassan Khalife. 2022. "An Enhanced Calculation Method of the Heat Rejection System of a Free-Piston Stirling Engine (FPSE) Operating on the Moon" Symmetry 14, no. 6: 1168. https://doi.org/10.3390/sym14061168

APA Style

Smirnov, S., Sinkevich, M., Antipov, Y., Tsarkov, I., Kupreev, S., & Khalife, H. (2022). An Enhanced Calculation Method of the Heat Rejection System of a Free-Piston Stirling Engine (FPSE) Operating on the Moon. Symmetry, 14(6), 1168. https://doi.org/10.3390/sym14061168

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop