Methods for operating solar-thermochemical processes

10801756 ยท 2020-10-13

Assignee

Inventors

Cpc classification

International classification

Abstract

Methods for controlling or operating solar thermochemical reactions process that maximize the two-step thermochemical energy cycle efficiency by a combination of pressure and temperature swing are disclosed.

Claims

1. A method of operating a thermochemical reactor system comprising: setting a thermal reduction temperature, pressure and oxide mass flow for an oxide in the thermochemical reactor; setting an oxidation temperature and oxygen partial pressure in the thermochemical reactor system for a selected feed stream; and operating the thermochemical reactor system with the set thermal reduction temperature, pressure and oxide mass flow for an oxide in the thermochemical reactor and the set oxidation temperature and oxygen partial pressure in the thermochemical reactor system; where an amount of heat to satisfy the oxidation temperature is equal to an amount of heat provided by the oxide mass flow and the thermal reduction temperature.

2. The method of claim 1, wherein the feed stream is carbon dioxide.

3. The method of claim 1, wherein the feed stream is water.

4. The method of claim 3, wherein the oxidation temperature is red is equal to a steam generation for water at the set oxygen partial pressure.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) Embodiments of the present invention are illustrated by way of example, and not limitation, in the figures of the accompanying drawings in which:

(2) FIG. 1 is a reactor schematic according to an embodiment of the disclosure.

(3) FIG. 2a shows the relationship between T.sub.iso(n.sub.w/h) and p.sub.O2.

(4) FIG. 2b shows the practical efficiency limits for isothermal two-step water-splitting cycles.

(5) FIG. 3 shows the efficiency .sub.R(p.sub.TR), for several values of T, .sub.R=0.5, and .sub.GG=0.8, according to an embodiment of the disclosure. In this non-isothermal example, the metal oxide is assumed to be cerium oxide (ceria).

(6) FIG. 4 shows calculated reactor efficiencies as a function of T. Solid lines are for .sub.GG=0.8, whereas dashed lines are for .sub.GG=0.6. Thick lines are for .sub.R=0, and thin ones for .sub.R=0.75. Lines of the same color correspond to the same p.sub.TR. The metal oxide in this example is ceria.

(7) FIG. 5 shows calculated reactor efficiencies (left y-axis) as a function of dT, for .sub.GG linearly increasing from .sub.GG (1773 K)=0.6 to .sub.GG (1273 K)=0.95 The metal oxide in this example is ceria.

(8) FIG. 6 shows .sub.R(T.sub.opt) as function of p.sub.TR, for several .sub.R and .sub.GG. The metal oxide in this example is ceria.

(9) FIG. 7 shows a flow chart of a method of determining the optimal temperature difference between oxidation and reduction for the control of a solar thermochemical process according to an embodiment of the disclosure.

DETAILED DESCRIPTION

(10) The present disclosure is directed to methods for conducting a two-step thermochemical cycle at maximum energy efficiency by controlling solar thermochemical reactions and designing and operating a reactor that drives the two reactions through an appropriate combination of pressure and temperature swings. The term temperature swings, as used herein, refers to the situation wherein the two different reactions (oxidation and reduction) are carried out at different temperatures so that the thermodynamic equilibrium state of each of is more favorable. The difference in temperature between the high temperature (oxide reduction) and low temperature (reoxidation) step is the temperature swing. In the following description, numerous details are set forth. It will be apparent, however, to one skilled in the art, that the present invention may be practiced without these specific details. In some instances, well-known methods and devices are shown in block diagram form, rather than in detail, to avoid obscuring the present invention. Reference throughout this specification to an embodiment means that a particular feature, structure, function, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase in an embodiment in various places throughout this specification are not necessarily referring to the same embodiment of the invention. Furthermore, the particular features, structures, functions, or characteristics may be combined in any suitable manner in one or more embodiments. For example, a first embodiment may be combined with a second embodiment anywhere the two embodiments are not mutually exclusive.

(11) The terms coupled and connected, along with their derivatives, may be used herein to describe structural relationships between components. It should be understood that these terms are not intended as synonyms for each other. Rather, in particular embodiments, connected may be used to indicate that two or more elements are in direct physical or electrical contact with each other. Coupled may be used to indicated that two or more elements are in either direct or indirect (with other intervening elements between them) physical or electrical contact with each other, and/or that the two or more elements co-operate or interact with each other (e.g., as in a cause and effect relationship).

(12) The reactive particles applicable to the systems and techniques described herein may generally be of any type known for thermochemical reactions that are further suitable for conveyance by the systems and techniques described herein. In an embodiment, the reactive particles may be metal oxides, for example, but not limited to binary, ternary, and quaternary metal oxides, doped or undoped spinels, perovskites, brownmillerites, or other particles including a material having a composition capable of cyclic redox reactions. In this disclosure, exemplary embodiments utilizing a metal oxide (MO.sub.x), such as ceria (CeO.sub.2), ferrites, manganites, cobaltites, perovskites and the like are disclosed, however, alternative embodiments may employ any known particle composition capable of similar cyclic redox reactions. Reactive particles applicable to the systems and techniques described herein may also vary in size significantly with smaller sizes having higher surface/volume ratios improving reaction rates, but potentially being more susceptible to sintering and/or melting. For one exemplary ceria particle embodiment, particle size is between about 5 m (microns) and 500 m (microns).

(13) Although the reactive particles are not consumed significantly with each reaction cycle in the exemplary embodiments described herein, one of skill in the art will note the systems and techniques described herein enable particle continuous addition and extraction and are therefore readily adaptable to embodiments where the reactive particles may be consumed (e.g., attritted or volatilized) and replenished. Reactive particles applicable to the systems and techniques described herein may be a solid media of homogenous or heterogeneous composition (e.g., carrier media coated with reactive media) and of various porosity.

(14) The two step process may be applied to various two-step thermochemical cycle processes, such as, but not limited to solar fuel production, energy storage, and air separations. In an embodiment, the two-step thermochemical cycle is for solar fuel production, and may be a water splitting process or a carbon dioxide splitting process. In water splitting, the reduced metal oxide is reoxidized with steam at a specific temperature to yield H.sub.2. Some steam may be unreacted. The newly oxidized material is separated from the steam/H.sub.2 gaseous environment and is then thermally reduced by raising the temperature, lowering the oxygen partial pressure, or as described herein a combination of both. The evolved oxygen is pumped or swept away and the material is then ready for reoxidation with steam after the temperature is lowered as prescribed. If material flow through the cycle is continuous, the heat released during cooling may be recuperated to help drive thermal reduction. The net result of the process is that a combination of heat and oxygen pressure manipulation split water into separate streams containing H.sub.2 and O.sub.2. In carbon dioxide splitting, the process is entirely analogous with CO.sub.2 filling the role of H.sub.2O, and CO and O2 being the products.

