US20210310699A1 - Solar energy collection system with symmetric wavy absorber pipe - Google Patents
Solar energy collection system with symmetric wavy absorber pipe Download PDFInfo
- Publication number
- US20210310699A1 US20210310699A1 US17/219,735 US202117219735A US2021310699A1 US 20210310699 A1 US20210310699 A1 US 20210310699A1 US 202117219735 A US202117219735 A US 202117219735A US 2021310699 A1 US2021310699 A1 US 2021310699A1
- Authority
- US
- United States
- Prior art keywords
- absorber pipe
- pipe
- absorber
- longitudinal axis
- central longitudinal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S10/00—Solar heat collectors using working fluids
- F24S10/40—Solar heat collectors using working fluids in absorbing elements surrounded by transparent enclosures, e.g. evacuated solar collectors
- F24S10/45—Solar heat collectors using working fluids in absorbing elements surrounded by transparent enclosures, e.g. evacuated solar collectors the enclosure being cylindrical
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S10/00—Solar heat collectors using working fluids
- F24S10/70—Solar heat collectors using working fluids the working fluids being conveyed through tubular absorbing conduits
- F24S10/74—Solar heat collectors using working fluids the working fluids being conveyed through tubular absorbing conduits the tubular conduits are not fixed to heat absorbing plates and are not touching each other
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S10/00—Solar heat collectors using working fluids
- F24S10/70—Solar heat collectors using working fluids the working fluids being conveyed through tubular absorbing conduits
- F24S10/75—Solar heat collectors using working fluids the working fluids being conveyed through tubular absorbing conduits with enlarged surfaces, e.g. with protrusions or corrugations
- F24S10/755—Solar heat collectors using working fluids the working fluids being conveyed through tubular absorbing conduits with enlarged surfaces, e.g. with protrusions or corrugations the conduits being otherwise bent, e.g. zig-zag
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S23/00—Arrangements for concentrating solar-rays for solar heat collectors
- F24S23/70—Arrangements for concentrating solar-rays for solar heat collectors with reflectors
- F24S23/71—Arrangements for concentrating solar-rays for solar heat collectors with reflectors with parabolic reflective surfaces
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S23/00—Arrangements for concentrating solar-rays for solar heat collectors
- F24S23/70—Arrangements for concentrating solar-rays for solar heat collectors with reflectors
- F24S23/74—Arrangements for concentrating solar-rays for solar heat collectors with reflectors with trough-shaped or cylindro-parabolic reflective surfaces
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24S—SOLAR HEAT COLLECTORS; SOLAR HEAT SYSTEMS
- F24S80/00—Details, accessories or component parts of solar heat collectors not provided for in groups F24S10/00-F24S70/00
- F24S80/50—Elements for transmitting incoming solar rays and preventing outgoing heat radiation; Transparent coverings
- F24S80/54—Elements for transmitting incoming solar rays and preventing outgoing heat radiation; Transparent coverings using evacuated elements
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/40—Solar thermal energy, e.g. solar towers
- Y02E10/44—Heat exchange systems
Definitions
- PTCs parabolic trough solar collectors
- Various approaches have been proposed to achieve desired thermal and optical characteristics of PTCs encompassing optimization and enhancement of collector geometry, optical and material properties, heat transfer fluids, and collector operation.
- PTCs have used straight absorber pipes to collect solar energy.
- Various implementations include a solar thermal energy collection system that comprises a solar energy concentrator, a heat transfer fluid, and an absorber pipe through which the heat transfer fluid flows for absorbing concentrated solar energy incident on the absorber pipe from the solar energy concentrator.
- the absorber pipe comprises a pipe wall that extends about and around a central longitudinal axis that extends between each end of the absorber pipe.
- the pipe wall comprises an inner surface and an outer surface.
- the inner surface defines a first contour that defines alternating peaks and troughs along a length of the absorber pipe, and the outer surface defines a second contour that defines alternating peaks and troughs along the length of the absorber pipe.
- the inner surface defines the entire flow path for the heat transfer fluid through the absorber pipe.
- the first contour as viewed through an axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis.
- the sinusoidal waves on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis, and at least a portion of the absorber pipe is disposed along a focal line of the solar energy concentrator.
- the second contour as viewed through the axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis.
- the sinusoidal waves of the second contour on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis.
- the pipe wall has a maximum radius as measured between the inner surface and the central longitudinal axis at each peak and a minimum radius as measured between the inner surface and the central longitudinal axis at each trough.
- the outer surface of the pipe wall has a second maximum radius as measured between the outer surface and the central longitudinal axis at each peak and a second minimum radius as measured between the outer surface and the central longitudinal axis at each trough.
- the peaks and troughs of the inner surface are radially aligned with the peaks and troughs of the outer surface.
- a thickness between the inner surface and the outer surface of the wall is constant along the length of the wall.
- the inner surface is shaped to induce a plurality of turbulent fluid flow regions along the length of the central longitudinal axis.
- troughs of the second contour define areas of concentrated solar energy.
- the solar thermal energy collection system is a parabolic trough solar energy absorber system.
- the solar energy concentrator comprises a parabolic reflector having a focal axis that is parallel to a central longitudinal axis of the parabolic reflector and is spaced apart a focal length from a vertex of the parabolic reflector.
- the parabolic trough solar energy absorber system further comprises a tube into which the absorber pipe is disposed. The outer surface of the absorber pipe is spaced apart from an inner surface of the tube, and a wall of the tube allows for the transmission of radiant energy therethrough.
- the tube is coupled to the parabolic reflector such that the central longitudinal axis of the absorber pipe is coincident with the focal axis.
- the solar thermal energy collection system further comprises a heat transfer fluid tube coupled to an entry end and an exit end of the absorber pipe.
- the heat transfer fluid tube and the absorber pipe form a closed loop through which the heat transfer fluid flows.
- the system further comprises a heat exchanger, wherein at least a portion of the heat transfer fluid tube is disposed inside the heat exchanger.
- the heat exchanger is a boiler tank.
- the system further comprises a fluid pump coupled to the heat transfer fluid tube, and the fluid pump causes the heat transfer fluid to flow through the heat transfer fluid tube and the absorber pipe.
- FIG. 1A shows an end view of a schematic diagram of a parabolic trough solar energy collector system comprising an absorber pipe according to one implementation.
- FIG. 1B shows an axial cross-sectional view along the C-C line of the parabolic trough solar energy collector system comprising the absorber pipe shown in FIG. 1A .
- FIG. 1C shows a schematic diagram of the system shown in FIG. 1A .
- FIG. 2A shows a perspective view of an absorber pipe according to another implementation.
- the absorber pipe shown in FIG. 2A is similar to the absorber pipe shown in FIG. 1A but has six peaks.
- FIG. 2B shows a partial axial cross-sectional view of the absorber pipe shown in FIG. 2A .
- FIG. 3A shows a local concentration ratio around a quarter of a conventional straight absorber pipe.
- FIG. 3B shows a local concentration ratio around a quarter of the absorber pipe, such as the absorber pipes shown in FIG. 1B or 2A .
- FIG. 4 shows a local concentration ratio profile around the absorber pipe in FIG. 1B .
- FIG. 5 is a convergence plot for an implementation of the absorber pipe, such as the absorber pipes shown in FIGS. 1B and 2A , that has seven waves and an amplitude of 0.006 m.
- FIG. 6A is a validation plot for the absorber pipe, such as the absorber pipes shown in FIGS. 1B and 2A .
- FIG. 6B is a validation plot for a complete PTC model with a conventional absorber pipe.
- FIG. 7A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipes shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 1 kg/s.
- FIG. 7B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 1 kg/s.
- FIG. 8A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 5 kg/s.
- FIG. 8B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 5 kg/s.
- FIG. 9A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 10 kg/s.
- FIG. 9B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , and amplitude for a mass flow rate of 10 kg/s.
- FIG. 10 shows streamlines at the last pipe segment cross section of an absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , for different numbers of waves and different amplitudes while the mass flow rate is 1 kg/s.
- FIG. 11 shows normalized exergy destruction and loss in an absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , as functions of amplitude and mass flow rate wherein the number of waves of the absorber pipe is 10.
- FIG. 12 shows velocity and temperature fields as functions of amplitude wherein the number of waves of the absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , is 10 and the mass flow rate is 10 kg/s.
- FIG. 13 shows fluid temperature distribution at the boundary with an absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A , in which the number of waves is ten and the amplitude is 0.003 m.
- FIG. 14 shows a pressure drop as a function of fluid mass flow rate for different numbers of waves and amplitudes through an absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A .
- FIG. 15 shows performance evaluation criteria as a function of fluid mass flow rate for different numbers of waves and amplitudes through an absorber pipe, such as the absorber pipe shown in FIGS. 1B and 2A .
- Various implementations include a solar thermal energy collection system comprising a solar energy concentrator, a heat transfer fluid, and an absorber pipe through which the heat transfer fluid flows for absorbing concentrated solar energy incident on the absorber pipe from the solar energy concentrator.
- the absorber pipe includes a pipe wall.
- the pipe wall of the absorber pipe extends along and about a central longitudinal axis extending through each end of the absorber pipe.
- the pipe wall has an inner surface and an outer surface.
- the inner surface has a first contour defining alternating peaks and troughs along a length of the absorber pipe.
- the outer surface has a second contour defining alternating peaks and troughs along the length of the absorber pipe.
- the inner surface defines the entire flow path for the heat transfer fluid through the absorber pipe.
- the first contour as viewed through an axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis.
- the sinusoidal waves on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis.
- At least a portion of the absorber pipe is disposed along a focal axis of the solar energy concentrator.
- the absorber pipe enhances thermal performance of the system by accelerating the heat transfer fluid across the trough portions of the first contour of the absorber pipe (or pipe throats), where most solar radiation incident on the absorber pipe is concentrated.
- the absorber pipe may be used for enhancing heat transfer in the heat transfer fluid flowing through the absorber pipe in the solar thermal energy collection system when the absorber pipe receives concentrated heat flux along a focal axis of the solar energy concentrator.
- solar energy concentrators include parabolic trough solar energy absorber systems, solar tower receivers, solar dish receivers, or any system comprising a heat exchanger that is subject to concentrated heat flux along a focal line.
- the solar thermal energy collection system is a parabolic trough solar energy absorber system.
- the solar energy concentrator comprises a parabolic reflector that has a focal axis that is parallel to a central longitudinal axis of the parabolic reflector and is spaced apart a focal length from a vertex of the parabolic reflector.
- the solar thermal energy collection system also includes a tube into which the absorber pipe is disposed. The outer surface of the absorber pipe is spaced apart from an inner surface of the tube, and a wall of the tube allows for the transmission of radiant energy therethrough.
- the tube is coupled to the parabolic reflector such that the central longitudinal axis of the absorber pipe is coincident with the focal axis.
- the absorber pipe described herein provides significant improvement in first law efficiency ( ⁇ c ), especially at mass flow rates ( ⁇ dot over (m) ⁇ f ) of 5 kg/s and greater due to the accelerated flow through the pipe throat portions, which in turn enhances the mixing and heat transfer between the absorber pipe and the heat transfer fluid.
- the second law efficiency ( ⁇ ex ) variations shift from an increasing to a decreasing trend with respect to amplitude (A) as mass flow rate ( ⁇ dot over (m) ⁇ f ) increases and the exergy destroyed due to pressure drop becomes evident.
- the absorber pipe performs better than the conventional straight absorber pipe as its peak performance evaluation criteria (PEC) approaches 1.25 despite its increased pressure drop across the pipe.
- PEC peak performance evaluation criteria
- the performance of the absorber pipe is not an improvement (PEC ⁇ 1) for all number of waves, or peaks (N), and mass flow rate ( ⁇ dot over (m) ⁇ f ) when the amplitude A>0.006 m due to an escalated pressure drop that overturns the heat transfer enhancement.
- PEC ⁇ 1 peak performance evaluation criteria
- N peak performance evaluation criteria
- mass flow rate ⁇ dot over (m) ⁇ f
- FIG. 1A shows a solar thermal energy concentration system 10 that includes an absorber pipe 100 , a tube 200 into which the absorber pipe 100 is disposed, a heat transfer fluid 300 , and a parabolic reflector 400 , according to one implementation.
- FIG. 1B illustrates an axial cross-sectional view of the system 10 as taken through the C-C line shown in FIG. 1A .
- the absorber pipe 100 in FIG. 1B includes a pipe wall 102 .
- the pipe wall 102 of the absorber pipe 100 extends along and about a central longitudinal axis 104 extending through each end 106 a, 106 b of the absorber pipe 100 .
- the pipe wall 102 has an inner surface 108 and an outer surface 110 .
- the inner surface 108 has a first contour defining alternating peaks 108 a and troughs 108 b along a length L of the absorber pipe 100 .
- the outer surface 110 has a second contour defining alternating peaks 110 a and troughs 110 b along the length L of the absorber pipe 100 .
- the inner surface 108 defines the entire flow path for the heat transfer fluid 300 through the absorber pipe 100 .
- the first contour of the inner surface 108 forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis 104 .
- the axial cross section refers to the cross section of the absorber pipe 100 as taken in a plane that includes the central longitudinal axis.
- the sinusoidal waves on each side of the central longitudinal axis 104 are symmetrical with respect to the central longitudinal axis 104 .
- the second contour of the outer surface 110 of the absorber pipe 100 forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis 104 .
- the sinusoidal waves of the second contour of the outer surface 110 on each side of the central longitudinal axis 104 are symmetrical with respect to the central longitudinal axis 104 .
- the pipe wall 102 has a maximum radius R max108 as measured between the inner surface 108 and the central longitudinal axis 104 at each peak 108 a and a minimum radius R min108 as measured between the inner surface 108 and the central longitudinal axis 104 at each trough 108 b.
- the peaks 108 a of the inner surface 108 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of the absorber pipe 100
- the troughs 108 b of the inner surface 108 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of the absorber pipe 100 .
- the outer surface 110 of the pipe wall 102 has a second maximum radius R max110 as measured between the outer surface 110 and the central longitudinal axis 104 at each peak 110 a and a second minimum radius R min110 as measured between the outer surface 110 and the central longitudinal axis 104 at each trough 110 b.
- the peaks 110 a of the outer surface 110 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of the absorber pipe 100
- the troughs 110 b of the outer surface 110 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of the absorber pipe 100 .
- the peaks 108 a and troughs 108 b of the inner surface 108 are also radially aligned with the peaks 110 a and troughs 110 b of the outer surface 110 .
- a thickness between the inner surface 108 and the outer surface 110 of the pipe wall 102 is constant along the length L of the pipe wall 102 .
- the first contour of the inner surface 108 is shaped to induce a plurality of turbulent fluid flow regions along the length L of the central longitudinal axis 104 .
- the troughs 110 b of the second contour of the outer surface 110 define areas of concentrated solar energy.
- the portions of the contour of the outer surface 110 that define the troughs emit a small amount of heat back toward themselves via radiative heat transfer due to the nature of sinusoidal waves.
- A, ⁇ ′′ b , L, L f , ⁇ dot over (m) ⁇ f , and r are amplitude, beam irradiance, absorber pipe length, focal length, fluid mass flow rate, and radius, respectively.
- ⁇ and ⁇ denote periodicity length and rim angle, whereas the subscripts ai, ao, ci, and co refer to the inner surface 108 of the absorber pipe 100 , the outer surface 110 of the absorber pipe 100 , the inner surface 202 of the tube 200 , and an outer surface 206 of the tube 200 , respectively.
- the first and second contours of the inner and outer surfaces 108 , 110 of the absorber pipe 100 is defined by the sinusoidal wavy curve given as
- ⁇ is defined as the ratio of L and the total number of waves N and z is the axial length.
- a ratio of the maximum radius R max108 to the minimum radius R min108 is at least 6:1 in some implementations. And, a ratio of the maximum radius R max110 to the minimum radius R min110 is at least 6:1 in some implementations.
- the outer surface 110 of the absorber pipe 100 is spaced apart from an inner surface 202 of the tube 200 , and a wall 204 of the tube 200 allows for the transmission of radiant energy therethrough.
- the tube 200 may include known tubes into which conventional straight absorber pipes are disposed, for example.
- the absorber pipe 100 is coupled within the tube 200 via one or more brackets that extend between the inner surface 202 of the tube 200 and the outer surface 110 of the absorber pipe 100 .
- the brackets may be known brackets currently used to couple conventional straight absorber pipes within tubes.
- the parabolic reflector 400 has a focal axis 404 that is parallel to (e.g., coincident with) a central longitudinal axis of the parabolic reflector 400 and is spaced apart a focal length L f from a vertex of the parabolic reflector 400 .
- the heat transfer fluid 300 is disposed within the absorber pipe 100 .
- the tube 200 is coupled to the parabolic reflector 400 such that the central longitudinal axis 104 of the absorber pipe 100 is coincident with the focal axis 404 .
- the tube is coupled to the parabolic reflector with brackets that extend between the tube and the reflector. Adjacent brackets are spaced apart from each other along the length of the tube.
- the brackets may be known brackets currently used to couple conventional tubes to reflectors.
- the parabolic reflector may be any known parabolic reflector.
- a vacuum condition is maintained between the outer surface 110 of the absorber pipe 100 and the inner surface 202 of the tube 200 .
- the vacuum condition is a condition in which the pressure is sufficiently low such that any convection through the annulus is negligible at the expected temperatures.
- the system 10 also includes a heat transfer fluid tube 500 coupled to the entry end 106 a and the exit end 106 b of the absorber pipe 100 such that the heat transfer fluid tube 500 and the absorber pipe 100 form a closed loop.
- the system 10 also includes a heat exchanger 600 .
- the heat exchanger is a boiler tank in one implementation. At least a portion of the heat transfer fluid tube 500 is disposed inside the heat exchanger 600 .
- the system 10 includes a fluid pump 700 coupled to the heat transfer fluid tube 500 . The fluid pump 700 causes the heat transfer fluid 300 to flow through the heat transfer fluid tube 500 and the absorber pipe 100 .
- FIG. 2A shows an absorber pipe 100 ′ according to another implementation.
- the absorber pipe 100 ′ has the same periodicity and amplitude as the absorber pipe 100 except that the absorber pipe 100 ′ has six peaks instead of the three peaks and is longer.
- FIG. 2B is a partial axial cross-sectional view of the absorber pipe 100 ′.
- the absorber pipe may have two or more peaks and one or more troughs.
- the heat transfer fluid 300 is a liquid that has high thermal conductivity, low viscosity, and a high boiling point.
- suitable heat transfer fluids include synthetic oils, such as silicone-based oils (e.g., SYLTHERMTM 800), mixtures of diphenyl oxide and biphenyl (e.g., THERMINOL® VP-1), or mixtures of diphenylethane and alkylated aromatics (e.g., DOWTHERMTM Q), and molten salt(s).
- the heat transfer fluid 300 used with this system may include any other known heat transfer fluid typically used with conventional straight absorber pipes.
- the absorber pipes 100 , 100 ′ are formed from non-corrosive materials having a high thermal conductivity such as that of nickel-based alloys or stainless steel, e.g., 20 W/m ⁇ K.
- copper is used as a foil or thin sublayer in between stainless steel and/or nickel-based alloys to prevent corrosion.
- suitable materials may include Ni—NiOx/Cu—Ni, CrN—CrxOy, TiNxOy, Ni/NiOx, Ag/Alumina, SS/Alumina, Cu/Alumina, and Ti/Silica.
- a foil is deposited on the outer surface of the absorber pipes 100 , 100 ′ for increasing solar energy absorption.
- the foil may include a suitable metal for increasing absorption by the pipe, such as Al, Cu, or SS.
- deposition methods for the foil may include sputtering and/or evaporation.
- the chart below provides examples of combinations of materials used for the absorber pipe and the foil along with deposition methods for the foil and the percentage of solar radiation absorbed by the pipe for the particular example combination.
- Absorber Deposition Materials Foil Method Absorption Ni—NiOx/Cu—Ni Al Sputtering 95% CrN—CrxOy Cu Sputtering 93% TiNxOy Cu Evaporation 95% Cu—Cr—C Cu Sputtering 93% Ni/NiOx Al Sputtering 95% Ag/Alumina Al Sputtering/ 95% Evaporation SS/Alumina SS Sputtering 94% Cu/Alumina Cu Sputtering 93% Ti/silica Al Sputtering 95%
- Models of the absorber pipes 100 , 100 ′ shown in FIGS. 1A and 2B and other absorber pipe models comprising varying numbers of peaks and amplitudes were analyzed using an experimentally validated 3D finite element method (FEM) coupled to a semi-analytical optical model.
- FEM 3D finite element method
- the absorber pipe models were assessed based on the first and second law efficiency as well as the performance evaluation criteria (PEC), which accounted for both heat transfer enhancement and pressure drop relative to a conventional straight absorber pipe design, as the wave amplitude, periodicity length, and fluid mass flow rate were varied.
- FEM 3D finite element method
- the computational domain was symmetric about the vertical centerline in the x ⁇ y (or r ⁇ ) plane, hence only half of the absorber system was modeled and simulated.
- the vacuum condition is a condition in which the pressure is sufficiently low such that any convection through the annulus is negligible at the expected temperatures.
- the absorber pipe boundary was shifted vertically by r ao and r ao (see FIG. 2B ) to ensure the mean solar concentration is the same as that on a conventional straight absorber pipe.
- the optical analysis aims to maximize the solar concentration ratio and minimize the errors associated with the tracking, geometric configurations, and operation.
- the optical component of the PTC was modeled by adapting an optical model based on the semi-finite formulation. Consequently, the concentrated heat flux distribution around the receiver was obtained quickly given the PTC geometry, optical properties, and incidence angle. Furthermore, the optical model has served as a reliable reference to verify several Monte Carlo Ray Tracing models proposed afterwards. The optical model is summarized herein for congruity and transparency of this study.
- G . s ′′ ⁇ - ⁇ m ⁇ m ⁇ ( 2 ⁇ G . b ′′ ⁇ ⁇ sin ⁇ s ) ⁇ u y ⁇ z ⁇ cos ⁇ ( ⁇ / 2 - ⁇ y ⁇ z ) ⁇ cos ⁇ ⁇ 0 ⁇ yz ⁇ cos ⁇ ( ⁇ / 2 ) ⁇ cos ⁇ ⁇ 0 ⁇ [ 1 - sin 2 ⁇ ⁇ 0 sin 2 ⁇ s ] 1 / 2 ⁇ cos ⁇ ⁇ 0 ⁇ d ⁇ ⁇ ( 2 )
- ⁇ ′′ b and s are direct normal irradiance and angular radius of the sun, respectively;
- u yz is the projected path length of the reflected central ray;
- ⁇ is the azimuth;
- ⁇ yz and ⁇ 0 are the projected skewness angle and its minimum value, respectively;
- ⁇ 0 tan ⁇ 1 (tan ⁇ i /cos ⁇ yz ) and ⁇ i is the incidence angle; and ⁇ yz is the reflected ray length.
- ⁇ ′′ s 0 for ⁇ 0 >s.
- the absorber pipe partially shades the reflector, thus the incoming sun ray is blocked by the receiver when
- the semi-analytical model described above was modified to account for the curvature of the symmetric wavy absorber pipe, such as the pipes shown in FIGS. 1B and 2A , in the axial direction. Subsequently, Eq. (2) was solved as a function of z and the resulting local concentration ratio (LCR) around the absorber pipe is displayed in FIG. 3B , along with the LCR for a conventional straight absorber pipe in FIG. 3A as a comparison.
- the LCR is the ratio between concentrated solar energy around an absorber and total direct beam solar radiation.
- FIG. 3A and 3B due to the symmetry about the vertical centerline.
- the LCR around the tube 200 in FIGS. 3A and 3B remains fixed as long as its geometry does not change.
- FIG. 4 also depicts the LCR profile around the outer surface of the absorber pipe surface of FIG. 1B in which the number of waves, or peaks, is three, as an example.
- FIG. 4 shows a continuous transition in the LCR between two peaks
- FIG. 3B depicts the distribution at three representative cross sections.
- FIG. 4 is a 3D distribution whereas FIG. 3B is a 2D distribution.
- u, p, ⁇ , ⁇ , c p , T, and Pr are velocity, pressure, dynamic viscosity, density, specific heat at constant pressure, temperature, and Prandtl number, respectively.
- S T The source term in Eq. (7), S T , is zero in the solid domain while
- ⁇ in Eq. (8) is the viscous stress tensor.
- v 2 ⁇ f model is employed to solve the turbulent flow from which ⁇ T results.
- v 2 ⁇ f model is similar to the standard k ⁇ model, but describes the anisotropy of the turbulence intensity in the turbulent boundary layer and non-local pressure-strain effects using two new equations in addition to the two equations for turbulence kinetic energy (k) and dissipation rate ( ⁇ ). The model has been demonstrated to outperform other turbulence models especially since it does not use wall functions.
- the boundary condition imposed at the outer surface of the tube is the sum of radiative and convective heat flux given by
- the convective heat transfer coefficient is computed as
- h _ 0 k c 2 ⁇ r c ⁇ o ⁇ ⁇ 0.3 + 0.62 ⁇ Re D 1 / 2 ⁇ P ⁇ r 1 / 3 ( 1 + ( 0.4 P ⁇ r ) 2 / 3 ) 1 / 4 ⁇ [ 1 + ( Re D 2 ⁇ 8 ⁇ 2 ⁇ 1 ⁇ 0 3 ) 5 / 8 ] 4 / 5 ⁇ , ( 10 )
- ⁇ c is the tube reflectance and ⁇ ′′ ci is the mutual irradiance between a point along the inner tube surface P ci and that on the outer surface of the absorber pipe P ao , and it can be expressed as
- ⁇ c is the tube absorptance
- n c and n a are normal vectors from P ci and P ao
- r is the distance between the two points
- J ao is the radiosity of the outer surface of the absorber pipe S ao .
- the integral (view factor) is evaluated by the hemicube method.
- Eq. (11) is imposed as the outer boundary condition of the absorber pipe by interchanging the subscripts ci and ao, and ⁇ ′′ ao is then the sum of mutual irradiance and the local concentrated irradiance ⁇ ′′ s,ao at the outer surface of absorber pipe.
- the local concentrated irradiance is obtained from the LCR profiles for the tube and absorber pipe given by the optical model described earlier.
- the profiles are projected onto the corresponding surfaces and interpolated. ⁇ ′′ s,co and ⁇ ′′ s,ao are therefore obtained as
- optical efficiency ⁇ opt ⁇ a ⁇ cl ⁇ c
- ⁇ represents the intercept factor accounting for absorber pipe shadowing, twisting and tracking error, geometry error, dirt on mirrors, dirt on HCE, and miscellaneous optical efficiency
- ⁇ cl and ⁇ c denote clean mirror reflectance and tube transmittance, respectively.
- the first law efficiency of a PTC, ⁇ c is defined as the ratio of useful heat ⁇ dot over (Q) ⁇ u and incoming solar energy,
- ⁇ dot over (X) ⁇ s is the solar exergy (exergy flow on the incoming undiluted solar irradiance) determined according to the Petela model as
- ⁇ dot over (X) ⁇ l ⁇ dot over (X) ⁇ l,opt + ⁇ dot over (X) ⁇ l, ⁇ dot over (Q) ⁇ , (20)
- ⁇ dot over (Q) ⁇ ′′ loss is the radiosity of the outer surface of the absorber pipe, i.e., heat leaving the absorber pipe via radiation as convection in the annulus is not considered.
- ⁇ dot over (X) ⁇ d ⁇ dot over (X) ⁇ d, ⁇ P + ⁇ dot over (X) ⁇ d, ⁇ dot over (Q) ⁇ a + ⁇ dot over (X) ⁇ d, ⁇ dot over (Q) ⁇ f , (23)
- the LCR profiles were obtained first by solving the optical model as a function of the absorber pipe geometry.
- the hybrid Streamline/Upwind Petrov-Galerkin finite element method was then employed to discretize and solve Eqs. (5)-(7) along with the boundary conditions.
- Adaptive damped Newton's method was used to update the variables and stepsize at each iteration, and the resulting systems of linear equations were solved using the Parallel Direct Sparse Solver developed to solve large sparse symmetric and asymmetric linear systems of equations on shared and distributed-memory architectures.
- the model was implemented in a commercial finite element package and MATLAB.
- the converged solution was obtained when the relative error in the norm of individual solution vectors (i.e., velocity, pressure, turbulence variables, elliptic function, and temperature) reached 10 ⁇ 3 .
- the relative residual was less than 10 ⁇ 14 in all cases.
- the model validation included two steps. First, only the symmetric wavy absorber pipe, such as shown in FIGS. 1B and 2A-2B , was simulated at different Reynolds numbers (Re), and corresponding average Nusselt numbers were compared against the results of a known wavy absorber pipe with a different contour than the contour described and shown herein, as shown in FIG. 6A .
- the complete PTC model with a conventional absorber pipe was then validated against experimental data with a conventional absorber pipe, and the first law efficiency was compared with the change in temperature in FIG. 6B .
- the average Nusselt number Nu from FIG. 6A is given by
- T b (z) is the bulk fluid temperature in the axial direction z.
- the problem included evaluating ⁇ c and ⁇ ext as functions of geometric parameters, namely N and A (thus ⁇ ), by means of exhaustive search under fixed length L of the tube and absorber pipe.
- This approach facilitated the quantification of the collector sensitivity to the considered parameter ranges as well as the determination of their physical implications.
- the considered ranges were 0.003 m ⁇ A ⁇ 0.01 m and 5 ⁇ N ⁇ 10.
- the mass flow rate was varied from 1 kg/s ⁇ dot over (m) ⁇ f ⁇ 10 kg/s while the inlet fluid temperature was fixed at 400K.
- the conventional straight absorber pipe was elongated to achieve the same absorber pipe volume as the modeled symmetric wavy absorber pipe and thus the material cost remained fixed.
- the absorber system design and operational parameters used in the simulations are listed in Table 1.
- Table 2 lists the grid convergence study results in which Nu and the maximum fluid temperature in the absorber system (T max ) are compared. T max depends on one nodal temperature and converges quicker as the mesh size increases, whereas Nu is computed based on multiple nodal values and exhibits a slower convergence. Mesh refinements continued until the relative error of both values became less than 1%.
- the maximum relative elongation of the straight absorber pipe was 3% of the length L defined for the symmetric wavy pipe, such as the pipes shown in FIGS. 1B and 2A .
- ⁇ c increases almost linearly along with A for all N.
- ⁇ c of the wavy absorber pipe is lower than that of the conventional pipe until A>0.006 m.
- mass flow rates e.g. 1 kg/s
- the heat transfer depends heavily on the presence of flow recirculation and deceleration in convex regions.
- ⁇ c of the conventional straight absorber pipe design decreases drastically as ⁇ dot over (m) ⁇ f increases according to FIGS. 8A-9B .
- the slope of ⁇ c becomes steeper as ⁇ dot over (m) ⁇ f increases for all N and the interval between the least and peak ⁇ c stretches.
- ⁇ ex of the conventional absorber pipe barely decrease with respect to increasing ⁇ dot over (m) ⁇ f
- the modeled symmetric wavy absorber pipe design yields ⁇ ex whose slope shifts gradually from positive to negative as ⁇ dot over (m) ⁇ f increases.
- Such a shift in ⁇ ex trend is investigated further based on the results in FIG. 11 .
- ⁇ tilde over (X) ⁇ d, ⁇ P is identified as the major factor causing the slope of ⁇ ex to shift from positive to negative with respect to increasing ⁇ dot over (m) ⁇ f as its contribution becomes evident only at higher ⁇ dot over (m) ⁇ f according to FIG. 11 .
- the exergy destroyed due to heat transfer from the sun to the absorber pipe surface, ⁇ tilde over (X) ⁇ d, ⁇ dot over (Q) ⁇ a remains nearly constant in all cases as the absorber pipe boundary is shifted vertically by r ai and r ao to ensure constant mean solar concentration.
- ⁇ tilde over (X) ⁇ d, ⁇ dot over (Q) ⁇ f is notable only at a low ⁇ dot over (m) ⁇ f and A.
- the pressure drop and heat transfer across the last segment of the modeled symmetric wavy absorber pipe are analyzed further based on the streamlines and temperature fields shown in FIG. 12 .
- FIG. 12 the convex area covered by the separated flow increases along with A and the effective pipe diameter dwindles thereby, causing a larger pressure drop.
- FIG. 13 depicts the fluid temperature distribution at the boundary with the absorber pipe whose hot spots are observed at the throats according to its local concentration profile.
- ⁇ c increases along with N and A and remains relatively high even at high mass flow rates unlike the conventional straight absorber pipe.
- the exergy destruction due to pressure drop ( ⁇ dot over (X) ⁇ d, ⁇ P ) is therefore pronounced at a high ⁇ dot over (m) ⁇ f (refer to FIG. 11 ), N, and A, causing ⁇ ex to decrease rapidly with respect to increasing N and A as shown in FIG. 9B .
- the physical implications of the qualitative trends in FIG. 15 can be discerned based on the results presented above, namely ⁇ c in FIGS. 7A-9B and ⁇ P in FIG. 14 .
- a symmetric wavy absorber pipe such as the pipes shown in FIGS. 1B and 2A , with a smaller N and A enhances heat transfer while diminishing the effect of pressure drop as compared to the conventional straight absorber pipe.
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Thermal Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Optical Elements Other Than Lenses (AREA)
Abstract
Description
- This application claims priority to U.S. Provisional Patent Application No. 63/002,834, filed Mar. 31, 2020, and entitled “Solar Energy Collection System with Symmetric Wavy Absorber Pipe,” the contents of which are herein incorporated by reference in their entirety.
- The promising future of solar energy as a prominent alternative to fossil fuels has resulted in substantial research and development of high-efficiency parabolic trough solar collectors (PTCs) with applications in solar thermal power plants and industrial processes such as drying, sterilization, and desalination. Various approaches have been proposed to achieve desired thermal and optical characteristics of PTCs encompassing optimization and enhancement of collector geometry, optical and material properties, heat transfer fluids, and collector operation. To date, PTCs have used straight absorber pipes to collect solar energy. However, there is a need to improve the thermal performance of PTCs.
- Various implementations include a solar thermal energy collection system that comprises a solar energy concentrator, a heat transfer fluid, and an absorber pipe through which the heat transfer fluid flows for absorbing concentrated solar energy incident on the absorber pipe from the solar energy concentrator. The absorber pipe comprises a pipe wall that extends about and around a central longitudinal axis that extends between each end of the absorber pipe. The pipe wall comprises an inner surface and an outer surface. The inner surface defines a first contour that defines alternating peaks and troughs along a length of the absorber pipe, and the outer surface defines a second contour that defines alternating peaks and troughs along the length of the absorber pipe. The inner surface defines the entire flow path for the heat transfer fluid through the absorber pipe. The first contour, as viewed through an axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis. The sinusoidal waves on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis, and at least a portion of the absorber pipe is disposed along a focal line of the solar energy concentrator.
- In some implementations, the second contour, as viewed through the axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis. The sinusoidal waves of the second contour on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis.
- In some implementations, the pipe wall has a maximum radius as measured between the inner surface and the central longitudinal axis at each peak and a minimum radius as measured between the inner surface and the central longitudinal axis at each trough.
- In some implementations, the outer surface of the pipe wall has a second maximum radius as measured between the outer surface and the central longitudinal axis at each peak and a second minimum radius as measured between the outer surface and the central longitudinal axis at each trough.
- In some implementations, the peaks and troughs of the inner surface are radially aligned with the peaks and troughs of the outer surface.
- In some implementations, a thickness between the inner surface and the outer surface of the wall is constant along the length of the wall.
- In some implementations, the inner surface is shaped to induce a plurality of turbulent fluid flow regions along the length of the central longitudinal axis.
- In some implementations, troughs of the second contour define areas of concentrated solar energy.
- In some implementations, the solar thermal energy collection system is a parabolic trough solar energy absorber system. The solar energy concentrator comprises a parabolic reflector having a focal axis that is parallel to a central longitudinal axis of the parabolic reflector and is spaced apart a focal length from a vertex of the parabolic reflector. The parabolic trough solar energy absorber system further comprises a tube into which the absorber pipe is disposed. The outer surface of the absorber pipe is spaced apart from an inner surface of the tube, and a wall of the tube allows for the transmission of radiant energy therethrough. The tube is coupled to the parabolic reflector such that the central longitudinal axis of the absorber pipe is coincident with the focal axis.
- In some implementations, a vacuum condition is maintained between the outer surface of the absorber pipe and the inner surface of the tube.
- In some implementations, the solar thermal energy collection system further comprises a heat transfer fluid tube coupled to an entry end and an exit end of the absorber pipe. The heat transfer fluid tube and the absorber pipe form a closed loop through which the heat transfer fluid flows.
- In some implementations, the system further comprises a heat exchanger, wherein at least a portion of the heat transfer fluid tube is disposed inside the heat exchanger.
- In some implementations, the heat exchanger is a boiler tank.
- In some implementations, the system further comprises a fluid pump coupled to the heat transfer fluid tube, and the fluid pump causes the heat transfer fluid to flow through the heat transfer fluid tube and the absorber pipe.
- Example features and implementations are disclosed in the accompanying drawings. However, the present disclosure is not limited to the precise arrangements and instrumentalities shown. Similar elements in different implementations are designated using the same reference numerals.
-
FIG. 1A shows an end view of a schematic diagram of a parabolic trough solar energy collector system comprising an absorber pipe according to one implementation. -
FIG. 1B shows an axial cross-sectional view along the C-C line of the parabolic trough solar energy collector system comprising the absorber pipe shown inFIG. 1A . -
FIG. 1C shows a schematic diagram of the system shown inFIG. 1A . -
FIG. 2A shows a perspective view of an absorber pipe according to another implementation. The absorber pipe shown inFIG. 2A is similar to the absorber pipe shown inFIG. 1A but has six peaks. -
FIG. 2B shows a partial axial cross-sectional view of the absorber pipe shown inFIG. 2A . -
FIG. 3A shows a local concentration ratio around a quarter of a conventional straight absorber pipe. -
FIG. 3B shows a local concentration ratio around a quarter of the absorber pipe, such as the absorber pipes shown inFIG. 1B or 2A . -
FIG. 4 shows a local concentration ratio profile around the absorber pipe inFIG. 1B . -
FIG. 5 is a convergence plot for an implementation of the absorber pipe, such as the absorber pipes shown inFIGS. 1B and 2A , that has seven waves and an amplitude of 0.006 m. -
FIG. 6A is a validation plot for the absorber pipe, such as the absorber pipes shown inFIGS. 1B and 2A . -
FIG. 6B is a validation plot for a complete PTC model with a conventional absorber pipe. -
FIG. 7A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipes shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 1 kg/s. -
FIG. 7B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 1 kg/s. -
FIG. 8A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 5 kg/s. -
FIG. 8B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 5 kg/s. -
FIG. 9A shows variations of first law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 10 kg/s. -
FIG. 9B shows variations of second law efficiency as functions of the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , and amplitude for a mass flow rate of 10 kg/s. -
FIG. 10 shows streamlines at the last pipe segment cross section of an absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , for different numbers of waves and different amplitudes while the mass flow rate is 1 kg/s. -
FIG. 11 shows normalized exergy destruction and loss in an absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , as functions of amplitude and mass flow rate wherein the number of waves of the absorber pipe is 10. -
FIG. 12 shows velocity and temperature fields as functions of amplitude wherein the number of waves of the absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , is 10 and the mass flow rate is 10 kg/s. -
FIG. 13 shows fluid temperature distribution at the boundary with an absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A , in which the number of waves is ten and the amplitude is 0.003 m. -
FIG. 14 shows a pressure drop as a function of fluid mass flow rate for different numbers of waves and amplitudes through an absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A . -
FIG. 15 shows performance evaluation criteria as a function of fluid mass flow rate for different numbers of waves and amplitudes through an absorber pipe, such as the absorber pipe shown inFIGS. 1B and 2A . - Various implementations include a solar thermal energy collection system comprising a solar energy concentrator, a heat transfer fluid, and an absorber pipe through which the heat transfer fluid flows for absorbing concentrated solar energy incident on the absorber pipe from the solar energy concentrator. The absorber pipe includes a pipe wall. The pipe wall of the absorber pipe extends along and about a central longitudinal axis extending through each end of the absorber pipe. The pipe wall has an inner surface and an outer surface. The inner surface has a first contour defining alternating peaks and troughs along a length of the absorber pipe. The outer surface has a second contour defining alternating peaks and troughs along the length of the absorber pipe. The inner surface defines the entire flow path for the heat transfer fluid through the absorber pipe. The first contour, as viewed through an axial cross section of the absorber pipe, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis. The sinusoidal waves on each side of the central longitudinal axis are symmetrical with respect to the central longitudinal axis. At least a portion of the absorber pipe is disposed along a focal axis of the solar energy concentrator. In some implementations, the absorber pipe enhances thermal performance of the system by accelerating the heat transfer fluid across the trough portions of the first contour of the absorber pipe (or pipe throats), where most solar radiation incident on the absorber pipe is concentrated.
- In various implementations, the absorber pipe may be used for enhancing heat transfer in the heat transfer fluid flowing through the absorber pipe in the solar thermal energy collection system when the absorber pipe receives concentrated heat flux along a focal axis of the solar energy concentrator. Examples of solar energy concentrators include parabolic trough solar energy absorber systems, solar tower receivers, solar dish receivers, or any system comprising a heat exchanger that is subject to concentrated heat flux along a focal line.
- For example, in certain implementations, the solar thermal energy collection system is a parabolic trough solar energy absorber system. The solar energy concentrator comprises a parabolic reflector that has a focal axis that is parallel to a central longitudinal axis of the parabolic reflector and is spaced apart a focal length from a vertex of the parabolic reflector. The solar thermal energy collection system also includes a tube into which the absorber pipe is disposed. The outer surface of the absorber pipe is spaced apart from an inner surface of the tube, and a wall of the tube allows for the transmission of radiant energy therethrough. The tube is coupled to the parabolic reflector such that the central longitudinal axis of the absorber pipe is coincident with the focal axis.
- In some implementations, the absorber pipe described herein provides significant improvement in first law efficiency (ηc), especially at mass flow rates ({dot over (m)}f) of 5 kg/s and greater due to the accelerated flow through the pipe throat portions, which in turn enhances the mixing and heat transfer between the absorber pipe and the heat transfer fluid. In addition, in some implementations, the second law efficiency (ηex) variations shift from an increasing to a decreasing trend with respect to amplitude (A) as mass flow rate ({dot over (m)}f) increases and the exergy destroyed due to pressure drop becomes evident. Furthermore, in some implementations, the absorber pipe performs better than the conventional straight absorber pipe as its peak performance evaluation criteria (PEC) approaches 1.25 despite its increased pressure drop across the pipe. However, the performance of the absorber pipe is not an improvement (PEC≤1) for all number of waves, or peaks (N), and mass flow rate ({dot over (m)}f) when the amplitude A>0.006 m due to an escalated pressure drop that overturns the heat transfer enhancement. Thus, as the number of peaks and amplitude decrease, the heat transfer performance increases while diminishing the effect of pressure drop as compared to the conventional straight absorber pipe.
- For example,
FIG. 1A shows a solar thermalenergy concentration system 10 that includes anabsorber pipe 100, atube 200 into which theabsorber pipe 100 is disposed, aheat transfer fluid 300, and aparabolic reflector 400, according to one implementation.FIG. 1B illustrates an axial cross-sectional view of thesystem 10 as taken through the C-C line shown inFIG. 1A . - The
absorber pipe 100 inFIG. 1B includes apipe wall 102. Thepipe wall 102 of theabsorber pipe 100 extends along and about a central longitudinal axis 104 extending through eachend 106 a, 106 b of theabsorber pipe 100. Thepipe wall 102 has an inner surface 108 and anouter surface 110. The inner surface 108 has a first contour defining alternating peaks 108 a andtroughs 108 b along a length L of theabsorber pipe 100. Theouter surface 110 has a second contour defining alternatingpeaks 110 a andtroughs 110 b along the length L of theabsorber pipe 100. The inner surface 108 defines the entire flow path for theheat transfer fluid 300 through theabsorber pipe 100. The first contour of the inner surface 108, as viewed through the axial cross section of theabsorber pipe 100, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis 104. The axial cross section refers to the cross section of theabsorber pipe 100 as taken in a plane that includes the central longitudinal axis. The sinusoidal waves on each side of the central longitudinal axis 104 are symmetrical with respect to the central longitudinal axis 104. In addition, the second contour of theouter surface 110 of theabsorber pipe 100, as viewed through the axial cross section of theabsorber pipe 100, forms sinusoidal waves on each side of and spaced apart from the central longitudinal axis 104. The sinusoidal waves of the second contour of theouter surface 110 on each side of the central longitudinal axis 104 are symmetrical with respect to the central longitudinal axis 104. - The
pipe wall 102 has a maximum radius Rmax108 as measured between the inner surface 108 and the central longitudinal axis 104 at each peak 108 a and a minimum radius Rmin108 as measured between the inner surface 108 and the central longitudinal axis 104 at eachtrough 108 b. The peaks 108 a of the inner surface 108 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of theabsorber pipe 100, and thetroughs 108 b of the inner surface 108 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of theabsorber pipe 100. - The
outer surface 110 of thepipe wall 102 has a second maximum radius Rmax110 as measured between theouter surface 110 and the central longitudinal axis 104 at each peak 110 a and a second minimum radius Rmin110 as measured between theouter surface 110 and the central longitudinal axis 104 at eachtrough 110 b. Thepeaks 110 a of theouter surface 110 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of theabsorber pipe 100, and thetroughs 110 b of theouter surface 110 are radially aligned and opposite each other with respect to the central longitudinal axis 104 of theabsorber pipe 100. - The peaks 108 a and
troughs 108 b of the inner surface 108 are also radially aligned with thepeaks 110 a andtroughs 110 b of theouter surface 110. In addition, a thickness between the inner surface 108 and theouter surface 110 of thepipe wall 102 is constant along the length L of thepipe wall 102. The first contour of the inner surface 108 is shaped to induce a plurality of turbulent fluid flow regions along the length L of the central longitudinal axis 104. And, thetroughs 110 b of the second contour of theouter surface 110 define areas of concentrated solar energy. In particular, the portions of the contour of theouter surface 110 that define the troughs emit a small amount of heat back toward themselves via radiative heat transfer due to the nature of sinusoidal waves. - In
FIGS. 1A and 1B , A, Ġ″b, L, Lf, {dot over (m)}f, and r are amplitude, beam irradiance, absorber pipe length, focal length, fluid mass flow rate, and radius, respectively. λ and ψ denote periodicity length and rim angle, whereas the subscripts ai, ao, ci, and co refer to the inner surface 108 of theabsorber pipe 100, theouter surface 110 of theabsorber pipe 100, theinner surface 202 of thetube 200, and anouter surface 206 of thetube 200, respectively. The first and second contours of the inner andouter surfaces 108, 110 of theabsorber pipe 100 is defined by the sinusoidal wavy curve given as -
- where λ is defined as the ratio of L and the total number of waves N and z is the axial length.
- In some implementations, a ratio of the maximum radius Rmax108 to the minimum radius Rmin108 is at least 6:1 in some implementations. And, a ratio of the maximum radius Rmax110 to the minimum radius Rmin110 is at least 6:1 in some implementations.
- The
outer surface 110 of theabsorber pipe 100 is spaced apart from aninner surface 202 of thetube 200, and awall 204 of thetube 200 allows for the transmission of radiant energy therethrough. Thetube 200 may include known tubes into which conventional straight absorber pipes are disposed, for example. Theabsorber pipe 100 is coupled within thetube 200 via one or more brackets that extend between theinner surface 202 of thetube 200 and theouter surface 110 of theabsorber pipe 100. The brackets may be known brackets currently used to couple conventional straight absorber pipes within tubes. Theparabolic reflector 400 has a focal axis 404 that is parallel to (e.g., coincident with) a central longitudinal axis of theparabolic reflector 400 and is spaced apart a focal length Lf from a vertex of theparabolic reflector 400. Theheat transfer fluid 300 is disposed within theabsorber pipe 100. Thetube 200 is coupled to theparabolic reflector 400 such that the central longitudinal axis 104 of theabsorber pipe 100 is coincident with the focal axis 404. For example, in some implementations, the tube is coupled to the parabolic reflector with brackets that extend between the tube and the reflector. Adjacent brackets are spaced apart from each other along the length of the tube. The brackets may be known brackets currently used to couple conventional tubes to reflectors. And, the parabolic reflector may be any known parabolic reflector. - In addition, a vacuum condition is maintained between the
outer surface 110 of theabsorber pipe 100 and theinner surface 202 of thetube 200. The vacuum condition is a condition in which the pressure is sufficiently low such that any convection through the annulus is negligible at the expected temperatures. - As shown in
FIG. 1C , thesystem 10 also includes a heattransfer fluid tube 500 coupled to the entry end 106 a and theexit end 106 b of theabsorber pipe 100 such that the heattransfer fluid tube 500 and theabsorber pipe 100 form a closed loop. Thesystem 10 also includes aheat exchanger 600. For example, the heat exchanger is a boiler tank in one implementation. At least a portion of the heattransfer fluid tube 500 is disposed inside theheat exchanger 600. In addition, thesystem 10 includes afluid pump 700 coupled to the heattransfer fluid tube 500. Thefluid pump 700 causes theheat transfer fluid 300 to flow through the heattransfer fluid tube 500 and theabsorber pipe 100. -
FIG. 2A shows anabsorber pipe 100′ according to another implementation. Theabsorber pipe 100′ has the same periodicity and amplitude as theabsorber pipe 100 except that theabsorber pipe 100′ has six peaks instead of the three peaks and is longer.FIG. 2B is a partial axial cross-sectional view of theabsorber pipe 100′. In other implementations, the absorber pipe may have two or more peaks and one or more troughs. - In some implementations, the
heat transfer fluid 300 is a liquid that has high thermal conductivity, low viscosity, and a high boiling point. For example, suitable heat transfer fluids include synthetic oils, such as silicone-based oils (e.g., SYLTHERM™ 800), mixtures of diphenyl oxide and biphenyl (e.g., THERMINOL® VP-1), or mixtures of diphenylethane and alkylated aromatics (e.g., DOWTHERM™ Q), and molten salt(s). Theheat transfer fluid 300 used with this system may include any other known heat transfer fluid typically used with conventional straight absorber pipes. - In some implementations, the
absorber pipes absorber pipes -
Absorber Deposition Materials Foil Method Absorption Ni—NiOx/Cu—Ni Al Sputtering 95% CrN—CrxOy Cu Sputtering 93% TiNxOy Cu Evaporation 95% Cu—Cr—C Cu Sputtering 93% Ni/NiOx Al Sputtering 95% Ag/Alumina Al Sputtering/ 95% Evaporation SS/Alumina SS Sputtering 94% Cu/Alumina Cu Sputtering 93% Ti/silica Al Sputtering 95% - Models of the
absorber pipes FIGS. 1A and 2B and other absorber pipe models comprising varying numbers of peaks and amplitudes were analyzed using an experimentally validated 3D finite element method (FEM) coupled to a semi-analytical optical model. The absorber pipe models were assessed based on the first and second law efficiency as well as the performance evaluation criteria (PEC), which accounted for both heat transfer enhancement and pressure drop relative to a conventional straight absorber pipe design, as the wave amplitude, periodicity length, and fluid mass flow rate were varied. - The following assumptions were made to simplify the modeling without losing significant accuracy:
- 1. The computational domain was symmetric about the vertical centerline in the x−y (or r−θ) plane, hence only half of the absorber system was modeled and simulated.
- 2. The annulus between the
absorber pipe tube 200 was in perfect vacuum—the convection was therefore neglected and the annulus was not included in the computational domain. In an actual system, the vacuum condition is a condition in which the pressure is sufficiently low such that any convection through the annulus is negligible at the expected temperatures. - 3. Radiation in participating media was neglected; instead, constant optical properties and tracking errors were assumed to compute the concentrated solar radiation reaching the absorber pipe and tube surfaces.
- 4. Adiabatic absorber pipe ends and negligible conduction through the support bracket.
- In addition, for modeling purposes, the absorber pipe boundary was shifted vertically by rao and rao (see
FIG. 2B ) to ensure the mean solar concentration is the same as that on a conventional straight absorber pipe. - The optical analysis aims to maximize the solar concentration ratio and minimize the errors associated with the tracking, geometric configurations, and operation. In this work, the optical component of the PTC was modeled by adapting an optical model based on the semi-finite formulation. Consequently, the concentrated heat flux distribution around the receiver was obtained quickly given the PTC geometry, optical properties, and incidence angle. Furthermore, the optical model has served as a reliable reference to verify several Monte Carlo Ray Tracing models proposed afterwards. The optical model is summarized herein for congruity and transparency of this study.
- The distribution of concentrated heat flux ({dot over (Q)}r″) around the receiver is given by
-
- for δ0≤s, where Ġ″b and s are direct normal irradiance and angular radius of the sun, respectively; uyz is the projected path length of the reflected central ray; ψ is the azimuth; δyz and δ0 are the projected skewness angle and its minimum value, respectively; ϕ0=tan−1(tan θi/cos δyz) and θi is the incidence angle; and ρyz is the reflected ray length. Note that Ġ″s=0 for δ0>s.
- The absorber pipe partially shades the reflector, thus the incoming sun ray is blocked by the receiver when
-
d=u yz sin(ψ−δyz)≤r ao, (3) - In such cases, the incoming rays strike the tube and absorber tube top surfaces and the irradiance is computed as
-
Ġ″ s =Ġ″ b cos θi cos(π−β). (4) - where β is the receiver azimuthal angle. Eq. (2) was solved numerically along with the shadowing condition in Eq. (3) using QUADPACK. The limits of integration and integrand were modified to account for the shadowing effect.
- The semi-analytical model described above was modified to account for the curvature of the symmetric wavy absorber pipe, such as the pipes shown in
FIGS. 1B and 2A , in the axial direction. Subsequently, Eq. (2) was solved as a function of z and the resulting local concentration ratio (LCR) around the absorber pipe is displayed inFIG. 3B , along with the LCR for a conventional straight absorber pipe inFIG. 3A as a comparison. The LCR is the ratio between concentrated solar energy around an absorber and total direct beam solar radiation. The LCR is defined as LCR=Ġ″s/Ġ″b, and only a quarter of the absorber pipe is featured inFIGS. 3A and 3B due to the symmetry about the vertical centerline. The LCR around thetube 200 inFIGS. 3A and 3B remains fixed as long as its geometry does not change.FIG. 3B depicts the LCR at three different axial lengths (z1, z2, and z3) illustrated in the diagram next to the plot. In both cases, the heat flux is lower when β=0 owing to the shade created by the absorber pipe, and the LCR eventually becomes relatively small as β decreases since only the direct beam on the receiver top persists.FIG. 4 also depicts the LCR profile around the outer surface of the absorber pipe surface ofFIG. 1B in which the number of waves, or peaks, is three, as an example.FIG. 4 shows a continuous transition in the LCR between two peaks, andFIG. 3B depicts the distribution at three representative cross sections.FIG. 4 is a 3D distribution whereasFIG. 3B is a 2D distribution. - The model was formulated based on mass, momentum and energy conservation laws given in the Cartesian tensor notation as
-
- where u, p, μ, ρ, cp, T, and Pr are velocity, pressure, dynamic viscosity, density, specific heat at constant pressure, temperature, and Prandtl number, respectively. The source term in Eq. (7), ST, is zero in the solid domain while
-
- in the fluid domain accounting for the work done by pressure changes and viscous dissipation, respectively. τ in Eq. (8) is the viscous stress tensor.
- The term in the parenthesis in Eq. (7) represents the conductive heat flux in Reynolds-averaged Navier-Stokes whose PrT is computed using Kays-Crawford heat transport turbulence model. In this model,
v 2−f model is employed to solve the turbulent flow from which μT results.v 2−f model is similar to the standard k−ε model, but describes the anisotropy of the turbulence intensity in the turbulent boundary layer and non-local pressure-strain effects using two new equations in addition to the two equations for turbulence kinetic energy (k) and dissipation rate (ε). The model has been demonstrated to outperform other turbulence models especially since it does not use wall functions. - The boundary condition imposed at the outer surface of the tube is the sum of radiative and convective heat flux given by
-
- where εc and σ are the tube emittance and Stefan-Boltzmann constant, respectively; Tsky=T0−8 in which T0 is the ambient temperature; and Ġ″s,co is the local concentrated irradiance discussed later. The convective heat transfer coefficient is computed as
-
- which is an averaged heat transfer correlation applied to forced convection on a cylinder in cross flow.
- Only radiation through the annulus (i.e., between absorber pipe and tube) is considered as mentioned earlier. Hence the inner tube boundary is defined as the following:
-
- where radiosity J is given by
-
J ci=ρc Ġ″ ci−εc σT 4. (12) - Here ρc is the tube reflectance and Ġ″ci is the mutual irradiance between a point along the inner tube surface Pci and that on the outer surface of the absorber pipe Pao, and it can be expressed as
-
- in which αc is the tube absorptance; nc and na are normal vectors from Pci and Pao; r is the distance between the two points; and Jao is the radiosity of the outer surface of the absorber pipe Sao. The integral (view factor) is evaluated by the hemicube method. Eq. (11) is imposed as the outer boundary condition of the absorber pipe by interchanging the subscripts ci and ao, and Ġ″ao is then the sum of mutual irradiance and the local concentrated irradiance Ġ″s,ao at the outer surface of absorber pipe.
- The local concentrated irradiance is obtained from the LCR profiles for the tube and absorber pipe given by the optical model described earlier. The profiles are projected onto the corresponding surfaces and interpolated. Ġ″s,co and Ġ″s,ao are therefore obtained as
-
Ġ″ s,co=γρcl Ġ″ bLCRco(x, y, z) (14) -
and -
Ġ″ s,ao=ηopt Ġ″ bLCRao(x, y, z) (15) - Here the optical efficiency ηopt=αaγρclτc where γ represents the intercept factor accounting for absorber pipe shadowing, twisting and tracking error, geometry error, dirt on mirrors, dirt on HCE, and miscellaneous optical efficiency; and ρcl and τc denote clean mirror reflectance and tube transmittance, respectively.
- The boundary conditions at the absorber pipe inlet included fully-developed average flow velocity evaluated at the inlet fluid temperature Tf,i since periodic flow conditions could not be imposed with temperature-dependent fluid properties; no-slip condition at the inner surface of the absorber pipe; symmetry about the vertical centerline in the x−y plane (∂T/∂n=0); and outflow condition at the absorber pipe outlet, i.e., μ(∂u/∂n)=0, ρ=0, and ∂T/∂n=0.
- The first law efficiency of a PTC, ηc, is defined as the ratio of useful heat {dot over (Q)}u and incoming solar energy,
-
- where the subscripts f, i and f, o denote fluid inlet and outlet while the aperture area Aap is determined as
-
- The second law efficiency (exergetic efficiency) is given by
-
- where {dot over (X)}s is the solar exergy (exergy flow on the incoming undiluted solar irradiance) determined according to the Petela model as
-
- in which Ts=5770K. The exergy loss {dot over (X)}l in Eq. (18) is given by
-
{dot over (X)} l ={dot over (X)} l,opt +{dot over (X)} l,{dot over (Q)}, (20) - where {dot over (X)}l,opt and {dot over (X)}l,{dot over (Q)} are energy losses associated with optical and thermal losses obtained as
-
- Here {dot over (Q)}″loss is the radiosity of the outer surface of the absorber pipe, i.e., heat leaving the absorber pipe via radiation as convection in the annulus is not considered.
- The destruction term in Eq. (18) can be expanded to
-
{dot over (X)} d ={dot over (X)} d,ΔP +{dot over (X)} d,{dot over (Q)}a +{dot over (X)} d,{dot over (Q)}f , (23) - in which the terms on the right represent exergy destruction due to pressure drop, heat transfer from the sun to the absorber pipe outer surface, and from the absorber pipe to the heat transfer fluid, respectively. The terms in Eq. (23) are computed as follows:
-
- The LCR profiles were obtained first by solving the optical model as a function of the absorber pipe geometry. The hybrid Streamline/Upwind Petrov-Galerkin finite element method was then employed to discretize and solve Eqs. (5)-(7) along with the boundary conditions. A combination of hexahedral and tetrahedral quadratic Lagrange elements (p=2) was employed to discretize fluid velocity and pressure while temperature was discretized with linear elements (p=1). Adaptive damped Newton's method was used to update the variables and stepsize at each iteration, and the resulting systems of linear equations were solved using the Parallel Direct Sparse Solver developed to solve large sparse symmetric and asymmetric linear systems of equations on shared and distributed-memory architectures. The model was implemented in a commercial finite element package and MATLAB.
- The converged solution was obtained when the relative error in the norm of individual solution vectors (i.e., velocity, pressure, turbulence variables, elliptic function, and temperature) reached 10−3. The relative residual was less than 10−14 in all cases.
FIG. 5 shows the convergence plot for N=7 and A=0.006 m case as an example. - The model validation included two steps. First, only the symmetric wavy absorber pipe, such as shown in
FIGS. 1B and 2A-2B , was simulated at different Reynolds numbers (Re), and corresponding average Nusselt numbers were compared against the results of a known wavy absorber pipe with a different contour than the contour described and shown herein, as shown inFIG. 6A . The complete PTC model with a conventional absorber pipe was then validated against experimental data with a conventional absorber pipe, and the first law efficiency was compared with the change in temperature inFIG. 6B . The average Nusselt numberNu fromFIG. 6A is given by -
- where l is the total absorber pipe boundary length and Tb(z) is the bulk fluid temperature in the axial direction z.
- Simulation and experiment data are in good qualitative and quantitative agreement in both cases according to
FIGS. 6A and 6B . The discrepancy in the observed and predicted absorber system efficiency at ΔT>300K inFIG. 6B is due to the absence of convection in the annulus which becomes predominant as the inlet fluid temperature increases. Nonetheless, all the other predictions are within the error margin and thus validate the model. - The problem included evaluating ηc and ηext as functions of geometric parameters, namely N and A (thus λ), by means of exhaustive search under fixed length L of the tube and absorber pipe. This approach facilitated the quantification of the collector sensitivity to the considered parameter ranges as well as the determination of their physical implications. The considered ranges were 0.003 m≤A≤0.01 m and 5≤N≤10. In addition, the mass flow rate was varied from 1 kg/s≤{dot over (m)}f≤10 kg/s while the inlet fluid temperature was fixed at 400K.
- The modeled symmetric wavy absorber pipe, such as the pipes shown in
FIGS. 1B and 2A , was compared against the conventional straight absorber pipe (i.e., A=0) under a fixed absorber pipe volume for a fair comparison. In other words, the conventional straight absorber pipe was elongated to achieve the same absorber pipe volume as the modeled symmetric wavy absorber pipe and thus the material cost remained fixed. The absorber system design and operational parameters used in the simulations are listed in Table 1. -
TABLE 1 Collector design and operational parameters [42]. Value Description (unit) d Absorber inner diameter 0.066 (m) d Absorber outer diameter 0.070 (m) d Cover inner diameter 0.105 (m) d Cover outer diameter 0.115 (m) Beam radiation 1000 (W/m2) L Collector length 1 (m) L Focal length 1.84 (m) P Annular pressure 0 (Pa) T0 Ambient temperature 293.15 (K) T Inlet fluid temperature 400 (K) U Ambient wind speed 0.2 (m/s) Absorber absorptance 0.95 Cover absorptance 0.023 Intercept factor 0.848 Absorber emittance 0.14 Cover emittance 0.86 Clean mirror reflectance 0.935 Envelope transmittance 0.95 ψ Rim angle 70° Absorber Material 321H SS Envelope material Pyrex glass Heat transfer fluid Syltherm 800 Selective coating Luz cermet indicates data missing or illegible when filed -
TABLE 2 Mesh refinement results. DOF Nu Difference Tmax Difference 91783 4.758 — 430.3 — 135191 4.918 3.2% 434.4 0.9% 501041 5.134 4.2% 437.4 0.7% 1695759 5.159 0.5% 438.9 0.3% - The heat transfer and flow characteristics of each design configuration were evaluated based on temperature fields, streamlines,
Nu , and ΔP. In particular, the absorber system performance was examined based on the performance evaluation criteria (PEC) obtained as follows: -
- where the
subscript 0 denotes the straight absorber pipe design. As its name suggests, PEC is useful in comparing the performance of novel absorber pipe designs against the traditional or previously developed ones while accounting for both pressure drop (pumping power) and thermal performance (heat transfer enhancement). - Table 2 lists the grid convergence study results in which
Nu and the maximum fluid temperature in the absorber system (Tmax) are compared. Tmax depends on one nodal temperature and converges quicker as the mesh size increases, whereasNu is computed based on multiple nodal values and exhibits a slower convergence. Mesh refinements continued until the relative error of both values became less than 1%. -
FIGS. 7A-9B depict the variations of ηc and ηex as functions of N and A for {dot over (m)}f=1, 5, and 10 kg/s. The dotted red horizontal lines indicate the corresponding efficiency obtained with the straight absorber pipe design (A=0), and only one red line is plotted in each graph for comparison purposes as the pipe elongation was minimal. The maximum relative elongation of the straight absorber pipe was 3% of the length L defined for the symmetric wavy pipe, such as the pipes shown inFIGS. 1B and 2A . InFIGS. 7A and 7B , ηc increases almost linearly along with A for all N. Although N=5 barely outperforms the conventional absorber pipe when A=0.01 m, the superior performance of N=7 is apparent when A=0.01 m. On the other hand, higher ηex is always achieved with the wavy absorber pipe for all N and A, and that of N=10 exhibits a nonlinear rise with respect to A. - According to
FIGS. 7A and 7B , ηc of the wavy absorber pipe is lower than that of the conventional pipe until A>0.006 m. At low mass flow rates (e.g., 1 kg/s), the heat transfer depends heavily on the presence of flow recirculation and deceleration in convex regions. When A=0.003 m, for instance, the heat transfer rate (or ηc) decreases as the convexity (or N) increases; hence, ηc of N=10 is lower than that of N=5 inFIG. 7A . However, the trend overturns as ηc of N=10 starts to surpass that of others as A increases.FIG. 10 displays streamlines at the last absorber pipe segment cross section for different N and A while {dot over (m)}f=1 kg/s, and in reference toFIGS. 7A and 7B , N=10 begins to achieve superior performance over the other designs as soon as it forms recirculation regions. The same behavior is observed with N=5 and N=7, both of which exhibit higher ηc than the conventional one once the flow starts to separate. - ηc of the conventional straight absorber pipe design decreases drastically as {dot over (m)}f increases according to
FIGS. 8A-9B . The slope of ηc becomes steeper as {dot over (m)}f increases for all N and the interval between the least and peak ηc stretches. The peak ηc for {dot over (m)}f=1 kg/s inFIG. 7A , for example, is 65% while that for {dot over (m)}f=10 kg/s is 67.14%. Such observations are examined further later in this section. Unlike ηc, ηex of the conventional absorber pipe barely decrease with respect to increasing {dot over (m)}f, whereas the modeled symmetric wavy absorber pipe design yields ηex whose slope shifts gradually from positive to negative as {dot over (m)}f increases. ηex is nearly independent of N and A and close to that of the conventional absorber pipe when {dot over (m)}f=5 kg/s, after which ηex becomes a decreasing function of A. Such a shift in ηex trend is investigated further based on the results inFIG. 11 . -
FIG. 11 shows the contribution of each absorber system component and energy transfer mechanism to the total exergy destruction and loss normalized as {tilde over (X)}={dot over (X)}i/{dot over (X)}s for N=10. The exergy destroyed due to pressure drop, {tilde over (X)}d,ΔP, increases along with A, especially when {dot over (m)}f=10 kg/s as the effective hydraulic diameter dwindles and a larger recirculation region is formed. {tilde over (X)}d,ΔP is identified as the major factor causing the slope of ηex to shift from positive to negative with respect to increasing {dot over (m)}f as its contribution becomes evident only at higher {dot over (m)}f according toFIG. 11 . The exergy destroyed due to heat transfer from the sun to the absorber pipe surface, {tilde over (X)}d,{dot over (Q)}a , remains nearly constant in all cases as the absorber pipe boundary is shifted vertically by rai and rao to ensure constant mean solar concentration. {tilde over (X)}d,{dot over (Q)}f is notable only at a low {dot over (m)}f and A. The pressure drop and heat transfer across the last segment of the modeled symmetric wavy absorber pipe are analyzed further based on the streamlines and temperature fields shown inFIG. 12 . - According to
FIG. 12 , the convex area covered by the separated flow increases along with A and the effective pipe diameter dwindles thereby, causing a larger pressure drop.FIG. 12 also demonstrates how A=0.01 m yields a higher mean fluid temperature and a more gradual temperature gradient near the absorber pipe boundary due to the accelerated flow and its resulting heat transfer enhancement and mixing. Higher convective heat transfer is achieved at narrower pipe throats (inlet and outlet of the convex region) as the fluid accelerates through them which are also subjected to higher local solar concentration (refer toFIG. 3B for the distribution of local concentration ratio in a symmetric wavy pipe, such as the pipes shown inFIGS. 1B and 2A ). Similar observations were made from a direct numerical simulation of forced convection in a wavy channel.FIG. 13 depicts the fluid temperature distribution at the boundary with the absorber pipe whose hot spots are observed at the throats according to its local concentration profile. As a result, ηc increases along with N and A and remains relatively high even at high mass flow rates unlike the conventional straight absorber pipe. -
FIG. 14 depicts the pressure drop (ΔP) as a function of mass flow rate for different N and A according to which the difference in ΔP between N=0 and the rest is remarkably high, especially at {dot over (m)}f=10 kg/s since ΔP is proportional to the square of the flow rate. The exergy destruction due to pressure drop ({dot over (X)}d,ΔP) is therefore pronounced at a high {dot over (m)}f (refer toFIG. 11 ), N, and A, causing ηex to decrease rapidly with respect to increasing N and A as shown inFIG. 9B . The peak ΔP is observed with N=10 and A=0.01 m and is approximately 20 times greater than that observed across the straight absorber pipe. - The PEC given by Eq. (28), which accounts for both heat transfer enhancement and pressure drop in the modeled symmetric wavy absorber pipe design relative to the straight absorber pipe, is plotted in
FIG. 15 against the fluid mass flow rate for different N and A. The physical implications of the qualitative trends inFIG. 15 can be discerned based on the results presented above, namely ηc inFIGS. 7A-9B and ΔP inFIG. 14 . The peak PEC is achieved with N=5 and A=0.003 m as the design exhibits the least pressure drop, and the PEC is above unity when the heat transfer enhancement outweighs the increased pressure drop at -
- and A=0.003 m for all N. The PEC increases along with {dot over (m)}f for all N as more useful heat is gained by the fluid. When A>0.006 m and {dot over (m)}f<5 kg/s, however, PEC≤1 due to an escalated pressure drop which offsets the heat transfer enhancement. From a thermal-hydraulic standpoint, N=5 and A=0.003 m is therefore the most preferred configuration in the design space explored herein.
- The performance of the symmetric wavy absorber pipe, such as the pipes shown in
FIGS. 1B and 2A , was evaluated herein by means of an experimentally validated 3D PTC model coupled to a semi-analytical optical model. Subsequently, the following observations and conclusions were deduced from this study: - 1. Significant improvement in ηc was observed with the modeled symmetric wavy absorber pipe design, especially at high {dot over (m)}f due to the accelerated flow across the pipe throat which in turn enhanced the mixing and heat transfer between the absorber pipe and the working fluid. For example, collector efficiency improves by 20% or more with the symmetrical wavy absorber pipe compared with using a conventional straight absorber pipe.
- 2. The variation of ηex shifted from an increasing to a decreasing trend with respect to A as {dot over (m)}f increased and the exergy destroyed due to pressure drop became evident accordingly.
- 3. The PEC verified the superior performance of the modeled symmetric wavy absorber pipe design over the conventional straight absorber pipe as the peak PEC≈1.25 despite its increased pressure drop across the pipe. However, PEC≤1 for all N and {dot over (m)}f when A>0.006 m due to an escalated pressure drop that overturned the heat transfer enhancement.
- In various implementations, a symmetric wavy absorber pipe, such as the pipes shown in
FIGS. 1B and 2A , with a smaller N and A enhances heat transfer while diminishing the effect of pressure drop as compared to the conventional straight absorber pipe. - A number of implementations have been described. Nevertheless, it will be understood that various modifications may be made without departing from the spirit and scope of the claims. Accordingly, other implementations are within the scope of the following claims.
- Certain terminology is used herein for convenience only and is not to be taken as a limitation on the present claims. In the drawings, the same reference numbers are employed for designating the same elements throughout the several figures. A number of examples are provided, nevertheless, it will be understood that various modifications can be made without departing from the spirit and scope of the disclosure herein. As used in the specification, and in the appended claims, the singular forms “a,” “an,” “the” include plural referents unless the context clearly dictates otherwise. The term “comprising” and variations thereof as used herein is used synonymously with the term “including” and variations thereof and are open, non-limiting terms. Although the terms “comprising” and “including” have been used herein to describe various implementations, the terms “consisting essentially of” and “consisting of” can be used in place of “comprising” and “including” to provide for more specific implementations and are also disclosed.
Claims (14)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/219,735 US11739984B2 (en) | 2020-03-31 | 2021-03-31 | Solar energy collection system with symmetric wavy absorber pipe |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US202063002834P | 2020-03-31 | 2020-03-31 | |
US17/219,735 US11739984B2 (en) | 2020-03-31 | 2021-03-31 | Solar energy collection system with symmetric wavy absorber pipe |
Publications (2)
Publication Number | Publication Date |
---|---|
US20210310699A1 true US20210310699A1 (en) | 2021-10-07 |
US11739984B2 US11739984B2 (en) | 2023-08-29 |
Family
ID=77921549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/219,735 Active 2041-07-25 US11739984B2 (en) | 2020-03-31 | 2021-03-31 | Solar energy collection system with symmetric wavy absorber pipe |
Country Status (1)
Country | Link |
---|---|
US (1) | US11739984B2 (en) |
Citations (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US1989999A (en) * | 1933-06-07 | 1935-02-05 | Niederle Max | Solar water heater |
US3227153A (en) * | 1963-09-04 | 1966-01-04 | American Mach & Foundry | Solar collector |
US3730229A (en) * | 1971-03-11 | 1973-05-01 | Turbotec Inc | Tubing unit with helically corrugated tube and method for making same |
US3954097A (en) * | 1974-09-16 | 1976-05-04 | Wilson Solar Kinetics | Solar heat collector |
US4007774A (en) * | 1975-09-23 | 1977-02-15 | Uop Inc. | Heat exchange apparatus and method of controlling fouling therein |
US4085491A (en) * | 1976-03-01 | 1978-04-25 | Mathwig Richard G | Solar energy heat exchanger and method for making same |
US4134392A (en) * | 1977-09-19 | 1979-01-16 | Spectrum Conversion, Inc. | Solar energy collection |
US4261333A (en) * | 1979-12-03 | 1981-04-14 | Stephen Guarnaschelli | Solar heat exchanger |
US4469088A (en) * | 1981-03-06 | 1984-09-04 | The Agency Of Industrial Science And Technology | Solar collector |
US4554908A (en) * | 1979-08-29 | 1985-11-26 | Alpha-Omega Development Inc. | Electromagnetic energy absorber |
US4558736A (en) * | 1982-10-21 | 1985-12-17 | Kabelmetal Electro Gesellschaft Mit Beschrankter Haftung | Heat exchanger, particularly for solar collectors and method of manufacture therefor |
US5181560A (en) * | 1990-10-17 | 1993-01-26 | Burn Mark N | Baffleless tube and shell heat exchanger having fluted tubes |
US5727585A (en) * | 1994-06-29 | 1998-03-17 | Daume; Jochen | Device for obtaining energy from sunlight with at least one solar collector |
US20020005275A1 (en) * | 1998-12-04 | 2002-01-17 | Beckett Gas. Inc. | Heat exchanger tube with integral restricting and turbulating structure |
US6637427B1 (en) * | 1999-04-22 | 2003-10-28 | Allan James Yeomans | Radiant energy absorbers |
US6877508B2 (en) * | 2002-11-22 | 2005-04-12 | The Boeing Company | Expansion bellows for use in solar molten salt piping and valves |
US7093619B2 (en) * | 2002-03-25 | 2006-08-22 | Sumitomo Electric Industries, Ltd. | Vacuum-insulated pipe |
US20080029243A1 (en) * | 2003-11-25 | 2008-02-07 | O'donnell Michael J | Heat exchanger tube with integral restricting and turbulating structure |
US20100126498A1 (en) * | 2006-12-07 | 2010-05-27 | Krones Ag | Apparatus for Generating Process Heat for a Packaging Arrangement |
US20100206518A1 (en) * | 2009-02-12 | 2010-08-19 | Patten James W | Corrugated heating conduit and method of using in thermal expansion and subsidence mitigation |
US20120031095A1 (en) * | 2009-01-08 | 2012-02-09 | Airlight Energy Ip Sa | Absorber pipe for the trough collector of a solar power plant |
US20130152917A1 (en) * | 2010-06-11 | 2013-06-20 | Centre National De La Recherche Scientifique | Solar receiver absorber and solar receiver comprising at least one such absorber |
US20150179910A1 (en) * | 2013-12-23 | 2015-06-25 | United Arab Emirates University | System For Converting Thermal Energy Into Electrical Energy |
US20150337811A1 (en) * | 2013-02-05 | 2015-11-26 | Zhongying Changjiang International New Energy Investment Co., Ltd. | Solar automatic heat collecting and equalizing tube, automatic heat equalizing trough-type module, solar-thermal complementary power generation system comprising the same, and power generation method using the same |
WO2016124207A1 (en) * | 2015-02-04 | 2016-08-11 | Demagh Yassine | Energy absorber provided with a tube having s-shaped curves for use in solar collectors |
US9624911B1 (en) * | 2012-10-26 | 2017-04-18 | Sunfolding, Llc | Fluidic solar actuator |
US20180187925A1 (en) * | 2010-12-27 | 2018-07-05 | Jeffrey Michael Citron | Trough shaped fresnel reflector solar concentrator |
US20200018938A1 (en) * | 2018-07-14 | 2020-01-16 | Donald Bennett Hilliard | Heat Concentrating Mirror |
US11319937B2 (en) * | 2019-10-16 | 2022-05-03 | Hpm Holding B.V. | Heat receiver for urban concentrated solar power |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2738167Y (en) | 2004-10-22 | 2005-11-02 | 张建城 | High temperature resisting and anti high pressure metal runner solar vacuum heat collecting tube |
CN200949937Y (en) | 2006-08-21 | 2007-09-19 | 河海大学 | Corrugated tube double-layer glass vacuum thermal-collecting tube |
EP2754977A1 (en) | 2013-01-11 | 2014-07-16 | Tre A S.r.l. | A solar collector having a corrugated tube |
CN105066476A (en) | 2015-08-12 | 2015-11-18 | 哈尔滨工业大学(威海) | Corrugated tube type solar heat absorber with functions of strengthening heat transferring and reducing heat deformation |
CN106016780A (en) | 2016-06-01 | 2016-10-12 | 哈尔滨工业大学(威海) | Asymmetric corrugated pipe type solar heat absorber with heat transfer strengthening capacity |
-
2021
- 2021-03-31 US US17/219,735 patent/US11739984B2/en active Active
Patent Citations (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US1989999A (en) * | 1933-06-07 | 1935-02-05 | Niederle Max | Solar water heater |
US3227153A (en) * | 1963-09-04 | 1966-01-04 | American Mach & Foundry | Solar collector |
US3730229A (en) * | 1971-03-11 | 1973-05-01 | Turbotec Inc | Tubing unit with helically corrugated tube and method for making same |
US3954097A (en) * | 1974-09-16 | 1976-05-04 | Wilson Solar Kinetics | Solar heat collector |
US4007774A (en) * | 1975-09-23 | 1977-02-15 | Uop Inc. | Heat exchange apparatus and method of controlling fouling therein |
US4085491A (en) * | 1976-03-01 | 1978-04-25 | Mathwig Richard G | Solar energy heat exchanger and method for making same |
US4134392A (en) * | 1977-09-19 | 1979-01-16 | Spectrum Conversion, Inc. | Solar energy collection |
US4554908A (en) * | 1979-08-29 | 1985-11-26 | Alpha-Omega Development Inc. | Electromagnetic energy absorber |
US4261333A (en) * | 1979-12-03 | 1981-04-14 | Stephen Guarnaschelli | Solar heat exchanger |
US4469088A (en) * | 1981-03-06 | 1984-09-04 | The Agency Of Industrial Science And Technology | Solar collector |
US4558736A (en) * | 1982-10-21 | 1985-12-17 | Kabelmetal Electro Gesellschaft Mit Beschrankter Haftung | Heat exchanger, particularly for solar collectors and method of manufacture therefor |
US5181560A (en) * | 1990-10-17 | 1993-01-26 | Burn Mark N | Baffleless tube and shell heat exchanger having fluted tubes |
US5727585A (en) * | 1994-06-29 | 1998-03-17 | Daume; Jochen | Device for obtaining energy from sunlight with at least one solar collector |
US20020005275A1 (en) * | 1998-12-04 | 2002-01-17 | Beckett Gas. Inc. | Heat exchanger tube with integral restricting and turbulating structure |
US6637427B1 (en) * | 1999-04-22 | 2003-10-28 | Allan James Yeomans | Radiant energy absorbers |
US7093619B2 (en) * | 2002-03-25 | 2006-08-22 | Sumitomo Electric Industries, Ltd. | Vacuum-insulated pipe |
US6877508B2 (en) * | 2002-11-22 | 2005-04-12 | The Boeing Company | Expansion bellows for use in solar molten salt piping and valves |
US20080029243A1 (en) * | 2003-11-25 | 2008-02-07 | O'donnell Michael J | Heat exchanger tube with integral restricting and turbulating structure |
US20100126498A1 (en) * | 2006-12-07 | 2010-05-27 | Krones Ag | Apparatus for Generating Process Heat for a Packaging Arrangement |
US20120031095A1 (en) * | 2009-01-08 | 2012-02-09 | Airlight Energy Ip Sa | Absorber pipe for the trough collector of a solar power plant |
US20100206518A1 (en) * | 2009-02-12 | 2010-08-19 | Patten James W | Corrugated heating conduit and method of using in thermal expansion and subsidence mitigation |
US20130152917A1 (en) * | 2010-06-11 | 2013-06-20 | Centre National De La Recherche Scientifique | Solar receiver absorber and solar receiver comprising at least one such absorber |
US20180187925A1 (en) * | 2010-12-27 | 2018-07-05 | Jeffrey Michael Citron | Trough shaped fresnel reflector solar concentrator |
US9624911B1 (en) * | 2012-10-26 | 2017-04-18 | Sunfolding, Llc | Fluidic solar actuator |
US20150337811A1 (en) * | 2013-02-05 | 2015-11-26 | Zhongying Changjiang International New Energy Investment Co., Ltd. | Solar automatic heat collecting and equalizing tube, automatic heat equalizing trough-type module, solar-thermal complementary power generation system comprising the same, and power generation method using the same |
US20150179910A1 (en) * | 2013-12-23 | 2015-06-25 | United Arab Emirates University | System For Converting Thermal Energy Into Electrical Energy |
WO2016124207A1 (en) * | 2015-02-04 | 2016-08-11 | Demagh Yassine | Energy absorber provided with a tube having s-shaped curves for use in solar collectors |
US20200018938A1 (en) * | 2018-07-14 | 2020-01-16 | Donald Bennett Hilliard | Heat Concentrating Mirror |
US11319937B2 (en) * | 2019-10-16 | 2022-05-03 | Hpm Holding B.V. | Heat receiver for urban concentrated solar power |
Non-Patent Citations (1)
Title |
---|
Machine translation of CN278167Y retrieved 03/22/2023 (Year: 2023) * |
Also Published As
Publication number | Publication date |
---|---|
US11739984B2 (en) | 2023-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Saravanan et al. | Performance assessment in V-trough solar water heater fitted with square and V-cut twisted tape inserts | |
Abed et al. | An extensive review of various technologies for enhancing the thermal and optical performances of parabolic trough collectors | |
Fan et al. | A comparative study on the performance of liquid flat-plate solar collector with a new V-corrugated absorber | |
Bellos et al. | Enhancing the performance of a linear Fresnel reflector using nanofluids and internal finned absorber | |
Bellos et al. | Optical and thermal analysis of a linear Fresnel reflector operating with thermal oil, molten salt and liquid sodium | |
Rodríguez-Sánchez et al. | Thermal design guidelines of solar power towers | |
Korres et al. | A new mini-CPC with a U-type evacuated tube under thermal and optical investigation | |
Osorio et al. | Performance analysis of parabolic trough collectors with double glass envelope | |
CN102084191B (en) | Trough collector for a solar power plant | |
Murugan et al. | Performance enhancement of centrally finned twist inserted solar collector using corrugated booster reflectors | |
El-Said et al. | Performance enhancement of a double pass solar air heater by using curved reflector: Experimental investigation | |
Gong et al. | Comparative study of heat transfer enhancement using different fins in semi-circular absorber tube for large-aperture trough solar concentrator | |
Mohamad et al. | Parabolic trough efficiency gain through use of a cavity absorber with a hot mirror | |
Kulahli et al. | Numerical simulation of a parabolic trough collector containing a novel parabolic reflector with varying focal length | |
Gouthamraj et al. | Design and analysis of rooftop linear fresnel reflector solar concentrator | |
Karwa et al. | Receiver shape optimization for maximizing medium temperature CPC collector efficiency | |
Saravanan et al. | Experimental studies on heat transfer and friction factor characteristics of twist inserted V-trough thermosyphon solar water heating system | |
Fang et al. | Effects of surface optical and radiative properties on the thermal performance of a solar cavity receiver | |
Reddy et al. | Estimation of heat losses due to wind effects from linear parabolic secondary reflector–receiver of solar LFR module | |
Yang et al. | 3D thermal-hydraulic analysis of a symmetric wavy parabolic trough absorber pipe | |
Chen et al. | Heat transfer enhancement for U-pipe evacuated tube solar absorber by high-emissivity coating on metal fin | |
Yehualashet et al. | Experimental study and numerical analysis of thermal performance of corrugated plate solar collector | |
Peng et al. | Evaluation of heat flux distribution uniformity around the receiver tube of parabolic trough collector based on six statistical and geometrical indices | |
Kumar et al. | A review of performance improvements in design features of liquid flat-plate solar collector | |
Kalidasan et al. | Linear cavity solar receivers: A review |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
AS | Assignment |
Owner name: THE FLORIDA STATE UNIVERSITY RESEARCH FOUNDATION, INC., FLORIDA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:YANG, SAM;ORDONEZ, JUAN CARLOS;REEL/FRAME:058826/0848 Effective date: 20220124 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |