Mechanism of heat transfer in granular media.
High power target systems require reliable high-temperature heat transfer media. Dense granular flow materials have potential to work as heat transfer media, especially at high temperature. In this chapter, a brief review of dense granular and their heat transfer properties is introduced, including basic concepts of heat transfer, thermal behavior of these materials, and factors that affect this behavior. The implementation of these materials as targets is addressed, where two targets were designed, constructed, and tested based on the concept of dense granular flow. The results of the application of a hopper flow-type target in 2.5 MW accelerator-driven system (ADS) and a chute flow-type target in the material irradiation facility will be presented.
- dense granular flow
- high power target
- prototype facility
In nuclear area, an accelerator is always helpful to generate particles uncommon or indiscrete in the natural world. Target materials are exposed to certain projectiles; as a result, secondary particles in the form of neutrons, photons, neutrino, and so on are generated.
For most cases, the higher particle flux outcome is more helpful. To achieve this, the target systems are designed to stand high beam intensity. Most of beam power deposits in the target material. The kinetic energy of the beam converts into internal energy and increase the temperature of target. This process also leads to radiation damage of target material and even target failure.
Lots of effort has been made to achieve high power target. But not until the 1990s, a specific megawatt target never came up. As examples, Weapons Neutron Research project in Los Alamos Neutron Science Center achieved 70 kW target with tungsten material, while Rutherford Appleton Laboratory reached 160 kW using 800 MeV protons and tungsten with tantalum cladding .
As mentioned above, cooling is essential to the high power target. Conventional cooling channels take too much space inside the target body and gradually get close to the limit. The MarkIII target of Swiss Spallation Neutron Source project in Paul Scherrer Institute (PSI) reaches 1.2 MW beam power based on such concept. Some projects adopt rotating target to even and lower the power density within the target body, such as European Spallation Source and China Spallation Neutron Source [1, 2, 3, 4].
For solid target, radiation material damage also causes swelling and creep, which might damage the cooling channel and cause cooling deterioration, as the displacement per atom (DPA) increases. One way to bypass the damage is using liquid metal as target material. Meanwhile, the flow of liquid metal removes the power deposit and dumps the heat in the heat exchanger. These properties make liquid metal target a promising solution of high power target. In PSI, MEGAwatt PIlot Experiment (MEGAPIE) achieved 0.8 MW with lead bismuth eutectic (LBE). Spallation Neutron Source in US and Japan Spallation Neutron Source also achieved 1.4 and 1 MW beam power with mercury .
Here we are going to present a third option of target cooling, which is based on the dense granular flow . Granular system makes solid materials flow like fluid and perform heat removing. Meanwhile, the candidate target materials are more various, which give better feasibility depending on the purpose of the facility.
2. Heat transfer media and heat transfer in granular material
2.1. Conventional heat transfer media
The continuous research of the thermal properties of materials is explained both by practical needs and fundamental science. Heat removal is always a crucial issue for energy conversion in the energy industry owing to increased levels of power requirement and energy efficiency. The search for materials that transfer heat well under high temperature and energy flux has covered all the industry area and never ends, and it is especially essential for design of the next generation of reactors and neutron facilities.
A medium’s transfer heat ability follows three basic heat transfer ways, conduction, radiation, and convection. Conduction dominates solids’ heat transfer behavior, while convection makes fluids absorb heat from a source and dump heat to a sink very efficiently. In field of nuclear energy, water is still the most important heat transfer material for cost and reliability reasons, while liquid metal such as sodium and gases such as helium take the rest portion. These fluid heat transfer media are capable of removing 102–105W/m2K from surface of power component and reaching up to the order of ~106W/m2K while boiling . At high temperature, some particular situations differ because of radiation. Radiation transfers heat from one surface to another through void following the Stephan-Boltzmann law, which indicates the radiation heat transfer rate is proportional to the cubic of temperature. Hence, if radiation happens effectively in the media, heat transfer can be more efficient at high temperature.
At high temperature, thermodynamics principles allow fluid’s thermal energy transfers more efficiently. But most practical heat transfer media failed under such temperature: Liquid water transforms to supercritical state over 385°C [7, 8]; steam reacts with common structural materials over 900°C , and about one-third of the zircaloy in the core create a sizable hydrogen bubble from the Zr-steam interactions ; liquid sodium or lithium, Na-K alloy, and mercury are all vaporized under 900°C; lead-bismuth eutectic, liquid lead’s corrosion-erosion effect grows significant at 480–550°C [11, 12]. What’s more, for nuclear usage, liquid metal loop needs multiple conditioning systems, which makes the system complicated and expensive, such as International Fusion Materials Irradiation Facility (IFMIF)  and eXperimental Accelerator-Driven System (XT-ADS) device . Gases seem to be another good choice, which has been testified in high-temperature gas-cooled reactors. But they have to be pressurized due to their low volume heat capacity .
2.2. Characteristic of granular flow
Granular material is a promising solution. Granular material is composed of discrete solid particles and gaseous environment where the particles immerse. Such material has been used for heating since ancient Rome by Lucretius . Modern research can be trace back to Maxwell’s era and got broad attentions in the twentieth century due to popularization of bed reactors . Nowadays, granular heat exchanger has been widely used in many industrial areas. They perform effective heating or cooling in so many industrial processes that indicate granular materials are realistic to be a promising heat transfer media. Packed bed and powder jet solution had been came up by research groups to solve the high power density removal issues of target [18, 19]. Pebble bed fusion reactors and fission first wall breeding layers also adopt granular material using the unique advantages in the high-temperature heat transfer [20, 21, 22, 23]. Flow of granular state enhances heat transport process, leading to better heat transfer performance and more realistic feasibility, especially for extreme internal heat source such as spallation reaction, compared to present solutions.
Its potential of heat transfer usage stems from flowability and thermal properties, which has been widely studied and practiced. The flowability allows granular materials to convey long distance by multiple methods. But granular materials’ flowability is quite different from fluids’. Although granular dynamics uses many hydrodynamics concepts, granular dynamics behaves quite different from fluid. Some of granular materials’ unique properties are critically helpful to be heat transfer media.
When granular materials flow out from a vessel under the action of gravity, the mass flow rate M is approximately independent of the head of material head H, just as an hourglass does. The most important research of this phenomenon was successively recorded by Beverloo et al., Aldin and Gunn, and Nedderman [24, 25, 26, 27, 28]. This phenomenon requires that the granular material must be not too smooth and small, while the bulk has a height not less than twice orifice diameter. A quantitative result of this flow gives a relationship of .
where m, ρs, D, dp, θ, and β represent mass flow rate, particle material density, orifice diameter, particle diameter, wall angle, and an angle related to friction of the material . This result implies the flow rate could be controlled by local valve.
As the immerging depth increase, the normal pressure of most fluids increases linearly and is proportional to the depth. But the normal pressure of granular system gradually reaches a constant. Different from fluids, inside the granular system and between the granules and the container walls, friction exists. The upward friction force exerted by the wall neutralizes the weight of grains. An important quantitative explanation given by Janssen is in accord with the experimental stress result .
Other granular properties also have special effects, which stem from the friction, dissipation, and size effect of granules. These effects lead to phenomenons in granular flow, such as the jamming and dilatancy properties prevent flow channel collapse ; while dissipation property prevents long-range disturbance propagation .
2.3. Thermal properties of granular flow
Granular materials’ thermal properties are decided by the components which make the effective properties vary largely. As high-temperature heat transfer media, high-melting-point metal alloy, ceramics, and other materials with good heat conductivity can be the solid particle candidate. There are relatively few gases performing heat transfer effectively, among which helium is representative. The large number of material combination makes granular media satisfy varieties of actual needs. For instance, the heat transfer media in traditional nuclear reactor need to minimize the neutron absorption and maximize the neutron moderation, while the heat transfer media had better to be the spallation material in a spallation target which needs high spallation neutron production.
Some effective thermophysical properties such as critical temperature, density, and specific heat can be estimated using simple calculations based on the properties of the initial components. But quantities such as heat transfer coefficient and effective heat conductivity are far more complicated than those. The following subsections will review critical thermal properties to get a better view of granular heat transfer media’s potential.
The mainstream viewpoint of heat transfer between granular material and a panel has two stages: heat transfers from the wall surface to the first granular layer next to the wall, and heat diffuses from the first layer to the inside of granular material [33, 34, 35, 36]. Sullivan and Sabersky  adopted concepts of fluid heat transfer such as Peclet and Nusselt number dealing with this process. An analytical result of temperature distribution can be derived from principles of heat transfer. This work was further extended by Spelt et al. [36, 37].
For both two stages, whatever discrete model or continuous model is applied, one should realize that the structure of granular material leads to more heat transfer mechanisms. As Yagi and Kunii reviewed, heat transfer in granular media happens in all three ways and their coupled effects . Most concerned mechanisms of heat transfer in granular media are listed in Table 1.
The first stage was observed by several early researchers, depicted as a sharp temperature drop at the wall [36, 39, 40]. In this stage, the maximum heat transfer is at the order of ~103W/m2K for granules’ diameter at the order of millimeter in air. Such order is determined by different researchers and methods [39, 41]. When the pressure drops to 0.1 Pa, the wall heat transfer coefficient decreases to the order of ~10 W/m2 K. Most significant heat transfer coefficient drop happens when pressure changes from ~1 kPa to 1 Pa, leading to the change of one and a half order.
In vacuum, the solid dominates the heat transfer rate, proved by the coefficient of bronze and glass, which are 30 and 5 W/m2 K, respectively. This result agrees well with the contact heat transfer analysis, and we will discuss this effect in the effective thermal conductivity of granular material in later section. In contrast, the gas phase decides the heat transfer in atmosphere, reaching the order of ~103W/m2K [35, 41]. It is believed that the porosity near the wall reaches almost 1 [42, 43, 44, 45], which means only little solid contact at the wall leads to a minor heat transfer mechanism only significant in vacuum.
Although most researches did not evaluate the specific limit of heat transfer rate, they did observe effect of various factors. Most of them have relationship with the second stage. In this stage the heat transport can be considered as heat diffusion in a continuous media with effective properties of granular material. The complexity of heat transfer mechanism makes multiple parameters dependence.
2.3.2. Effect of material constituent
Granular material component combination is the most principle factor of heat transfer properties. In the past century, a number of various granules’ heat transfer properties have been studied. A collection of experimental data is shown in Table 2. It shows gas phase plays a decisive role in granular materials’ heat transfer property, while the solid particle geometry plays a role as important as its inherent properties. Large surface-volume ratio makes the temperature of fluid and solid phase nearly uniform [46, 47].
In energy area, most heat transfer is performed under pressure. Pressurization enhanced heat transfer in two aspects: increase boiling temperature and volume heat capacity. High-melting-point solid granules need none of these improvements and hence may significantly reduce engineering difficulty, risk and cost. The pressure concept in granular system also has two aspects: the gas pressure and particle compress pressure.
Early stage of ETC of granular material focuses on experimental determination. Back to 1926, Aberdeen and Laby studied effective thermal conductivity (ETC) of monox powder when immersed in air, in carbon dioxide, and in hydrogen at various pressures has been determined. They found a linear relation between the conductivity of the powder and the logarithm of the pressure of the gas in which it is immersed. The result fits in
where k is effective heat transfer coefficient and kf is fluid heat transfer coefficient and
2.3.3. Effect of temperature
Microscope heat transfer mechanism decides most condensed matters’ thermal conductivity deteriorates as temperature increases while gases are contrary. In a granular material, solid phase has more significant heat conductivity than gas phase. On the other hand, the packed solid particle structure is point contact for most cases; hence, the gas voids play an important role in the effective heat conduct which increases as temperature rises. Radiation also significantly improves effective heat conductivity of granular material under high temperature. Solid surfaces and voids allow radiation to perform efficient heat transfer, while such mechanism does not exist in bulk solid or pure gas.
Yagi and Kunii  determined ETC for several materials, in the temperature range of 100–800°C. The result shows an accelerated increase of ETC as the temperature increases. Most materials’ ETC doubled or tripled for the test range. They proposed a model considering the effect of radiation . As Botterill published in 1989, the radiation is proportional to T3, σ, and the character length represented by Dp; hence a radiation term of ETC falls in the form of
where σ is Stefan-Boltzmann coefficient and χ is coefficient depending on the specific model. A comparative research of various radiation models and experimental results was also given by Botterill .
In a word, unlike traditional heat transfer material, granular material could perform better under high temperature; as previous study of gas cooling reactor, heat conductivity of graphite pebble bed exceeds pure graphite for temperature higher than 1400°C . Further experimental results were reported in [38, 49, 50].
2.3.4. Effect of gas motion
As is well known, forced convection enhanced heat transfer properties. Packed bed is motionless granules with gas flows through. Researches of these systems explain the motion effect of gas phase, which also implies the same effect in granular flow.
Bunnell et al. studied a wall-heated packed bed with different superficial gas velocity, measured the inner temperature distribution, and calculated ETC through these results. They found effective thermal conductivity dependence on mass velocity of the gas, and the gas phase and solid phase temperatures are nearly uniform .
Solid particle motion enhanced cold and hot particles mixing, hence increasing the heat transfer rate. A lot of research about packed bed also showed that only fluid flow rate increase also leads to enhancement of ETC [46, 50, 52, 53, 54, 55].
2.3.5. Models of ETC
To get a quick evaluation of granular materials, researchers had come up with various models, from simple to complex. Three most important models are as follows:
This model is probably the earliest method for granular ETC evaluation. It represents a series model that considers the system as a parallel-serial thermal resistance system  which got its popularity in the early age.
Yagi and Kunii :
in which , , . hrs, hrv, lp, ls, and lv represent radiation heat transfer coefficient of solid to solid, void to void, effective length of centers, thermal conduction, and fluid film at the contact point. They considered effect of fluid flow using Ranz fluid mixing equation and shows the particle size effect. The result fits temperature effect well .
Zehner and Schlunder's model:
in which B is a deformation factor. This model represents a significant improvement of contact term which applies for various material properties.
2.3.6. Total heat transfer coefficient
Schlunder gave explanation of the combination of wall heat transfer and ETC as a serial of heat resistance as follows :
Because the heat diffusive process inside the granules varies as time changes, total heat transfer coefficient also depends on time. Most measured total heat transfer coefficient is time averaged value:
In three situations, (1) t goes to 0, (2) bed heat capacity is infinity, and (3) the bed is under perfect mixing, αETC reaches a maximum and hence derives the maximum heat transfer coefficient, i.e., the wall heat transfer coefficient. A typical time dependence of αETC based on Fourier theory is given by
The total heat transfer depends on time, hence given by.
Through the experiments of various components’ combination, Ernst showed α decreases as time increases, which is less significant for lower pressure and larger granules . Such time-dependent phenomenon has been observed by other researches [36, 56]. The total heat transfer coefficient fell in the range of 50–200 W/m2·K which is also a typical range of most facility situation .
3. Implementation of granular flow target
With the concept of dense granular flow cooling, the target can be designed according to different demands. Two different types of granular flow target have been developed for China Initiative Accelerator-Driven System (CIADS) project and China material irradiation facility project. For both designs helium was selected as the cover gas, concerning its high thermal conductivity and safety. In principle, pressure was set lower than 1 atm to prevent leakage and help target-beam line coupling. For each case, gas pressure was variable parameter. The difference of the two targets includes structure, target material, beam energy, spectrum, and so on.
By the end of 2017, two prototype facilities were constructed and tested with electrical heating, up to about 300°C. The beam coupling experiments are planned to be performed later, depending on the progress of radiation shielding hall.
3.1. ADS target
For CIADS, the target needs to prove high neutron yield inside a subcritical reactor core. To replenish the neutrons in a subcritical power reactor, the target-beam power is very high. The primary design chose proton beam of 500 MeV @ 5 mA and beam power of 2.5 MW.
Meanwhile, the space for the target installation is very limited to lower the neutron leakage and increase the effective neutron multiplication factor (keff) of the reactor. Usually the installation space is in vacancy of several removed fuel components. In CIADS, the diameter of the target tube needs to be under 30 cm. To fit this space, a thin and long vertical tube was applied as target tube. The dense granular flow went downward within. The beam line installed coaxially in the target tube sent protons deep into the target granules and reacted at the same level of the core center. Figure 1 shows a schematic view of CIADS, which adopts a loop target section.
Tungsten was selected as the target material. Minor iron was added to optimize the mechanical properties. The alloy was made into 1 mm sphere. Tungsten generates abundant neutrons in spallation reaction which is widely adopted in neutron resources. With the basic idea of dense granular flow, the tungsten also plays a role as coolant. Considering the beam heat load and material’s temperature endurance, the circulation mass flow is required about 200 kg/s. Compared to the fluidized dilute granular flow developed in the Rutherford Appleton Laboratory, we believed that the dense flow is more stable and lowers the requirement of gas system, for this specific case.
To sustain the operation, the system also included conveying, cooling, conditioning, monitoring, cleaning, and other affiliate systems. To lift the spheres back into the target, a chain conveyor is employed. The spheres inside flow vertically under gravity. When the spheres flow to the bottom of the loop, the elevator lifts them back to the top, and the spheres reenter into the target. The elevator actively supplies the energy to sustain the circulation; hence the reliability and feasibility should be carefully considered. Other confines came from the ADS application, such as airtight and compact. Due to the grain size and density, pneumatic system does not fit the system well. Both conventional mechanical conveyer and a new developed magnetic conveyer are chosen as potential candidate. The heat of spheres will be removed in a panel heat exchanger. Through the steel shell of the heat transfer panel, the heat is transferred to the counter water flow. The width of parallel flow channel formed by the panels is large enough to avoid arching and clogging. Figure 2 shows an electron beam coupling test facility constructed in 2015. Without conditioning and some other minor systems, the facility is simplified to test the circulation and heat loading effect. As shown in Figure 3, the coupling section is about 1:5 to the full scale.
The full-scale testing was setup later in 2017, as shown in Figure 4. The model facilities were put into practice as a preliminary design verification by Institute of Modern Physics, Lanzhou. Other than the components of previous system, sieving and cover gas systems were installed. Conveying system was updated to a magnet conveyer, and heat exchange test was elevated to about 250°C.
3.2. Target for material irradiation facility
Material irradiation facility also requires intensive neutron source to test material tolerance of radiation damage. The China material irradiation facility aims at a small-scale and low-cost facility. To maximize the radiation effect, the facility design took the advantage of forward neutrons. The system adopted deuterium beam of 20 MeV@5 mA and beryllium as the target material. Due to 20 MeV deuterium only has a very low penetration at the order of millimeter; there is hardly effective neutrons leak out, with the previous configuration. Hence, a new structure based on chute flow was adopted.
To demonstrate the heat load effect of this chute flow target, this progress is simulated by our GPU code . Figure 5 showed one of the simulation results of heat deposition in a chute flow. The result showed the 20 MeV beam particles only penetrate several millimeters into the flow. According to the design and simulation, mass flow at about 2–3 m/s is sufficient to remove the heat. The parameters of feeding rate, chute slope, beam spot position, and flow channel width system could be well adjusted to required reaction region. After its first test, target section of the facility in Figure 2 was replaced by a slope chute. Figure 6 showed the test loop facility and targets applied in the system. An acrylic chute was chosen to observe the flow status, and a stainless steel chute was installed for beam coupling effect of the facility. The experimental result is showed lower right.
The target design is first tested on remodeled small-scale circulation loop for ADS target in 2015. The facility is also radiated by electron beam. Temperature of the flow increased slightly, which showed that the system heat removal is workable.
Later in 2017, a full-scale prototype facility finished its major construction, as shown in Figure 7. Due to the radiation hazards, instead of actual beam, system used laser and electric heating to test temperature endurance and heat removal effect. 30 kW electric heating offered major overall heating power, while laser provided a high power density up to 3 MW/cm2. Granular flow reached an average temperature more than 300°C and maximum temperature more than 500°C. High-temperature testing was performed over 200 h, among which the longest continuous run lasts more than 100 h.
To test the neutron performance of this target material, a static granular target as shown in Figure 8 was installed in the beamline of a 20 MeV linear accelerator. The neutron yield data is tested through proton beam and the same 1 mm beryllium in static vacuum chamber. Five stacks of activated foils are located around the chamber. The measurement of the foils deduced a total neutron flux of 2.14×1013/cm2·mA·s.
The unique property of dense granular flow makes it possible to transfer heat effectively at high temperature. The major mechanism and factors of granular were reviewed. Using the granular material as heat transfer media and target material, two different configurations came up depending on the application. The prototype testing facilities were constructed. More designs and applications based on the granular heat transfer media are to be developed hopefully.
The major work mentioned above is achieved by the staff members from spallation target department and advance nuclear material department of Institute of Modern Physics of the CAS.
Bauer G. Overview on spallation target design concepts and related materials issues. Journal of Nuclear Materials. 2010; 398(1-3):19-27
Buligins L, Thomsen K, Lielausis O, et al. Internal geometry and coolant choices for solid high power neutron spallation targets. Nuclear Instruments & Methods In Physics Research Section A-Accelerators Spectrometers Detectors And Associated Equipment. 2014; 761:58-68
Hao J, Chen Q, Xu Y, et al. Target thickness optimization design of a spallation neutron source target cooling system. Applied Thermal Engineering. 2013; 61(2):641-648
Bauer Go OECD. Towards reliable high power spallation target design. 2003. pp. 81-93
Yang L, Zhan W. New concept for ADS spallation target: Gravity-driven dense granular flow target. SCIENCE CHINA Technological Sciences. 2015; 58(10):1705-1711
Incropera F. Fundamentals of Heat and Mass Transfer. John Wiley & Sons; 2011
Lehr J. Introduction. In: Nuclear Energy Encyclopedia. Hobeken, NJ, USA: John Wiley & Sons, Inc; 2011. pp. xi-xii
Wolf J. Supercritical water reactor. In: Nuclear Energy Encyclopedia. Hobeken, NJ, USA : John Wiley & Sons, Inc.; 2011. pp. 305-308
Suzuki M, Kawasaki S. Development of computer code PRECIP-II for calculation of Zr-steam reaction: Comparison of calculation with experiments in temperature transient. Journal of Nuclear Science and Technology. 1980; 17(4):291-300
Marques JG. Safety of nuclear fission reactors: Learning from accidents. In: Nuclear Energy Encyclopedia. Hobeken, NJ, USA: John Wiley & Sons, Inc.; 2011. pp. 127-149
Gromov B et al. Use of lead-bismuth coolant in nuclear reactors and accelerator-driven systems. Nuclear Engineering and Design. 1997; 173(1-3):207-217
Hill R, Grandy C, Khalil H. Generation-IV sodium-cooled fast reactors (SFR). In: Nuclear Energy Encyclopedia. John Wiley & Sons, Inc.; 2011. pp. 353-364
Martone M. IFMIF: International Fusion Materials Irradiation Facility Conceptual Design Activity: Final report. ENEA, Frascati, Italy: Dipt. Energia; 1997
De Bruyn D et al. From Myrrha to XT-ADS: The design evolution of an experimental ads system. AccApp'07; 2007. pp. 848-854
Taylor P. Generation-IV gas-cooled fast reactor. In: Nuclear Energy Encyclopedia. John Wiley & Sons, Inc.; 2011. pp. 349-351
Faraday M. On a peculiar class of acoustical figures; and on certain forms assumed by groups of particles upon vibrating elastic surfaces. Philosophical Transactions of the Royal Society of London. 1831; 121:299-340
Duran J. Sands, Powders, and Grains—An Introduction to the Physics of Granular Materials. 2000
Davies T et al. The production and anatomy of a tungsten powder jet. Powder Tech-nology. 2010; 201(3):296-300
Magistris M. Radiological considerations on multi-MW targets. Part II: After-heat and temperature distribution in packed tantalum spheres. Nuclear Instruments & Methods in Physics Research Section A—Accelerators Spectrometers Detectors and Associated Equipment. 2005; 545(3):823-829
Colominas S et al. Octalithium plumbate as breeding blanket ceramic: Neutronic performances, synthesis and partial characterization. Fusion Engineering and Design. 2012; 87(5-6):482-485
Donne M, Sordon G. Heat-transfer in pebble beds for fusion blankets. Fusion Technology. 1990; 17(4):597-635
Rousseau P, du Toit C, Landman W. Validation of a transient thermal-fluid systems CFD model for a packed bed high temperature gas-cooled nuclear reactor. Nuclear Engineering and Design. 2006; 236(5-6):555-564
Reimann J et al. Measurements of the Thermal Conductivity of Compressed Beryllium Pebble Beds. 2005
Herbst T. Shear phenomena in granular random packings. Soil Science. 1964; 98(4):280
Wieghardt K. Experiments in granular flow. Annual Review of Fluid Mechanics. 1975; 7:89-114
Beverloo W, Leniger H, Vandevelde J. The flow of granular solids through orifices. Chemical Engineering Science. 1961; 15(3-4):260
Aldin N, Gunn D. The flow of non-cohesive solids through orifices. Chemical Engineering Science. 1984; 39(1):121-127
Nedderman R. Statics and Kinematics of Granular Materials. Cambridge University Press; 2005
Rao K, Nott P, Sundaresan S. An Introduction to Granular Flow. New York: Cambridge University Press; 2008
Gutiérrez G et al. Silo collapse: An experimental study. In: Traffic and Granular Flow’07. Springer; 2009. pp. 517-523
Hall J, Richart F. Dissipation of Elastic Wave Energy in Granular Soils. 1963
Broughton J, Kubie J. Heat-transfer mechanism as applied to flowing granular media. International Journal of Heat and Mass Transfer. 1976; 19(2):232-233
Balakrishnan A, Pei D. Heat transfer in gas-solid packed bed systems. 1. A critical review. Industrial & Engineering Chemistry Process Design and Development. 1979; 18(1):30-40
Schlunder E. Heat transfer to packed and stirred beds from the surface of immersed bodies. Chemical Engineering and Processing. 1983; 18:31-53s
Sullivan W, Sabersky R. Heat transfer to flowing granular media. International Journal of Heat and Mass Transfer. 1975; 18(1):97-107
Spelt J, Brennen C, Sabersky R. Heat transfer to flowing granular material. International Journal of Heat and Mass Transfer. 1982; 25(6):791-796
Yagi S, Kunii D. Studies on effective thermal conductivities in packed beds. AICHE Journal. 1957; 3(3):373-381
Ernst R. Wärmeübergang an Wärmeaustauschern im Moving Bed. Chemie Ingenieur Technik. 1960; 32(1):17-22
Schlunder E. Wärmeübergang an bewegte Kugelschüttungen bei kurzfristigem Kontakt. Chemie Ingenieur Technik. 1971; 43(11):651-654
Wunschmann J, Schlunder E. Heat transfer from heated plates to stagnant and agitated beds of spherical shaped granules under normal pressure and vacuum. In: Proc. 5th Int. Heat Transfer Conf.; Tokyo. 1974
Martin H. Low Peclet number particle-to-fluid heat and mass transfer in packed beds. Chemical Engineering Science. 1978; 33(7):913-919
Du Toit C. Radial variation in porosity in annular packed beds. Nuclear Engineering and Design. 2008; 238(11):3073-3079
de Klerk A. Voidage variation in packed beds at small column to particle diameter ratio. AIChE Journal. 2003; 49(8):2022-2029
Gnielinski V. Heat and mass-transfer in packed-beds. Chemie Ingenieur Technik. 1980; 52(3):228-236
Sissom L, Jackson T. Heat exchange in fluid-dense particle moving beds. Journal of Heat Transfer. 1967; 89(1):1
Chukhano Z. Heat and mass transfer between gas and granular material. 2. International Journal of Heat and Mass Transfer. 1970; 13(12):1805
Aberdeen J, Laby T. Conduction of heat through powders and its dependence on the pressure and conductivity of the gaseous phase. Proceedings of the Royal Society of London Series a-Containing Papers of a Mathematical and Physical Character. 1926; 113:459-447
Kannuluik W, Martin L. Conduction of heat in powders. Proceedings of the Royal Society of London Series A-Containing Papers of a Mathematical and Physical Character. 1933; 141(843):144-158
Sokolov V, Yablokova M. Thermal-conductivity of a stationary granular bed with upward gas-liquid flow. Journal of Applied Chemistry of the USSR. 1983; 56(3):551-553
Gelperin J, Kagan A. Study of heat transfer in gas flow in tubes packed with granular materials. Chemical Engineering Progress. 1966; 62(8):88
Zumbrunnen D, Viskanta R, Incropera F. Heat-transfer through granular beds at high-temperature. Warme Und Stoffubertragung-Thermo and Fluid Dynamics. 1984; 18(4):221-226
Jaguaribe E, Beasley D. Modeling of the effective thermal conductivity and diffusivity of a packed bed with stagnant fluid. International Journal of Heat and Mass Transfer. 1984; 27(3):399-407
Einav I et al. 10,000—A reason to study granular heat convection. In: Yu A et al., editors. Powders and Grains 2013; 2013. pp. 38-45
Botterill J et al. The effect of gas and solids thermal properties on the rate of heat transfer to gas-fluidized beds. In: Proc. Int. Symp. on Fluidization; Eindhoven. Amsterdam: Netherlands University Press; 1967
Baird M et al. Heat transfer to a moving packed bed of nickel pellets. Canadian Journal of Chemical Engineering. 2008; 86(2):142-150
Tian Y, Qi J, Lai J, et al. A heterogeneous CPU-GPU implementation for discrete elements simulation with multiple GPUs. In: 2013 International Joint Conference on Awareness Science and Technology and Ubi-Media Computing (i CAST-UMEDIA); 2013. pp. 547-552
Roberts I. Determination of the vertical and lateral pressures of granular substances. Proceedings of the Royal Society of London. 1883; 36(228-231):225-240
Yagi S, Kunii D, Wakao N. Studies on axial effective thermal conductivities in packed beds. AIChE Journal. 1960; 6(4):543-546