(15) The disclosed methods determine the operating conditions and control parameters for operating a two-step thermochemical cycle at a point wherein the largest fraction of solar energy is converted to chemical energy.

(16) The methods determine the preferred combination of temperature swing and pressure swing between the reduction and oxidation steps for a cyclic two-step thermochemical process for a given set of the reduction temperature, and the effectiveness of recuperation of heat from the gaseous and solid streams. The preferred combination is that which provides the greatest chemical energy output per unit solar energy input. The thermal reduction temperature will typically be at or near an upper limit determined by materials or system limitations while the recuperation will be variable based on design and cost considerations for example; foregoing recuperation is one possible configuration. The methods are applicable to any solid phase metal-oxide system and any reactor embodiment that includes the option of both temperature and pressure swing.

(17) The methods for given reactive material include the following steps: 1. Determine the thermal reduction temperature based on materials limitations or other constraints, e.g. optimizing thermal receiver performance. 2. Determine the operating pressure (oxygen partial pressure) in the thermal reduction step. 3. Determine the solid-solid and gas-gas recuperation efficiencies. 4. Calculate the optimum temperature for the oxidation step. 5. Iterate/repeat to develop efficiency map as a function of reduction temperature, temperature swing, pressure swing, and recuperation efficiencies. 6. For water splitting, the preferred oxidation temperature is that wherein the oxide heating and steam generation heat requirements are equal or approximately equal. 7. Performance/efficiency mapping can be embedded into larger system model including the solar field, receiver, etc. to determine preferred thermal reduction temperature for overall system efficiency. 8. Performance/efficiency mapping can be embedded into larger system cost models to determine preferred configuration of thermal recuperation, pump sizes, etc. that gives the largest amount of stored chemical energy per unit of capital cost.

(18) The methods for operating the two-step thermochemical cycle process are further disclosed by referencing a two-step solar thermochemical reactive metal oxide water splitting cycle discussed below.

(19) Reactions (1) and (2) below, generically describe a two-step thermochemical cycle for H.sub.2 production via water splitting, based on a reactive metal oxide (MO.sub.x). Reaction (1) is an endothermic thermal reduction of the oxide carried out at a temperature T.sub.TR and pressure P.sub.TR. Reaction (2) is the mildly exothermic reoxidation of the reduced oxide with H.sub.2O, at temperature T.sub.WS, which yields H.sub.2 and restores the oxide to its initial state. In carbon dioxide splitting, the process can also be used to produce CO from CO.sub.2 in reaction (2). The sum of the two reactions is heat-driven H.sub.2O splitting, described by reaction (3):

(20) 1 TR - WS MO x - WS .fwdarw. 1 TR - WS MO x - TR + 1 2 O 2 ( 1 ) thermal reduction at T TR , P TR , 1 TR - WS MO x - TR + H 2 O .fwdarw. 1 TR - WS MO x - WS + H 2 ( 2 ) H 2 O splitting at T WS , H 2 O .fwdarw. H 2 + 1 2 O 2 ( 3 ) H 2 O thermolysis .

(21) Here, .sub.TR and .sub.WS are the extents of reduction of the oxide following the thermal reduction and water splitting steps. Their difference, =.sub.TR.sub.WS, is the reversible oxygen capacity, as realized in the cycle. The oxygen partial pressure for the reduction reaction is p.sub.O2.

(22) At any temperature the Gibbs free energies of the above reactions are related:
G.sub.1=G.sub.3G.sub.2(4),
meaning that reactions (1) and (2) are thermodynamically spontaneous only in distinct and non-intersecting regions, i.e. G.sub.1(T.sub.TR), G.sub.2 (T.sub.WS)<0 only for T=T.sub.TRT.sub.WS>0. If the cycle is carried out in these temperature regions, it is driven entirely by thermal energy, i.e. a primary heat source. However, cyclically heating and cooling the oxide between T.sub.TR and T.sub.WS opens the possibility for thermal losses that may render H.sub.2 production inefficient in practice.

(23) At temperatures where reactions (1) and (2) are not spontaneous, additional energy must be provided to the system (generally in some form other than heat, i.e. a secondary energy source) in order to drive the reactions towards the same endpoints. The minimum amount of thermal energy, Q.sub.min that must be supplied to carry out the cycle in this general case is:

(24) Q min ( T TR , T WS ) = H 1 ( T TR ) + ( 1 - .Math. R ) F R C p T + G 1 ( T TR ) 1 + G 2 ( T WS ) 2 . ( 5 )

(25) The 1.sup.st term on the right-hand side is simply the reduction endotherm. The 2.sup.nd term accounts for thermal cycling, where C.sub.p is the molar specific heat capacity of the reactive oxide, F.sub.R is the molar fraction of the solid that is reactive (as opposed to inert), and .sub.R is the effectiveness of solid-solid heat recovery. The 3.sup.rd and 4.sup.th terms generically (the process is not specified) account for the scenario(s) where reactions (1) and/or (2) are carried out at unfavorable temperatures, i.e. when G.sub.1(T.sub.TR)>0 and/or G.sub.2(T.sub.WS)>0, and represent the minimum work must be performed to drive reaction (1) and/or (2). This work is determined by G.sub.1 and/or G.sub.2 at the respective conditions. Coefficients .sub.1 and .sub.2 account for losses in converting heat to the required amount of work. Additionally, they implicitly include the possibility that waste heat (i.e. heat of reoxidation and unrecovered oxide sensible heat) can be used to provide part of the work in the 3.sup.rd and 4.sup.th term. (See also discussion of eq. 11.) The limiting case (Q.sub.min=H.sub.1) applies to a system operating in the favorable temperature regimes with ideal heat recovery (.sub.R=1).

