The urgent need to mitigate global warming requires a swift reduction in CO2 emissions from fossil fuel utilization. As such, promising technologies that can reduce reliance on fossil fuels and accelerate the shift toward renewable energy are highly desirable. The intermittent nature of renewable electricity from wind and solar requires considerable compensation from other energy sources, limiting their shares in the power grid. Therefore, energy storage systems, which reduce the temporal and spatial imbalances between electricity generation and consumption, will play a critical role in accelerating the renewable transition.
To store different forms of energy, various energy storage approaches have been proposed such as pumped hydro, batteries, compressed air, flywheels, and thermal energy storage (TES). Pumped hydro comprises the vast majority of the current storage capacity. However, due to the need for large reservoirs separated by elevation differences, pumped hydro is typically limited to mountainous areas, as opposed to population centers, and requires high water usage. With practical limits on increasing pumped hydro capacity, other energy storage approaches are needed to fully address the capacity required for a complete transition to renewable energy. Battery storage systems can directly store and supply electricity. However, the cost of electrode materials and the need for extensive electrical equipment make them very expensive. As such, cheaper technologies are necessary to spur 100% renewable adoption. Thermal energy storage (TES), with low-cost and direct uptake of thermal energy, can be complementary to battery storage.
Among these, TES is considered to be one of the most cost-effective approaches to overcoming the intermittency of concentrated solar power. In addition, TES can directly utilize the widely available heat resources from industrial processes such as steel mills, glass furnaces, and thermal power plants. By 2022, TES has the second-highest installed capacity of 234 GWh and is expected to reach 800 GWh by 2030.
TES technology can be further categorized into three different types, i.e., sensible thermal energy storage (STES), latent thermal energy storage (LTES), and thermochemical energy storage (TCES). At present, the only commercially available TES technology is molten salt-based STES. However, the solidification issue at low temperatures and the instability at high temperatures are the main barriers for the molten salt-based STES. The operating temperature window is typically limited to 200˜600° C., leading to a small energy density of 36˜180 kJ/kg. LTES technology relies on the heat uptake and release of phase change materials (PCMs) during their physical phase change, which can achieve a 2˜3 times larger energy density within a narrower temperature window. To date, medium/high-temperature LTES has yet to be implemented due to various technical challenges such as the high-capital cost, poor thermal conductivity, and equipment corrosion by the liquid PCMs.
Various TES materials such as molten-salt, inert particles, concrete, and redox materials have been investigated, each with its advantages and drawbacks. Molten-salt sensible-heat systems are more mature and allow excellent heat transfer. However, salt induced corrosion and limited stability at high temperatures can limit their selection and increase the cost. It is also not mechanically practical to utilize the latent heat of the salts' solid-liquid phase transition. Concrete-based TES systems store the sensible heat from high-quality steam with a low-cost storage medium. However, as with other systems that do not involve phase/chemical changes, the energy density is limited (−0.85 kJ/kg/° C.). Concerns with cracking of the concrete also limit its operating temperature range. Redox-active particles can absorb a significant amount of energy as they release lattice oxygen during energy storage. However, they are often limited by slow redox kinetics, large temperature swings and/or high cost. “High performance” redox-oxides such as Co3O4 are prohibitively expensive. Carbonates of Ca, Sr, or Ba can store energy through the decomposition to their oxides and CO2. However, they tend to deactivate or require very high operating temperatures (>900° C.). Moreover, the need for integration with a CO2 source or storage adds to the system complexity.
In summary, existing TES materials face various limitations, such as high cost, low energy storage density, large temperature swings, high operating temperatures, and/or complex design for long storage durations. As such, significant advancements in TES materials and systems are necessary to increase the power grid robustness, flexibility, and affordability, while maximizing the adoption of renewable energy. These needs and other needs are satisfied by the present disclosure.
In accordance with the purpose(s) of the present disclosure, as embodied and broadly described herein, the disclosure, in one aspect, relates to an oxide-molten salt (OMS) composite material that addresses the shortcomings of existing TES technologies by offering excellent energy storage capacity at low cost. The OMS composite material can be provided in at least three forms, including a porous, perovskite-oxide particle as a matrix having pores filled with unary salt or salt mixtures; an external shell of perovskite-oxide particles and an internal core of salt or a salt mixture; and/or a porous, perovskite-oxide particle as a matrix filled with unary salt or salt mixtures, wherein the whole particle is enclosed with an oxide outer shell. Further in this aspect, the outer shell can be dense but either mixed ionic-electronically conductive or porous.
In another aspect, the disclosed methods and systems are quick and simple to build and can be used for grid-scale energy storage by retrofitting existing steam turbine-based power plants and concentrated solar plants with TCES systems utilizing the disclosed OMS particles. The disclosed OMS particles have a high energy density and a small charge-discharge temperature window. In a further aspect, the disclosed methods, systems, and materials are advantageous over current energy storage technology using batteries since the disclosed OMS composite material poses no danger to the environment.
Other systems, methods, features, and advantages of the present disclosure will be or become apparent to one with skill in the art upon examination of the following drawings and detailed description. It is intended that all such additional systems, methods, features, and advantages be included within this description, be within the scope of the present disclosure, and be protected by the accompanying claims. In addition, all optional and preferred features and modifications of the described embodiments are usable in all aspects of the disclosure taught herein. Furthermore, the individual features of the dependent claims, as well as all optional and preferred features and modifications of the described embodiments are combinable and interchangeable with one another.
Many aspects of the present disclosure can be better understood with reference to the following drawings. The components in the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating the principles of the present disclosure. Moreover, in the drawings, like reference numerals designate corresponding parts throughout the several views.
Additional advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or can be learned by practice of the invention. The advantages of the invention will be realized and attained by means of the elements and combinations particularly pointed out in the appended claims. It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention, as claimed.
In one aspect, disclosed herein is an oxide molten salt (OMS) composite material including a plurality of perovskite-oxide particles in a molten sort. The OMS composite can take several forms including, but not limited to, redox active perovskite-oxide particles having a porous matrix, wherein a plurality of pores in the porous matrix are filled with the molten salt. In some aspects, the OMS composite material can be in the form of a core-shell particle, wherein the particle has an outer shell that includes an oxide, such as, for example a second perovskite oxide, Al2O3, or any combination thereof. In any of these aspects, porosity of the redox active perovskite-oxide can be increased by ball milling. In an aspect, and without wishing to be bound by theory, the redox active perovskite-oxide stores sensible heat and thermochemical heat, while the molten salt stores the sensible heat and latent heat. Further in this aspect, heat can be stored by the disclosed materials via a temperature-swing mode and/or a pressure swing mode such as, for example, a swing between a low temperature and/or an oxygen-depleted atmosphere and a high temperature and/or an oxygen-rich atmosphere. Alternatively, the swing can be between a low temperature and/or a carbon dioxide depleted atmosphere and a high temperature and/or a carbon dioxide rich atmosphere.
In an alternative aspect, redox active perovskite-oxide particles can form a shell surrounding an internal core, wherein the internal core includes the molten salt. In any of these aspects, the molten salt can be a unary salt or a mixture of two or more salts. In an aspect, the molten salt can be selected from a chloride salt, a carbonate salt, a fluoride salt, a sulfate salt, a phosphate salt, a nitrate salt, a tetrafluoroborate salt, or any combination thereof. In another aspect, the chloride salt can be CaCl2, NaCl, KCl, MgCl2, LiCl, ZnCl2, CsCl, SrCl2, CrO3, or any combination thereof. In still another aspect, the carbonate salt can be Li2CO3, Na2CO3, K2CO3, or any combination thereof. In any of these aspects, the molten salt can have a melting temperature of from about 400° C. to about 1000° C., or of about 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, or about 1000° C., or a combination of any of the foregoing values, or a range encompassing any of the foregoing values. In one aspect, the molten salt can make up from about 10 wt % to about 90 wt % of the OMS composite material, or can be about 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, or about 90 wt % of the OMS composite material, or a combination of any of the foregoing values, or a range encompassing any of the foregoing values.
In another aspect, the redox active perovskite-oxide particles can have a perovskite structure, a brownmillerite structure, a double perovskite structure, or any combination thereof. In a further aspect, the perovskite structure can be selected from CaxA1-xMnyB1-yO3-δ, LaxA1-xFeyB1-yO3-δ, SrxA1-xFeyB1-yO3-δ, or a combination thereof; wherein A is selected from the group Ca, Sr, Ba, La, or a combination thereof; and wherein B is selected from Fe, Mn, Al, Ti, Mg, Ce, Co, Cr, or a combination thereof. In an alternative aspect, the double perovskite structure can be CaxA2-xMnAlO5-δ; wherein A is selected from Sr, Ba, or any combination thereof. In any of these aspects, x and y, if present, can independently be from about 0 to about 1, or can be about 0, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, or about 1, or a combination of any of the foregoing values, or a range encompassing any of the foregoing values.
Without wishing to be bound by theory, the disclosed system stores energy via (1) latent heat due to the phase transition of the molten salt, (2) chemical energy from the oxide redox reactions, and (3) the sensible heat of the composite.
In one aspect, the redox active perovskite-oxide particles in any of the disclosed forms can have the formula SrxA1-xFeyB1-yO3-δ; wherein A comprises Ba, Ca, La, or any combination thereof; wherein B comprises Mg, Mn, Ti, or any combination thereof; wherein x is from 0 to 1; wherein y is from 0.125 to 1; and wherein δ is from 0 to 0.5. In an aspect 6 can vary based on whether a system including the particles is in an endothermic charge state or an exothermic discharge state and δ>1 indicates a nonstoichiometric amount of oxygen is present at a given time. In a further aspect, x can be from about 0.125 to about 0.375 and y can be from about 0.125 to about 0.25. In some aspects, 6 can be from about 0 to about 0.25.
In any of these aspects, the redox active perovskite-oxide particles can be selected from SrFeO3-δ, Sr0.75Ba0.25Fe0.875Mg0.125O3-δ, SrFe0.875Mn0.125O3-δ, Sr0.875Ba0.125FeO3-δ, Sr0.75Ba0.25FeO3-δ, Sr0.625Ba0.375FeO3-δ, Sr0.625Ba0.375Fe0.875Mg0.125O3-δ, Sr0.875Ba0.125Fe0.875Mn0.125O3-δ, Sr0.875Ca0.125Fe0.875Mn0.125O3-δ, Sr0.25Ba0.75Fe0.875Mg0.125O3-δ, SrFe0.75Ti0.25O3-δ, Sr0.125Ba0.875Fe0.875Ti0.125O3-δ, Sr0.625Ba0.375Fe0.875Mn0.125O3-δ, Sr0.75Ba0.25Fe0.875Mn0.125O3-δ, Sr0.875Ba0.125Fe0.625Ti0.375O3-δ, Sr0.75Ca0.25Fe0.75Ti0.25O3-δ, Sr0.75Ca0.25Fe0.75Mn0.25O3-δ, Sr0.25Ba0.75Fe0.875Ti0.125O3-δ, Sr0.875Ca0.125Fe0.625Ti0.375O3-δ, Sr0.75Ca0.25Fe0.625Mn0.375O3-δ, Sr0.75La0.25Fe0.375Mn0.62O3-δ, Sr0.875La0.125Fe0.75Mn0.25O3-δ, Sr0.875Ca0.125Fe0.875Mg0.125O3-δ, Sr0.875La0.125Fe0.625Mn0.375O3-δ, Sr0.125Ca0.875Fe0.25Mn0.75O3-δ, Sr0.75La0.25Fe0.75Mg0.25O3-δ, Sr0.25Ca0.75Fe0.25Mn0.75O3-δ, Sr0.75Ba0.25Fe0.625Ti0.375O3-δ, Sr0.75Ba0.25Fe0.5Ti0.5O3-δ, Sr0.25Ca0.75Fe0.375Mn0.625O3-δ, Sr0.625Ca0.375Fe0.625Ti0.375O3-δ, Sr0.375Ca0.625Fe0.125Mn0.875O3-δ, Sr0.375Ca0.625Fe0.5Mn0.5O3-δ, Sr0.625Ca0.375Fe0.5Ti0.5O3-δ, Sr0.25Ca0.75Fe0.625Ti0.375O3-δ, Sr0.5Ca0.5Fe0.625Ti0.375O3-δ, Sr0.25Ca0.75Fe0.625Mn0.375O3-δ, BaFe0.625Mn0.375O3-δ, Sr0.375Ca0.625Fe0.75Ti0.25O3-δ, Sr0.125Ca0.875Fe0.75Mn0.25O3-δ, Sr0.25Ca0.75Fe0.75Mn0.25O3-δ, SrFe0.375Mn0.625O3-δ, Sr0.125Ca0.875Fe0.25Ti0.75O3-δ, Sr0.125Ca0.875Fe0.625Ti0.375O3-δ, Sr0.5La0.5Fe0.125Mn0.875O3-δ, or any combination thereof. In a further aspect, the redox active perovskite-oxide particles can be Sr0.125Ca0.875Fe0.25Mn0.75O3-δ or Sr0.375Ca0.625Fe0.125Mn0.875O3-δ.
In one aspect, the OMS composite material exhibits an oxygen release capacity of from about 0.1 to about 5 wt % of the redox active perovskite-oxide particles, or at least about 2 wt %, or at least about 2.5 wt %. In a further aspect, oxygen release capacity can be measured under a first redox condition of 400° C./0.2 atm O2 and a second redox condition of 800° C./0.01 atm O2. In still another aspect, the OMS composite material exhibits a chemical energy storage density of at least about 80 kJ/kg ABO3 under an isobaric air condition between 400° C. and 800° C., and exhibits a chemical energy storage density of at least about 150 kJ/kg ABO3 under a temperature and oxygen partial pressure swing between 400° C./0.2 atm O2 and 1100° C./0.01 atm O2.
Also disclosed herein is a system for thermochemical energy storage and release, the system including (a) an endothermic charge stage, wherein the OMS composite material of any one of claims 1-17 has a δ of 0 at a beginning time point of the endothermic charge stage, and wherein the OMS composite material is exposed to heat, thereby releasing molecular oxygen such that δ increases to greater than 0 at an ending time point of the endothermic charge stage; and (b) an exothermic discharge stage, wherein the OMS composite material has a δ greater than 0 at a beginning time point of the exothermic discharge stage, wherein molecular oxygen combines with the OMS composite material such that δ decreases to 0 and heat is released at an ending time point of the exothermic discharge stage.
In one aspect, the endothermic charge stage can operate with an oxygen partial pressure of from about 0.01 atm to about 5 atm, or of about 0.01, 0.1, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, or about 5 atm, or a combination of any of the foregoing values, or a range encompassing any of the foregoing values. In another aspect, the endothermic charge stage and the exothermic discharge stage operate with a temperature difference of from about 50° C. to about 500° C., or from about 100° C. to about 200° C.
In one example, as illustrated in
In one aspect, heat is supplied in the endothermic charge stage by solar power, heat radiated by an industrial process such as, for example, steel milling or glass melting, a thermal power plant, or any combination thereof. In an aspect, heat released from the exothermic discharge stage can be used to create steam to power steam turbines. Also disclosed are power plants including the disclosed systems.
In one aspect, and without wishing to be bound by theory, the high tunability of the TES material allows it to operate under a wider range of conditions, as may be suited for specific heat sources and/or heat demands with high efficiency. In an aspect, the temperature difference between the charging and discharging steps can be from about 50° C. to about 500° C. with a difference of 250° C. or less, or can be between about 100° C. and about 200° C. In another aspect, the oxygen release between the charging and discharging steps can be from about 0.1 to about 5 wt % of the perovskite material, or from about 1 wt. % to about 5 wt. %, about 2 to about 3 wt. %, or about 3 to about 5 wt. %. For the preferable perovskite-based materials, the variation in non-stoichiometry between the charging and discharging steps can be from about 0.1 to about 0.5, from about 0.1 to about 0.2, or from about 0.3 to about 0.5. The timing, temperature, and flowrate of air or oxygen-rich gas for the exothermic oxygen regeneration can be adjusted to maintain heat release within a narrow range of temperature (e.g., 200° C. or 150° C., or about 100° C., 50° C., or 25° C.). In still another aspect, oxygen vacancies can also be formed during contact with transient higher temperature source (e.g., 900-1200° C., or 800° C.-1100° C., or temperatures more than 1200° C.), where after the heat or reaction during regeneration is recovered at a lower operating temperature as is suitable for the heat use such as from about 350° C. to about 700° C., or from about 400 to about 500° C., or about 550° C. In still another aspect, the oxygen-rich stream can have a partial pressure of from about 0.01 atm to about 5 atm of O2 and can be an oxygen-rich flue gas containing CO2, H2O, N2, and O2 along with other combustion byproducts. In this embodiment, “oxygen-rich” means a stream with a sufficient oxygen pressure to react with the TES material to fill oxygen vacancies. The stream can contain less, the same or similar, or more oxygen than air, and can be air. In an aspect, the oxygen-rich stream can have an oxygen molar composition of from about 1% to about 100%, or from about 5% to about 95%, or from about 3% to about 21%, about 5% to about 10%, about 2.5% to about 20%, or from about 10% to about 20%. In one embodiment, the stream would be air (with an oxygen content of ˜20 mol %). In another embodiment, enriched or purified oxygen (e.g., from about 90 to about 100% or from about 95% to about 99%) would be contacted to limit the loss of energy to sensible heat of oxygen diluents. In one aspect, such a “purified” stream would be generated by the TES material during the heat uptake mode either recovered directly or after condensation of a steam sweep gas.
In any of these aspects, the oxide can provided as a perovskite, brownmillerite, double perovskite, or any combination thereof. For example, a perovskite oxide can have a formula selected from the group consisting of CaxA1-xMnyB1-yO3-δ, LaxA1-xFeyB1-yO3-δ, SrxA1-xFeyB1-yO3-δ, and a combination thereof. Further in this aspect, A can be selected from the group Ca, Sr, Ba, La, or a combination thereof. In another aspect, B can be independently selected from the group Fe, Mn, Al, Ti, Mg, Ce, Co, Cr, or a combination thereof. In one aspect, a double perovskite oxide can have a formula selected from the group consisting of CaxA2-xMnAlO5-δ. Further in this aspect, A can be selected from Sr, Ba, or any combination thereof. In one aspect, the range of x can be between 0-1, e.g. about 0, about 0.125, about 0.25, about 0.375, about 0.5, about 0.625, about 0.75, about 0.875. or about 1. In another aspect, the range of y can be 0-1, e.g. about 0, about 0.05, about 0.1, about 0.125, about 0.25, about 0.375, about 0.5, about 0.625, about 0.75, about 0.875. or about 1.
In an aspect, the melting temperatures of salt or salt mixture can be within the range of from about 400 to about 1000° C. In a further aspect, the melting temperatures of salt mixtures can a single temperature or a temperature range. In one aspect, the suitable ranges of melting temperature can be from about 50 to about 200° C. or from about 50 to about 100° C. or can be less than about 50° C. In any of these aspects, the salt can be a chloride salt, carbonate salt, fluoride salt, sulfate salt, phosphate salt, nitrate salt, tetrafluoroborate salt, or a combination thereof. In a further example, the chloride salt can be selected from CaCl2, NaCl, KCl, MgCl2, LiCl, ZnCS2, CsCl, SrCl2, CrCl3, or any combination thereof. In another aspect, the carbonate salt can be selected from Li2CO3, Na2Co3, K2O3, or any combination thereof. In still another aspect, the salt content in the OMS composite particle can be from about 10 to about 90 wt %, or from about 30 to about 50 wt. %, or from about 70 to about 90 wt. %.
Additional exemplary oxides and salts are provided in Table 1:
Preferred oxide and salt combinations are provided in Table 2:
Many modifications and other embodiments disclosed herein will come to mind to one skilled in the art to which the disclosed compositions and methods pertain having the benefit of the teachings presented in the foregoing descriptions and the associated drawings. Therefore, it is to be understood that the disclosures are not to be limited to the specific embodiments disclosed and that modifications and other embodiments are intended to be included within the scope of the appended claims. The skilled artisan will recognize many variants and adaptations of the aspects described herein. These variants and adaptations are intended to be included in the teachings of this disclosure and to be encompassed by the claims herein.
Although specific terms are employed herein, they are used in a generic and descriptive sense only and not for purposes of limitation.
As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosure.
Any recited method can be carried out in the order of events recited or in any other order that is logically possible. That is, unless otherwise expressly stated, it is in no way intended that any method or aspect set forth herein be construed as requiring that its steps be performed in a specific order. Accordingly, where a method claim does not specifically state in the claims or descriptions that the steps are to be limited to a specific order, it is no way intended that an order be inferred, in any respect. This holds for any possible non-express basis for interpretation, including matters of logic with respect to arrangement of steps or operational flow, plain meaning derived from grammatical organization or punctuation, or the number or type of aspects described in the specification.
All publications mentioned herein are incorporated herein by reference to disclose and describe the methods and/or materials in connection with which the publications are cited. The publications discussed herein are provided solely for their disclosure prior to the filing date of the present application. Nothing herein is to be construed as an admission that the present invention is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided herein can be different from the actual publication dates, which can require independent confirmation.
While aspects of the present disclosure can be described and claimed in a particular statutory class, such as the system statutory class, this is for convenience only and one of skill in the art will understand that each aspect of the present disclosure can be described and claimed in any statutory class.
It is also to be understood that the terminology used herein is for the purpose of describing particular aspects only and is not intended to be limiting. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the disclosed compositions and methods belong. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the specification and relevant art and should not be interpreted in an idealized or overly formal sense unless expressly defined herein.
Prior to describing the various aspects of the present disclosure, the following definitions are provided and should be used unless otherwise indicated. Additional terms may be defined elsewhere in the present disclosure.
As used herein, “comprising” is to be interpreted as specifying the presence of the stated features, integers, steps, or components as referred to, but does not preclude the presence or addition of one or more features, integers, steps, or components, or groups thereof. Moreover, each of the terms “by,” “comprising,” “comprises,” “comprised of,” “including,” “includes,” “included,” “involving,” “involves,” “involved,” and “such as” are used in their open, non-limiting sense and may be used interchangeably. Further, the term “comprising” is intended to include examples and aspects encompassed by the terms “consisting essentially of” and “consisting of.” Similarly, the term “consisting essentially of” is intended to include examples encompassed by the term “consisting of.
As used in the specification and the appended claims, the singular forms “a,” “an” and “the” include plural referents unless the context clearly dictates otherwise. Thus, for example, reference to “a perovskite oxide particle,” “a molten salt,” or “a dopant,” include, but are not limited to, mixtures or combinations of two or more such perovskite oxide particles, molten salts, or dopants, and the like.
It should be noted that ratios, concentrations, amounts, and other numerical data can be expressed herein in a range format. It will be further understood that the endpoints of each of the ranges are significant both in relation to the other endpoint, and independently of the other endpoint. It is also understood that there are a number of values disclosed herein, and that each value is also herein disclosed as “about” that particular value in addition to the value itself. For example, if the value “10” is disclosed, then “about 10” is also disclosed. Ranges can be expressed herein as from “about” one particular value, and/or to “about” another particular value. Similarly, when values are expressed as approximations, by use of the antecedent “about,” it will be understood that the particular value forms a further aspect. For example, if the value “about 10” is disclosed, then “10” is also disclosed.
When a range is expressed, a further aspect includes from the one particular value and/or to the other particular value. For example, where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also included in the disclosure, e.g. the phrase “x to y” includes the range from ‘x’ to ‘y’ as well as the range greater than ‘x’ and less than ‘y.’ The range can also be expressed as an upper limit, e.g. ‘about x, y, z, or less’ and should be interpreted to include the specific ranges of ‘about x’, ‘about y’, and ‘about z’ as well as the ranges of ‘less than x’, less than y’, and ‘less than z’. Likewise, the phrase ‘about x, y, z, or greater’ should be interpreted to include the specific ranges of ‘about x,’ ‘about y,’ and ‘about z’ as well as the ranges of ‘greater than x,’ greater than y,’ and ‘greater than z.’ In addition, the phrase “about ‘x’ to ‘y’”, where ‘x’ and ‘y’ are numerical values, includes “about ‘x’ to about ‘y’”.
It is to be understood that such a range format is used for convenience and brevity, and thus, should be interpreted in a flexible manner to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly recited. To illustrate, a numerical range of “about 0.1% to 5%” should be interpreted to include not only the explicitly recited values of about 0.1% to about 5%, but also include individual values (e.g., about 1%, about 2%, about 3%, and about 4%) and the sub-ranges (e.g., about 0.5% to about 1.1%; about 5% to about 2.4%; about 0.5% to about 3.2%, and about 0.5% to about 4.4%, and other possible sub-ranges) within the indicated range.
As used herein, the terms “about,” “approximate,” “at or about,” and “substantially” mean that the amount or value in question can be the exact value or a value that provides equivalent results or effects as recited in the claims or taught herein. That is, it is understood that amounts, sizes, formulations, parameters, and other quantities and characteristics are not and need not be exact, but may be approximate and/or larger or smaller, as desired, reflecting tolerances, conversion factors, rounding off, measurement error and the like, and other factors known to those of skill in the art such that equivalent results or effects are obtained. In some circumstances, the value that provides equivalent results or effects cannot be reasonably determined. In such cases, it is generally understood, as used herein, that “about” and “at or about” mean the nominal value indicated ±10% variation unless otherwise indicated or inferred. In general, an amount, size, formulation, parameter or other quantity or characteristic is “about,” “approximate,” or “at or about” whether or not expressly stated to be such. It is understood that where “about,” “approximate,” or “at or about” is used before a quantitative value, the parameter also includes the specific quantitative value itself, unless specifically stated otherwise.
As used herein, the term “effective amount” refers to an amount that is sufficient to achieve the desired modification of a physical property of the composition or material. For example, an “effective amount” of an OMS composite refers to an amount that is sufficient to achieve the desired improvement in the property modulated by the formulation component, e.g. achieving the desired level of energy storage for a given time within a given temperature range. The specific level in terms of wt % in a composition required as an effective amount will depend upon a variety of factors including the arrangement of the various parts of the composite material, amount and type of shell around the OMS core if a shell is present, operation temperature of a TCES device using the OMS composite, and identities of perovskite oxide including any dopants as well as the particular molten salt used in the OMS.
As used herein, the terms “optional” or “optionally” means that the subsequently described event or circumstance can or cannot occur, and that the description includes instances where said event or circumstance occurs and instances where it does not.
Unless otherwise specified, temperatures referred to herein are based on atmospheric pressure (i.e. one atmosphere).
Now having described the aspects of the present disclosure, in general, the following Examples describe some additional aspects of the present disclosure. While aspects of the present disclosure are described in connection with the following examples and the corresponding text and figures, there is no intent to limit aspects of the present disclosure to this description. On the contrary, the intent is to cover all alternatives, modifications, and equivalents included within the spirit and scope of the present disclosure.
The following examples are put forth so as to provide those of ordinary skill in the art with a complete disclosure and description of how the compounds, compositions, articles, devices and/or methods claimed herein are made and evaluated, and are intended to be purely exemplary of the disclosure and are not intended to limit the scope of what the inventors regard as their disclosure. Efforts have been made to ensure accuracy with respect to numbers (e.g., amounts, temperature, etc.), but some errors and deviations should be accounted for. Unless indicated otherwise, parts are parts by weight, temperature is in ° C. or is at ambient temperature, and pressure is at or near atmospheric.
The structural and compositional flexibility of perovskite oxides, and their complex yet tunable redox properties offer unique optimization opportunities for thermochemical energy storage (TCES). To improve upon the relatively inefficient, empirical-based approaches, the present study reports a high-throughput combinatorial approach for accelerated development and optimization of perovskite oxides for TCES.
Over the past decade, TCES, which stores thermal heat via a reversible chemical reaction, has become an important research direction owing to its promise in wide operating temperature ranges and high energy storage density. Moreover, since the energy is stored in a chemical form, TCES is capable of long-term, long-distance, or even seasonal energy storage, which is difficult for STES and LTES. Various materials have been proposed for TCES. Among them, redox-active metal oxides have shown the greatest promise due to their ability to operate at high temperatures and without gas storage. Recent modeling studies indicate that up to 55% round trip efficiency can be realized with air as the only reactant. In terms of TCES materials, perovskite oxides have been frequently investigated owing to their unique compositional/structural flexibility, high activity, and cyclic stability.
Many perovskites have been reported for high-temperature TCES. They generally belong to the doped CaMnO3-family, the La-family, or the Sr-family. Doped CaMnO3-based perovskites are the most investigated high-temperature TCES materials. Although CaMnO3-δ itself is unstable under an oxygen-lean environment at high temperatures (e.g., >1000° C.), doping of A-site and B-site cations have shown to improve both the stability and thermodynamic properties. Examples include CaTi0.2Mn0.8O3-δ, CaAl0.2Mn0.8O3-δ, CaFe0.1Mn0.9O3-δ, CaFe0.3Mn0.9O3-δ, CaCr0.05Mn0.05O3-δ, Ca0.95Sr0.05MnO3-δ, Ca0.9Sr0.1MnO3-δ, CaCo0.05Mn0.95O3-δ, CaCu0.2Mn0.8O3-δ, and CaMg0.1Mn0.9O3-δ. The maximum energy density reported to date was 571 kJ/kg ABO3 from CaMn0.95Co0.05O3-δ under a redox cycle between 500° C./0.2 atm O2 and 1000° C./10−5 atm O2. La-based perovskites exhibit high reaction enthalpies but require high operating temperatures (>1000° C.) without doping. Cobalt-doping of LaxSr1-xFeO3-δ and LaxSr1-xMnO3-δ can effectively lower the operating temperatures and enhance the activity. La0.3Sr0.7Co0.9Mn0.1O3-δ demonstrated the largest chemical energy density of 250 kJ/kg ABO3 among the La-family perovskites under a redox cycle between 200° C./0.9 atm O2 and 1000° C./0.001 atm O2. SrFeO3-based perovskites present another promising option for “low-temperature” (400-650° C.) operations. The energy density of undoped SrFeO3 was 81.7 kJ/kg ABO3 under a redox cycle between 400° C./0.2 atm O2 and 1000° C./Ar. The effects of various dopants on the energy storage performance of SrFeO3 have also been investigated. For instance, Mn doping can significantly improve the redox enthalpy and cyclic stability of SrFeO3, Cu substitution contributes to increasing the redox capacity of SrFeO3.
Oxygen capacity, reaction enthalpy, and operating temperature range are the three most critical criteria for the selection of TCES materials. The high flexibility of the perovskite structure and the complex yet tunable properties of perovskite oxides offer unique properties for TCES optimization. However, previous research has mainly relied on an empirical-based approach. There is an urgent need for efficient screening and optimization of perovskite oxides for TCES applications.
The first principle density function theory (DFT) has shown promise in predicting the redox properties of perovskite candidates for TCES. Previous studies showed that the oxygen capacity is closely related to the oxygen vacancy formation energy (Ev). Computationally predicting the reaction enthalpy of perovskite oxides has also been attempted. Most of these studies started from a defect-free structure and directly calculated the reduction enthalpy from perfect perovskite phases to brownmillerite phases. However, oxygen defects often exist in perovskite oxides, even in their oxidized form. High-temperature TCES operations may also exceed the order-disorder transition temperature, and the brownmillerite phase may not form. Therefore, oxygen vacancy concentration and structure as a function of oxygen vacancy levels need to be carefully considered. Additionally, although a few high-throughput computational studies investigated perovskites in the context of chemical looping, no previous research has systematically computed and optimized the perovskite materials for TCES. Relevant experimental data are also limited.
Using high-throughput DFT calculations, this study demonstrates an accelerated approach to developing perovskite oxides for TCES. Over 2000 A/B-site doped SrFeO3 were computed to identify promising TECS materials with satisfactory oxygen capacity and reaction enthalpy. Experimentally, 61 materials were synthesized, 45 samples with pure perovskite phases were tested for oxygen capacities, and 20 promising candidates were chosen for measuring the standard oxidation enthalpy. A chemical energy storage density of up to 157 kJ/kg ABO3 (˜766 kJ/kg ABO3 in total) was achieved by a moderate pressure swing (400° C./0.2 atm O2 and 1100° C./0.01 atm O2). These experimental and simulation data were used to develop effective optimization criteria, which are shown to be highly satisfactory in predicting the standard oxidation enthalpy (within 25.4% deviation) and oxygen capacity (a correlation coefficient of 0.64). Besides accelerating the optimization of TCES materials, the methods developed in this study can be applied in related fields such as chemical looping and catalysis.
Standard reaction enthalpy and redox oxygen capacity are critical for TCES systems since their product determines the TCES density. Redox oxygen capacity is closely related to the operating temperature. As such, multiple criteria are required for choosing a suitable candidate for TCES. At a specific temperature, ΔG needs to be within a suitable range to facilitate oxygen release and uptake. Since a small ΔG and a large ΔH are desirable for TCES, they were both used as the descriptors when analyzing the results from the high-throughput DFT calculations of the 2,003 SrxA1-xFeyB1-yO3-δ candidates.
The selection of TCES materials is also dependent upon the operating temperature range. A 400-800° C. temperature range, which is readily compatible with thermal power plants and many industrial processes, was chosen in this study. At the upper (800° C.) and the lower (400° C.) limits, the suitable range of ΔG for triggering the oxygen release can be directly calculated from Eq. (3).
Most studies on TCES materials expose the samples to very low oxygen concentrations (e.g., ≤10−5 atm) during the energy storage step. Although this would increase the measured energy storage capacity, it is unlikely to be feasible from a practical standpoint, especially for using packed beds. Therefore, this study adopts a pressure swing between PO2=0.01 (energy storage) and PO2=0.2 atm (energy release). This corresponds to ΔG ranges of 0.047-0.134 eV at 400° C., and 0.074-0.213 eV at 800° C. Since a previous investigation indicated a DFT simulation error of −0.3 eV to 0.5 eV, potentially suitable materials should meet the following requirements:
Past experience indicates that oxygen non-stoichiometry (δ) in TCES materials is usually less than 0.3 under a moderate oxygen partial pressure of >0.001 atm. Therefore, the initial screening focused on a δ range of 0-0.25, i.e., the average between ΔG0-0.125 and ΔG0.125-0.25 calculated from DFT was compared to the aforementioned criteria.
775 samples are potentially suitable based on this criterion. Most samples containing Ti, Sm, La, and Cu were removed from consideration due to undesirable ΔG. Ti, Sm, and La dopants tend to increase the temperatures for oxygen release, leading to an excessively large ΔG. In contrast, Cu dopant has the opposite impact on ΔG.
With respect to the other descriptor, i.e., ΔH,
Besides ΔG and ΔH, the cost of the material is also important. Therefore, Ni, Co, Cu, Sm, and Y were excluded from experimental investigations. K-doped materials were also excluded due to the potential corrosion issues. 209 candidates were sorted by ΔH (δ=0-0.25) normalized by molecular weight, 61 materials were synthesized, and 45 samples with negligible phase impurities were tested for their redox oxygen capacity. A complete list of tested samples can be seen in Table 3 below.
The thermochemical energy storage of the oxides can also be triggered by CO2 pressure swing. A CO2-rich gas stream, such as the flue gas from power plants or other industrial processes, is flowed into the reactor during the heat discharging period. The mixed metal oxides will react with the CO2 to form the carbonate and release a significant amount of heat at the same time. The O2 release will offset a certain amount of heat release, but the overall reaction is still highly exothermic. Taking Sr0.8Ba0.2Fe0.4Co0.6O3-δ as an example, the carbonation reaction of perovskite with CO2 is described as follows:
During the heat charging mode, high-temperature steam and O2 flow into the reactor. The carbonate will decompose and absorb a significant amount of heat. At the same time, the perovskite will be regenerated so that the sintering can be avoided. The decarbonation reaction and the regeneration reaction can be described as follows:
The gas flow downstream contains concentrated CO2, which is desirable for CO2 capture and utilization. Based on this reaction, CO2 can be combined with thermal energy storage.
The oxygen capacities under different conditions are summarized in Table.1. Under the redox condition of 400° C./0.2 atm O2 and 800° C./0.01 atm O2, the highest oxygen capacity was 2.11 wt. % for Sr0.75Ba0.25Fe0.875Mg0.125O3-δ, while Sr0.5La0.5Fe0.125Mn0.85O3-δ show negligible oxygen release. Ba-doped or Ba/Mg-co-doped SrFeO3-δ exhibited considerable oxygen capacities at low temperatures and high oxygen partial pressures. Compared to SrFeO3-δ, Sr0.75Ba0.25Fe0.875Mg0.125O3-δ and Sr0.875Ba0.125FeO3-δ favored oxygen release at higher oxygen pressure. The oxygen absorption/desorption kinetics also increased with the Ba dopant, which would be beneficial to the power density. However, when Ba doping was larger than 0.375 (e.g., Sr0.5Ba0.5FeO3-δ), phase segregation of BaFeO3 and SrFeO3 was observed through XRD. Interestingly, phase segregation at high Ba levels can be inhibited by substituting Fe with small amounts of Mg and Ti in the B-site, e.g., Sr0.25Ba0.75Fe0.875Mg0.125O3-δ and Sr0.125Ba0.875Fe0.875Ti0.125O3-δ. Even though the oxygen capacities were slightly lower, they all demonstrated sufficient oxygen capacity under an isothermal condition at 400° C. even though 0.01 atm O2 was used as the lower bound for oxygen partial pressure, which is three orders of magnitude higher than typical testing conditions. In general, Ba-based perovskites exhibited superior performance in O2 redox capacity within the 400-800° C. range, which is highly compatible with steam-based power plants.
SrxCa1-xFeyMn1-yO3-δ perovskites were another group of materials showing promising TCES performance. They exhibited the highest changes in δ owing to their smaller molar masses. Mn-doping tended to increase the temperature window for oxygen release, a considerable amount of oxygen was released in the range of 800-1100° C. Sr0.125Ca0.875Fe0.25Mn0.75O3-δ, Sr0.25Ca0.75Fe0.375Mn0.625O3-δ, and Sr0.375Ca0.625Fe0.125Mn0.875O3-δ all exhibited oxygen releases >2.8 wt % under the redox condition between 400° C./0.2 atm O2 and 1100° C./Ar. In particular, Sr0.375Ca0.625Fe0.125Mn0.875O3-δ and Sr0.125Ca0.875Fe0.25Mn0.75O3-δ had an oxygen capacity >2.1 wt % even with a redox condition between 400° C./0.2 atm O2 and 1100° C./0.01 atm O2. This moderate oxygen atmosphere is helpful for the flexible operation of industrial TCES applications.
Based on the oxygen capacity under the redox condition between 400° C./0.2 atm O2 and 1100° C./Ar, 20 materials with oxygen releases >1 wt. % were selected to measure their reaction enthalpies.
The effects of Mn-dopants varied with A-site dopant types. For the Sr—Ca perovskites, Mn doping increased the oxygen release temperature but concurrently led to a higher ΔH. For instance, for the Sr0.25Ca0.75FeyMn1-yO3-δ materials and the Sr0.375Ca0.625FeyMn1-yO3-δ materials, ΔH significantly increased with the increase in the substitution of Mn in the B-site. However, a small quantity of Mn-doping in Sr0.75Ba0.25FeO3-δ significantly affected ΔH as well as the oxygen release performance. Among all the tested materials, Sr0.125Ca0.575Fe0.25Mn0.75O3-δ displayed the highest ΔH of 160 kJ/mol 0 with Ab of 0.3 under a redox condition between 400° C./0.2 atm O2 and 1100° C./Ar.
The chemical energy storage densities (ΔHchem) of the screened materials are summarized in
As shown in
Although the initial screening criteria were sufficient for predicting the ΔH, they led to non-negligible deviation in predicting oxygen capacities. Based upon the experimental data, alternative screening criteria were tested to improve the effectiveness in predicting the oxygen capacity of perovskite oxides. Due to the limit to the supercell size, DFT calculations have limited resolutions for δ change. In the current calculation, the minimum δ change is 0.125. Therefore, the experimentally measured oxygen capacity was correlated with the DFT calculated ΔG at different δ variation ranges to determine the most effective descriptor for the oxygen capacity. Correlation coefficients, defined in Eq.(1), were used to describe the relationship between oxygen capacity and ΔG within different δ ranges.
A smaller but non-negative ΔG indicates the thermodynamic favorability of oxygen release from the perovskite. Therefore, oxygen capacity and ΔG should present a negative correlation. Since some DFT calculated ΔG were less than zero and are likely to be inaccurate, the correlation coefficients were also calculated by resetting the negative ΔG to 0. The oxygen capacities under the isothermal redox condition between 0.2 atm O2 and Ar were used for the calculation. As shown in
The current study reports an effective high-throughput combinatorial approach for accelerated development and optimization of perovskite oxides for thermochemical energy storage. Based on Density Functional Theory (DFT) calculation results of more than two thousand A-site and/or B-site substituted SrFeO3-δ perovskites at various oxygen non-stoichiometry levels (δ), 61 promising TCES candidates were selected and prepared. Nearly all the samples formed perovskite phases as their main compositions. Of these, 45 samples showed negligible phase impurities and were thoroughly investigated experimentally, including detailed measurements of redox thermodynamic parameters of 20 samples. Many of these samples showed promising TCES properties. Sr0.875Ba0.125FeO3-δ exhibited a chemical energy storage density of 85 kJ/kg ABO3 under an isobaric air condition between 400° C. and 800° C., and Sr0.125Ca0.875Fe0.25Mn0.75O3-δ achieved an energy density of 157 kJ/kg ABO3 under a temperature and oxygen partial pressure swing between 400° C./0.2 atm O2 and 1100° C./0.01 atm O2. The experimental results supported the effectiveness of the high-throughput predictions for oxygen capacity and standard oxidation enthalpy. An improved set of screening criteria was also developed based on the experimental and DFT results. This improved approach, with an average deviation of 25.4% for predicting the standard oxidation enthalpy and a correlation coefficient of 0.64 for predicting the oxygen capacity, would be highly useful for accelerated development of redox-active perovskite oxides.
First-principles DFT simulations were executed using the Vienna ab initio Simulation package (VASP) software. Here the frozen-core all-electron projector augmented wave (PAW) model and Perdew-Burke-Ernzerhof (PBE) functional were utilized. The energy cutoff was set to 450 eV, and the force and energy converged with the criteria of 0.01 eV Å−1 and 10−5 eV, respectively. Gaussian smearing was employed with a width of 0.1 eV for optimization simulation. 1×1×1 and 1×2×2 k-points were used for the 2×2×2 SrxA1-xFeyB1-yO3-δ perovskite supercells and brownmillerite structures, respectively. Based on previous studies, DFT+U method was utilized for d orbitals of Fe, Co, Cu, Mn, Ni, and Ti with Ueff=4, 3.4, 4, 3.9, 6, and 3, respectively. Only FM phase magnetic ordering was considered for all the doped structures owing to the negligible effect of magnetic ordering on oxygen vacancy formation and migration. The initial spin moment for Fe, Co, Mn, Ni were set to 4, 5, 5, 5, respectively.
A Monte Carlo special quasi-random structures (MCSQS) method was applied to determine the positions of all A- and B-site dopants and oxygen vacancies to approach randomly disordered structures. The zero-point energy (ZPE) was calculated using the Phonopy code. For an optimized SrxA1-xFeyB1-yO3-δ crystal structure, its dynamical matrix of force constants was obtained by using the forces from DFT calculations on a 2×2×2 supercell. The enthalpy of O2 was computed using the CBS-QB3 method in Gaussian 16. A more detailed description of the simulation approach can be found in a previous publication.
Given an increment of 0.125 for both dopant fraction and nonstoichiometry, 2401 perovskite models and 9604 different conditions were constructed, including various A-site cations (Ca, K, Y, Ba, La or Sm) and/or B-site cations (Co, Cu, Mn, Mg, Ni or Ti). Some unstable perovskite structures were excluded by the preliminary screening with two criteria. Firstly, the compositions should be charge neutral, leading to the removal of 168 materials from the screening list. Then, a modified tolerance factor was used to further eliminate high-distortion candidates, leaving 2003 materials for high-throughput calculations of ΔG, ΔH, and ΔS.
All the perovskite oxides were synthesized using a solid-state method. In a typical synthesis of SrxBa1-xFeyMn1-yO3-δ, stoichiometric amounts of SrCO3, BaO, Fe2O3 and MnO2 were weighed and put in a 5 mL PTFE vial. Then, 2 mm ZrO2 beads were added into the vial with a mass ratio of 5:1. 2 mL ethanol (>99 vol % purity) was further added to the mixture to prevent powders from sticking to the vial walls. Four vials with four different materials were housed in a stainless-steel sample jar and then balled milled with 1200 RPM for 3 h or 24 h. Based on the XRD patterns, excessive ball-milling time will cause minor impurities in some materials. The resulting wet mixture in each vial was dried at 95° C. for 0.5 h and 130° C. for 15 min to remove ethanol. The powder mixture was separated from ZrO2 beads and then fired at 1000° C. in a muffle furnace for 10 h to obtain the perovskite structure. The heating and cooling rates were all set to be 3° C./min. Finally, the perovskite samples were sieved to two desired particle size ranges, i.e., 180-250 μm for oxygen release experiments and 0-180 μm for XRD and reaction enthalpy experiments. 61 different materials were synthesized in this study. Precursors used were SrCO3 (>99.9%), CaCO3 (>99.9%), BaO (>98%), La2O3 (>99%), Fe2O3 (>99%), MnO2 (>99%), TiO2 (>99%) and MgO (>99%).
Crystal structures of the samples were determined on an Empyrean PANalytical XRD with Cu-Kα radiation (λ=1.5406 Å) operating at 45 kV and 40 mA. The scan was conducted from 2θ of 15° to 80° with a step size of 0.0262° and a hold time of 0.2 s for each step. The XRD phases were identified using Highscore Plus software. 45 samples were pure perovskite phases or contained negligibly small phase impurities, but the rest of the 16 samples contained notable impurities.
Two types of experiments were carried out in a thermogravimetric analyzer (TGA, TA-Instrument, SDT Q650). The oxygen capacity experiments were used as the first-step screening for the above-mentioned 40 samples. In a typical experiment, 40-50 mg samples with a particle size of 180-250 μm were loaded into an Al2O3 crucible with a 6.5 mm I.D. and then placed in the TGA. The total flow rate was maintained at 200 mL/min, and three oxygen concentrations were tested, i.e., 0.2 atm, 0.01 atm, and Ar (Airgas UHP 5.0 grade). The oxygen concentration was varied by mixing pure oxygen (Airgas extra dry grade O2) with Ar. The oxygen partial pressure in the Ar flow, which was approximately 5×105 atm, was measured by an oxygen analyzer (Setnag). The sample was increased from room temperature to 180° C. under a 0.2 atm O2 atmosphere and then was held for 10 min to eliminate moisture. After that, the temperature was ramped to 700° C. and followed by a 10 min isothermal step as a pretreatment. Then, the temperature was decreased to 400° C. and held for 15 min to obtain the first particle weight mi. Another two points were measured at 800° C. and 1100° C., respectively. Afterward, the temperature was ramped back to 400° C., and the oxygen concentration was changed to 1%. The next cycle of the temperature test was started using the same program. After measuring the particle weights at different temperatures and oxygen concentrations, oxygen capacities were obtained under different conditions. The oxygen release kinetics of the Ba-doped or Ba/Mg co-doped SrFeO3 were also tested.
Twenty samples with an oxygen capacity >1 wt. % were then tested to obtain the oxygen nonstoichiometry and reduction enthalpy. In a typical experiment, ˜30 mg samples with a particle size <180 μm were tested in the TGA. The total flow rate was also maintained at 200 ml/min, and seven different oxygen concentrations were tested, i.e., 0.8, 0.2, 0.05, 0.01, 0.003, 0.0005 atm O2, and Ar. The 0.003 atm and 0.0005 atm oxygen concentrations were realized by mixing 1.048 vol % O2 calibration gas (balance Ar) with pure Ar. The temperature was first increased to 1100° C. under the 80% O2 atmosphere and then held for a specific time until the particle weight change reached equilibrium. The temperature was then sequentially decreased in a stepwise mode by 400° C. with an increment of 100° C. Afterward, the temperature was ramped back to 1100° C., and then the gas atmosphere was switched to 20% O2. The next cycle of the temperature test was started under a different oxygen concentration using the same program.
A Ca2AlMnO5/ binary chloride (CaCl2/NaCl) composite material was prepared by the following procedures. Ca2AlMnO5 was first synthesized via a modified Pechini method and then was balled mill with CaCl2 and ZrO2 beads. The resulting powder mixtures were rinsed with ethanol to remove CaCl2. Afterward, the remaining wet particles were dried in the air to obtain the high-porosity Ca2AlMnO5. CaCl2 and NaCl powders were mixed with a molar ratio of 52:48 and then calcined above the melting point (˜500° C.) to form the binary chloride molten salts. High-porosity Ca2AlMnO5 and molten salt powders were mixed with a specific mass ratio and then calcined at 550° C. to form the composite material.
The BET results revealed that the BJH pore volume increased from 0.0002 to 0.14 cm3/g after the ball milling step. As shown in
Other perovskite/metal oxides and eutectic salt mixtures were prepared by mixing in a mortar. A small portion of the particles was separated for TGA experiments, and the rest of the powders were pelletized in a 6 mm I.D. pellet dies at 10 mPa for 3 min with a hydraulic press. This was followed by heating up to 50° C. above the melting points of the salts in a muffle furnace. This step was to ensure the uniform mixing of the salt and perovskite within the pellets by melting the salt. The pellets were heated up by first the ramping rate 3° C./min and then 1° C./min during the last 150° C. to reach a temperature 50° C. above the melting point of the salt mixture, which was kept for 1 hour. This was followed by cooling it down to room temperature. After this step, the pellets were crashed into fine particles for XRD analysis.
The compatibility results based on the XRD patters were summarized in Table 5. 13 composite materials were found to be compatible among the 26 samples screened. The XRD patterns of the compatible samples were presented in
The performance of the composite material was tested in a thermogravimetric analyzer (TGA) under a cyclic operation between 400˜600° C. The latent heat of the molten salt was obtained via a differential scanning calorimeter (DSC) method. The heat flow was first calibrated using a sapphire with a determined heat capacity. As shown in
The exhaust gas emitted by furnaces in chemical industries carries substantial waste heat in the form of a reducing gas composition. This composition, often rich in reducing agents like H2 and CO, holds valuable potential for utilization in thermochemical energy storage applications. By harnessing this waste heat and incorporating it into thermochemical energy storage systems, industries have the opportunity to boost their energy efficiency. Therefore, an experiment under a reducing gas atmosphere was performed to achieve high thermochemical energy through the reduction of perovskite in the composite materials. Among the tested materials with La-based perovskite oxides, La0.6Sr0.4FeO3-δ was the only one that showed a redox capacity in the regarding operating conditions (575-725° C.). Therefore, La0.6Sr0.4FeO3-δ:Li2MoO4 was selected with the weight ratio of 0.6:0.4 as a potential candidate for this experiment. After the sample was completely reduced under 20% H2, it was exposed to 20% O2 for reoxidation. The thermochemical energy storage from the peak during oxidation was determined as the oxidation occurs suddenly, so the amount of released heat can be found from the DSC method. The plots related to these measurements are presented in
Long term stability test was performed for the sample of La0.5Sr0.2FeO3-δ:LiF—NaF—CaF2 as it was one of the most promising ones. The sample was operated during 100 redox cycles under the conditions of 510/20% O2-660° C./Ar. The results of 100 cycles was summarized in
In the long-run test, porous perovskite oxide particles were used in the composite materials. The porosity of the perovskites was increased by the above-mentioned salt grinding method. The porous structure of perovskite oxide can be seen from the SEM image (
The compatibility between Ca2AlMnO5 and binary chloride molten salt was tested using Nano-computed tomography (Nano-CT). The composite materials were first packed with 10 tons of pressure into a pellet of 6 mm O.D. and 10 mm height. Then, the pellet was calcined at 550° C. in a tube furnace for one hour. As shown in
DFT was used to conduct high-throughput screening of A-site and/or B-site substituted SrFeO3 based perovskites with 2401 different compositions. Oxygen capacity, reaction enthalpy as well as material price were considered to pre-screen the suitable candidates for TCES. Then, some perovskite candidates were synthesized using the solid-state method. The oxygen releases under different redox conditions were measured in a TGA. The non-stoichiometry of the perovskite oxide was further obtained by varying the oxygen concentrations and temperatures. The corresponding reaction enthalpy was calculated via Van′t Hoff equation.
As aforementioned, standard reaction enthalpy and redox oxygen capacity are critical for TCES systems since their product determines the overall TCES density. Redox oxygen capacity is closely related to the operating temperature. As such, multiple criteria are required for choosing a suitable candidate for TCES. Gibbs free energy change ΔG has been demonstrated as an adequate descriptor for redox oxygen capacity based on previous publications. At a specific temperature ΔG needs to be within a suitable range to facilitate oxygen release and uptake, while a higher standard reaction enthalpy ΔH for a similar oxygen capacity can lead to a higher energy storage density. Thus, both ΔG and ΔH were used as the descriptors for the high-throughput DFT calculations on the 2,003 SrxA1-xFeyB1-yOδ candidates.
The simulation results also revealed a less significant correlation between ΔH and temperature, which has been reported in many previous studies. The ΔH heatmaps of the screened perovskite candidates are shown in
In terms of the operating conditions, various materials can be developed depending on the temperatures of heat sources and oxygen partial in the working fluid. 400-800° C., which is widely compatible with thermal power plants and many other industrial processes, was chosen in this study for material design and experimental validation. At the upper (800° C.) and the lower (400° C.) limits, the suitable range of ΔG for triggering the oxygen release can be directly calculated from the pressure swing using the following equation:
Even though most of the previous studies on material design treated the samples in the atmosphere with considerably low oxygen concentrations, an oxygen concentration closer to air is more feasible and energy-efficient for TCES. A pressure swing between PO2=0.01 and PO2=0.2 atm was utilized in this study when considering both the heat charging and discharging processes. This corresponds to ΔG ranges of 0.047-0.134 eV at 400° C., and 0.074-0.213 eV at 800° C. Based on a previous investigation, the DFT simulation error is −0.3 eV to 0.5 eV. Therefore, the suitable materials should meet the requirements for ΔGDFT as set out previously.
Previous experiments indicate that δ dramatically varies among different TCES materials but is usually less than 0.3 under a moderate oxygen partial pressure >0.001 atm. As such, the initial focus was on the δ range between 0-0.25 and used the average ΔG results of ΔG0-0.125 and ΔG0.125-0.25 for the pre-screening step.
As shown in
The redox oxygen capacities of selected materials under different conditions is listed in Table.4. 45% of the predicted materials released >0.15 wt % O2 with the pressure swing between 0.2 atm O2 to 0.01 atm O2 at 400° C., and that number increased to 75% at 800° C. Considering the simulation errors, the computed ΔG can be used a reasonable reference for pre-screening materials with good oxygen capacity.
In the redox condition of 400° C./0.2 atm O2 and 800° C./Ar, the highest oxygen capacity was 2.19 wt. % for Sr0.625Ba0.375FeO3-δ, while the lowest one was 0.18 wt. % for CaFeO3-δ. Sr0.75Ba0.25FeO3-δ has the second largest oxygen release capacity by 2.04 wt. %, which might indicate that higher Ba doping in the A-site led to a higher oxygen capacity at lower temperatures when B-site only contained Fe. This is consistent with Bush et al.'s study on Sr1-xBaxFeO3-δ for air separation. However, when Ba doping was larger than 0.375 (e.g., Sr0.5Ba0.5FeO3), phase segregation of BaFeO3 and SrFeO3 in the XRD patterns was observed, which might be due to the crystal distortion caused by the large atom of Ba. Pure phases could still be restored for high Ba dopants by substituting Fe with small amounts of Mg and Ti in the B-site, such as Sr0.25Ba0.75Fe0.875Mg0.125O3-δ and Sr0.125Ba0.875Fe0.875Ti0.125 O3-δ. Even though the oxygen capacities were slightly lower, they all demonstrated sufficient oxygen capacity under an isothermal condition at 400° C., even with a redox condition of 0.2 atm/0.01 atm O2. The general trend of Ba-based perovskites exhibited superior performance in O2 redox capacity within the medium-high temperature range, which was promising for TCES applications that require medium-high temperatures, e.g., the recovery of industrial heat. Additionally, when the temperature was ramped from 400° C. to 800° C. at 0.2 atm O2, oxygen release of most of the Ba-containing perovskites (from Sr0.625Ba0.375FeO3-δ to Sr0.75Ba0.25Fe0.875Mn0.125O3-δ) were all larger than 1 wt %. In other words, the energy storage of Ba-containing perovskites can be triggered with only a temperature swing in the air.
SrxCa1-xFeyMn1-yO3 perovskites were another group of materials showing promising performance for TCES applications, which had the highest changes in δ owing to their smaller molar masses. As Mn-doping tended to increase the temperature window for oxygen release, a considerable amount of oxygen was released in the range of 800-1100° C. Sr0.375Ca0.625Fe0.125Mn0.875O3, Sr0.25Ca0.75Fe0.75Mn0.25O3, and Sr0.375Ca0.625Fe0.125Mn0.875O3 all exhibited oxygen releases>2.5 wt. % under the redox condition between 400° C./0.2 atm O2 and 1100° C./Ar. In particular, Sr0.375Ca0.625Fe0.125Mn0.875O3 had an oxygen capacity >2 wt. % even with a redox condition between 400° C./0.2 atm O2 and 1100° C./0.01 atm O2. This moderate oxygen atmosphere would be a potential benefit for industrial TCES applications with flexible operation.
Among all the tested materials, A-site dopant types and concentrations appeared to have more significant effects on the redox capacity than B-sites. For instance, for perovskite co-doped with Sr—Ca in the A-site, higher Mn-doped materials released more substantial amounts of oxygen than Mg and Ti-doped ones. However, Mn-doping resulted in lower oxygen capacities in those materials co-doped with Sr—Ba. Additionally, increasing the doping amount of Mg and Ti adversely affected phase purity or redox performance, which was another factor that required consideration.
Based on the oxygen capacity under the redox condition between 400° C./0.2 atm O2 and 800° C./Ar, 15 good materials with high oxygen release (>1 wt. %) and 2 unfavorable materials (<1 wt. %) with low oxygen release were selected to measure their reaction enthalpies. ΔH of all tested materials varied from 20 to 150 kJ/mol O, while δ almost covered the whole range of 0 to 0.5. In general, ΔH increased with the oxygen vacancy, which is consistent with most reported data of SrFeO3-based perovskite. It is clear that the dopant types and concentrations significantly changed ΔH, but some general trends can be captured from the tested materials. Firstly, a trade-off between the oxygen release and the reaction enthalpy can be easily seen. As mentioned above, Ba doping tended to reduce the temperature window for oxygen release, but the corresponding ΔH was also lower compared to La-doped, Ca-doped, or un-doped samples. Mg-doping stabilized the SrxBa1-xFeO3-δ at high-Ba doping conditions and facilitated the oxygen release but negatively affected the ΔH. Additionally, δ of both Sr0.625Ba0.375FeO3-6 and Sr0.75Ba0.25FeO3-5 reached 0.45 at 800° C./0.01 atm O2 condition, so their energy storage density would be limited by the insufficient oxidation. In other words, they might be suitable for TCES at an even lower temperature (<400° C.), which would be benefited from further oxidation and a smaller 5.
In contrast, Mn doping enhanced the temperature to release oxygen but concurrently led to a higher ΔH. For instance, for the three Sr0.25Ca0.75FeyMn1-yO3-δ materials, ΔH significantly increased by about 100 kJ/mol 0 when the substitution of Mn in the B-site increased from 37.5% to 75%. Moreover, a small quantity of Mn-doping in Sr0.75Ba0.25FeO3-δ significantly improved the oxygen release performance but barely affected ΔH, leading to an elevated energy density at the same condition. Among all the tested materials, Sr0.375Ca0.625Fe0.125Mn0.875O3-δ displayed the highest ΔH of 160 kJ/mol 0 with Ab of 0.3 under a redox condition between 400° C./0.2 atm O2 and 1000° C./Ar.
The chemical energy storage densities (ΔHchem) of screened materials are summarized in Table.3. The redox reaction of Ba-doped samples can be easily triggered by the air. A higher temperature or a lower oxygen partial pressure only leads to a minor increase in ΔHchem, especially for Sr0.25Ba0.75Fe0.575Mg0.125O3-δ. Sr0.625Ba0.375FeO3-δ exhibited the largest ΔHchem of 67 kJ/kgABO3 under a isobaric air condition between 400° C. and 800° C. In contrast, SrxCa1-xFeyMg1-yO3-δ barely released O2 and absorbed heat at low temperatures, but ΔHchem was significantly improved when the temperature was above 800° C. Sr0.375Ca0.625Fe0.125Mn0.875O3-δ achieved a ΔHchem of 162 kJ/kgABO3 under a redox condition between 400° C./0.2 atm O2 and 800° C./0.01 atm O2. The therotical energy densities of composite materials were also calculated for each individual perovskite. The composite material was assumed to contain 50 wt % of Li/Na binary eutatic carbonate salt (44 wt. % Li2CO3-56% K2CO3) and 50 wt % of perovskite oxide. The redox condition was assumed to be under 400° C./0.2 atm O2 and 800° C./Ar or 400° C./0.2 atm O2 and 1100° C./Ar. For the carbonate salts, the melting temperature is about 496° C., the heat of fusion is about 380 kJ/kg, and the specific heat capacity is about 2.36 kJ/kg-K. The specific heat capacity of the perovskite was assumed to be 0.8 kJ/kg-K. The composite materials under specific conditions were all considerably high.
Chemical energy storage densities of screened materials are provided in Table 7:
It should be emphasized that the above-described embodiments of the present disclosure are merely possible examples of implementations set forth for a clear understanding of the principles of the disclosure. Many variations and modifications may be made to the above-described embodiment(s) without departing substantially from the spirit and principles of the disclosure. All such modifications and variations are intended to be included herein within the scope of this disclosure and protected by the following claims.
This application claims the benefit of U.S. Provisional Application No. 63/489,242 filed on Mar. 9, 2023, which is incorporated herein by reference in its entirety.
This invention was made with government support under grant number CBET1923468 awarded by the National Science Foundation. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63489242 | Mar 2023 | US |