{'fraction_non_alphanumeric': 0.030160003775900317, 'fraction_numerical': 0.012366073535658658, 'mean_word_length': 4.511119781506047, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A Water Balloon as an Innovative Energy Storage Medium

Soft rubbery materials are capable of withstanding large deformation, and stretched rubber contracts when heated. Additionally, rubber balloons exhibit non-monotonic pressure–volume curves. These unique properties have inspired numerous ingenious inventions based on rubber balloons. To the authors’ knowledge, however, it is surprising that these properties have not inspired any study that exploits the elasticity of rubber balloons for energy storage. Motivated by these, this study examines the performance of water balloons as energy storage media. In each experiment, a single water balloon is implemented using a flat membrane, and it is subject to repeated inflation, heating, deflation, and cooling. Inflating the balloon deposits energy into it. The heating simulates the recycling of waste heat. The balloon delivers work during its deflation. Finally, the cooling completes the energy-storage cycle. The performance is evaluated in terms of the balloon’s transferred energies, efficiencies, and service life. Simple as it is, a water balloon is actually an impressively efficient energy storage medium. The efficiency is 85–90% when a water balloon stores and releases energy at room temperature. Recycling waste heat can boost a balloon’s efficiency beyond 100%, provided that the cost of the heat is negligible so that the heat is not taken as part of the input energy. However, heating shortens the service life of a balloon and reduces the total energy it can accommodate. By running fatigue tests on balloons, this study reveals the trade-off between a water balloon’s efficiency and its longevity. These results shall serve as a useful guide for implementing balloon-based mechanical devices not limited to energy-storage applications.


Introduction
Soft rubber membranes possess many attractive properties for implementing mechanical applications. The capability of withstanding large deformation is perhaps the most unique one. This property enables a rubber membrane to form a balloon and enclose a fluid with several hundred times its volume. In fact, a balloon formed by a flat membrane is known to possess a non-monotonic relation between its pressure P and its volume V [1][2][3][4]. The non-monotonicity of the pressure-volume (P-V) curve incurs snap-through and snapback instabilities [5][6][7]. These properties enable a rubber membrane to exert forces with a large displacement volume and/or a fast response. In published works, the properties have been exploited to implement pumps [8][9][10][11] for flow control, catheters [12][13][14][15][16] for angioplasty [17], artificial muscles for robots [18][19][20], harvesters for renewable energy [21][22][23], and elastocaloric heat pumps for thermal energy management [24]. When a pressurized fluid inflates a rubber membrane into a balloon, the membrane stores the work done by the fluid as its elastic potential energy. The energy can be withdrawn for utility by deflating the balloon and releasing the pressurized fluid. The inflation and deflation of a balloon are analogous to the charge and discharge of a battery. Motivated by these, this study develops water balloon energy storage (WBES), in which a water balloon functions as an energy storage medium to accommodate the energy transported by subcooled liquid water. A first scientific question here is the long-term evolution of the energies associated with the repeated inflation and deflation of a water balloon. Additionally, a second scientific effect for energy storage incurs a trade-off between the balloon's longevity and efficiency of energy storage.
This study tests the performance of a single water balloon as an energy storage medium. In each experiment, a balloon is subject to cyclic mechanical and thermal loads. The performance shall be quantified in terms of the balloon's transferred energies, efficiencies, and service life. Instead of using compressed air, this study adopts subcooled liquid water as the working fluid for its higher density and its incompressibility. For fluid streams with the same volume flow rate, the one with a higher density can deliver a larger inertia force for work production. On the other hand, the incompressibility enables precise and repeatable prescription of the balloon's volume that is immune to the disturbance associated with any temperature fluctuation. A single balloon accommodates only a limited amount of energy. For practical utility, the capacity can be scaled up by implementing a large number of balloons in arrays. Working with multiple interconnected balloons requires suitable control of the scavenging between the balloons [59]. For simplicity, this study focuses on the behavior of one single balloon and saves the multi-balloon problem for a future study. In what follows, Section 2 reports the operation and performance evaluation of WBES, and Section 3 presents all results and analyses, followed by the conclusion in Section 4.

Materials and Methods
This section reports the operation of WBES and the methods for testing its performance. In what follows, Section 2.1 proposes the energy-storage cycle of WBES, its synchronization with daily power demand, and the candidate power sources for its operation. Subsequently, Section 2.2 presents the experimental setup for testing the performance of WBES, and Section 2.3 defines the parameters for performance evaluation.

Energy-Storage Cycle
Much like other energy storage media, a water balloon operates in a cycle to store and supply energy. Figure 1 illustrates the processes of the cycle, in which a balloon is inflated with cold water, heated, deflated when it is warm, and cooled when it is flat. The inflation process deposits energy into the balloon, and the balloon's membrane transforms the work done by the water into its elastic potential energy. The heating process simulates the recycling of waste heat by heating the balloon to boost its pressure. This can make up for a substantial proportion of the pressure drop due to the membrane's relaxation. Subsequently, the balloon releases pressurized water to propel a water turbine (not shown) and produce work during the deflation process. Finally, the cooling process completes the cycle by bringing the balloon's membrane back to the initial (i.e., pre-inflation) temperature. The processes in Figure 1 can be readily synchronized with the daily demand for power in cities. The demand is typically higher during the day and lower at night. As such, the balloon should be inflated at night and heated with solar heat or other waste heat during the day. Then the balloon can release its pressurized water for power generation when energy is most needed. Based on this, a water balloon repeats the cycle in Figure 1 once every day.
WBES accommodates the energy of the pressurized water, and the water may be pumped from different sources. One intuitive option is to pump the water mechanically with the electricity during the off-peak hours or from a renewable source. Alternatively, a WBES system can be integrated with a rainwater-harvesting (RWH) system [60][61][62]. This can be implemented, for example, by installing the WBES system in the basement of a building and connecting it either directly to the gutter through the downpipes or a header tank right below the roof [62]. In that case, WBES stores the gravitational energy that the building retains in the rain. This integrates WBES and the rainwater harvester into an eco-friendly energy harvesting system.

Experimental Setup
For each experiment in this study, a balloon is prepared by clamping a membrane (model 8166NS, by Cranberry International Sdn Bhd, Shah Alam, Malaysia) around the end of a circular tube, as shown in Figure 2. The material of the membrane is natural rubber vulcanized with zinc oxide, sulfur, ZDEC, ZDBC, antioxidant, 20% Aerodisp W7520, peppermint oil, and spearmint oil according to the manufacturer. The membrane itself is originally flat and has a uniform thickness of h = 216 ± 4 µm. It is clamped around the tube with a worm-drive clamp, and a rubber seal is sandwiched between the clamp and the membrane to prevent leakage. The radius of the circular tube is 8.3 mm, which is the balloon's base radius r b . Once assembled, a balloon is degassed in a water tank and kept there to be connected with other mechanical parts for the experiment. (Color online) A schematic illustration of an inflated balloon with pressure P, volume V, and base radius r b in this study. The balloon is prepared by clamping a flat piece of rubber membrane around the end of an aluminum tube and injecting de-ionized water from inside the tube.
Once prepared, a balloon is tested with the setup in Figure 3. The setup consists of a volume control unit, a temperature control unit, and a signal processing unit. The volume control unit applies the mechanical load to the balloon and executes the inflation and the deflation processes in Figure 1. The components are a cylinder module (MBR40-200, Chanto Air Hydraulics Co., Ltd., Taichung City, Taiwan), a linear motor (model CLMS-1010, by Chieftek Precision, Tainan City, Taiwan), a motor controller (model Will1-B-3, by Chieftek Precision, Tainan City, Taiwan), and a pressure transducer (PT, model PT-VP-2P-G4-H-S1-C-24V, by Atlantis Inc., Taipei City, Taiwan). Prior to each experiment, the cylinder module is degassed in a water tank and connected to a degassed balloon in the tank. Then the piston rod of the cylinder module is connected to the linear motor on a workbench. Subsequently, the inflation process is executed by having the cylinder module extrude water into the balloon. This stores energy in the balloon. During the deflation process, conversely, the cylinder module withdraws the water from the balloon. This simulates the situation when the balloon releases its pressurized water for work production. In both processes, the linear motor actuates the cylinder to control its displacement and speed. The controller drives the motor and keeps track of its displacement, and the motion is prescribed through an operating software (CPC GUI version 0.5.20, by Chieftek Precision, Tainan City, Taiwan). The motor's displacement is monitored to deduce the balloon's volume, and the pressure transducer reports the balloon's pressure. With the volume and pressure, the energies associated with the inflation and deflation processes can be obtained through numerical integration. A volume control unit such as the one in Figure 3 is certainly not the only way to apply mechanical load to a balloon. Alternatively, the mechanical load can be applied by prescribing the water's pressure instead of its volume (flow rate). Practically, such pressurecontrolled mechanical load can be implemented, for example, by connecting the balloon to the header tank of a rainwater harvesting system, see Section 2.1. However, balloons are known to possess non-monotonic P-V curves. Where the pressure and the volume do not increase with one another (a P-V curve's descending segment), the pressure-controlled mechanical load incurs the snap-through instability during inflation and the snap-back instability during deflation. These instabilities, in return, prohibit access to the balloon's pressure and volume at least over the descending segment. The pressure-controlled mechanical load thus makes a complete P-V curve unavailable and evaluating transferred energies impossible. In contrast, the volume control unit in Figure 3 prescribes the water's volume (flow rate) and measures the balloon's pressure. The resulting mechanical load is volume-controlled and incurs neither the snap-through nor the snap-back instability, thus making the entire P-V curve accessible through direct measurement. Because of these, the volume-controlled mechanical load is preferred and has been adopted in this study.
The temperature control unit applies the thermal load to the balloon and manages the heating and the cooling processes of Figure 1. This unit consists of the thermal bath, the water pump, the cold-water reservoir, and the thermocouple (TC) in Figure 3. The thermal bath is a small aluminum container with ten embedded heaters (not shown). It is filled with tap water at room temperature prior to each experiment. The balloon is immersed in the bath throughout the experiment. The temperature of the bath is monitored through the thermocouple. During the heating process, the heaters are turned on to raise the temperature of the water in the bath. This heats up the balloon. During the cooling process, the water pump delivers water at room temperature from the reservoir to chill the (hot) water in the bath. On the side wall of the bath, an overflow orifice is implemented to guide excess water back to the reservoir, thereby keeping a constant water level in the bath.
The signal processing unit in Figure 3 serves two purposes. The first purpose is to record the pressure of the balloon, the temperature of the thermal bath, and the displacement of the linear motor. The second purpose is to control the temperature of the thermal bath. These functions are implemented with a data-acquisition system (DAQ, model cDAQ-9178, NI-9203, NI-9403, NI-9212, and NI-9361, by National Instruments, Austin, TX,USA).
In experiments, the signals of the pressure transducer, the thermocouple, and the linear motor are recorded by the DAQ with a 25 Hz sampling frequency. For the heating process, the DAQ turns the heaters on in the beginning. Then it keeps track of the thermocouple measurement to shut the heating power when the prescribed thermal load temperature T t is reached. For the cooling process, the DAQ controls the water pump through an electromagnetic relay. No data are recorded during the cooling process because the balloon is flat during that process.
Following Figure 1, this study tests the response of a water balloon when it is subject to cyclic mechanical and thermal loads. These experiments require the collaboration of all units in Figure 3. To illustrate the experimental conditions, Figure 4 shows the mechanical load in (a) and the corresponding pressure of a balloon in (c). Starting at t = 0, the mechanical load is applied by injecting 9 mL of de-ionized water at 0.126 mL/s into the balloon. The corresponding speed of the linear motor is 0.01 mm/s, and the process takes 72 s. The heating process follows the inflation to start at t = 72 s, lasts for 180 s, and ends at t = 252 s. As mentioned, the thermal bath is heated during the process until the prescribed thermal load temperature T t is reached. If T t is reached within 180 s (say, 120 s), the balloon is left to relax for the rest of the time (i.e., 60 s). In this study, thermal loads are applied with T t = 30 • C, 40 • C, 60 • C, and 80 • C.  Each heated balloon is deflated to deliver pressurized water for power generation. This is the deflation process. In Figure 4a, the deflation starts at t = 252 s when the cylinder withdraws the water from the balloon at 0.126 mL/s. This is the same flow rate for the inflation. After the deflation, cold water is pumped from the reservoir into the thermal bath to cool the bath. The temperature of the reservoir is close to the room temperature, about 23 • C. The process continues for 3 min, which is sufficient to reduce the bath's temperature back to the room temperature.
Experiments without thermal load are also conducted to reveal the effects of heating. The setup in Figure 3 is used. However, the temperature control unit is not activated, and the balloon is only subject to cyclic mechanical load when it is immersed in the thermal bath filled with water at room temperature. For these experiments, Figure 4 shows the profile of the mechanical load in (b) and the corresponding pressure of the balloon in (d).
In the beginning, each balloon is inflated with 9 mL of water at 0.126 mL/s, the same flow rate as the experiments with the thermal load. The process takes 72 s. Subsequently, an inflated balloon is allowed to relax for 120 s until t = 192 s. Then the balloon is deflated with the same flow rate for its inflation. The deflation takes another 72 s and terminates at t = 264 s. A balloon without heating requires no cooling, and it is left to idle for only 30 s after its deflation and before its next inflation.
Each balloon in this study possesses a non-monotonic relation between its pressure and volume. As a reminder, the balloon is prepared from a clamped circular membrane that is originally flat and uniformly thick. For such a balloon, the membrane's curvature and (nonlinear) elasticity underlie the pressure-volume relation. The membrane's elastic modulus increases relatively slowly [63] when the balloon starts to grow in size and the stretch is relatively small. In that case, the influence of the elasticity upon the pressure is relatively minor, and the evolution of the pressure roughly follows that of the membrane's curvature. For the balloons in Figure 4, this causes the pressure increase between t = 0 s and t = 13 s and the pressure drop afterwards. Further inflating the balloons raises the membranes' elastic moduli more rapidly [63], and the influence of the elasticity becomes more significant. In Figure 4c,d, the growing influence of the elasticity makes the pressure decrease at more and more gentle rates. In fact, the pressure would start to rise if the balloons were further inflated beyond 9 mL after t = 72 s. Based on the discussion thus far, the competition between the membranes' curvature and elasticity causes the non-monotonic pressure evolution during the balloon's inflation.
The competition between the membrane's curvature and elasticity persists during the balloon's deflation. As such, the pressure also evolves non-monotonically as it returns to zero. In Figure 4c, the pressure evolution during the deflation resembles the time-reversed (i.e., from t = 72 s back to t = 0) pressure evolution of the inflation, and the pressure is slightly lower because of the membrane's relaxation. In Figure 4d, the balloon without heating relaxes much more, and the pressure drops significantly once the deflation starts. The pressure drop continues until t = 230 s when the membrane's curvature dominates again. Then the pressure gently rise to a maximum at t = 253 s and quickly returns to zero at t = 264 s.
The non-monotonic relation between the pressure and the volume is an interesting feature of the balloons. In Figure 4, the curve segments with opposite trends of pressure and volume evolutions cause the non-monotonicity by creating the descending segment in a balloon's in a balloon's P-V curve. The descending segment incurs the snap-through and snap-back instabilities when the mechanical load is pressure-controlled. These instabilities are attractive features of a balloon as they enable the invention of devices with fast responses [19,20,24]. For energy storage, the descending segment is also attractive because operating a balloon over this segment enables the balloon to be charged more and more easily during inflation. Additionally, it also makes the balloon deliver more and more momentum as it discharges its energy during the deflation. In this study, 9 mL is chosen as the maximum volume of each balloon to cover the descending segment of the P-V curve. Indeed, the purpose is to demonstrate the utility of these features to the advantage of energy storage.

Performance Indicators
In this study, a water balloon's performance is evaluated in terms of its service life, the energy it stores during its inflation, the energy it delivers during its deflation, its roundtrip efficiency in each cycle, and its average efficiency throughout its service life. All these parameters are available from the P-V curves of balloons. In what follows, the P-V curve of an inflation process shall be referred to as an inflation curve and that of a deflation process a deflation curve.
In each experiment, a balloon is tested until it fails. The balloon's service life is the maximum number of energy-storage cycle N max that the balloon survives, and it can be read off from the P-V curves without any numerical data analysis. As an example, Figure 5 shows the P-V curves of two consecutive cycles during which a balloon fails. In both plots, the arrows indicate the directions in which the balloon's state evolves. An arrow pointing to the upper right means that the curve segment next to it corresponds to an inflation process. Conversely, an arrow pointing to the lower left indicates that the balloon deflates along the curve segment nearby.
The failure of a balloon is recognized when the balloon's pressure starts to turn negative during a deflation. In Figure 5a, the balloon behaves as it has been in all previous cycles, and its pressure remains positive. The 64th cycle is the last one without a negative pressure and hence the one right before the balloon fails. In the 65th cycle, however, the pressure drops to zero as the balloon deflates but should still possess 2 mL of water. The 2 mL is apparently lost, and withdrawing that water ends up creating a suction in the balloon. This results in negative pressure, which is a clear sign of the balloon's failure. Because the balloon operates without failure for 64 cycles, the service life of that balloon is N max = 64 cycles.  Besides their service lives, the balloons are also examined in terms of the energies they store and release and their efficiencies of energy storage. During the j-th cycle, the energies that a balloon stores (W ij ) and releases (W oj ) are obtained by numerically integrating the pressure with respect to the volume as where V max = 9 mL, the first subscript specifies whether a parameter pertains to inflation (i) or deflation (o), and the second subscript j specifies the number of the cycle to which the parameters pertain. Throughout its service life, a balloon stores W i and delivers W o of energies in total. These total energies are obtained by summing over W ij and W oj as While W ij and W oj clearly reflect the evolution of a balloon's capacity, they are not useful for comparing the evolutions of the energies between different balloons. This is because two balloons typically have different pressures at the same volume since the beginning of the experiments, and the energies depend on the pressures according to Equation (1). To enable the comparison, W ij and W oj are normalized by W i1 , the energy for the balloon's first inflation, Through W * ij and W * oj , the evolutions of W ij and W oj can be directly compared among different balloons. Finally, the round-trip efficiency η j of the j-th cycle is defined as and the average efficiency η of a balloon throughout its service life is The definitions of the efficiencies in (4) and (5) essentially ignore the cost of the heat for boosting a balloon's pressure. This is because low-grade waste heat is meant to be used, and the cost of the heat is negligible. As mentioned in Section 1, the temperature of the heat is only tens of degrees Celsius above the room temperature. Such waste heat is indeed freely available from the solar radiation in summer or the condenser of an air-conditioning system. Consequently, an efficiency beyond 100% is possible.

Results and Discussion
This study includes 15 experiments, each testing a single water balloon as an energy storage medium. The key parameters are presented in Table 1, where T m is the average and σ T is the standard deviation of the thermal bath temperature during the deflation process. All other parameters are defined in Section 2. The experiments are organized by their thermal load temperatures T t and their chronological orders. No thermal load is applied in experiments 1-3, and their temperatures T t and T m are specified as room temperatures (RT). In the following discussions, each balloon shall be referred to by the identity of its experiment and each thermal load by its temperature T t . Accordingly, two or more thermal loads shall be referred to as the same thermal load if they are applied with the same temperatures T t ; otherwise, they are referred to as different thermal loads. Table 1. The parameters of all experiments in this study. For each balloon, T t is the thermal load temperature of its heating process, T m and σ T are the average (T m ) and the standard deviation (σ T ) of the thermal bath temperature during the deflation process, N max is the service life, W i and W o are the total energies it absorbs (W i ) and delivers (W o ) throughout its service life, W i1 is the energy it absorbs during its first inflation, and η is the average efficiency. According to Table 1, temperature overshoot (i.e., T m > T t ) is typical in the experiments. The overshoot occurs because raising the thermal bath's temperature to T t requires the heaters (and perhaps the aluminum container as well) to be warmer than T t . As such, the heaters continue to heat up the water even after their power is turned off. The overshoot could have been minimized if the heating were managed with a temperature controller. According to Table 1, however, T m of the same T t are fairly close. Additionally, σ T is less than 1.2 • C for balloons #4-#15 and less than 0.6 • C for 10 of the 12 experiments. These facts suggest that the overshoot is repeatable across different experiments. Because of this, no temperature controller is adopted for the experiments of this study.
The content of Table 1 provides an overview of the experiments. The relation between the temperatures T t , T m and σ T has just been reported, and the rest of the parameters are discussed in subsequent sections. In what follows, Section 3.1 investigates the factors that influence a balloon's service life, and Section 3.2 presents the experimental observations through the P-V curves of two selected balloons. Subsequently, Section 3.3 discusses the evolution of transferred energies, and Section 3.4 examines the efficiencies of WBES. Finally, Section 3.5 reveals an interesting scientific discovery about the permanent deformation of a flat membrane based on the experimental observations.

Service Life
According to Table 1, the balloons survive tens to hundreds of cycles depending on the thermal load. The balloons do possess different service lives even if they are subject to the same thermal loads. Overall, however, heating accelerates a balloon's failure. This is manifest from the correspondence between a higher T t and a smaller N max : without any thermal load, a balloon survives several hundred cycles. With T t = 30 • C, N max falls below 200 cycles. Further raising T t to 40 • C leads to an N max around 100 cycles, and T t = 60 • C reduces N max to 50-80 cycles. Finally, with T t = 80 • C, the balloons fail in 50 cycles. As mentioned earlier, stretching a rubber induces crystallization in the rubber, and the crystallites strengthen the rubber. However, heating retards the crystallization and weakens the rubber. These effects accelerate crack growth in the rubber and eventually shorten the service life of a balloon.
The service lives of balloons in Table 1 may seem too short to adopt WBES for practical utility. For the other energy storage systems mentioned in Section 1, indeed, the service lives are on the order of several years at least. For WBES without heating, however, the membrane needs replacement every 12-15 months if the energy-storage cycle repeats once a day as proposed in Section 2.1. When the thermal load is applied, the membrane needs replacement at least twice a year. Regarding this, it should be noted that the service lives in Table 1 are meant to reveal how N max varies with the thermal load (i.e., T t ). The values of N max are not the focus here because practically, N max can be increased by modifying the balloon. For example, the balloons can be prepared using membranes with a suitable filler [64]. Alternatively, the membranes can be tailored so that they are thicker wherever a uniformly thick membrane is more likely to start failing. In this study, however, the thinnest membranes available to the authors are purposefully chosen to expedite the fatigue failure. This restricts the durations of the experiments to hours or days instead of weeks, thus revealing the relation between N max and T t within an affordable time frame.

Evolution of P-V Curves
The P-V curves of balloons evolve as the loads repeat. As an illustration, the P-V curves of four selected cycles of balloon #2 (no heat) and balloon #14 (heated, 80 • C) are plotted in Figure 6a   The results in Figure 6a,b reveal several important facts about water balloon energy storage. First of all, heating alters a balloon's response to the mechanical load. This is manifest in the fact that balloon #2 and balloon #14 deflate differently ever since their first cycles. To contrast the responses of the balloons, Figure 6 overlays the P-V curves of their first cycles in cand those of their 11th cycles in d. For both cycles (and the cycles in between), heating helps balloon #14 maintain a higher pressure than that of balloon #2 at the same volume during its deflation. This reduces the area between the inflation and deflation curves and recovers a large proportion of the input energy that would otherwise be wasted.
Second, heating shortens a balloon's service life, as discussed in Section 3.1. In Figure 6a, balloon #2 does not respond with a negative pressure (i.e., fail) until its 501st cycle. In Figure 6b, however, balloon #14 fails during its 49th cycle. These results suggest that heating can effectively reduce the service life of a balloon by one order of magnitude. This is undesirable for an energy storage system. Despite the differences just mentioned, the P-V curves of balloon #2 and balloon #14 do have common features. To illustrate these features, Figure 6e replots three inflation curves of Figure 6a, and Figure 6f does three deflation curves of Figure 6b. In Figure 6a,b, the deflation process never returns to where the inflation of the same cycle starts. Instead, each cycle leaves a dead volume below which the balloon possesses negligible pressure and stores no energy. In fact, each cycle contributes a small increment to the dead volume. As such, the dead volume increases as the cycle repeats. This is manifest from the inflation curves in Figure 6e and the deflation curves in Figure 6f. In a later cycle, balloon #2 ( Figure 6e) requires a larger volume to start building up its pressure. On the other hand, the pressure of balloon #14 (Figure 6f) vanishes at larger and larger volumes as the cycle repeats. By raising the dead volume, the repetition of the cycle shifts the inflation and deflation curves of one cycle to the right of those in the previous cycles. Besides the rightward shift, the P-V curves also migrate downwards. That is, the inflation curve of each cycle lies below those of all previous cycles: in Figure 6e, the inflation curve of the 11th cycle lies below that of the first cycle, and the inflation curve of the 301st cycle is beneath those of both the first and the 11th cycles. This is also the case for each deflation curve according to Figure 6f. These features are common to all balloons in this study. On a pressure-volume diagram like Figure 6, therefore, the responses of the balloons only evolve towards the lower right. Accordingly, a balloon accommodates less and less energy as the cycle repeats.

Transferred Energy
For each balloon in this study, Table 1 reports the total energies that the balloon absorbs (W i ) and delivers (W o ) throughout its service life. Depending on the temperature T t , W i and W o vary between several Joules to tens of Joules. These energies may seem small enough to render WBES completely useless. However, the energies can be readily scaled up in many different ways. For example, thicker and stiffer membranes can be adopted to implement balloons of the same base radius. Alternatively, the same membrane can be used to implement an array of balloons of a smaller base radius. With these changes, a WBES system would possess a higher pressure and can store more energy with the same amount of water.
Heating reduces W i and W o , as the data in Table 1 manifest. To find the cause of the decrease, one should examine how the absorbed (W ij ) and delivered (W oj ) energies evolve and accumulate. As mentioned in Section 2.2, W ij and W oj are scaled by W i1 so that their evolutions can be compared between different balloons in terms of W * ij and W * oj . The energy W i1 is chosen for the scaling because it is based on the data before a balloon is heated for the first time. For different balloons, W i1 conveys the difference between their assemblies without any effect of the thermal load. Scaling W ij and W oj by W i1 thus preserves the effects of heating upon these energies. As shown in Table 1, W i1 varies between 0.174 J and 0.209 J, and the variation appears to be independent of T t . Recall that each balloon in this study is assembled by clamping a membrane around the end of a tube. The stress for clamping a membrane is not precisely controlled, and the membrane of each newly assembled balloon warps slightly and differently. The different extents of warping then differentiate the pressures of the balloons and cause the ±9% variation of W i1 .
In terms of W * ij and W * oj , Figure 7 compares the evolutions of the transferred energies between balloons subject to different thermal loads. For each thermal load, the data pertain to the balloon with the longest service life, c.f. Table 1. Only these data are plotted to avoid cluttering. According to Figure 7, balloons accommodate less and less energy as the cycle repeats. This is consistent with the observations from Figure 6, in which P-V curves evolve towards the lower right. For energy storage, the decays in W * ij and (especially) W * oj indicate a capacity reduction, which is undesirable. Two factors possibly underlie the decays, both are related to the Mullins effect [65,66]. First of all, the cyclic mechanical load increases the balloon's dead volume. Second, the inflation very likely damages the membrane on the microscopic scale. As the damage accumulates, it weakens the membrane and reduces the balloon's pressure. As discussed in Figure 6, these factors shift the P-V curves to the lower right as the cycle repeats. Admittedly, the microscopic damage is speculated but has not been examined in this study. However, the speculation seems consistent with the experimental observations. Inflating a balloon stretches its membrane, and the damage should presumably increase with the stretch. Accordingly, the apex of a balloon is the balloon's most vulnerable spot as it endures the largest stretch [2,3,67]. In experiments, indeed, the neighborhood of the apex is where all balloons are observed to leak at the end of their service lives. The results in Figure 7 also reveals how heating reduces W i and W o through its impact over W * ij and W * oj . According to Figure 7a, a higher T t causes a faster decay of W * ij . For example, W * ij stays above 0.7 until the 370th cycle when no thermal load is applied. When the balloon is heated with T t = 30 • C, W * ij decays faster and drops below 0.7 in the 132nd cycle. Raising T t to 40 • C makes W * ij reach 0.7 in the 100th cycle. With T t = 60 • C, W * ij = 0.7 occurs in 49 cycles. With 80 • C, finally, W * ij reaches 0.7 in the 35th cycle. In addition to accelerating the decay of W * ij , heating also shortens the service life. Because heating makes a balloon absorb less energy in every cycle and survive fewer cycles, a lower W i at a higher T t is a natural consequence.
Heating reduces W o in a similar way to how it influences W i , although it changes W * ij and W * oj differently. According to Figure 7b, heating raises W * oj by several percent during the first 20 cycles. This reflects the situation in Figure 6c,d, where the heating recovers more of the input energy. However, heating can reduce the service life up to an order depending on T t , thereby imposing a larger sanction on the accumulation of W o . Much like W i , W o is decreased by heating because of the balloon's shortened service life.
According to Figure 7, heating seems to influence the transferred energies differently during different stages of a balloon's service life. For both W * ij and W * oj , heating does not seem to exacerbate the decay significantly during the first 20 cycles. This is manifest from the similar slopes of the data in each plot during these cycles. Accordingly, the decay of W * ij and W * oj are independent of the thermal load during the early stage of a balloon's service life. Subsequently, the presence of a thermal load leads to faster decays in energies, and a higher T t makes the decays occur earlier. These effects are more noticeable from Figure 7b. Without being heated, balloon #3 shows no clear sign of significant slope change in its data throughout its service life. For the heated balloons, the faster decay starts around the 80th cycle for T t = 30 • C, the 60th cycle for T t = 40 • C, the 30th cycle for T t = 60 • C, and the 20th cycle for T t = 80 • C. Similar trends of decays can also be found in Figure 7a for W * ij . The softening of a balloon's membrane underlies the faster decays of W * ij and W * oj . The softening is expected to result from the accumulation of the microscopic damage in the membrane speculated previously. To keep a balloon's energy-storage capacity, one should adopt a rubber with a better resistance to crack growth. Alternatively, one can tailor the thickness profile of the membrane to re-distribute the stress and stretch more evenly. These modifications shall strengthen the membrane, mitigate crack growth, and retard the decays of the energies.

Efficiency
The high efficiency is the most attractive feature of WBES. According to Table 1, the average efficiency η varies between 81% and 100%. With waste heat as low-grade as 30 • C, the efficiency can be pushed beyond 90%. The round-trip efficiency η i underlies the attractive values of η, and Figure 8 shows the evolutions of η i for balloons subject to different thermal loads. For each balloon, η i increases drastically in the second cycle. Then it persists until close to the end of the balloon's service life. Without the help of heating, the efficiency varies between 85%-90%. Heating improves the efficiency of WBES, and the improvement depends on the temperature T t . With T t = 30 • C, η i stays around 90%. Raising T t to 40 • C increases η i by 5%. Further raising T t to 60 • C raises the efficiency by another 5%, pushing η i close to 100%. In that case, the water balloon recovers almost all the work that the water does to inflate it. Finally, adopting T t = 80 • C enables the balloon to achieve η i beyond 100%. These results suggest that a water balloon can effectively utilize low-grade waste heat to boost its efficiency. To avoid any confusion, it should be emphasized that the efficiencies do not account for the cost of the low-grade heat. With the heat being freely available, an efficiency beyond 100% becomes possible and does not violate the conservation of energy.

Permanent Deformation
As discussed with Figure 6, the cyclic mechanical load creates a dead volume that increases monotonically as the cycle repeats. From an engineering perspective, the dead volume is undesirable as it keeps no pressure, stores no energy, and degrades the balloon's performance. Scientifically, however, it is interesting to note that the dead volume seems to approach the volume of a hemisphere with the same base radius as that of the balloons. This is most clearly shown by comparing the balloon's base radius r b to its average radius of curvature R a corresponding to the dead volume V d . The radius R a is the radius of curvature of a spherical cap with base radius r b and volume V d . The parameters R a , r b and V d are related by In Figure 9, the ratio of R a to r b is plotted against the number of cycles for most of the balloons in Figures 7 and 8a. To avoid cluttering, the curve for balloon #14 (T t = 80 • C) is not plotted because it lies right between those of balloon #8 and #10. For each balloon in Figure 9, R a exhibits an overall decreasing trend with some occasional fluctuations. Additionally, R a decreases at more and more gentle rates as the cycle repeats. The heated balloons fail relatively early before their R a /r b ratios have a chance to exhibit a clear lower bound. For balloon #3, which is not heated, R a ends up becoming 1.003r b . Based on these, r b seems to be the lower limit that R a 'strives' to approach over the several hundred cycles. For spherical caps with the same base radius, a hemisphere is the one with the smallest radius of curvature. The observation apparently suggests that a pinned circular membrane evolves to approximate a hemisphere and minimize its radius of curvature when it is subject to cyclic mechanical load.

Conclusions
This study proposes the use of water balloons as an energy storage medium and tests their performances experimentally. In each experiment, a water balloon is subject to cyclic mechanical and thermal loads, and its response is examined in terms of their pressures, transferred energies, efficiencies, and service lives. With the results already presented, a first contribution of this study is the discovery of the excellent efficiencies that a water balloon achieves as an energy storage medium. A second contribution is the demonstration of recycling low-grade waste heat with a water balloon. Finally, a third contribution is the disclosure of the trade-off between the efficiency and the service life of a balloon. Candidate solutions are also proposed to enlarge the capacity and extend the service life of a balloon when it serves as an energy storage medium. These results shall serve as a useful guide for developing balloon-based thermo-mechanical devices in a broad range of applications including but not limited to the storage of energy.

Data Availability Statement:
The data that support the findings of this study are available within the article. Further details are available from the corresponding author upon reasonable request.