(26) To understand the implications of eq. 5, consider the situation where T=0 (ITWS), at a temperature intermediate to those of the thermodynamically favorable regimes. In this case, no energy is required for thermal cycling (the 2nd term is zero), but work must be added to both reactions (1) and (2). As T.sub.TR=T.sub.WS change in value, the 2nd term remains zero and G.sub.1 and G.sub.2 increase and decrease in opposition to one another. As discussed below, the thermodynamic limits of this special case can be determined without any knowledge of the reactor design or properties of the working materials.

(27) More interestingly, from an application perspective, if T.sub.TR is increased and T.sub.WS is decreased (TSWS), the thermal cycling energy requirement (2nd term) becomes nonzero and grows, but the work requirements embodied by G.sub.1 and G.sub.1 decrease and are eventually eliminated. Thus, eq. 5 suggests that for a given material (defining the thermodynamics) and set of reactor and process designs (defining the various efficiencies, practical and physical limitations, parasitic losses, etc.) there should be such T.sub.TR and T.sub.WS, that minimize the sum of the last three terms in eq. 5, and therefore Q.sub.min.

(28) Determining the conditions under which a two-step reactor operates at maximum efficiency requires the knowledge of the thermodynamic properties of the working material, as well as a reactor model. Specific thermodynamic properties of reactive oxides are not generally known, and are typically determined by extensive experimentation. Generalizations can therefore be made, but there is no universal solution.

(29) For the purpose of demonstrating efficiency calculations, we use assume CeO.sub.2 as the active oxide, a well characterized candidate material for solar thermochemical water splitting. The thermodynamics of CeO.sub.2 are obtained from the work of Zinkevich et al. (M. Zinkevich, D. Djurovic, F. Aldinger Thermodynamic Modelling of the Cerium-Oxygen System Solid State Ionics, 2006; 177, 989.). These authors performed a comprehensive, critical review of the thermodynamic literature concerning cerium oxides and applied a Calphad analysis to derive models covering all relevant phases (including both liquids and solids). This comprehensive approach covers the entire temperature range relevant to our analysis, and is necessary to accurately describe non-stoichiometric phases, such as the ceria fluorite phase (CeO.sub.2-) of concern here, which persists from =0 to 0.3 at 1773 K. This model also generally predicts lower values of (p,T) than models based solely on the experiments by Panlener (R. J. Panlener, R. N. Blumenthal, J. E. Gamier A Thermodynamic Study of Nonstoichiometric Cerium Dioxide Journal of Physics and Chemistry of Solids 1975, 36, 1213.), used in almost all previous studies. This difference leads to lower overall efficiency predictions for the comprehensive approach, but the conclusions of the analysis remain unchanged by the choice of model. The difference between ceria reduction models is illustrated in FIG. 3.

(30) FIG. 1 is a schematic of a reactor according to an embodiment of the disclosure. As can be seen in FIG. 1, heat exchangers and their effectiveness (.sub.R and .sub.GG) are indicated in the oxide and steam flows. The input solar heat power {dot over (Q)}.sub.TH is used to add sensible heat to the ceria (Q.sub.SH) and to partially reduce it (Q.sub.TR). If the heat from the H.sub.2 production reaction (Q.sub.MOX) and the unrecovered sensible heat of the oxide (Q.sub.SH,L), is insufficient to heat steam to T.sub.WS (Q.sub.H2O), part of {dot over (Q)}.sub.TH is used. For part of T.sub.H simplicity and because it is of minor importance, heat recovery from the O.sub.2 stream is not shown, but it is included in the calculations (cf. eq. 11).

(31) The reactor is assumed to accomplish four primary unit operations, common to many thermochemical reactor designs. These are a thermal reduction, solid-solid heat exchange in a thermal recuperator, H.sub.2 production (water splitting), and steam pre-heating. During operation, concentrated solar radiation heats and thermally reduces the reactive oxide in a thermal reduction chamber. The oxide then moves through the recuperator (entering at the hot inlet) and then into an H.sub.2 production chamber, where it is exposed to a pre-heated steam flow in a countercurrent arrangement, producing H.sub.2. The reoxidized material is then brought back to the reduction chamber, via the recuperator, where heat is exchanged between the two oxide flows.

(32) The reactor heat-to-H.sub.2 efficiency or reactor efficiency, .sub.R, is defined as:

(33) R = n . H 2 HHV H 2 Q . A , ( 6 )
where {dot over (n)}.sub.H.sub.2 is the hydrogen molar production rate, HHV.sub.H2 is its higher heating value, and {dot over (Q)}.sub.A is the solar power at the reactor aperture. Efficiency calculations incorporate an important assumption: Only solar primary energy is used for the entire operation of the reactor. This includes all heat needs, such as thermal reduction, ceria and feedstock heating, and heat equivalents of mechanical work, such as vacuum pumping, compression, oxide moving, etc. Equally importantly, all major losses and inefficiencies are also included, (e.g. black body radiation through the aperture at a concentration ratio C.sub.A=3000, conversion efficiency from heat to mechanical work, etc.) Radiation losses are included in the reactor efficiency for completeness and for ease of comparison with other work. The sensible heat of the product gasses is used to preheat steam (with a gas-gas heat recovery effectiveness .sub.GG), and to generate power. Conduction losses through insulated reactor walls are considered negligible in a large device. The sensible heat in the O.sub.2 product is comparatively small, and its omission or inclusion has little effect on the efficiency calculations.

(34) The H.sub.2 molar production rate in eq. 6 can be expressed in terms of the heat input power ({dot over (Q)}.sub.TH) and the heat required for production of 1 mol H.sub.2 (Q.sub.mol):

(35) n . H 2 = Q . TH Q mol . ( 7 )

(36) After losses to aperture intercept (A=0.95) and thermal re-radiation (P.sub.rad), {dot over (Q)}.sub.TH can be expressed as:
{dot over (Q)}.sub.TH=A*{dot over (Q)}.sub.AP.sub.rad(8),
whereas Q.sub.mol is:
Q.sub.mol=Q.sub.TR+Q.sub.SH+Q.sub.AUX(9)

(37) Here the individual terms (roughly ordered by decreasing temperature) correspond to those in eq. 5 as follows: Q.sub.TR=H.sub.r(CeO.sub.2) is the thermal reduction endotherm. The energy required for heating the oxide (sensible heat) from T.sub.WS to T.sub.TR (assuming F.sub.R=1) is:

(38) Q SH = C p T ( 1 - .Math. R ) , ( 10 )
where the molar heat capacity of CeO.sub.2 is C.sub.p80 J/mol K..sup.22 Finally, Q.sub.AUX encompasses the heat equivalents of other, auxiliary, energy requirements:
Q.sub.AUX(Q.sub.H.sub.2.sub.O+Q.sub.pump+Q.sub.mech+Q.sub.sep)(Q.sub.MOX+Q.sub.SH,L+Q.sub.O.sub.2)(11).

(39) Here Q.sub.H2O is the energy required to heat steam by T.sub.I/O=T.sub.WST.sub.0, i.e. from ambient temperature (T.sub.0) to T.sub.WS, and it includes preheating by hot product streams. The heat equivalents of the pumping of products (in both chambers) and mechanical and separation work are Q.sub.pump, Q.sub.mech, and Q.sub.sep, respectively. The negative terms represent the waste heat available from the product gasses, mainly the H.sub.2OH.sub.2 mix, which consists of the heat released at T.sub.WS in the reoxidation reaction, Q.sub.MOX=H.sub.rH.sub.cH2.sup.0, and the unrecovered sensible heat of the oxide, Q.sub.SH,L. Steam, in the fuel production chamber, acts as both a reactant (oxidizer) and a coolant. The sensible heat in the oxygen exhaust is Q.sub.O2.

(40) Importantly, Q.sub.AUX is forced to be non-negative, i.e. it is set to zero when the waste heat exceeds the first three terms in eq. 11, since heat at T.sub.WS cannot contribute to either Q.sub.TR or Q.sub.SH. The quantities in eq. 11 are heat equivalents, so conversion efficiency terms are included where applicable, such as the conversion of heat to mechanical or pump work. An efficiency of 10% was used for heat-to-pump work and for the oxide moving work. Thermal reduction and water splitting are assumed to end in their thermodynamic equilibrium states, i.e. kinetic limitations are not considered.

(41) The appeal of ITWS lies in the perceived simplification of reactor design and operation, as it eliminates the need for solid-solid heat recovery and, depending on the design, the frequent temperature cycling of reactor components. Coincidentally, this special case lends itself to straightforward theoretical analysis. To begin, we use well-known relationships for each of the reactions (1), (2), and (3):
G=RT ln K and G=HTS(12) and (13),
where R is the gas constant. The equilibrium constants for reactions (1) and (2) depend on the reactant and product activities:

(42) K 1 = ( a O 2 ) 1 / 2 ( a MO x - TR ) 1 / ( a MO x - WS ) 1 / and K 2 = ( a H 2 ) ( a MO x - WS ) 1 / ( a H 2 O ) ( a MO x - TR ) 1 / ( 14 ) and ( 15 ) .

(43) At all relevant operating pressures, the gas activities can be expressed as partial pressures:

(44) a O 2 = p O 2 and a H 2 a H 2 O = p H 2 p H 2 O , ( 16 ) and ( 17 )
where p.sub.O.sub.2 is measured relative to standard pressure. Substituting eqs. 12, 13, 16 and 17 into eq. 4 gives:

(45) - RT iso ln ( p O 2 1 / 2 ) - RT iso ln ( a MO x - TR 1 / a MO x - WS 1 / ) = H 3 - T iso S 3 + RT iso ln ( p H 2 p H 2 O ) ++ RT iso ln ( a MO x - WS 1 / a MO x - TR 1 / ) . ( 18 )

(46) Here T.sub.iso, is the isothermal operating temperature. Solving for T.sub.iso is facilitated by the exact cancellation of the oxide terms on the left and right side of eq. 18, giving:

(47) T iso = - H 3 R ( ln ( p O 2 1 / 2 ) - ln ( p H 2 O p H 2 ) ) - S 3 . ( 19 )

(48) FIG. 1a shows the relationship between T.sub.iso(n.sub.w/h) and p.sub.O2. For the water splitting reaction constant thermodynamic values of

(49) 0 H 3 = 250.8 kJ mol and S 3 = 57.35 J mol - K
are assumed. These are the values at 1673 K and best represent the temperature range of practical interest. Adopting values for a different temperature introduces very small differences in the results. The fading rectangle in FIG. 2a roughly indicates a reasonable operating parameter space for a practical reactor and shows that ITWS largely falls outside of it.

(50) FIG. 2b shows practical efficiency limits for ITWS. Efficiency curves are plotted for n.sub.w/h<10.sup.5. The lack of results below certain T.sub.iso for the higher p.sub.O2 values indicates that ITWS is not possible for n.sub.w/h<10.sup.5. The efficiency scale in FIG. 2b was chosen for ease of comparison with FIGS. 3, 4, and 5.

(51) Referring again to FIG. 2a, the relationship between T.sub.iso and n.sub.w/h, for several values of p.sub.O2 (n.sub.w/h=n.sub.H2O/n.sub.H2 resulting from the oxidation step (ideal gas behavior is assumed so p.sub.H2O/p.sub.H2 and n.sub.H2O/n.sub.H2 are used interchangeably) is shown. It is evident from FIG. 2a that relaxing one parameter results in the restriction of one or both of the other two. For example, decreasing T.sub.iso requires either a p.sub.O2 decrease or an n.sub.w/h increase, or both. Importantly, FIG. 1a strictly applies to an ideal isothermal cycle. A real isothermal cycle would require even more stringent operating conditions than an ideal one (e.g. a higher T.sub.iso at the same values of p.sub.O2 and n.sub.w/h), for example due to kinetic limitations.

(52) It is important to understand that the results in FIG. 1a are general in the sense that a cycle represented by a valid combination of the three parameters can be performed by multiple possible materials which must satisfy eq. 4 at T.sub.iso. Conversely, this requirement means that an isothermal cycle, defined by a point in the graph, cannot be realized by arbitrarily chosen materials. The high temperatures (in this case T.sub.TR=T.sub.WS=T.sub.iso), low p.sub.O2, and high n.sub.w/h values indicated in FIG. 1a, raise questions regarding the ultimate feasibility of ITWS. To understand the implications of these results and to outline a realistic space of operating parameters, we consider the practical limitations regarding T.sub.TR, p.sub.O2, and n.sub.w/h.

(53) Increasing T.sub.iso corresponds to less strict p.sub.O2 (higher) and n.sub.w/h (lower) requirements, and isothermal temperatures as high as 2173 K have been considered. However, radiation losses through the reactor aperture, as well as oxide sublimation and reactivity with reactor structures, limit T.sub.TR to no more (and possibly substantially less) than 1773 K in devices of practical relevance. To appreciate the challenge of ITWS under the extreme conditions considered in the literature, it is instructive to note that at 2173 K, ceria has a vapor pressure p.sub.CeO29.3 Pa, leading to a swift and irreversible oxide loss via sublimation, as observed experimentally by Abanades et al.

(54) It is assumed that low p.sub.O2 is achieved by pumping, i.e. lowering of p.sub.TR, the total pressure in this step (therefore p.sub.O2=p.sub.TR). As shown in a previous analysis, the heat equivalent of pump work is not a major contributor to the total energy requirement, but the lowest p.sub.TR is limited by other factors, such as oxygen volumetric flow and entering the molecular flow regime, to no less than 1 Pa.

(55) The alternative, isothermal inert gas sweeping, was examined by Bader et al., who showed that, even under best-case conditions, the amount of required N.sub.2 by far exceeds the amount of the H.sub.2 product (n.sub.N2/n.sub.H2700). In addition to requiring an N.sub.2 purification plant, the only manner of somewhat efficient ITWS was found to require exceptionally high levels of heat recovery (>95%) between the incoming and outgoing N.sub.2 gasat T.sub.iso. Finally, the vast majority of the products are the inert-oxygen mix and steam, not H.sub.2. It can thus be concluded that sweeping is not an option for ITWS.

(56) The above practical limitations regarding T.sub.TR and p.sub.TR give context to the results in FIG. 1a: For example, even for very low p.sub.TR (e.g. 1 Pa), and the highest operationally relevant T.sub.iso of 1773 K, large amounts of steam must be provided for a very low H.sub.2 yield (n.sub.h/w=n.sub.H2/n.sub.H2O<<1); i.e. the majority of the reactor product would then be unreacted steam, not H.sub.2. In short, ITWS requires impractically high thermal reduction temperatures, or exceptionally low thermal reduction pressures, or results in an exceptionally low H.sub.2 fraction in the output stream. As the derivation of eq. 19 shows, these results and associated limitations are a consequence of water thermodynamics and the basic relationship for the Gibbs free energy of the reduction and oxidation reactions (eq. 4), and can therefore not be circumvented by either redox material choice or innovative reactor design.

(57) For a deeper insight into the implications of the tradeoffs associated with ITWS, the low H.sub.2 fraction or high n.sub.w/h values can be viewed in the context of separation work, i.e. the work that must be performed to separate H.sub.2 from the H.sub.2-poor output stream. Separation work is directly related to efficiency, and practical ITWS efficiency limits, as defined in eq. 6, can be estimated by including some well-established efficiencies of the constituent processes. Re-radiation losses are given in eq. 8 (T.sub.iso.sup.4). Separation work (2.sup.nd law) depends on T.sub.iso, n.sub.h/w, and the final H.sub.2 purity, assumed here to be a modest x.sub.H2=99.9%. The practical separation efficiency (i.e. theoretical 2.sup.nd law work vs. actual work) is generally 15%, albeit not at the high temperatures considered here. The heat-to-power efficiency of a Rankine cycle in concentrated solar power plants, necessary to perform pump and separation work, is at best 40%. Finally, neglecting other work and heat requirements, and keeping the assumption for the heat-to-pumping efficiency (10%) from Section 2.2, plotted in FIG. 1b shows the practical limits for ITWS efficiency (h.sub.R).

(58) The results show that the ITWS .sub.R values are low, even at high T.sub.iso and low p.sub.O2. The limits are almost independent of T.sub.iso and depend very weakly on p.sub.O2. In light of the outstanding operating conditions shown in FIG. 1a, low efficiencies should not be surprising. The relative insensitivity of the results on T.sub.iso and p.sub.O2 may, on the other hand, be initially unexpected. This, however, is a manifestation of the mutual dependence of the operating conditions (p.sub.O2, n.sub.w/h, and T.sub.iso), which are dictated solely by water thermodynamics as expressed in eq. 19. Simply put, the results in FIG. 1b show that, whether one in practice accepts higher re-radiation losses (to increase iso), or chooses to invest work in pumping (to lower p.sub.O2) or in H.sub.2 separation (to allow higher n.sub.w/h), the overall ITWS efficiency limit is roughly the same, in what can be viewed as a thermodynamic zero-sum game.

(59) The above results describe the conditions required for ITWS and indicate the practical difficulties in realizing it, including inherent limitations on efficiency. In this section, the effects of operational parameters on efficiency in the general case where T0 are examiner. In this case, both material and reactor assumptions, as discussed above are included. p.sub.TR is limited to values between 1 Pa and 1 kPalow enough to meaningfully facilitate thermal reduction, but not too low to be entirely unfeasible in a reactor in the field. Likewise, T.sub.TR 1773 K is set.

(60) The above results describe the conditions required for ITWS and indicate the practical difficulties in realizing it, including inherent limitations on efficiency. In the following paragraphs, the effects of operational parameters on efficiency in the general case where T0 are examined. In this case, both material and reactor assumptions as detailed above are included. p.sub.TR is limited to values between 1 Pa and 1 kPalow enough to meaningfully facilitate thermal reduction, but not too low to be entirely unfeasible in a reactor in the field. Likewise, we set T.sub.TR 1773 K.

(61) FIG. 3 shows the efficiency .sub.R(p.sub.TR), for several values of T, .sub.R=0.5, and .sub.GG=0.8. In FIG. 3, solid lines show the calculated reactor efficiency as function of p.sub.TR, for several values of T (left y-axis). Note the nil efficiency for ITWS above p.sub.TR75 Pa, and the initial efficiency increase, followed by a decrease, as T increases from 0 K to 400 K. Dashed lines, on the right y-axis, show the extent of thermal reduction, .sub.TR (assuming T.sub.TR=1773 K). The black dashed line is for the comprehensive model used in this work (based on Zinkevich et al.), whereas the white line is based on Panlener only.

(62) As can be seen in FIG. 3, two important features are immediately evident in the plot. Firstly, ITWS (T=0) yields the lowest efficiency for every p.sub.TR in the range. Secondly, for p.sub.TR75 Pa and above, ITWS efficiency is zero, whereas TSWS (T>0) for the same p.sub.TR yields a positive efficiency.

(63) The latter is a special case of a more general result, also evident in FIG. 3: for low values of T, the efficiency is zero above some critical p.sub.TR. This stems from the presence of oxygen in steam as a result of partial dissociation, characterized by a well-known pressure, p.sub.O2,d (e.g. p.sub.O2,d(1773 K)75 Pa) For some values of p.sub.TR and T.sub.WS, the equilibrium oxygen pressure above the oxide, p.sub.O.sub.2.sub.,CeO.sub.2-TR (T.sub.WS), exceeds p.sub.O2,d(T.sub.WS), leading to further oxide reduction (not reoxidation) in steam, and yielding no H.sub.2 product.

(64) T.sub.opt, defined as the value of T(p.sub.TR, .sub.R, .sub.GG) for which Q.sub.min (eq. 5) is the smallest, and efficiency (.sub.R) is the highest, can now be determined. To visually introduce the T.sub.opt concept, the efficiency curves in FIG. 3, from yellow to black, at a constant p.sub.TR can be followed. A trend is evident: .sub.R(T) initially increases and then decreases with T. Focusing on this trend, FIG. 4 shows the efficiency as function of T, .sub.R(T), for .sub.R=0 and .sub.R=0.75, with two values of .sub.GG (0.8 and 0.6) and several values of p.sub.TR. The chosen values for .sub.R and .sub.GG reflect the awareness that heat recovery at high temperature is challenging, and a high effectiveness, however desirable, may not be possible in practice.

(65) A peak efficiency, .sub.R(T.sub.opt), exists for any of the combinations of p.sub.TR, .sub.R, and .sub.GG in FIG. 4. For the chosen range of parameters, T.sub.opt varies roughly between 180 K and 440 K. Furthermore, for most regimes, ITWS efficiency is several-fold lower than that at T.sub.opt, and it also yields the lowest efficiency for all .sub.R and .sub.GG. This is especially pronounced for .sub.GG=0.6. Other trends are also apparent. For example, p.sub.TR decrease and .sub.R decrease generally lead to a T.sub.opt decrease, whereas .sub.GG decrease leads to a T.sub.opt increase. Efficiency generally increases with decreasing p.sub.TR, and increasing .sub.R and .sub.GG. Conversely, the lowest efficiencies would be realized at high (e.g. atmospheric) p.sub.TR, in the absence of heat recovery.

(66) Some of the efficiency curves in FIG. 4 coincide. First, below T.sub.opt(.sub.R=0), .sub.R(.sub.R=0)=.sub.R(.sub.R=0.75). Second, above T.sub.opt(.sub.GG=0.6), .sub.R(.sub.GG=0.6)=.sub.R(.sub.GG=0.8). Finally, curves for p.sub.TR and .sub.R=0 coincide with curves for 10*p.sub.TR and .sub.R=0.75, below T.sub.opt(p.sub.TR, .sub.R=0). The first two coincidences result from the balance of energy requirements in the cycle, including those that lead to low isothermal efficiency, and can be understood by considering the interplay of factors that determine T.sub.opt. First, we recall that the unrecovered sensible heat of the reduced oxide and the heat of the reoxidation reaction are used to preheat the steam feedstock to near T.sub.WS, via the hot H.sub.2OH.sub.2 exhaust (FIG. 1). Furthermore, oxide reduction (Q.sub.TR), oxide heating (Q.sub.SH), and steam heating (Q.sub.H2O) are by far the three largest energy requirements in the system (depending on T). Of these, Q.sub.SH and Q.sub.H2O can be considered parasitic, since they do not directly contribute to H.sub.2 production (whereas Q.sub.TR does, with an efficiency given by HHV.sub.H2/Q.sub.TR). Therefore, at the point where Q.sub.SH and Q.sub.H2O are roughly balanced, the least amount of direct solar input is required for steam heating or the least amount of the oxide sensible heat has to be rejected, thus leading to peak efficiency.

(67) The coinciding efficiencies in FIG. 4 can now be explained as follows: First, below T.sub.opt, H.sub.H2O is the largest energy requirement and it requires direct solar input (in addition to preheating). Specifically, the region of coincidence, where .sub.R(.sub.R=0)=.sub.R(.sub.R=0.75), corresponds to a situation in which direct solar input is required for steam heating irrespective of .sub.R. More generally, below T.sub.opt(.sub.R=0), Q.sub.H2O is sufficiently large to make efficiency independent on .sub.R. The situation corresponds to a low and eventually zero 2.sup.nd right-hand term in eq. 5 (at T=0), and a large 4.sup.th term.

(68) In light of this, the reasons behind the very low efficiency for ITWS are straightforward: While no energy is needed to heat the oxide after the water splitting reaction (Q.sub.SH=0 in FIG. 1), Q.sub.H2O is at its maximum, even with extensive preheating by the H.sub.2OH.sub.2 output. This is because, at T=0, both n.sub.w/h and T.sub.I/O are at their maximum values. Since Q.sub.SH,L=0, virtually all of Q.sub.H2O requires direct solar input. This dominant role of Q.sub.H2O in ITWS (but not for T>T.sub.opt) also strongly cautions against omitting this energy requirement from efficiency calculations, as is sometimes done.

(69) Second, for T>T.sub.opt, oxide heating is the largest energy requirement and it requires direct solar input (to the extent that .sub.R<1). In this region, Q.sub.H2O is small, mainly because, at low T.sub.WS, comparatively little steam is needed to oxidize ceria back to equilibrium (n.sub.w/h is small), but also because T.sub.I/O is smaller, compared to ITWS. Therefore, Q.sub.SH/Q.sub.H2O>1. In fact, at T>T.sub.opt, Q.sub.AUX<0, necessitating heat rejection. Alternatively, this high quality waste heat can be used for other purposes, even though it cannot directly increase reactor efficiency. This corresponds to a small 4.sup.th right-hand term in eq. 5 (compared to the 2.sup.nd term).

(70) Lastly, coinciding efficiency curves for .sub.R=0 and .sub.R=0.75 reflect the particulars of ceria thermodynamics. As it happens, below T.sub.opt(p.sub.TR), decreasing p.sub.TR by a factor of 10 is equivalent, in efficiency terms, to increasing .sub.R from 0 to 0.75.

(71) The high .sub.GG required for efficient ITWS and low T operation warrants some further consideration. High levels of gas-gas heat recovery (>97%) are attainable at temperatures up to roughly 923 K using stainless steel recuperators. At higher temperatures, creep and corrosion limitations require the use of nickel alloys (up to 1273 K), and no recuperators operating at 1773 K or higher have been reported. Since in ITWS Q.sub.H2O is the highest, the role of gas-gas heat exchange is critical. For ITWS at p.sub.TR=1 Pa, T.sub.iso=1773 K and T.sub.I/O=1475 K, nearly 100 mol of steam must be heated per mole of produced H.sub.2 (FIG. 1a). The associated energy requirement considerably exceeds the chemical energy content of the H.sub.2 product, even for high .sub.GG.

(72) The above further underscores the challenges associated with ITWS: If an exceptionally high .sub.GG cannot be achieved in practice, the efficiencies for ITWS would be even lower than those shown in FIG. 4. On the other hand, as T increases (i.e. T.sub.WS decreases), higher .sub.GG values become more realistic. To account for this, we calculate efficiency in a more flexible fashion, by using a variable .sub.GG, such that it is realistic at both high and low T.sub.WS (FIG. 5).

(73) FIG. 5 shows calculated reactor efficiencies (left y-axis) as a function of T, for .sub.GG linearly increasing from .sub.GG (1773 K)=0.6 to .sub.GG (1273 K)=0.95. This roughly corresponds to a situation where heat from the output steam-H.sub.2 mixture is rejected above 1273 K. Solid lines are for .sub.R=0, dashed lines are for .sub.R=0.75. The dash-dot line shows .sub.GG on the right y-axis.

(74) The results in FIG. 5 give insight into what might be expected in practice, and show an even larger efficiency advantage of TSWS over ITWS. Although this analysis is for a pumped reactor, it is worth repeating at this point that efficiently recovering heat from a sweep gas would face the same obstacles as in the case of steam, with the added construction difficulty that it must be performed at T.sub.TR, irrespective of T. However, operating at T.sub.opt and under best-case conditions, would significantly decrease sweep gas requirements compared to ITWS, thus much decreasing, possibly even eliminating the need for heat recovery from the inert-oxygen exhaust and associated hardware complications. An inert gas purification plant would still be necessary.

(75) Solid-solid heat recovery being also challenging, we have consistently included the limiting case of .sub.R=0. Understanding that .sub.R=0.75 is near the upper end of what may be possible in practice, plots such as that in FIG. 5 can be helpful for roughly estimating maximum reactor efficiency: for every p.sub.TR, it lies approximately between the peaks of the corresponding solid and dashed lines. The results in FIGS. 4 and 5 show that for ceria and moderate values of the operating parameters, such as p.sub.TR>10 Pa, .sub.R<0.5, and .sub.GG<0.7, efficiency ranges roughly from 10 to 25%, whereas the T.sub.opt range is roughly between 250 K and 400 K.

(76) It should be noted at this point that the assumption that all reactions end in their thermodynamic equilibrium states is more important in ITWS, when is small and the system operates near equilibrium, than for T>0, when is comparatively large. Therefore, any deviations from this assumption, as is certain to be the case in practice, would further disfavor ITWS.

(77) In addition to maximizing efficiency for a given reactor/material system, knowing T.sub.opt can guide reactor design.

(78) FIG. 6 shows calculated reactor efficiencies at T.sub.opt (ranging from 130 K to 460 K) as function of p.sub.TR. The three groups of lines correspond to three values of .sub.R. Within the groups, solid lines correspond to .sub.GG=0.9, dashed to .sub.GG=0.7, and dash-dot lines to .sub.GG=0.5. It is evident that, for a ceria-based reactor operating at T.sub.opt, .sub.R influences efficiency more than .sub.GG does. For example, increasing .sub.R from 0 to 0.5 allows the same efficiency to be achieved at a roughly 10 times higher p.sub.TR. This type of plot can be used to help evaluate the cost of achieving a certain p.sub.TR, .sub.R, or .sub.GG vs. corresponding efficiency benefits.

SUMMARY

(79) As seen in FIG. 1a, a reactor producing H.sub.2 from H.sub.2O in an isothermal two-step cycle, must operate at high thermal reduction temperatures, very low thermal reduction pressures, or at an exceptionally high n.sub.w/h, or some combination of the three. High operating temperatures place extraordinary demands on reactor materials, and require unrealistically high concentration ratios from the solar collector and concentrator system. Low thermal reduction pressures require very large receivers and pumpsan undertaking of possibly prohibitive cost. Finally, high values of n.sub.w/h make water near ambient temperature, not H.sub.2, the main reactor product. This necessitates either high temperature H.sub.2 separation, or the addition of a high-throughput, high-efficiency steam heat recovery system, which preheats the input water while cooling the output steam-H.sub.2 mix. Each of these factors complicate reactor and plant design and compromise economics: Therefore, when the entire balance of plant is considered, ITWS makes matters more difficult, not easier, as one might hope.

(80) Further, in addition to posing extraordinary design and operational demands, ITWS offers no efficiency payoff. On the contrary, it appears to be the most inefficient fashion of producing H.sub.2 from H.sub.2O in a two-step cycle (FIGS. 3 and 4). Even in the most favorable case, without solid-solid heat recovery, ITWS with a high level of steam-steam heat recovery (.sub.GG=0.8), is less efficient than TSWS at T.sub.opt for a far more plausible .sub.GG=0.6 (FIG. 4). This disadvantage only widens for .sub.R>0 or when a realistic T.sub.WS-dependent .sub.GG is included (FIGS. 4 and 5).

(81) Recalling that the question in solar fuel production is not one of feasibility, but of efficient solar utilization and minimization of the product cost, no ITWS advantages are evident in our analysis. Rather, requiring that T.sub.TR=T.sub.WS seems to be an unnecessary and counterproductive limitation. On the other hand, TSWS at T.sub.opt maximizes solar resource utilization, and the associated low steam requirement simplifies plant design and operation.

(82) The very existence of T.sub.opt may seem peculiar if one thinks of thermochemical reactors as engines that reverse combustion, i.e. where heat is the input and fuel (chemical work) is the output. It may initially appear most plausible that maximizing T would also maximize efficiency. In the case of perfect heat recovery (.sub.R=.sub.GG=1), this would be correct. It is because of the reality of non-ideal heat recovery that T.sub.opt exists at all.

(83) Some general implications regarding reactor operation follow from the above results. Under all conditions, the slopes of the efficiency curves are shallower for T>T.sub.opt than for T<T.sub.opt. Furthermore, the reasons for suboptimal efficiency are different in these two regions. For T>T.sub.opt, the oxide heating requirement results in the production of high quality waste heat (at T.sub.WS), which can be used elsewhere in the plant, even if not for H.sub.2 production directly. For T<T.sub.opt, however, efficiency decreases because of the steam heating requirement, with waste heat available at low temperature (following steam-steam heat recovery) and of little value. This suggests that it may be prudent to err on the side of T>T.sub.opt, rather than the opposite, in order to achieve higher average efficiency under the variable environmental conditions present in practice. Even though the feasibility of thermochemical fuel production using ceria as a working oxide has been demonstrated, a consensus exists regarding the need for material improvements. Most notable of them is the need for a material that provides a higher at a higher p.sub.TR and lower T.sub.TR than, for example, ceria, yet with similar kinetics and stability. In addition to the possibility of less demanding operation (at a lower required T.sub.TR), a higher would enable higher yields, so that less oxide must be heated per unit produced H.sub.2 or fuel in general. Importantly, this increase must not be offset by an equal increase in C.sub.p or more precisely C.sub.pT (cf. eq. 10 and compare, for example, C.sub.p for ceria with that of LaMnO.sub.3 or La.sub.1-xSr.sub.xMnO.sub.3). A material that reduces more easily than ceria (larger under identical T.sub.TR and p.sub.TR) will generally also be more difficult to reoxidize, assuming similar and temperature independent values for S. This assumption is justified as S is largely a function of the evolution of oxygen into the gas phase. Operationally, this suggests that advanced materials may require a lower T.sub.WS (larger T) for reoxidation to achieve the same n.sub.w/h2. In other words, an insufficiently large /C.sub.p increase could be offset by an increase in the 2.sup.nd right-hand term in eq. 5, and therefore decrease efficiency.

(84) Understanding this, it is to be expected that optimal operation with advanced materials is likely to involve an increase of T.sub.opt, not its decrease in the direction of ITWS. Precise details would additionally depend on reactor design and various internal efficiencies such as .sub.R, .sub.GG, etc. This also follows from general thermodynamic consideration regarding the maximum theoretical efficiency of a two-step thermochemical process. Because of strict thermodynamic limitations, ITWS is unlikely to benefit from the use of advanced materials.

(85) While some key design and operating parameters of a thermochemical reactor for two-step H.sub.2 production present obvious tradeoffs between difficulty and efficiency, T is not one of them. For example, the higher the degree of heat recovery (.sub.R or .sub.GG), which is increasingly difficult to accomplish, the higher the efficiency. One must therefore find a balance that minimizes the levelized cost of the H.sub.2 product. The choice of T is easier to make: operating difficulty being largely independent of T, it would always be advantageous to operate at T for which efficiency is the highest.

(86) Finally, it should be noted that even though much of this analysis is specific to CeO.sub.2 as a reactive oxide, it can be applied to any material for which sufficient thermodynamic information is available. With appropriate small adaptations (mainly heat capacity and dissociation coefficients), the analysis can also be applied to the solar-thermochemical production of CO from CO.sub.2, with largely the same conclusions.

(87) Lastly, the conditions (T.sub.iso, p.sub.O2 and n.sub.w/h) required for isothermal two-step thermochemical water splitting (T.sub.TR=T.sub.WS) can be determined based on water thermodynamics and are highly mutually dependent: choosing values for any two defines the third. This analysis shows that isothermal water splitting is impractical, being a choice between high thermal reduction temperatures, very low oxygen partial pressures for thermal reduction, or exceptionally high steam requirements (i.e. high separation work)or some combination of the three.

(88) Isothermal water splitting is substantially less efficient than the same process at T>0. This is true even in the complete absence of solid-solid heat recovery in the latter case, and assuming a high steam heat recovery effectiveness at high temperature. The low efficiency of ITWS is primarily a result of the exceptionally high energy requirement for steam heating in the water splitting step of the cycle.

(89) Given a specific reactor/material combination, an optimal T=T.sub.TRT.sub.WS can be found to maximize efficiency. For reasonable values of process parameters in a ceria-based cycle, this T.sub.opt ranges roughly between 250 K and 400 K, and is expected to increase in well-designed advanced materials.

(90) A combination of pressure and temperature swing, rather than either individually, is by a wide margin the most efficient mode of operation of a two-step cycle thermochemical reactor for hydrogen production. Efficiency being of paramount importance for the practical application of this technology, temperature and pressure swing reactors appear to be the most promising direction for future research and development.

(91) FIG. 7 shows a flow chart of a method of determining the optimal temperature difference between oxidation and reduction for the control of a solar thermochemical process according to an embodiment of the disclosure. In an embodiment, methods are disclosed for determining oxidation and reduction operating temperatures. It has been unexpectedly disclosed that by setting the thermal reduction temperature, pressure and oxide mass flow, and then determining the oxidation temperature where the oxide heating and steam generation heat requirements are equal or approximately equal, peak efficiency of the cycle is attained.

(92) As so described, the artisan will appreciate that with many independent design parameters, a design protocol for the reactors described herein may include first determining thermal reduction properties of the reactive particles (e.g., reaction kinetics at a chosen reduction temperature and oxygen partial pressure). A solar concentrator power, overall particle mass flow rate through the reduction chamber, and feed stock pumping speed may be independently set to establish the desired reduction temperature, pressure, and residence time within the reactor. Similarly, the production temperature, and production residence time may be set by appropriate component sizing and adjusting the flow of the reactant fluid (two more independent parameters).

(93) It is to be understood that the above description is illustrative, and not restrictive. For example, while flow diagrams in the figures show a particular order of operations performed by certain embodiments of the invention, it should be understood that such order is not required (e.g., alternative embodiments may perform the operations in a different order, combine certain operations, overlap certain operations, etc.). Furthermore, many other embodiments will be apparent to those of skill in the art upon reading and understanding the above description. Although the present invention has been described with reference to specific exemplary embodiments, it will be recognized that the invention is not limited to the embodiments described, but can be practiced with modification and alteration within the spirit and scope of the appended claims. The scope of the invention should, therefore, be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled.