<<

University of Kentucky UKnowledge

Theses and Dissertations--Biosystems and Agricultural Engineering Biosystems and Agricultural Engineering

2021

MODELING ENERGY FLOWS IN FLOATING IN- RACEWAYS UTILIZING SOLAR POWER BACK-UP

Bo Smith University of Kentucky, [email protected] Author ORCID Identifier: https://orcid.org/0000-0003-1142-1823 Digital Object Identifier: https://doi.org/10.13023/etd.2021.129

Right click to open a feedback form in a new tab to let us know how this document benefits ou.y

Recommended Citation Smith, Bo, "MODELING ENERGY FLOWS IN FLOATING IN-POND RACEWAYS UTILIZING SOLAR POWER BACK-UP" (2021). Theses and Dissertations--Biosystems and Agricultural Engineering. 78. https://uknowledge.uky.edu/bae_etds/78

This Master's Thesis is brought to you for free and open access by the Biosystems and Agricultural Engineering at UKnowledge. It has been accepted for inclusion in Theses and Dissertations--Biosystems and Agricultural Engineering by an authorized administrator of UKnowledge. For more information, please contact [email protected]. STUDENT AGREEMENT:

I represent that my thesis or dissertation and abstract are my original work. Proper attribution has been given to all outside sources. I understand that I am solely responsible for obtaining any needed copyright permissions. I have obtained needed written permission statement(s) from the owner(s) of each third-party copyrighted matter to be included in my work, allowing electronic distribution (if such use is not permitted by the fair use doctrine) which will be submitted to UKnowledge as Additional File.

I hereby grant to The University of Kentucky and its agents the irrevocable, non-exclusive, and royalty-free license to archive and make accessible my work in whole or in part in all forms of media, now or hereafter known. I agree that the document mentioned above may be made available immediately for worldwide access unless an embargo applies.

I retain all other ownership rights to the copyright of my work. I also retain the right to use in future works (such as articles or books) all or part of my work. I understand that I am free to register the copyright to my work.

REVIEW, APPROVAL AND ACCEPTANCE

The document mentioned above has been reviewed and accepted by the student’s advisor, on behalf of the advisory committee, and by the Director of Graduate Studies (DGS), on behalf of the program; we verify that this is the final, approved version of the student’s thesis including all changes required by the advisory committee. The undersigned agree to abide by the statements above.

Bo Smith, Student

Dr. Joseph Dvorak, Major Professor

Dr. Donald Colliver, Director of Graduate Studies

MODELING ENERGY FLOWS IN FLOATING IN-POND RACEWAYS UTILIZING SOLAR POWER BACK-UP

______

THESIS ______A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Biosystems and Agricultural Engineering in the Colleges of Agriculture and Engineering at the University of Kentucky

By Bo H. Smith Lexington, Kentucky Director: Dr. Joseph Dvorak, Professor of Biosystems and Agricultural Engineering Lexington, Kentucky 2021

Copyright © Bo H. Smith 2021 https://orcid.org/0000-0003-1142-1823

ABSTRACT OF THESIS

MODELING ENERGY FLOWS IN FLOATING IN-POND RACEWAYS UTILIZING SOLAR POWER BACK-UP

The In-pond (IPR) is a novel option for production , depending on water moving devices to provide constant flow. Device failure may result in catastrophic fish loss, requiring power backup systems to mitigate risk in case of power outages. Because these systems must be dependable and many suitable locations are remote, off-grid solar photovoltaic (PV) systems with battery storage have been considered since they eliminate need for utility power. Such systems can be hard to size and expensive. This study modeled system requirements using an energy balance to determine whether systems could withstand varying loads based on climatological conditions. Sizing was iterative, with battery storage and panel size increasing until the model predicted continuous power was provided year-round. This study found failure events were clustered over multiple days in winter. Therefore, it determined undersized systems were suitable if there was no stocking in these months. Further work found an integrated generator backup system would decrease necessary system size. Likewise, substitution of continuous motor loadings with variable speed motors operated based on need may further decrease system demand. The presented modelling approach has broad implications for feasibility of IPR systems, providing reduced startup costs and possibilities for greater implementation of this novel technology.

KEYWORDS: Aquaculture; Distributed Power Generation; Energy Modeling; Photovoltaics; Renewable Energy Technologies; Solar Power

Bo H. Smith

04/20/2021 Date

MODELING ENERGY FLOWS IN FLOATING IN-POND RACEWAYS UTILIZING SOLAR POWER BACK-UP

By

Bo H. Smith

Dr. Joseph Dvorak Director of Thesis

Dr. Donald Colliver Director of Graduate Studies

04/20/2021 Date

DEDICATION

To the people of Hopkins County and to the rural people of Kentucky and the United States at large, who shape my worldview and bring inspiration to my life daily. I had the good fortune to come of age as a rural person, and have seen first-hand that, though we may be but little, we are fierce. In all that I do, I remember my roots.

ACKNOWLEDGMENTS

The following thesis, while an individual work, benefited from the insights and direction of several people. First, my Thesis Chair, Dr. Joseph Dvorak, exemplifies the high-quality scholarship to which I aspire. Starting as an undergraduate and moving into my graduate school career, he trusted in me and my abilities as a researcher and granted me the opportunity to thrive in my learning. To say that I am grateful would be an understatement. Next, I wish to thank the complete Thesis Committee: Dr. Tiffany

Messer, Dr. Donald Colliver, and Dr. Kenneth Semmens. Each individual provided insights that guided and challenged my thinking, substantially improving the finished product.

In addition to the technical and instrumental assistance above, I received equally important assistance from family and friends. My partner, Harper Oglesby, provided on- going support throughout the thesis process, as well as technical assistance critical for completing the project in a timely manner. I also thank my parents, Gary and Melissa

Smith, for providing significant support and motivation to me throughout my studies, without which, this would not have been possible. Finally, I wish to thank the many friends I brought into my inner circle along the way toward the completion of this degree.

From my tight circle of former roommates to the executive teams of Alpha Epsilon and the BAE Graduate League of Students, I have accrued a large support network, which I have leaned on heavily in my time at UK.

In addition, I would like to thank all of those from both the University of

Kentucky and Kentucky State University who pitched in where they were needed to assist me in the completion of this project. Namely, I thank Cynthia Arnold, James Brown, John

Kelso, and Dawson Armstrong for providing their timely and instructive comments and

iii evaluation at every stage of the thesis process, allowing me to complete this project on schedule.

iv

TABLE OF CONTENTS

ACKNOWLEDGMENTS ...... iii

LIST OF FIGURES ...... ix

LIST OF ADDITIONAL FILES ...... xi

CHAPTER 1. INTRODUCTION ...... 1 1.1 Introduction ...... 1 1.2 Literature Review ...... 2 1.2.1 Floating IPR Systems ...... 2 1.2.2 Determining Device Performance ...... 5 1.2.3 Decision Support Tools ...... 7 1.2.4 Off‐Grid Solar PV Modeling ...... 8 1.2.5 Backup Generator Systems...... 10 1.3 Objectives ...... 11 1.4 Thesis Outline ...... 11

CHAPTER 2. METHODOLOGY ...... 13 2.1 Site Description ...... 13 2.2 Raceway Flow Data Collection ...... 13 2.3 System Performance Characteristics ...... 16 2.4 Modeling Parameters ...... 18 2.4.1 Available Dissolved Oxygen ...... 18 2.4.2 Daily Energy Demand ...... 20 2.4.3 Solar and Climate Data ...... 21 2.4.4 Tools of the Trade: SAM and PVWatts ...... 22 2.4.4.1 PV Watts ...... 22 2.4.4.2 SAM ...... 23 2.4.5 Battery System Sizing and Efficiency ...... 23 2.5 Iterative System Sizing Models ...... 25 2.5.1 Off‐Grid Solar PV Plus Battery (No Generator) ...... 25 2.5.2 Off‐Grid Solar PV Plus Battery (With Generator) ...... 26 2.6 Conceptualizing a Variable Speed Propeller System ...... 29

CHAPTER 3. RESULTS AND DISCUSSION ...... 31 3.1 IPR System Sizing Tool ...... 31 3.2 Estimating Fish Carrying Capacity ...... 32 3.3 Off‐Grid Solar Recommendations...... 32

v

3.3.1 Off‐Grid Solar PV Plus Battery (No Generator) ...... 33 3.3.2 Off‐Grid Solar PV Plus Battery (With Generator) ...... 38 3.4 Implemented Variable Speed Propeller System ...... 47

CHAPTER 4. CONCLUSIONS AND FUTURE WORK ...... 51 4.1 Conclusions ...... 51 4.2 Future Work ...... 52

APPENDICES ...... 54 [APPENDIX 1. IPR DECISION SUPPORT TOOL] ...... 54 [APPENDIX 2. SOLAR MODEL CODE] ...... 66 [APPENDIX 3. GENERATOR CODE] ...... 73

REFERENCES ...... 81

VITA ...... 90

vi

LIST OF TABLES

Table 2.1 A summary of findings for device performance characterization...... 18 Table 3.1 Carrying capacity for channel (kg) in a raceway assuming DO concentrations of 4 mg l-1 or higher in the pond with raceway effluent DO of 3 mg l-1 at 12.9°C...... 32 Table 3.2 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Lexington, KY...... 33 Table 3.3 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Louisville, KY...... 34 Table 3.4 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Jackson, KY...... 34 Table 3.5 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Paducah, KY...... 34 Table 3.6 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Cincinnati, OH/Covington, KY...... 35 Table 3.7 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Bowling Green, KY...... 35 Table 3.8 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in London, KY...... 35 Table 3.9 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Lexington, KY which use Fewer Than 100L per Year of Propane...... 39 Table 3.10 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Louisville, KY which use Fewer Than 100L per Year of Propane...... 39 Table 3.11 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Jackson, KY which use Fewer Than 100L per Year of Propane. ..39 Table 3.12 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Paducah, KY which use Fewer Than 100L per Year of Propane. .40 Table 3.13 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Cincinnati, OH/Covington, KY which use Fewer Than 100L per Year of Propane...... 40

vii

Table 3.14 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Bowling Green, KY which use Fewer Than 100L per Year of Propane...... 40 Table 3.15 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in London, KY which use Fewer Than 100L per Year of Propane. ..41 Table 3.16 Comparison of system sizes for Solar PV plus Battery Systems with and without Backup Generators in Lexington, KY...... 41 Table 3.17 Estimated worst-case generator fuel consumption values by Equation 8. These values were found by assuming a worst-case operation of 8760 operation hours for year- round operations and 5880 hours for seasonal operation...... 42

viii

LIST OF FIGURES

Figure 1.1. Floating raceways described by Hartleb (2004)...... 4 Figure 2.1 The site location at the Harold R. Benson Research and Demonstration Farm in Frankfort, KY ...... 13 Figure 2.2 Flow testing as completed within the Floating In-pond Raceway ...... 15 Figure 2.3 The three installed flow generation devices, clockwise from top: a 1 Hp airlift, 3/4 Hp Propeller aerator, and the 1/15 Hp Slow-Rotating Paddlewheel ...... 17 Figure 2.2 A map of ASOS site locations from around the state of Kentucky utilized in this research...... 21 Figure 3.1 The application to size In-pond Raceway production systems ...... 31 Figure 3.2 Available battery capacity for a system utilizing 10kW Panels and a 40 kWh Battery to power a 0.25 kW load at the Lexington site location...... 37 Figure 3.3. Number of hours at different battery state of charge ((a) 0% to 50%, (b) 40% to 80%, and (c) 70% to 100%) over the 22-year dataset for the 10kW Panel and 40 kWh Battery while supporting a 0.25 kW load in Lexington, KY...... 38 Figure 3.4 System fuel use trends for a 0.25kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries...... 43 Figure 3.5 System fuel use trends for a 0.25kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario of no Solar Panels and no batteries...... 43 Figure 3.6 System fuel use trends for a 0.5kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries...... 44 Figure 3.7 System fuel use trends for a 0.5kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario of no Solar Panels and no batteries...... 44 Figure 3.8 System fuel use trends for a 0.75kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries...... 45

ix

Figure 3.9 System fuel use trends for a 0.75kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario of no Solar Panels and no batteries...... 45 Figure 3.10 System fuel use trends for a 1-kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries...... 46 Figure 3.11 System fuel use trends for a 1-kW load in Lexington, KY for Year-Round Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario of no Solar Panels and no batteries...... 46 Figure 3.12 The industrial control box as installed in Frankfort, KY ...... 48 Figure 3.13 The industrial control system as installed, including all controls and the 1.5kVA transformer...... 49

x

LIST OF ADDITIONAL FILES

Supplemental Data. SAM Application Outputs ...... ZIP 238,973 KB Supplemental Data. Solar Failure Tables ...... XML 327 KB Supplemental Data. Generator Table Outputs ...... ZIP 642 KB Supplemental Data. PLC Send-Receive Code ...... PDF 62 KB

xi

CHAPTER 1. INTRODUCTION

1.1 Introduction

Recently, interest in aquaculture has grown among small producers (FAO, 2020; National

Marine Service, 2020). One newer method of production is the In-pond Raceway

(IPR), which relies on an aerator device to maintain sufficient flow of water for fish confined in a flowing water system. Installing and operating IPR systems can be expensive.

Studies have found total construction and implementation costs for a raceway system to be in excess of $113,000 (Brown, Chappell, & Boyd, 2011), with the operational cost of energy for constant aeration accounting for as much as 10% of total production costs

(Fullerton, 2016). Since these devices are required to operate continuously, power back-up systems are also an essential added cost in the case of a power failure (Masser, 2012).

Because of this expense and the remote nature of many farm , there is interest in operating independent of the electrical grid, and solar energy systems pose a novel and sustainable option. These systems also enable the use of ponds in more rural locations where the expense of a connection to the electrical grid can be exorbitant and power failures are often more frequent or continue for longer durations (Mukherjee, Nateghi, & Hastak,

2018). However, starting costs for a PV system with battery storage can be high depending on the scale needed to meet capacity of the battery storage and size of the PV array.

To address these concerns, this thesis research was completed to model energy flows and consumption by water moving devices, with the intent of developing a series of model equations to assist in sizing of cost-effective PV plus battery backup systems, as well as examining novel technologies to reduce operational power costs for IPR systems.

1

1.2 Literature Review

1.2.1 Floating IPR Systems

Commercial aquaculture production is typically performed in large, specialized facilities such as split-ponds or fixed raceways (Fornshell, Hinshaw, & Tidwell, 2012; Tucker &

Hargreaves, 2012). Recently, there has been a push toward greater intensification of production (Brune, Schwartz, Eversole, Collier, & Schwedler, 2003; Brune, Tucker,

Massingill, & Chappell, 2012) and facilities which can utilize existing ponds (Brown,

Chappell, & Boyd, 2011; Füllner, Gottschalk, & Pfeifer, 2007). IPRs show great promise, allowing for more efficient fish production in existing freshwater ponds, though at a greater net energy input than the traditional systems (Brown, Hanson, Chappell, Boyd, & Wilson,

Jr., 2014) As world land and water resources remain static, food source demands have increased (Gao, 2012; Fukase & Martin, 2017), and renewable energy has become cheaper and more abundant (Ellabban, Abu-Rub, & Blaabjerg, 2014). Therefore, energy independent systems are becoming more relevant and may be a useful solution for existing enterprises.

Another important advantage of floating IPR systems is that they are able to utilize existing water resources. This stands to benefit producers in Kentucky especially, as the state has more than 2700 lakes and reservoirs, more than a third of which are larger than 10 acres and 46 inches of precipitation annually to fill them (Kentucky Geological Survey, 2020).

Floating raceway technology may facilitate production of fish in areas which have not previously been considered, offering producers within the state the chance to diversify their income in unique and exciting ways.

2

However, floating IPR systems are not without concern. The main driving factor for producers tends to be profit (Shang & Costa-Pierce, 1983), and IPR systems are limited in their production capabilities by factors such as waste production and their ability to assimilate into the pond environment (Brown, Chappell, & Hanson, 2010). These issues may cause IPR systems to be less profitable than comparable systems (Kumar, et al., 2018), with investigators at Auburn University documenting that cost of production for catfish grown in floating raceways at a small scale was slightly higher than pond culture

(Bernardez, 1995; Masser, 2012). Brown et al. (2011) found total construction and implementation costs for a large raceway system to be in excess of $113,000, while Hartleb

(2004) reported the cost of a flowing water system with three much smaller floating IPR units designed by Superior Raceways, LLC (Superior Raceway Systems, 2018) to be nearer to $10,000.

If the water body is large enough, it is possible to take advantage of the scalable nature of floating raceways by adding additional units. However, there is a large requirement of energy to each raceway in these systems, as water must constantly be flowing into the system (Masser, 2012). While moving water is required for waste exchange, the most demanding requirement for a continuous flow of water is to ensure sufficient oxygen is available to the fish. A recent study reported that energy used for aeration is around 10% of total operating cost (Fullerton, 2016). However, within newly developed, floating in- pond systems, there is some uncertainty about the most effective methods of water movement in terms of energy usage and water moving efficiency.

3

Figure 1.1. Floating raceways described by Hartleb (2004).

Oxygen demand by a crop of fish and availability of dissolved oxygen in the water flow

are critical elements in determining the minimum safe conditions for fish production

(Kramer, 1987). Unfortunately, these parameters are highly variable. Oxygen availability

is set by the level of dissolved oxygen in the water. This can vary based on growth,

which itself is impacted by temperature, weather, sunlight, pH, and turbidity (Round, 1981;

Wurts & Durborow, 1992; Boyd & Tucker, 1998). There are significant variations of

dissolved oxygen depending on rates of algae growth and respiration, with the highest natural concentrations occurring in colder months and the lowest concentrations occurring just before dawn in late summer (Boyd, Romaire, & Johnston, 1978; Ouyang, Nkedi-Kizza,

Wu, Shinde, & Huang, 2006). The amount of oxygen required by the fish varies with number of fish, fish size, fish species, and water temperature (Masser, 2012). Certain events such as feeding can cause a temporary surge in activity and a resulting increase in

4

oxygen demand. To maintain growth rates, dissolved oxygen should be higher than 3.5 mg

l-1 for fish to feed well (Andrews, Murai, & Gibbons, 1973; Andrews & Matsuda, 1975;

Torrans, 2005), and significant fish mortality may result once concentrations fall below 2

mg l-1 (Lai-Fa & Boyd, 1988).

1.2.2 Determining Device Performance

Much work has been done to determine device performance in larger systems, with the

assumption that findings would scale with smaller systems. In split pond aquaculture, slow rotational paddlewheels (SRP) were found to be the most efficient water moving devices

(Brown & Tucker, 2013). These results were confirmed in a later study (Brown, Tucker,

& Rutland, 2016), and it was determined that paddlewheels are the most efficient water moving devices in split pond aquaculture in terms of supplied flow per unit power when specifically compared to commercially available paddle aerators (Brown & Tucker, 2014).

Not compared, however, were two popular methods of water movement in small-scale production: the propeller aerator (Rivara, Tetrault, & Patricio, 2002) and an airlift device powered by regenerative blower (Huang, et al., 2016; Parker & Suttle, 1987).

Measuring flow within these structures has also traditionally been an issue. In large-scale production systems for algae, computational fluid dynamics (CFD) modeling has been used to understand mixing of flow in a and report its flow field graphically (Huang, et al., 2015). However, complicated CFD modeling is often not the most effective method for producers measuring flow in small scale production facilities, as they often lack the financial resources and technical expertise. Previous studies have shown that measurements taken by an Acoustic Doppler Velocimeter (ADV), a device which 5

monitors the relative shift in frequency of sound waves to estimate water velocity, within

the raceway environment to be trustworthy (Viadero, Jr., Rumberg, Gray, Tierney, &

Semmens, 2006; Fizer, Gray, & Semmens, 2013). These measurements also supported use

of the ADV in future studies, showing that measuring flow within the raceway is possible

with such a device and would be reliable for a more detailed engineering analysis.

There is a question, however, of which protocol to follow in flow measurements with the

ADV. One option is the protocol for use of an ADV device in stream gauging (Rehmel,

2007). In this method, the USGS standard for stream gauging is followed with a SonTek

ADV in place of a standard flowmeter. The study found that ADV measurements were not

statistically dissimilar to the standard and had reduced error in more turbulent flows,

meaning an ADV can be used to measure discharge in a greater variety of streams and

channels than the standard, given that USGS protocol is followed.

Besides flow measurement, systems must also be effective at improving water quality.

Typically, this has been gauged by measuring the system’s ability to meet nighttime

dissolved oxygen demand (Boyd, Romaire, & Johnston, Predicting Early Morning

Dissolved Oxygen Concentrations in Channel Catfish Ponds, 1978). More recently, work

has been done to clarify the difference between organic matter and organic carbon loadings

to set appropriate budgets (Brown, Boyd, & Chappell, 2015). This work concluded that

estimates of fish production based on feed rate can be useful, but there must be additional

measures taken to ensure nutrient loading in the water stays at a safe level in this system

of production. Work has also been done to quantify the effects of waste settling and

removal in raceway systems (Cripps & Bergheim, 2000), finding that this distribution of

solid wastes will deplete water available oxygen concentrations by increasing rates of

6 microbial respiration if not accounted for and/or removed. Due to these considerations, water quality must be monitored continuously to determine effects of water moving devices on the raceway environment.

As with other industries, most opportunities for improvement in aquaculture will depend on some degree of automation. Both flow measurements and water quality measurements will be necessary for any future system which depends on automation (Chen, Sung, & Lin,

2015). Adding these measurements will result in already energy expensive systems becoming more financially taxing and requiring greater technical support. Such a complex system is useful in large-scale facilities, but any model for smaller producers would likely need to be less involved, possibly relying on only a few important variables to control a small range of outputs. In addition, these systems may require decision support tools for system design, especially in the sizing of backup power systems.

1.2.3 Decision Support Tools

Demand for practical decision-making tools will increase as there is greater reliance on technology in agriculture. Advances in agriculture soon will manifest a greater need for decision support systems, which can import large sums of available data behind an easy- to-navigate user interface (Zhai, Martínez, Beltran, & Martínez, 2020). A number of examples of such systems exist in modern agriculture: crop-field suitability (Mbugwa,

Prager, & Krall, 2015), expected economic benefit of grain storage practices (Dvorak,

Shockley, Mason, & McNeill, 2018), pest populations (Rupnik, et al., 2019), managing irrigation (Navarro-Hellín, Martínez-del-Rincon, Domingo-Miguel , Soto-Valles, &

Torres-Sánchez, 2016), environmental threat from the changing climate (Han, Ines, & 7

Baethgen, 2017), and estimating level of food security (Ferjani, Mann, & Zimmermann,

2018), among others. The challenge for this study was to follow the examples of previous decision selection tools and of relevant power models to build a sizing system for novel sustainable energy systems used in IPR aquaculture.

Modeling to size system infrastructure in aquaculture is itself a novel technique. Computer modeling in the aquaculture industry typically focuses on topics such as spatial distribution and carrying capacity (Filgueira, Grant, & Strand, 2014), economic viability and ecological impact (Nobre, Musango, de Wit, & Ferreira, 2009), or to estimate disease prevalence (de

Blas, Muniesa, Vallejo, & Ruiz-Zarzuela, 2020). There has been research on the necessary levels of aeration and oxygen supply (Andrews, Murai, & Gibbons, 1973; Andrews &

Matsuda, 1975; Boyd, Romaire, & Johnston, Predicting Early Morning Dissolved Oxygen

Concentrations in Channel Catfish Ponds, 1978; Lai-Fa & Boyd, 1988) and the costs of power to run aeration devices are included in operational budgets (Brown, Chappell, &

Boyd, 2011; Fullerton, 2016). However, to date, research into IPR systems has not directly investigated optimization of the power and energy requirements associated with this aeration.

1.2.4 Off-Grid Solar PV Modeling

Once considered expensive and impractical, solar PV systems are expected to greatly decrease in cost in the coming decade (United States Energy Information Administration,

2020). Studies completed in the Philippines now exhibit solar PV systems are useful in aquaculture (Hendarti, Wangidjaja, & Septiafani, 2018) and may reduce lifetime system cost as compared to gas-powered generators (Hendarti & Septiafani, 2020). Thus, the burden of cost has shifted away from solar PV systems and fallen more on systems storing 8

produced energy. Lead-acid batteries are traditionally the most popular method of energy

storage, but lithium-ion based batteries are becoming more popular because of their higher

energy density (Chang, 2019). As batteries are highly dependent on environmental

conditions, they often need to be oversized, increasing expense. Thus, it is important to

correctly size the battery banks for the systems and eliminate unnecessary costs in order to

achieve the full benefit of a solar PV system (Comello & Reichelstein, 2019).

In the more common grid-tied system, solar energy is converted to a useful form and

supplied to the electrical grid (Fernandez-Infantes, Contreras, & Bernal-Agustin, 2006;

Yang, Li, Zhao, & He, 2010). In a non-grid-tied system, this energy is instead carried to a

battery array for storage where it can be accessed as necessary by connected electrical

devices (Mohanty, Sharma, Gujar, Kolhe, & Azmi, 2016). These systems risk insufficient

storage, as weather conditions are variable and if solar radiation falls below the minimum

threshold to meet energy demands, it is possible battery storage will be depleted and the

system will not function until the batteries can be recharged.

In the past, energy modeling has been used to size solar PV systems. Because of the recent

push to utilize more renewable energy sources, research into energy storage and generation

in solar PV systems is plentiful; however, they were often viewed in terms of economic

optimization (Hesse, et al., 2017) or with a focus on carbon offset for residential systems

(Celik, Muneer, & Clarke, 2008). Those studies which did focus primarily on sizing for

optimum storage capacity tended to focus on applications in which the load was a part of

a grid-tied system (Ru, Kleissel, & Martinez, 2013) or in combination with a number of

other methods of energy generation (Fathima & Palanisamy, 2015). These projects

generally assume the loads are variable and controllable with only minor penalties if the

9

available energy is insufficient. Because of our focus on a constant off-grid agricultural

load, where continuous operation was much more important than power supplied to the grid or carbon offset, this thesis research required a unique approach.

1.2.5 Backup Generator Systems

Off-grid solar PV systems are a promising technology for rural electrification, but studies have determined the economics of such systems are often not feasible on their own

(Baurzhan & Jenkins, 2016; Irfan, Zhao, Ahmad, & Rehman, 2019). Often, researchers examine such systems in combination with other renewable energy sources such as wind or hydroelectric sources (Krishna & Kumar, 2015). These hybrid renewable energy sources have been shown to effectively meet demand (Rahil, Gammon, & Brown, 2018) and to reduce total carbon emissions when compared to fossil fuel sources (Hossain, Mekhelif, &

Olatomiwa, 2017). However, when compared to a backup generator system, the renewable sources are still not as beneficial at maintaining a constant electrical load (Madziga, Rahil,

& Mansoor, 2018).

Backup generator systems alone are often used to provide a more dependable micro-grid for energy production (Dong, et al., 2018). Significant research has been done on combined systems using PV with backup generator power (Vivas, De las Heras, Segura,

& Andújar, 2018; Bajpai & Dash, 2012; Abedi, Alimardani, Gharehpetian, Riahy, &

Hosseinian, 2012), but most work is done considering variable loading instead of constant, critical loads (Muceka, Kukeera, Alokore, Noara, & Groh, 2018). Since the loads associated with IPR aquaculture are constant and critical, there exists a need to examine similar Hybrid PV-plus-Generator systems used for these purposes.

10

1.3 Objectives

The overall objective of this thesis was to examine the energy systems of floating In-pond

Raceways to determine whether these systems can be fitted with new methods of energy backup systems, including off-grid solar with and without backup generator systems. In addition, this thesis considers the addition of a new, variable load motor in relationship to these solar energy systems to further examine methods of reduction for daily energy demand.

Four primary objectives of this research are shown below and act as a basis followed heretofore by the ensuing thesis document:

1. Developing a model which accounts for energy flows in floating IPR systems

utilizing three different water moving devices: an airlift pump, propeller aerator

motor, and a slow-rotating paddlewheel.

2. Developing a model to provide size range estimates for off-grid solar PV plus

battery systems based on hourly climatological data at seven testing locations in

the state of Kentucky.

3. Examining the effect to sizes of off-grid PV plus battery systems following the

addition of backup generators for a range of system loads.

4. Implementing a full-scale variable speed propeller aerator system and determining

the effect such a system when used in place of a constant motor loading.

1.4 Thesis Outline

Chapter 1: Introduction. Contains an overview of floating IPR aquaculture, a review of literature regarding the modeling of energy flows in similar aquacultural systems, listed objectives of the thesis research, and outlines the thesis contents.

11

Chapter 2: Methodology. Details the modeling basis and framework for estimation of energy flows in floating IPR systems, protocol for sizing off-grid solar PV plus battery systems, outlines the process for estimating modified solar PV and battery sizes with the inclusion of generator backup charging, and conceptualizes the operation of a variable speed propeller aerator motor system.

Chapter 3: Results and Discussion. Summarizes the results for estimated motor flow efficiency, maximum fish weight supported, average sized solar PV plus battery systems both including and excluding additional generator backup charging and details the variable speed motor and solar energy systems as installed at the Frankfort, KY testing location.

Chapter 4: Conclusions and Future Work. Details the main findings of this research and provides examples of future work to build off these established concepts.

Chapter 5: References Cited. Lists all works cited within the thesis research.

12

CHAPTER 2. METHODOLOGY

2.1 Site Description

To meet the objectives of the thesis research, studies were done on three assembled

Superior Raceways Model 11000 Floating In-pond Raceway systems at the Kentucky State

University Harold R. Benson Research and Demonstration Farm facilities in Frankfort,

KY. Each raceway had a length of 40 ft, with a width of 8 ft and a depth of 6 ft at its deepest

point (Superior Raceway Systems, 2018). These raceways sit in a 1.6-acre pond at the research farm location (Figure 2.1) and recirculate water from this pond environment by means of water mover devices which are constantly in motion.

Figure 2.1 The site location at the Harold R. Benson Research and Demonstration Farm

in Frankfort, KY

2.2 Raceway Flow Data Collection

To monitor flow provided by each raceway device, researchers used a SonTek FlowTracker

2 Acoustic Doppler Velocimeter (ADV). Measurements were taken in discharge mode

within the raceway channels because the raceways were determined to be a continuous,

13

straight reach of uniform area with little impact from vegetation and obstructions (Corbett,

1943; Rantz, 1982). Flow velocity could not accurately be measured at the head of the

raceway due to the turbulent mixing of water near the flow generation device, and it also

could not be measured too near the end of the raceway close to a constriction in flow.

Because of this, measurements took place inside the raceway 24 feet from the raceway inlet

end-wall, as this was approximately 2/3 of the raceway length and the internal flow profile

was well developed by this point and not impacted by entrance and exit effects (Gordon,

McMahon, Finlayson, Gippel, & Nathan, 2004)

For the use of the ADV, there were established standards of measurement which were

adhered to (ISO, 2007; WMO, 2010). The width of the raceway was divided into a series of verticals at which depth and velocity were recorded. The optimal number of these vertical stations was 25, with 20 stations as a minimum (Rehmel, 2007). But because these stations must be a minimum of 5 cm apart from each other (SonTek, 2019) and the

raceways had a total width of only 8 feet, 20 stations were deemed to be sufficient for an

accurate discharge measurement. Additionally, two stations logged the raceway sides as

bank locations where no velocity measurements were taken.

To do this, a standard soft tapeline was installed across the raceway 24 feet (7.3 m) downstream from the end-wall closest to the water movement device. The first measurement was taken at the left wall of the raceway, with its location along the tapeline recorded as a bank measurement. From there, flow velocity measurements were taken every 4.75” (12.07 cm) horizontally, minus the reading at the point 47.5” from the left wall

of the raceway which was replaced with a reading at 48 in. To find average discharge, it is

typically measurements at 80% and 20% of depth which are necessary (Missouri State

14

University OEWRI, 2007). However, as the generated flow profiles were more turbulent

than flows typically measured by the FlowTracker2, it was recommended measurements

take place at 20%, 60%, and 80% of depth. At each station, the horizontal distance was recorded, and depth was measured by lowering the wading rod into the raceway until the probe contacted the raceway bottom and recording the depth measurement given by the metered tick-marks on the outer edge of the wading rod. Once the flow measurements at the various depths were recorded for each station and the entire width of the raceway was spanned, the location of the right-side wall was recorded, and testing was completed for that raceway.

Figure 2.2 Flow testing as completed within the Floating In-pond Raceway

The SonTek FlowTracker 2 ADV has a built-in system of data storage in which each

velocity measurement is sorted by horizontal location and vertical depth and then

15

integrated to find average discharge over the span. Each measurement of velocity was

broken into X and Y components, with X being downstream and Y being cross-sectional.

For this study, the most important measurement was the X component of velocity, as the downstream velocity was the component which generated water recirculation in the raceway environment.

2.3 System Performance Characteristics

A survey of previous work determined the three most plausible flow generation devices for use in IPR production systems. Airlifts are a traditional approach to providing flow for In-

Pond Raceways (Brown, Chappell, & Boyd, 2011; Parker & Suttle, 1987), with Slow-

Rotating Paddlewheels and Propeller Aerators gaining in prevalence in larger split pond systems and in select IPR environments (Brown, Tucker, & Rutland, 2016; Brune, Kirk, &

Eversole, 2004; Brune, Schwartz, Eversole, Collier, & Schwedler, 2003; Leavitt, et al.,

2008). An example of each of these devices were constructed for testing in this research.

The airlift was assembled by crew at the University of Kentucky Agricultural Machine

Research Lab and utilized a 1 Hp regenerative air pump. The propeller was a ¾ Hp unit purchased online from Kasco Marine, Inc (Kasco Marine, 2021). The Paddlewheel was a custom-built 1/15th Hp unit, rotating at 5.5 rotations per minute and consisting of six blades, each with a surface area of 5.72 square feet (Nguyen, Sines, Smith, & Whitlock, 2019).

Each device was installed at the head of its own respective raceway at the Frankfort, KY test location.

16

Figure 2.3 The three installed flow generation devices, clockwise from top: a 1 Hp airlift,

3/4 Hp Propeller aerator, and the 1/15 Hp Slow-Rotating Paddlewheel

Flow was then measured through the raceway for each device, alongside its electrical current draw. These flow measurements were taken on the SonTek FlowTracker 2

(SonTek/Xylem, Inc., San Diego, California, USA). Meanwhile, the current draw and power usage of the three devices was measured and a flow per watt of electrical power was calculated (Table 2.1).

17

Table 2.1 A summary of findings for device performance characterization.

Flow Power Flow

Produced Required Efficiency

(LPM) (Watts) (LPM Watt-1)

Airlift 10,900 1193 9.10

Propeller 11,000 800 13.8

Paddlewheel 4300 166 25.9

2.4 Modeling Parameters

2.4.1 Available Dissolved Oxygen

Boyd’s equation for nighttime dissolved oxygen consumption in catfish was used to predict the minimum required mass flow of oxygen through the raceway in terms of milligrams of oxygen per gram of fish per hour (Boyd, Romaire, & Johnston, Predicting Early Morning

Dissolved Oxygen Concentrations in Channel Catfish Ponds, 1978).

𝑂10... .. . (1)

Where:

W = Average mass per fish (g)

T = Water temperature (°C)

-1 -1 O = Oxygen demand (mgoxygen gfish hour )

∗∗ 𝑂 (2)

Where:

N = Number of fish in raceway

18

-1 Opopulation = Oxygen demand for entire fish population (mgoxygen min )

Boyd’s equations above are intended to solve for oxygen demand in milligrams per minute

within a raceway stocked with channel catfish. Channel catfish have been shown to be

tolerant of low oxygen environments in the past (Andrews, Murai, & Gibbons, 1973);

however, in most applications it is recommended that dissolved oxygen concentration

should remain above 3.5 mg l-1 for fish to feed well (Andrews, Murai, & Gibbons, 1973;

Andrews & Matsuda, 1975; Torrans, 2005), and above 2 mg l-1 for fish to survive (Lai-Fa

& Boyd, 1988). Therefore, the equation can be extended to solve for a required flow in

liters per minute which would provide for a desired outlet dissolved oxygen concentration falling within the recommended safe range (Eq. 3). To solve for this desired flow within the raceway, it was assumed that effluent dissolved oxygen was set constant at a minimum value of 3 mg l-1 so as to avoid unsafe conditions. Flow was determined at a range of inlet

dissolved oxygen contents (4 mg l-1, 5 mg l-1, 7 mg l-1, and 9 mg l-1), which represented the

general dissolved oxygen concentration levels in the pond. Although it was not considered

in this study, device aeration can be important, as it adds oxygen when water flowing into

the raceway may be oxygen depleted. These events may happen many times over the course

of a season, but their occurrence is often difficult to predict in intensive systems and may

depend more on complete nutritive budgets than stocked fish consumption (Boyd, Torrans,

& Tucker, 2018). For the purpose of energy flow modeling, it was assumed that

oxygenation by the water circulation device would be negligible.

Given these assumptions, values for required flow were produced and later used to estimate

a maximum carrying capacity in kg.

19

𝑄 (3)

Where:

Q = Required flow of water through the raceway (l min-1)

-1 Oin = Dissolved Oxygen measured at the raceway inlet (mg l )

-1 Oadded = Dissolved Oxygen added by the water moving device, 0 in this study (mg l )

-1 Oout = Desired Dissolved Oxygen at the raceway outlet (mg l )

2.4.2 Daily Energy Demand

Given device performance and required minimum flow from the above equations, it is

possible to estimate the minimum energy required to meet oxygen demand within the

raceway (Eq. 4).

𝐸 ∗ 24ℎ (4)

Where:

Q = Required flow (l min-1)

n = Device flow efficiency (l min-1 W-1)

E = Device Energy Demand (Wh)

Because there are grave effects on fish welfare once the water moving device ceases

operation, the Energy Demand found by the above operation must be applied as a constant

load. When sizing off-grid PV systems, it is typical to assume energy loads will be variable

and dependent on sophisticated mathematical modeling to determine the average daily loading (Mandelli, Merlo, & Colombo, 2016; Narayan, et al., 2020). Assuming energy

demand to be a constant load simplifies the process of determining average daily loading

and streamlines the sizing process for off-grid PV systems. On the other hand, this constant

load is also critical, and the system must ensure it can always provide that level of power. 20

2.4.3 Solar and Climate Data

Necessary to the project purposes was the daily weather data at the testing location. As the

main goal was to determine solar panel sizing for charging a battery backup system, the

main factors were knowing the solar irradiance at the location and the daily maximum air

temperatures. A survey of available online data found that the National Solar Radiation

Database (NSRDB) provided the most useful data source, providing accurate estimates for

solar energy historically provided at a given time and location anywhere in the United

States (NREL, 2019). This source provided data for both irradiance and meteorological

conditions by hour over the period from 1998 to 2019, primarily collected from airports.

For this study, we examined conditions from seven Automated Surface Observing System

(ASOS) sites located in the state of Kentucky: Lexington, Louisville, Jackson, Paducah,

Cincinnati/Covington, Bowling Green, and London.

Figure 2.2 A map of ASOS site locations from around the state of Kentucky utilized in this research.

21

2.4.4 Tools of the Trade: SAM and PVWatts

2.4.4.1 PV Watts

The National Renewable Energy Lab (NREL) has developed equations to relate available

daily solar radiation, nominal nameplate panel power ratings, and the total energy produced

by the system (Dobos, 2014). They embed these in their online solar performance

estimation tool, PVWatts, which can be used to estimate energy production for a given

solar power system for one hour of operation. (Dobos, 2014). This PVWatts model is not intended to perform discrete analysis of an individual cell within a solar array, meaning it does not break the total nameplate power into a set of power values for an integer number of cells. Instead, the nameplate power for the entire solar array is assumed to be the total maximum power output. This value is then subjected to a default efficiency value of 96% to predict power supplied from the inverter on the solar collector unit.

𝑃 𝑃 1 𝛾𝑇 𝑇 (5)

Where:

Pdc = DC power supplied by PV array (kW)

-2 Itr = Plane of Array Irradiance (W m )

Pdc0 = PV nameplate capacity (kW)

γ = Temperature coefficient (% °C-1)

Tcell = Temperature of PV cell at a specified location (°C)

Tref = Standard reference temperature (25°C)

22

2.4.4.2 SAM

The NREL has also developed a techno-economic model called the System Advisor Model

(SAM) to assist in decision making for work done with renewable energy systems (Blair,

et al., 2018). This model is typically used to assist with the sizing and development of grid-

tied energy systems and has a variety of photovoltaic energy models embedded, though it

can also be used to estimate system battery charging for off-grid systems. This study used

the PVWatts model for Distributed Residential owners (Eq. 5) embedded within SAM to

estimate solar power production per hour of a variety of system sizes.

To do this, hourly weather and irradiance data was imported from the NSRDB from each of the selected Kentucky ASOS station sites over the period from the start of January 1998

to the end of December 2019. These hourly weather files also imported geographic

information important to the estimation of solar power output, including the site location’s

latitude, which was also used in the system sizing tab of the model inputs page. For each

weather file, the hourly solar power production was estimated for systems of the following

sizes: 2 kW, 3 kW, 4 kW, 5 kW, 10 kW, 15kW, and 20 kW. Each sized system was assumed

to be oriented due south, with a tilt equal to the latitude of that location. SAM default values

were used for other system parameters.

2.4.5 Battery System Sizing and Efficiency

The energy required to support electric loading must be produced from a battery pack in

our application, making the PV system an off-grid or non-grid-tied system. Battery packs

will have inefficiencies, which depend on weather and battery composition. Lead-acid and

lithium-ion batteries differ in their various inefficiencies due to differences in material

properties and temperature factors (Chang, 2019). Additionally, there will be inefficiencies

23

in charging these batteries that result from the battery’s inability to store all the provided energy, the inability of the charge controller to provide all available energy to the batteries

for them to store, and the solar inverter’s inability to convert all available energy into useful

energy to operate loads.

Typically, charging efficiency depends on the battery’s State of Charge (SOC). For a lead-

acid battery charging is 91% efficient under 84% SOC. However, this efficiency drops to

55% with higher SOC (Stevens & Corey, 1996). Our model is concerned with complete

battery depletion events rather than precisely capturing efficiency, so this variation based

on SOC was not critical, and we focused on the overall charging efficiency across the entire

range. To accommodate for the SOC effect on charge efficiency, the weighted average over

all SOC values, which was 85%, was used as the battery charging efficiency. Standard lead

acid batteries suffer considerable loss of operating life when the state of charge falls below

20% and battery manufacturers advise not discharging beyond this point (Trojan Battery

Company, 2019). Therefore, our model assumed the minimum acceptable state of charge

was 20%. High discharge rates also negatively affect the capacity of a battery and can result in obtaining less energy than expected. However, the battery in this system must be able to

continuously discharge overnight at a constant rate. This requires a ratio of battery capacity and loads that provides a discharge time of at least eight hours, and the acceptable systems displayed in the results all have discharge times of at least 60 hours. With these discharge rates, battery capacity is stable and not dependent on discharge rate (EnerSys, 2016), so

the effect of discharge rate on capacity was not included in the model.

24

2.5 Iterative System Sizing Models

2.5.1 Off-Grid Solar PV Plus Battery (No Generator)

Once the system sizes and efficiencies were estimated, an energy analysis function was

formulated to simulate the simultaneous charging and discharging of the battery to

determine if the energy requirements would be met. An hourly energy balance was created

using the SAM tool estimates to determine energy produced by the system and combined

with required energy loads and battery storage capacity at the hour of measurement (Eq.

6).

⎧ 𝐸 𝐸 ∗𝜂 ∗𝜂 , 𝑖𝑓 𝐸 ∗𝜂 ⎪ 𝐸 (6) ⎨ ⎪ 𝐸 𝐸 ∗𝜂 , 𝑖𝑓 𝐸 ∗𝜂 ⎩

Where:

Ei = Initial battery storage (kWh)

Eo = Remaining battery storage (kWh)

Ed= Energy demand required to operate the flow device for one hour (kWh)

Es = Energy supplied by the solar PV system (kWh)

ηch = Charge Controller Efficiency

ηi = Inverter Efficiency

ηbatt = Battery Storage Efficiency, only affects charging into the battery

Efficiency values for the charge controller, ηch, inverter, ηi, and battery storage, ηbatt, were set at 0.97, 0.93, and 0.85, respectively. These were obtained from the specification sheets of commercially available components in this size range – a Cotek SP1000-148 inverter

(Cotek Electronic Industries, Taoyuan City, Taiwan), a FLEXmax60 charge controller

25

(Outback Power, Phoenix, Arizona, United States) and lead-acid batteries (discussed in

section 2.7). This energy balance was used with the historical solar radiation and

temperature data to estimate the battery’s SOC on an hourly basis. If at any point the battery reached a SOC of 20% or lower, it was recorded as an energy depletion event. These events were summed up over the dataset to find the total estimated number of energy depletion events for each system. A successful system was defined as one with no events in which battery capacity ever fell below a 20% state of charge in the 22-year dataset.

The testing process compared system outputs against continuous system loads of 0.25 kW,

0.5 kW, 0.75 kW, and 1 kW, each value representing a corresponding unique system, which were substituted into Equation 6 above as the required hourly energy demand. These loads were applied to each of the sized systems from SAM, each system in conjunction with battery storage capacities of 5 kWh, 10 kWh, 15 kWh, 20 kWh, 30 kWh, 40 kWh, 50 kWh, and 60 kWh.

The worst-case scenario for energy depletion was winter operation with short, cloudy days followed by snowfall that temporarily covers the panel. Because of this, two models were produced: one to account for year-round operation, and one which ceases operation in the winter months (not operating from November to February). Given that the associated risks of energy depletion were so high for IPR aquaculture, this seasonal model was chosen to

be most representative of likely operations, though year-round operation is still a

possibility.

2.5.2 Off-Grid Solar PV Plus Battery (With Generator)

A supplemental energy analysis function, based on prior energy analysis for solely off-grid

systems, was formulated to simulate the effect of charging by a small generator system

26

once SOC fell below a value of 20%. For this analysis, a generator capable of 2 kW output

was used as a reference. It was assumed to shut off once SOC had exceeded a value of

85%. An hourly energy balance was created using the SAM tool estimates to determine

energy produced by the system in tandem with generator charging and combined with

required energy loads and battery storage capacity at the hour of measurement (Eq. 7).

⎧ 𝐸 𝐸 ∗𝜂 𝐸 ∗𝜂 , 𝑖𝑓 𝐸 ∗𝜂 ⎪ 𝐸 (7) ⎨ ⎪ 𝐸 𝐸 ∗𝜂 𝐸 , 𝑖𝑓 𝐸 ∗𝜂 ⎩

Where:

Ei = Initial battery storage (kWh)

Eo = Remaining battery storage (kWh)

Ed= Energy demand required to operate the flow device for one hour (kWh)

Es = Energy supplied by the solar PV system (kWh)

Egen= Energy supplied by supplemental generator system

ηch = Charge Controller Efficiency

ηi = Inverter Efficiency

ηbatt = Battery Storage Efficiency, only affects charging into the battery

Efficiency values for the charge controller, ηch, inverter, ηi, and battery storage, ηbatt, were again set at 0.97, 0.93, and 0.85, respectively. These were obtained from the specification sheets of commercially available components in this size range – a Cotek SP1000-148 inverter (Cotek Electronic Industries, Taoyuan City, Taiwan), a FLEXmax60 charge controller (Outback Power, Phoenix, Arizona, United States) and lead-acid batteries.

This energy balance was used with the historical solar radiation and temperature data to estimate the battery’s SOC on an hourly basis. If at any point the battery reached a SOC of 27

20% or lower, the generator output was added to the system to simulate supplemental

backup charging. The total number of these supplemental charging events were summed

up over the dataset to find estimate the number of hours of operation for the 2kW generator

system over the 22-year time period. Using, hours of operation, the fuel consumption of

prospective generator systems was determined to be 0.297 gallons per watt-hour for

gasoline (Honda Motor Company, 2019) and 0.164 gallons per watt-hour for propane

(Ericson & Olis, 2019).

This approach was deemed satisfactory to estimate system sizing based on generator

supplied energy. However, at the worst-case loading, it was known that generator fuel

consumption would not hold at a constant value, and therefore, a supplemental analysis

was necessary to determine the variable fuel consumption at various supported loads. Using

reported fuel consumption values for a Honda EU2200i generator modified to run on

propane (Genconnex, 2019), a linear regression was performed to generate an equation of

best fit to estimate propane consumption values in liters hour-1 for loads of 0.25kW, 0.5kW,

0.75kW, and 1kW.

𝑄0.7012𝑃0.3303 (8)

Where:

Q = Fuel Consumption Rate (Liters Hour-1)

P = Generator Load (kW)

As with the prior analysis, system loads of 0.25 kW, 0.5 kW, 0.75 kW, and 1 kW were utilized. Battery storage capacities that were modeled were 10 kWh, 15 kWh, 20 kWh, 30 kWh, 40 kWh, and 50 kWh, and 60 kWh. These capacities ensured battery capacity was

28

large enough to provide at least 10 hours of discharge at even the largest load level.

Discharge levels that are higher than one tenth of capacity can have a detrimental effect on

battery capacity and battery life, so the 10kWh limit ensured we were not at risk of lower

capacity levels and damaging battery life (EnerSys, 2016). This analysis was performed to

find the total number of hours the generator backup was in operation over the 22-year data set at the seven testing locations.

2.6 Conceptualizing a Variable Speed Propeller System

Prior work in this research focused on constant motor loads to move water within the floating IPR systems. However, loads used in off-grid solar energy systems are often variable (Mandelli, Merlo, & Colombo, 2016; Narayan, et al., 2020). The constant load assumption will cause the off-grid solar energy systems to be significantly larger than their sustained loads, as the constant system demands cannot be easily supported in the winter while solar energy production is at its lowest. This effect was theorized to be alleviated minorly by seasonal operation, and more so by the addition of a generator power backup system, but there does exist a third option: to lower energy demand by utilizing a variable load system based on system oxygen demand.

Presently, the designed flow generation systems provide constant flow of water through the raceways year-round. However, in the winter months, dissolved oxygen concentrations are higher (Ouyang, Nkedi-Kizza, Wu, Shinde, & Huang, 2006) and metabolic rates for both fish and microbial communities will decrease. Thus, the need for high rates of flow may be reduced. This coincides with the time of lessened solar potential in the state of Kentucky, meaning demand may be lowered at the time in which

29

the least energy is available. In typical backup systems, batteries and panels must be

oversized to provide storage to allow operation in these months. However, reducing loads

in relation to need will allow for more efficient operation and possibly decrease off-grid

solar system sizes further.

Work has been done with propeller devices in floating upwellers supporting shellfish seed in marine environments (Rivara, Tetrault, & Patricio, 2002), including variable speed systems, though they still depend on manual alteration of flow values (Bearon

Manufacturing, 2021). This introduces potential human error to the system and does not guarantee system operation based on need. The task for this phase of the project was to

develop an automated system which adjusts motor speeds based on measured dissolved

oxygen values within the raceway environment. These values must be used to scale motor

speeds with the intent of producing only enough flow to meet demand, allowing for more energy efficient operation of water moving devices without requiring manual operation

by producers.

30

CHAPTER 3. RESULTS AND DISCUSSION

3.1 IPR System Sizing Tool

Of the final tools produced for modeling these systems, one was a decision support

application for producers built with MATLAB R2020a’s App Designer studio and

packaged for individual download. This decision support tool combines all of the previous

work done to size IPR systems into an easy to navigate dashboard with inputs of system device type, estimated flow within the raceway, and costs of fuel associated with backup systems if generators were used, and outputs the minimum required system motor power in Horsepower, minimum generator sizing in Watts with related fuel consumption per day for both gasoline (Honda Motor Company, 2019) and propane (Ericson & Olis, 2019), and the corresponding estimates for carrying capacity at varying levels of pond dissolved oxygen concentration (Fig. 3.1).

Figure 3.1 The application to size In-pond Raceway production systems

31

3.2 Estimating Fish Carrying Capacity

Using Boyd’s equation, it was possible to predict both a required flow generated by the water circulation device and carrying capacity for channel catfish at the estimate ranges of dissolved oxygen content. Through site analysis at the Frankfort, KY testing location, it was determined 4 mg l-1 would be the average minimum condition for pond dissolved

oxygen content. The following table summarizes expected flow in liters per minute as well

as carrying capacity in kilograms for each water circulation device type at the specified

electric loads in kW. Estimated weight of fish supported was given assuming base loading

at an average water temperature of 12.9°C, though the values of estimated fish supported

may be within a margin of ±10% given meteorological data for the tested locations and

time of year.

Table 3.1 Carrying capacity for channel catfish (kg) in a raceway assuming DO concentrations of 4 mg l-1 or higher in the pond with raceway effluent DO of 3 mg l-1 at 12.9°C.

Device Type Load (kW) Paddlewheel Propeller Airlift Flow (lpm) 4900 2600 1900 0.25 Fish Supported (kg) 2300 1200 880 Flow (lpm) 9500 4900 3400 0.50 Fish Supported (kg) 4400 2300 1600 Flow (lpm) 14000 7600 4900 0.75 Fish Supported (kg) 6700 3500 2300 Flow (lpm) 19000 10000 6800 1.00 Fish Supported (kg) 8900 4700 3200

3.3 Off-Grid Solar Recommendations

Economics of Solar PV backup are not stationary, with prices falling every year (Lorenz,

Pinner, & Seitz, 2008; Feldman, et al., 2012). Even so, costs of such systems can be prohibitively high. In supporting a constant load, there is a tradeoff between panel array size and battery capacity. Increasing panel array size or increasing battery capacity will

32

both increase the likelihood that a system would be able to support the continuous load during a temporary period of lower solar radiation. To account for variation in economics between the costs of batteries and panels, tables were provided with multiple recommendations at each load level that illustrate this tradeoff. These tables, produced for

systems both with and without additional generator charging and at seven different

locations, are shown in the following sections.

3.3.1 Off-Grid Solar PV Plus Battery (No Generator)

With the panel (up to 20kW) and battery (up to 60 kWh) sizes tested, it was not possible to

support constant loads greater than or equal to 1 kW. Though seasonal operation would

still not support 1 kW loads, the necessary system sizes to support loads between 0.25 and

0.75 kW were greatly reduced when the system was only operated between March and

October. The following tables display several options for the smallest acceptable system

size for many loads at several locations. In these cases, the acceptable panel size could be

decreased by increasing the battery capacity. Given the rapid changes in battery and solar

component prices as well as location characteristics such as limited space for solar panels

or protected storage for batteries, it is not possible to tell producers which combination of

panel and battery size would definitively be best for their installation.

Table 3.2 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Lexington, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 10 40 3 40 0.25 15 30 0.25 4 30 20 20 10 15 ** ** 10 50 0.5 0.5 15 40 ** ** 20 60 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

33

Table 3.3 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Louisville, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 3 50 0.25 ** ** 0.25 4 40 10 20 ** ** 15 60 0.5 0.5 20 40 ** ** ** ** 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

Table 3.4 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Jackson, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 10 40 3 40 0.25 20 50 0.25 4 30 10 15 ** ** 10 50 0.5 0.5 15 40 ** ** ** ** 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

Table 3.5 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Paducah, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 5 60 3 50 0.25 10 40 0.25 4 30 15 30 10 20 10 60 0.5 ** ** 0.5 15 50 20 40 ** ** ** ** 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

34

Table 3.6 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Cincinnati, OH/Covington, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 3 50 ** ** 4 40 0.25 0.25 5 30 10 20 ** ** 10 60 0.5 0.5 15 40 ** ** 20 60 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

Table 3.7 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in Bowling Green, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 3 40 0.25 ** ** 0.25 5 30 10 20 10 60 0.5 ** ** 0.5 15 50 20 40 ** ** ** ** 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

Table 3.8 Smallest acceptable system sizes for both seasonally (March to October) and Year-Round loads in London, KY. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 10 60 3 40 0.25 15 50 0.25 4 30 20 40 10 15 10 50 0.5 ** ** 0.5 15 40 20 30 ** ** 20 60 0.75 0.75

**Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

35

The results of this study show that there is variation in the sizes of recommended systems

based on both the geographic location within the state and whether the systems were run year-round or seasonally (from March to October). Further, the recommended system sizes

in all locations were large in comparison to the loads they were able to sustain. Costs of

such systems will be high, showing that seasonal operation may be preferable for aspiring

producers. In many of the cases, it was not possible to identify an acceptable system that required a 20 kW or less solar panel array or 60 kWh or less of battery capacity. In several locations, seasonal operation was the only way non-grid-tied solar energy systems would successfully sustain even a 0.25 kW load.

With the suggested systems, the battery was almost always operating at a high state of charge (Figure 3.2). There were periodic decreases in SOC, and these occurred on a yearly cycle during the winter months. However, annual decreases often only lowered values to a minimum of 75% SOC. It was rarely when a certain set of weather and light conditions created a sustained depletion, and the battery approached its minimum level.

36

Figure 3.2 Available battery capacity for a system utilizing 10kW Panels and a 40 kWh

Battery to power a 0.25 kW load at the Lexington site location.

The battery state of charge distribution (Figure 3.3) also confirmed significant depletion events were exceedingly rare. Out of a total of 192,720 hours, 188,640 hours or ~ 97% had a battery state of over 85%. This pattern of shallow discharges should provide excellent battery life as depth of discharge is negatively correlated with the number of cycles in lead- acid battery life. Only 124 hours or 0.06% were at a battery SOC of less than 50%. The full battery capacity was only needed on a few occasions during this 22-year weather dataset, but such a large battery pack was required to ensure operation during these extreme events. Long-term effects of a changing climate are unknown and were not considered in this model. However, the efficacy of solar panel PV systems is known to decrease with

rising temperatures (Carbon Disclosure Project, 2009) and a shift in global climate would cause variability in expected solar irradiance (Solaun & Cerdá, 2019), meaning extreme weather events may be more common in the future and would further necessitate the use of the large battery packs recommended by the model equation.

37

Figure 3.3. Number of hours at different battery state of charge ((a) 0% to 50%, (b) 40% to 80%, and (c) 70% to 100%) over the 22-year dataset for the 10kW Panel and 40 kWh

Battery while supporting a 0.25 kW load in Lexington, KY.

3.3.2 Off-Grid Solar PV Plus Battery (With Generator)

The addition of backup generator systems greatly altered the feasibility of the off-grid PV systems. A passable system was defined as systems which use under 100L per year of a defined fuel source, in this case propane, system recommendation sizes greatly decreased for every monitored location. Likewise, every monitored location was found to be able to support loads up to 1kW with such systems and seasonally operated loads of 0.25kW could be sustained by any tested combination of PV-Plus-Generator configurations at all sites.

This is illustrated in the tables below.

38

Table 3.9 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Lexington, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 30 0.25 3 15 0.25 2 10 4 10 5 60 4 30 0.5 10 30 0.5 5 20 15 15 10 10 10 60 10 30 0.75 15 40 0.75 15 20 20 30 20 15 20 60 10 50 1.0 1.0 15 40 20 30

Table 3.10 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Louisville, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 30 0.25 3 15 0.25 2 10 5 10 5 60 4 30 0.5 10 30 0.5 5 20 15 15 10 10 15 40 10 30 0.75 20 0.75 15 20 30 20 15 20 60 10 50 1.0 1.0 15 40 20 30

Table 3.11 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Jackson, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 50 0.25 3 15 0.25 2 10 5 10 10 30 4 30 0.5 15 20 0.5 5 20 20 10 10 10 15 40 10 30 0.75 0.75 15 20 20 30 20 15 20 60 10 60 1.0 1.0 15 40 20 30

39

Table 3.12 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Paducah, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 30 0.25 3 15 0.25 2 10 5 10 5 50 4 30 0.5 10 30 0.5 5 20 15 15 10 10 10 60 10 30 0.75 15 40 0.75 15 15 20 30 20 10 20 60 10 50 1.0 1.0 15 30 20 20

Table 3.13 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Cincinnati, OH/Covington, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 60 0.25 3 15 0.25 2 10 5 10 10 30 4 40 0.5 15 20 0.5 5 20 20 15 10 10 15 40 10 30 0.75 0.75 15 20 20 30 20 15 10 60 1.0 ** ** 1.0 15 40 20 30 **Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

Table 3.14 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in Bowling Green, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 30 0.25 3 15 0.25 2 10 5 10 5 60 4 30 0.5 10 30 0.5 5 20 15 15 10 10 15 40 10 30 0.75 0.75 15 20 20 30 20 15 20 60 10 50 1.0 1.0 15 30 20 20

40

Table 3.15 Smallest acceptable system sizes for Solar PV plus Battery Systems with Backup Generators in London, KY which use Fewer Than 100L per Year of Propane. Year Round Seasonal (March to October) System Size System Size Battery Battery Load Panel Size Capacity Load Panel Size Capacity (kW) (kW) (kWh) (kW) (kW) (kWh) 2 30 0.25 3 15 0.25 2 10 5 10 5 60 4 30 0.5 10 30 0.5 5 20 15 15 10 10 10 60 10 30 0.75 15 40 0.75 15 15 20 30 20 10 20 60 10 50 1.0 1.0 15 40 20 30

The benefits of incorporating propane generator backup systems are numerous for potential producers. It allows for the demonstrated decrease in necessary PV plus battery system size, which will decrease the cost of such systems and stands to improve their levels of implementation, shown in table 3.16 below. They also lower the risk of a fish-kill event during periods of extreme weather and provide security for the investments of farmers.

They even have reduced emissions when compared to gasoline generators of a similar size, with 100L of gasoline producing 36% more CO2 into the atmosphere than the same quantity of propane (U.S. EIA, 17).

Table 3.16 Comparison of system sizes for Solar PV plus Battery Systems with and without Backup Generators in Lexington, KY. Year Round Seasonal (March to October) Without Generator With Generator Without Generator With Generator Battery Battery Battery Battery Load Panel Size Capacity Panel Size Capacity Load Panel Size Capacity Panel Size Capacity (kW) (kW) (kWh) (kW) (kWh) (kW) (kW) (kWh) (kW) (kWh) 10 40 2 30 3 40 0.25 15 30 3 15 0.25 4 30 2 10 20 20 5 10 10 15 5 60 10 50 4 30 0.5 ** ** 10 30 0.5 15 40 5 20 15 15 10 10 10 60 20 60 10 30 0.75 ** ** 15 40 0.75 15 20 20 30 20 15 20 60 10 50 1.0 ** ** 1.0 ** ** 15 40 20 30 **Acceptable system would require panel size > 20 kW or battery capacity > 60 kWh

41

Another benefit is the expansion of choice granted to a producer by the generator backup

system. Presently, systems which run the specified load while utilizing under 100L of

propane per year are selected as the recommended system sizes. There are other systems,

however, which may utilize slightly over 100L per year of fuel but greatly decrease the

solar PV plus battery system size required to sustain the critical motor load, shown in

Figures 3.4 through 3.11 below. Based on these findings, producers may make the decision

to build systems slightly outside of the recommended table of sizes and fall within a

reasonable level of fuel usage per year. However, producers wishing to size systems which

fall outside of the recommended ranges listed above should exhibit restraint. This is because systems which are significantly undersized may consume levels of fuel which are

unsustainable over year-round operation, with a worst-case operation mimicking that of a

system with 0kW panels and a 0kWh battery. This condition is shown in the table below.

Table 3.17 Estimated worst-case generator propane consumption values by Equation 8. These values were found by assuming a worst-case operation of 8760 operation hours for year-round operations and 5880 hours for seasonal operation.

Year-Round Seasonal Fuel Load Fuel Consumption Fuel Consumption Consumption (kW) (L/hr) (L/yr) (L/yr)

0.25 0.51 4429 2973

0.50 0.68 5965 4004

0.75 0.86 7500 5034

1.0 1.03 9036 6065

42

Figure 3.4 System fuel use trends for a 0.25kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries.

Figure 3.5 System fuel use trends for a 0.25kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario

of no Solar Panels and no batteries.

43

Figure 3.6 System fuel use trends for a 0.5kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries.

Figure 3.7 System fuel use trends for a 0.5kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario

of no Solar Panels and no batteries.

44

Figure 3.8 System fuel use trends for a 0.75kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries.

Figure 3.9 System fuel use trends for a 0.75kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario

of no Solar Panels and no batteries.

45

Figure 3.10 System fuel use trends for a 1-kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), including the worst-case fuel use scenario of no Solar Panels and no batteries.

Figure 3.11 System fuel use trends for a 1-kW load in Lexington, KY for Year-Round

Operation (left) and Seasonal Operation (right), without the worst-case fuel use scenario of no Solar Panels and no batteries.

The above figures show the trends in fuel use for a propane backup generator system in combination with an off-grid PV plus battery system. Areas of the diagram which are gray are areas in which zero liters per year of fuel are used by the backup generators in Figures

3.5, 3.7, 3.9, and 3.11, and sizes estimated to be near-zero in Figures 3.4, 3.6, 3.8, and 3.10.

Fuel use in these figures is shown in hues of dark blue where fuel use is lowest, and bright

46 yellow where fuel use is highest. As a general trend, smaller loads may subsist on lesser fuel usage, with certain configurations running loads of up to 0.75kW able to be sustained with no generator backup at all. As battery and panel size combinations decrease in size, fuel usage is expected to increase rapidly toward the maximum fuel usage expected in systems running entirely off of generator power. Thus, producers should take caution when attempting to size smaller systems than those recommended in the presented tables.

3.4 Implemented Variable Speed Propeller System

Work done for the variable speed propeller system in this research focused on design and implementation of such a system. To achieve this, a Power House Vari-Speed propeller motor with a maximum power rating of 1 Hp was used as the water moving device at the head of Raceway 2 at the site location. This motor has been used in similar operations for shellfish farming (Bearon Manufacturing, 2021) but depended on manual inputs from producers by use of a designed control panel. For the purposes of the research project, this system must autonomously accept dissolved oxygen sensor measurements from the existing AM-2300 sensor system.

This was done by incorporating a Click Programmable Logic Controller (PLC) device in series with the AM-2300 controller to accept DO sensor inputs from the raceway environment. These sensor values were then used as system inputs for the built-in PID loop control function. The PID loop control in Click PLC uses its own autotuning operations to determine the optimum values for PID variables to scale an output, an output in this case that was motor input frequency. This input frequency was carried by RS-485 communication to a WEG CFW300 Variable Frequency Drive (VFD), which conveyed this frequency value to the Vari-Speed motor. Because the AM-2300 continuously

47

monitored the DO concentrations in the raceway environment, this process was repeated

continuously with the intent of holding the DO concentration at the raceway exit that was

as near to the set point of 7.5 mg/l as possible. This program was written in the Click PLC

Programming Software window.

Figure 3.12 The industrial control box as installed in Frankfort, KY

48

This research ran into significant blocks in implementation which had to be overcome before the system could be considered operational. The first major issue was that the electrical system available at the site location produced too much noise in the input voltage signal for the WEG CFW300 VFD to operate properly. This caused the VFD to sense a ground fault error and shut down power supplied to the Vari-Speed motor. Upon further examination, it was determined that an isolation transformer would be required to reduce noise from the input signal, which was alleviated by the inclusion of an Emerson SolaHD

1.5 kVA isolating transformer.

Figure 3.13 The industrial control system as installed, including all controls and the

1.5kVA transformer.

49

There were also significant challenges in establishing RS-485 Communication between the

PLC and VFD device. Communication between the devices was established using

MODBUS protocol, a device address of 1, and the offset value of 400,000 as mandated by the Click PLC communications help file (CLICK PLCs, 2021). Parameters 222 through

228 on the WEG VFD were set to values which allowed for remote operation by serial inputs to establish which variables would be altered by the Click PLC serial

communications. This operation was tested by the implementation of a rocker switch on

VFD digital input 4 to alternate between local and remote operation. If the device

successfully switched from local to remote control, parameter 680 would change value, showing remote control was enabled. Serial port status could likewise be determined by monitoring Parameter 316 on the VFD. A test code was developed in the Click PLC

Programming window to determine functionality of the RS-485 communications.

50

CHAPTER 4. CONCLUSIONS AND FUTURE WORK

4.1 Conclusions

In-pond Raceways require a continuous supply of power to operate water moving devices,

which are necessary to provide dissolved oxygen for fish production. The power

requirements varied based on the water moving device. Slow moving paddle wheels were

the most efficient and the traditional airlifts were the least efficient in terms of amount of

flow provided per unit of power required. However, it is important to note that the

traditional airlift unit provides additional aeration into the raceway environment, which

may be beneficial in more oxygen depleted environments.

Testing for recommended energy systems determined the required solar system sizes to

support off-grid solar power to operate an IPR all-year were very large for even small

continuous electrical loads. Since major energy depletion events were clustered during winter months, the system sizes are smaller and larger loads can be considered if raceways

are only operated between March and October. Overall, the batteries experience very low

depths of discharge in the average cycle; therefore, battery life should be excellent in these

systems.

The inclusion of a generator backup system improved the relationship between necessary

sized solar energy systems and sustained loads. Systems up to 1kW were supported on

panels under 20kW with batteries under 60kWh while consuming under 100L of propane

per year for generator backup, making the hybrid Off-Grid Solar and Generator energy

systems more feasible for producers in aquaculture than the previously sized solely Off-

Grid Solar PV plus Battery systems. This addresses the issue of seasonal interruption and

51

eliminated the need to drastically oversize a solar PV plus battery system to accommodate even small loads.

4.2 Future Work

This work provides a new view of energy flows within floating IPR systems; however, more work is needed to greater improve feasibility for small-scale producers. Off-grid solar

systems can be improved by operating seasonally or in tandem with generator backup.

Likewise, it was conceptualized a variable speed load based on demand for oxygen in the

raceway environment could reduce energy demand overall, thus decreasing necessary

backup system size. Though it has been theorized, the efficacy of such a system is not yet

known, especially in comparison to existing and modelled system variations.

Presently, the variable speed propeller system modifies motor speed based on an input of dissolved oxygen concentration at the raceway’s exit. Before testing the control systems for operation, the flow rate and energy demand should be measured with the motor operating at different levels. To test for system performance, it is necessary to design an experiment which would alter dissolved oxygen concentrations in the raceway to mimic the behavior of fish and microbial consumption of algae occurring within the water system.

Fish and microbial respiration, as well as the growth and decomposition of algae, will deplete usable oxygen in the raceways as they respire, but these factors can be cumbersome to monitor in a testing environment and pose significant risk for loss of fish if the system is not appropriately calibrated. Instead, a designed experiment should utilize sodium sulfite as a method to remove dissolved oxygen from the raceway, similar to the method devised to test for aeration potential of different devices for aerating pond systems used in rearing of fish (Boyd C. E., 2015). This test can be done by mixing concentrations of sodium sulfite

52 and feeding them into the raceway with the variable speed propeller motor to monitor system response.

Control systems often exhibit behavior where there exists a high number of oscillations

(Thornhill & Hägglund, 1997). One of the main causes for this oscillation is lag in the system (Ang, Chong, & Li, 2005) – the system adjusts the input, but the output does not change in a timely manner. If speed increases to the maximum or drops to the minimum, it will take the time required for the water to circulate before one can see the final system response. The main concern of the testing with concentrations of sodium sulfite would be to stabilize the oscillations of such a system in response to a spike in oxygen demand, such as after a heavy feeding in a fully stocked raceway. Researchers can intentionally create the worst-case lag scenario by moving the DO probes from normal water to the O-depleted water rapidly and waiting to switch the sensor back into the normal water after an elapsed lag time as found in the earlier experiment. It must stabilize after this shock to be considered fully operational.

Once these tests for operation are completed, it would be necessary to run the configured system over an entire growing season while also running a separate raceway with a similar propeller motor system comparable to what a producer would be able to buy off market and install in their own systems. This comparison would show whether or not there is a marked difference in energy demand to the system when a variable speed propeller motor is used. An economic analysis of all system variations for off-grid solar PV plus battery backup - seasonal operation, including generator backup, and variable speed motor systems

- could then be completed to determine which, if any, of the systems pose the greatest benefit to small producers in terms of cost benefit.

53

APPENDICES

[APPENDIX 1. IPR DECISION SUPPORT TOOL] classdef Solar_Tool < matlab.apps.AppBase

% Properties that correspond to app components properties (Access = public) UIFigure matlab.ui.Figure PriceperkWhforElectricityinYourAreaEditField matlab.ui.control.NumericEditField PriceperkWhforElectricityinYourAreaLabel matlab.ui.control.Label PriceofPropaneinYourAreaEditField matlab.ui.control.NumericEditField PriceofPropaneinYourAreaEditFieldLabel matlab.ui.control.Label PriceofGasinYourAreaEditField matlab.ui.control.NumericEditField PriceofGasinYourAreaEditFieldLabel matlab.ui.control.Label RacewayFlowGPMEditField matlab.ui.control.NumericEditField RacewayFlowGPMEditFieldLabel matlab.ui.control.Label SystemTypeListBox matlab.ui.control.ListBox SystemTypeListBoxLabel matlab.ui.control.Label WeightSupportedLabel matlab.ui.control.Label PondDOLabel matlab.ui.control.Label MotorHp matlab.ui.control.Label MinMotorSizeHpLabel matlab.ui.control.Label Pleasenoteadditionalaerationmayberequiredunder3mglLabel matlab.ui.control.Label kWhCostLabel matlab.ui.control.Label kWh matlab.ui.control.Label PropCostLabel matlab.ui.control.Label GasCostLabel matlab.ui.control.Label DailyCostofElectricityLabel matlab.ui.control.Label kWhofElectricityUsedperDayLabel matlab.ui.control.Label DailyCostofPropaneLabel matlab.ui.control.Label DailyCostofGasolineLabel matlab.ui.control.Label lbsLabel_4 matlab.ui.control.Label lbsLabel_3 matlab.ui.control.Label lbsLabel_2 matlab.ui.control.Label lbsLabel matlab.ui.control.Label mglLabel_4 matlab.ui.control.Label mglLabel_3 matlab.ui.control.Label mglLabel_2 matlab.ui.control.Label mglLabel matlab.ui.control.Label

54

EstimatedWeightofCatfishSupportedLabel matlab.ui.control.Label TogetstartedwiththisapplicationfirstseLabel matlab.ui.control.Label PropaneUsed matlab.ui.control.Label GallonsofPropaneUsedperDayLabel matlab.ui.control.Label SystemSizingLabel matlab.ui.control.Label Gasoline matlab.ui.control.Label Generator matlab.ui.control.Label GallonsofGasUsedperDayLabel matlab.ui.control.Label MinGeneratorSizeWLabel matlab.ui.control.Label InPondRacewaySystemSizingToolLabel matlab.ui.control.Label end

properties (Access=private) T W N9 N7 N5 N3 system hr Q I V FileName num SolarI Tvar P P1 n Qr Preq Preq1 AvgS AvgT Gas Propane Value Value1 Value2 GasPrice PropPrice ElectricPrice

55

end methods (Access = private) end

% Callbacks that handle component events methods (Access = private)

% Code that executes after component creation function startupFcn(app) app.T=12.9; %Average temperature at Frankfort location (in Celsius) app.W=454; %Assumption for average individual fish weight in grams app.I=1.36; app.Q=1141.3*3.785; app.V=122.4; %Measured constant voltage app.P=app.I*app.V; %Power required for device app.P1=(app.P*1.341)/1000; %Converted from kW to Hp app.n=(app.Q)/app.P; %Effective efficiency for device in lpm/kW app.GasPrice=2.08; app.PropPrice=2.07; app.ElectricPrice=0.08; app.hr=24; end

% Value changed function: SystemTypeListBox function SystemTypeListBoxValueChanged(app, event) app.Value = app.SystemTypeListBox.Value; app.system=app.Value; if app.system=="Paddlewheel" app.Q=(1141.3*3.785); %Measured device flow, converted to lpm app.I=1.36; %Measured electrical current draw to device in Amps elseif app.system=="Propeller" app.Q=(2914*3.785); %Measured device flow, converted to lpm app.I=6.54; %Measured electrical current draw to device in Amps else app.Q=(2868.6*3.785); %Measured device flow, converted to lpm app.I=9.75; %Measured electrical current draw to device in Amps end app.P=app.I*app.V; %Power required for device app.P1=(app.P*1.341)/1000; %Converted from kW to Hp app.n=(app.Q)/app.P; %Effective efficiency for device in lpm/kW if app.Value2>0

56

app.RacewayFlowGPMEditFieldValueChanged() end end

% Value changed function: RacewayFlowGPMEditField function RacewayFlowGPMEditFieldValueChanged(app, event) app.Value2 = app.RacewayFlowGPMEditField.Value; app.Qr =(app.Value2*3.785); app.N9=BoydReverse(app.T, app.W, app.Qr, 9, 3); app.N7=BoydReverse(app.T, app.W, app.Qr, 7, 3); app.N5=BoydReverse(app.T, app.W, app.Qr, 5, 3); app.N3=BoydReverse(app.T, app.W, app.Qr, 4, 3); app.Preq=(app.Qr/app.n); %Required power to provide required flow app.Preq1=((app.Qr/app.n)*1.341)/1000; %Converted from kW to Hp app.BatStorage= BatBackup(app.Qr, app.n, app.hr, app.Tvar); %System battery storage in lead acid and lithium app.MotorHp.Text=num2str(app.Preq1,'%.2f'); app.lbsLabel.Text=strcat(num2str(app.N9,'%.0f'),' lbs'); app.lbsLabel_2.Text=strcat(num2str(app.N7,'%.0f'),' lbs'); app.lbsLabel_3.Text=strcat(num2str(app.N5,'%.0f'),' lbs'); app.lbsLabel_4.Text=strcat(num2str(app.N3,'%.0f'),' lbs'); app.Generator.Text=num2str(app.Preq,'%.0f'); app.Gas=((app.Preq/1000)*.296875)*app.hr; %Fuel Consumption assumed from Honda EU2200i generator app.GasPrice=app.Gas*(app.PriceofGasinYourAreaEditField.Value); app.GasCostLabel.Text=strcat('$',num2str(app.GasPrice, '%.2f')); app.Gasoline.Text=num2str(app.Gas,'%.1f'); app.Propane=((app.Preq/1000)*.163755)*app.hr; %Fuel consumption taken from NREL fuel type paper app.PropPrice=app.Propane*app.PriceofPropaneinYourAreaEditField.V alue; app.PropCostLabel.Text=strcat('$',num2str(app.PropPrice,'%.2f')); app.PropaneUsed.Text=num2str(app.Propane,'%.1f'); app.kWh.Text=num2str((app.Preq/1000)*app.hr,'%.0f'); app.ElectricPrice=(app.Preq/1000)*app.hr*app.PriceperkWhforElectr icityinYourAreaEditField.Value; app.kWhCostLabel.Text=strcat('$',num2str(app.ElectricPrice,'%.2f' )); end end

% Component initialization methods (Access = private)

57

% Create UIFigure and components function createComponents(app)

% Create UIFigure and hide until all components are created app.UIFigure = uifigure('Visible', 'off'); app.UIFigure.Color = [0.9412 0.9412 0.9412]; app.UIFigure.Position = [100 100 634 563]; app.UIFigure.Name = 'MATLAB App';

% Create InPondRacewaySystemSizingToolLabel app.InPondRacewaySystemSizingToolLabel = uilabel(app.UIFigure); app.InPondRacewaySystemSizingToolLabel.HorizontalAlignment = 'center'; app.InPondRacewaySystemSizingToolLabel.FontSize = 24; app.InPondRacewaySystemSizingToolLabel.Position = [112.5 520 411 30]; app.InPondRacewaySystemSizingToolLabel.Text = 'In‐Pond Raceway System Sizing Tool';

% Create MinGeneratorSizeWLabel app.MinGeneratorSizeWLabel = uilabel(app.UIFigure); app.MinGeneratorSizeWLabel.Position = [332 409 135 22]; app.MinGeneratorSizeWLabel.Text = 'Min. Generator Size (W)';

% Create GallonsofGasUsedperDayLabel app.GallonsofGasUsedperDayLabel = uilabel(app.UIFigure); app.GallonsofGasUsedperDayLabel.Position = [332 384 162 22]; app.GallonsofGasUsedperDayLabel.Text = 'Gallons of Gas Used per Day';

% Create Generator app.Generator = uilabel(app.UIFigure); app.Generator.Position = [547 411 32 22]; app.Generator.Text = '0000';

% Create Gasoline app.Gasoline = uilabel(app.UIFigure); app.Gasoline.Position = [547 385 66 22]; app.Gasoline.Text = '0000';

58

% Create SystemSizingLabel app.SystemSizingLabel = uilabel(app.UIFigure); app.SystemSizingLabel.FontWeight = 'bold'; app.SystemSizingLabel.Position = [415 479 87 22]; app.SystemSizingLabel.Text = 'System Sizing';

% Create GallonsofPropaneUsedperDayLabel app.GallonsofPropaneUsedperDayLabel = uilabel(app.UIFigure); app.GallonsofPropaneUsedperDayLabel.Position = [332 334 185 22]; app.GallonsofPropaneUsedperDayLabel.Text = 'Gallons of Propane Used per Day';

% Create PropaneUsed app.PropaneUsed = uilabel(app.UIFigure); app.PropaneUsed.Position = [547 334 56 22]; app.PropaneUsed.Text = '0000';

% Create TogetstartedwiththisapplicationfirstseLabel app.TogetstartedwiththisapplicationfirstseLabel = uilabel(app.UIFigure); app.TogetstartedwiththisapplicationfirstseLabel.FontSize = 11; app.TogetstartedwiththisapplicationfirstseLabel.Position = [16 394 284 60]; app.TogetstartedwiththisapplicationfirstseLabel.Text = {'To get started with this application, first select'; 'system type from the dropdown list. Next, you may use '; 'the prepopulated values for gas and electricity '; 'or enter your own. Finally, choose target Raceway Flow, '; 'type in your target value and press the Enter key.'};

% Create EstimatedWeightofCatfishSupportedLabel app.EstimatedWeightofCatfishSupportedLabel = uilabel(app.UIFigure); app.EstimatedWeightofCatfishSupportedLabel.FontWeight = 'bold'; app.EstimatedWeightofCatfishSupportedLabel.Position = [344 213 230 22]; app.EstimatedWeightofCatfishSupportedLabel.Text = 'Estimated Weight of Catfish Supported';

% Create mglLabel app.mglLabel = uilabel(app.UIFigure);

59 app.mglLabel.Position = [343 164 38 22]; app.mglLabel.Text = '9 mg/l';

% Create mglLabel_2 app.mglLabel_2 = uilabel(app.UIFigure); app.mglLabel_2.Position = [343 133 38 22]; app.mglLabel_2.Text = '7 mg/l';

% Create mglLabel_3 app.mglLabel_3 = uilabel(app.UIFigure); app.mglLabel_3.Position = [343 99 38 22]; app.mglLabel_3.Text = '5 mg/l';

% Create mglLabel_4 app.mglLabel_4 = uilabel(app.UIFigure); app.mglLabel_4.Position = [343 67 38 22]; app.mglLabel_4.Text = '4 mg/l';

% Create lbsLabel app.lbsLabel = uilabel(app.UIFigure); app.lbsLabel.Position = [497 164 84 22]; app.lbsLabel.Text = '0000 lbs';

% Create lbsLabel_2 app.lbsLabel_2 = uilabel(app.UIFigure); app.lbsLabel_2.Position = [497 133 84 22]; app.lbsLabel_2.Text = '0000 lbs';

% Create lbsLabel_3 app.lbsLabel_3 = uilabel(app.UIFigure); app.lbsLabel_3.Position = [497 99 84 22]; app.lbsLabel_3.Text = '0000 lbs';

% Create lbsLabel_4 app.lbsLabel_4 = uilabel(app.UIFigure); app.lbsLabel_4.Position = [497 67 84 22]; app.lbsLabel_4.Text = '0000 lbs';

% Create DailyCostofGasolineLabel

60

app.DailyCostofGasolineLabel = uilabel(app.UIFigure); app.DailyCostofGasolineLabel.Position = [333 359 124 22]; app.DailyCostofGasolineLabel.Text = 'Daily Cost of Gasoline';

% Create DailyCostofPropaneLabel app.DailyCostofPropaneLabel = uilabel(app.UIFigure); app.DailyCostofPropaneLabel.Position = [333 309 122 22]; app.DailyCostofPropaneLabel.Text = 'Daily Cost of Propane';

% Create kWhofElectricityUsedperDayLabel app.kWhofElectricityUsedperDayLabel = uilabel(app.UIFigure); app.kWhofElectricityUsedperDayLabel.Position = [331 269 174 22]; app.kWhofElectricityUsedperDayLabel.Text = 'kWh of Electricity Used per Day';

% Create DailyCostofElectricityLabel app.DailyCostofElectricityLabel = uilabel(app.UIFigure); app.DailyCostofElectricityLabel.Position = [332 245 128 22]; app.DailyCostofElectricityLabel.Text = 'Daily Cost of Electricity';

% Create GasCostLabel app.GasCostLabel = uilabel(app.UIFigure); app.GasCostLabel.Position = [545 359 36 22]; app.GasCostLabel.Text = '$0.00';

% Create PropCostLabel app.PropCostLabel = uilabel(app.UIFigure); app.PropCostLabel.Position = [545 309 36 22]; app.PropCostLabel.Text = '$0.00';

% Create kWh app.kWh = uilabel(app.UIFigure); app.kWh.Position = [546 268 64 22]; app.kWh.Text = '0000';

% Create kWhCostLabel app.kWhCostLabel = uilabel(app.UIFigure); app.kWhCostLabel.Position = [544 245 66 22]; app.kWhCostLabel.Text = '$0.00';

61

% Create Pleasenoteadditionalaerationmayberequiredunder3mglLabel app.Pleasenoteadditionalaerationmayberequiredunder3mglLabel = uilabel(app.UIFigure); app.Pleasenoteadditionalaerationmayberequiredunder3mglLabel.FontS ize = 11; app.Pleasenoteadditionalaerationmayberequiredunder3mglLabel.Posit ion = [315 12 288 48]; app.Pleasenoteadditionalaerationmayberequiredunder3mglLabel.Text = {'*Please note, estimates of weight supported assume that'; 'effluent DO does not fall below 3 mg/l. Additional aeration'; 'may be required for ponds with DO concentrations 3 mg/l'; 'or less.'};

% Create MinMotorSizeHpLabel app.MinMotorSizeHpLabel = uilabel(app.UIFigure); app.MinMotorSizeHpLabel.Position = [331 450 116 22]; app.MinMotorSizeHpLabel.Text = 'Min. Motor Size (Hp)';

% Create MotorHp app.MotorHp = uilabel(app.UIFigure); app.MotorHp.Position = [549 450 29 22]; app.MotorHp.Text = '0.00';

% Create PondDOLabel app.PondDOLabel = uilabel(app.UIFigure); app.PondDOLabel.FontWeight = 'bold'; app.PondDOLabel.Position = [333 185 57 22]; app.PondDOLabel.Text = 'Pond DO';

% Create WeightSupportedLabel app.WeightSupportedLabel = uilabel(app.UIFigure); app.WeightSupportedLabel.FontWeight = 'bold'; app.WeightSupportedLabel.Position = [472 185 109 22]; app.WeightSupportedLabel.Text = 'Weight Supported';

% Create SystemTypeListBoxLabel app.SystemTypeListBoxLabel = uilabel(app.UIFigure); app.SystemTypeListBoxLabel.Position = [16 342 116 22]; app.SystemTypeListBoxLabel.Text = 'System Type';

62

% Create SystemTypeListBox app.SystemTypeListBox = uilistbox(app.UIFigure); app.SystemTypeListBox.Items = {'Paddlewheel', 'Propeller', 'Airlift', ''}; app.SystemTypeListBox.ValueChangedFcn = createCallbackFcn(app, @SystemTypeListBoxValueChanged, true); app.SystemTypeListBox.Position = [131 290 125 74]; app.SystemTypeListBox.Value = {};

% Create RacewayFlowGPMEditFieldLabel app.RacewayFlowGPMEditFieldLabel = uilabel(app.UIFigure); app.RacewayFlowGPMEditFieldLabel.Position = [16 115 132 22]; app.RacewayFlowGPMEditFieldLabel.Text = 'Raceway Flow (GPM)';

% Create RacewayFlowGPMEditField app.RacewayFlowGPMEditField = uieditfield(app.UIFigure, 'numeric'); app.RacewayFlowGPMEditField.RoundFractionalValues = 'on'; app.RacewayFlowGPMEditField.ValueDisplayFormat = '%.0f'; app.RacewayFlowGPMEditField.ValueChangedFcn = createCallbackFcn(app, @RacewayFlowGPMEditFieldValueChanged, true); app.RacewayFlowGPMEditField.Position = [147 113 109 22];

% Create PriceofGasinYourAreaEditFieldLabel app.PriceofGasinYourAreaEditFieldLabel = uilabel(app.UIFigure); app.PriceofGasinYourAreaEditFieldLabel.Position = [16 234 140 22]; app.PriceofGasinYourAreaEditFieldLabel.Text = 'Price of Gas in Your Area';

% Create PriceofGasinYourAreaEditField app.PriceofGasinYourAreaEditField = uieditfield(app.UIFigure, 'numeric'); app.PriceofGasinYourAreaEditField.Position = [171 234 85 22]; app.PriceofGasinYourAreaEditField.Value = 2.08;

63

% Create PriceofPropaneinYourAreaEditFieldLabel app.PriceofPropaneinYourAreaEditFieldLabel = uilabel(app.UIFigure); app.PriceofPropaneinYourAreaEditFieldLabel.Position = [16 200 164 22]; app.PriceofPropaneinYourAreaEditFieldLabel.Text = 'Price of Propane in Your Area';

% Create PriceofPropaneinYourAreaEditField app.PriceofPropaneinYourAreaEditField = uieditfield(app.UIFigure, 'numeric'); app.PriceofPropaneinYourAreaEditField.Position = [195 200 61 22]; app.PriceofPropaneinYourAreaEditField.Value = 2.07;

% Create PriceperkWhforElectricityinYourAreaLabel app.PriceperkWhforElectricityinYourAreaLabel = uilabel(app.UIFigure); app.PriceperkWhforElectricityinYourAreaLabel.Position = [16 154 126 28]; app.PriceperkWhforElectricityinYourAreaLabel.Text = {'Price per kWh for '; 'Electricity in Your Area'};

% Create PriceperkWhforElectricityinYourAreaEditField app.PriceperkWhforElectricityinYourAreaEditField = uieditfield(app.UIFigure, 'numeric'); app.PriceperkWhforElectricityinYourAreaEditField.Position = [157 160 99 22]; app.PriceperkWhforElectricityinYourAreaEditField.Value = 0.08;

% Show the figure after all components are created app.UIFigure.Visible = 'on'; end end

% App creation and deletion methods (Access = public)

64

% Construct app function app = Solar_Tool

% Create UIFigure and components createComponents(app)

% Register the app with App Designer registerApp(app, app.UIFigure)

% Execute the startup function runStartupFcn(app, @startupFcn)

if nargout == 0 clear app end end

% Code that executes before app deletion function delete(app)

% Delete UIFigure when app is deleted delete(app.UIFigure) end end end

65

[APPENDIX 2. SOLAR MODEL CODE]

% This model is intended to estimate Battery storage for a specified % solar energy system using an energy balance equation based off of % the law of conservation of energy.

% This system has inputs of estimated system load, estimate system battery % capacity, and SAM produced system output data. This model will use a % power balance to determine the presence of insufficient energy events, % and provide a number of such events expected per year.

clc; clear; close all;

% Simulation Variables Eout= [0.25 0.5 0.75 1]; %System load in kW

% The panel and battery sizes must be used together. The first battery size % will be used with the first panel size. The second with the second and so % on.

E= [5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... ]; %Initial Battery Storage estimate in kWh

panel = [2, 2, 2, 2, 2, 2, 2, 2 ... 3, 3, 3, 3, 3, 3, 3, 3 ... 4, 4, 4, 4, 4, 4, 4, 4 ... 5, 5, 5, 5, 5, 5, 5, 5 ... 10, 10, 10, 10, 10, 10, 10, 10 ... 15, 15, 15, 15, 15, 15, 15, 15 ... 20, 20, 20, 20, 20, 20, 20, 20 ...

66

]; % Panel Sizes

location = ["Bowling Green", "Cincinnati", "Jackson", "Lexington", "London","Louisville", "Paducah"];

StartMonth = 11; % Seasonal: Ignore starting on the 1st day of this month EndMonth = 2; % Seasonal: Ignore stopping on the last day of this month minDOD = .2; invEff = 0.93; % Inverter Efficiency. (From DC to AC) charConEff = 0.97; % Charge Controller Efficiency (From PV to DC) DoDlimit = 0.2; % This is the minimum acceptable Depth of Discharge

% Battery Efficiency battCharEff = 0.85; % Efficiency of adding charge to the battery

% This model assumes constant charging and discharging rates for each hour. % If the power input from solar is greater than the power output to the % load, then the system is charging. This efficiency is used. The % efficiency only applies to the power difference as this is what goes to % the battery. % If the power input from solar is less than the power output to the % load, then the system is discharging. This value does not apply in % discharging.

% Initialize Data Variables - These are cell variables timeStamp = cell(size(location,2),size(panel,2)); Ein = cell(size(location,2),size(panel,2));

%Import the SAM Data for pidx = 1:size(panel,2) for li = 1:size(location,2) [timeStamp{li, pidx},Ein{li, pidx}] = ... importSAMfile(join([location(li), "_", panel(pidx), "kW"],"")); %Array of system production values in kW end

67

end

Ebat = cell(size(location,2),size(panel,2),size(Eout,2)); fail = cell(size(location,2),size(panel,2),size(Eout,2)); EbatSeasonal = cell(size(location,2),size(panel,2),size(Eout,2)); failSeasonal = cell(size(location,2),size(panel,2),size(Eout,2)); runningSeason = cell(size(location,2),size(panel,2),size(Eout,2));

% These outputs are tables per location eventPerYear = cell(size(location,2),1); events = cell(size(location,2),1); eventPerYearSeasonal = cell(size(location,2),1); eventsSeasonal = cell(size(location,2),1);

% Process results for li = 1:size(location,2) % Initialize Output tables for each location eventPerYear{li} = zeros(size(Eout,2),size(panel,2)); events{li} = zeros(size(Eout,2),size(panel,2)); eventPerYearSeasonal{li} = zeros(size(Eout,2),size(panel,2)); eventsSeasonal{li} = zeros(size(Eout,2),size(panel,2)); for pidx = 1:size(panel,2) for Ei = 1:size(Eout,2) Ebat{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); %An array of zeros to be populated later fail{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); Ebat{li, pidx, Ei}(1)=E(pidx); %This sets the initial battery value as being fully charged

EbatSeasonal{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); %An array of zeros to be populated later failSeasonal{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); runningSeason{li, pidx, Ei}=true(1,numel(Ein{li, pidx})); EbatSeasonal{li, pidx, Ei}(1)=E(pidx); %This sets the initial battery value as being fully charged

for i=1:(numel(Ein{li, pidx})-1) % Normal

68

if Ein{li, pidx}(i)*charConEff>Eout(Ei)/invEff % charging Ebat{li, pidx, Ei}(i+1)=Ebat{li, pidx, Ei}(i)+(Ein{li, pidx}(i)*charConEff- Eout(Ei)/invEff)*battCharEff; else % discharging Ebat{li, pidx, Ei}(i+1)=Ebat{li, pidx, Ei}(i)+Ein{li, pidx}(i)*charConEff-Eout(Ei)/invEff; end

if Ebat{li, pidx, Ei}(i+1)>E(pidx) Ebat{li, pidx, Ei}(i+1)=E(pidx); end if Ebat{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) fail{li, pidx, Ei}(i+1)=1; Ebat{li, pidx, Ei}(i+1)=E(pidx); end % Seasonal cur_m = month(timeStamp{li, pidx}(i)); if StartMonth > EndMonth % Loop around January if ((cur_m < StartMonth) && (cur_m > EndMonth)) % Operating Period runningSeason{li, pidx, Ei}(i) = true; if i~=1 % Matlab doesn't loop like python. Can't check initial state. if runningSeason{li, pidx, Ei}(i) && not(runningSeason{li, pidx, Ei}(i-1)) % Restarting with full battery EbatSeasonal{li, pidx, Ei}(i)=E(pidx); end end if Ein{li, pidx}(i)*charConEff>Eout(Ei)/invEff % charging EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+(Ein{li, pidx}(i)*charConEff-Eout(Ei)/invEff)*battCharEff; else % discharging

69

EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+Ein{li, pidx}(i)*charConEff-Eout(Ei)/invEff; end

if EbatSeasonal{li, pidx, Ei}(i+1)>E(pidx) EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx); end if EbatSeasonal{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) failSeasonal{li, pidx, Ei}(i+1)=1; EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx); end else % Not operating runningSeason{li, pidx, Ei}(i) = false; end else % In the same calendar year if ((cur_m < StartMonth) || (cur_m > EndMonth)) % Operating Period runningSeason{li, pidx, Ei}(i) = true; if i~=1 % Matlab doesn't loop like python. Can't check initial state. if runningSeason{li, pidx, Ei}(i) && not(runningSeason{li, pidx, Ei}(i-1)) % Restarting with full battery EbatSeasonal{li, pidx, Ei}(i)=E(pidx); end end EbatSeasonal{li, pidx, Ei}(i+1)=(EbatSeasonal{li, pidx, Ei}(i)-Eout(Ei))+Ein{li, pidx}(i);

if EbatSeasonal{li, pidx, Ei}(i+1)>E(pidx) EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx);

70

end if EbatSeasonal{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) failSeasonal{li, pidx, Ei}(i+1)=1; EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx); end else % Not operating runningSeason{li, pidx, Ei}(i) = false; end end end

eventPerYear{li}(Ei, pidx)=(sum(fail{li, pidx, Ei})/21); events{li}(Ei, pidx) = sum(fail{li, pidx, Ei});

eventPerYearSeasonal{li}(Ei, pidx)=(sum(failSeasonal{li, pidx, Ei})/21); eventsSeasonal{li}(Ei, pidx) = sum(failSeasonal{li, pidx, Ei}); end end end

tableheaderYear = ["Insufficient Energy Events All Year", "Batt Size (kWh) [E]", "panel size (kw)", Eout+" kW Load"]'; tableheaderSeason = [["Insufficient Energy Events Seasonal: Month "+StartMonth+" to "+EndMonth], "Batt Size (kWh) [E]", "panel size (kw)", Eout+" kW Load"]'; tableSeasonStart = size(tableheaderYear,1)+2;

for li = 1:size(location,2) filename = 'failureTables.xlsx'; % All year data writematrix(tableheaderYear,filename,'Sheet', location(li),'Range','A1'); writematrix(E,filename,'Sheet', location(li),'Range','B2'); writematrix(panel,filename,'Sheet', location(li),'Range','B3');

71

writematrix(events{li},filename,'Sheet', location(li),'Range','B4');

writematrix(tableheaderSeason,filename,'Sheet', location(li),'Range',["A"+tableSeasonStart]); writematrix(E,filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+1)]); writematrix(panel,filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+2)]); writematrix(eventsSeasonal{li},filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+3)]); end

% The outputs will be in four tables. The tables are cell arrays with each % corresponding to a location. They are in order of the location list. % For each location, the rows are loads (Eout), and the columns are the % panel/battery system sizes (panel+E). The first column will be for the % first value in both the panel and battery size (E) array. The second will % be for the second set of values in each array (E(2) and panel(2), and so % on.

72

[APPENDIX 3. GENERATOR CODE]

% This model incorporates solar, battery, and generator backup.

% This model is intended to estimate Battery storage for a specified % solar energy system using an energy balance equation based off of % the law of conservation of energy.

% This system has inputs of estimated system load, estimate system battery % capacity, and SAM produced system output data. This model will use a % power balance to determine the battery State of Charge. If the battery % reaches a minimum value, a generator is used to restore power to a % minimum level.

clc; clear; close all; % Simulation Variables

Eout= [0.25 0.5 0.75 1]; %System load in kW % The panel and battery sizes must be used together. The first battery size % will be used with the first panel size. The second with the second and so % on.

E= [5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... 5, 10, 15, 20, 30, 40, 50, 60 ... ]; %Initial Battery Storage estimate in kWh

panel = [2, 2, 2, 2, 2, 2, 2, 2 ... 3, 3, 3, 3, 3, 3, 3, 3 ...

73

4, 4, 4, 4, 4, 4, 4, 4 ... 5, 5, 5, 5, 5, 5, 5, 5 ... 10, 10, 10, 10, 10, 10, 10, 10 ... 15, 15, 15, 15, 15, 15, 15, 15 ... 20, 20, 20, 20, 20, 20, 20, 20 ... ]; % Panel Sizes

location = ["Bowling Green", "Cincinnati", "Jackson", "Lexington", "London", "Louisville", "Paducah"]; StartMonth = 11; % Seasonal: Ignore starting on the 1st day of this month EndMonth = 2; % Seasonal: Ignore stopping on the last day of this month minDOD = .2; invEff = 0.93; % Inverter Efficiency. (From DC to AC) charConEff = 0.97; % Charge Controller Efficiency (From PV to DC) DoDlimit = 0.2; % This is the minimum acceptable Depth of Discharge genOutput = 2; % The output of the generator in kW genShutOffLevel = 0.85; % The generator turns on at the DoDlimit and runs until this level genOn = false; % Flag to determine if the generator is on. genOnSeason = false; % Flag to determine if the generator is on (Seasonal)

% Battery Efficiency battCharEff = 0.85; % Efficiency of adding charge to the battery

% This model assumes constant charging and discharging rates for each hour. % If the power input from solar is greater than the power output to the % load, then the system is charging. This efficiency is used. The % efficiency only applies to the power difference as this is what goes to % the battery. % If the power input from solar is less than the power output to the % load, then the system is discharging. This value does not apply in % discharging.

% Initialize Data Variables - These are cell variables timeStamp = cell(size(location,2),size(panel,2));

74

Ein = cell(size(location,2),size(panel,2));

%Import the SAM Data for pidx = 1:size(panel,2) for li = 1:size(location,2) [timeStamp{li, pidx},Ein{li, pidx}] = ... importSAMfile(join([location(li), "_", panel(pidx), "kW"],"")); %Array of system production values in kW end end

Ebat = cell(size(location,2),size(panel,2),size(Eout,2)); genRun = cell(size(location,2),size(panel,2),size(Eout,2)); EbatSeasonal = cell(size(location,2),size(panel,2),size(Eout,2)); genRunSeasonal = cell(size(location,2),size(panel,2),size(Eout,2)); runningSeason = cell(size(location,2),size(panel,2),size(Eout,2)); % These outputs are tables per location genHoursPerYear = cell(size(location,2),1); genHours = cell(size(location,2),1); genHoursPerYearSeasonal = cell(size(location,2),1); genHoursSeasonal = cell(size(location,2),1);

% Process results for li = 1:size(location,2) % Initialize Output tables for each location genHoursPerYear{li} = zeros(size(Eout,2),size(panel,2)); genHours{li} = zeros(size(Eout,2),size(panel,2)); genHoursPerYearSeasonal{li} = zeros(size(Eout,2),size(panel,2)); genHoursSeasonal{li} = zeros(size(Eout,2),size(panel,2)); for pidx = 1:size(panel,2) for Ei = 1:size(Eout,2) Ebat{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); %An array of zeros to be populated later genRun{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); Ebat{li, pidx, Ei}(1)=E(pidx); %This sets the initial battery value as being fully charged

75

EbatSeasonal{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); %An array of zeros to be populated later genRunSeasonal{li, pidx, Ei}=zeros(1,numel(Ein{li, pidx})); runningSeason{li, pidx, Ei}=true(1,numel(Ein{li, pidx})); EbatSeasonal{li, pidx, Ei}(1)=E(pidx); %This sets the initial battery value as being fully charged

genOn = false; % Ensure Generator starts off. genOnSeason = false; % Ensure Generator starts off.

for i=1:(numel(Ein{li, pidx})-1) % Normal if (Ein{li, pidx}(i)+genOn*genOutput)*charConEff>Eout(Ei)/invEff % charging - some power into battery and affected by % battery charging efficiency Ebat{li, pidx, Ei}(i+1)=Ebat{li, pidx, Ei}(i)+((Ein{li, pidx}(i)+genOn*genOutput)*charConEff- Eout(Ei)/invEff)*battCharEff; else % discharging - all power direct to load Ebat{li, pidx, Ei}(i+1)=Ebat{li, pidx, Ei}(i)+(Ein{li, pidx}(i)+genOn*genOutput)*charConEff- Eout(Ei)/invEff; end

if Ebat{li, pidx, Ei}(i+1)>E(pidx) % Cannot be more than 100% charged. Ebat{li, pidx, Ei}(i+1)=E(pidx); end if Ebat{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) % Detect when minimum is reached genOn = true; % Turn on generator end if genOn == true genRun{li, pidx, Ei}(i+1)=true; % Record that the generator is on if Ebat{li, pidx, Ei}(i+1)>= genShutOffLevel*E(pidx) % We have reached shutoff level genOn = false; % Stop Generator end

76

end % Seasonal cur_m = month(timeStamp{li, pidx}(i)); if StartMonth > EndMonth % Loop around January if ((cur_m < StartMonth) && (cur_m > EndMonth)) % Operating Period runningSeason{li, pidx, Ei}(i) = true; if i~=1 % Matlab doesn't loop like python. Can't check initial state. if runningSeason{li, pidx, Ei}(i) && not(runningSeason{li, pidx, Ei}(i-1)) % Restarting with full battery EbatSeasonal{li, pidx, Ei}(i)=E(pidx); end end if (Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff>Eout(Ei)/invEff % charging EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+((Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff- Eout(Ei)/invEff)*battCharEff; else % discharging EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+(Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff-Eout(Ei)/invEff; end

if EbatSeasonal{li, pidx, Ei}(i+1)>E(pidx) % Cannot be more than 100% charged. EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx); end if EbatSeasonal{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) % Detect when minimum is reached genOnSeason = true; % Turn on generator end if genOnSeason == true

77

genRunSeasonal{li, pidx, Ei}(i+1)=true; % Record that the generator is on if EbatSeasonal{li, pidx, Ei}(i+1)>= genShutOffLevel*E(pidx) % We have reached shutoff level genOnSeason = false; % Stop Generator end end else % Not operating runningSeason{li, pidx, Ei}(i) = false; end else % In the same calendar year if ((cur_m < StartMonth) || (cur_m > EndMonth)) % Operating Period runningSeason{li, pidx, Ei}(i) = true; if i~=1 % Matlab doesn't loop like python. Can't check initial state. if runningSeason{li, pidx, Ei}(i) && not(runningSeason{li, pidx, Ei}(i-1)) % Restarting with full battery EbatSeasonal{li, pidx, Ei}(i)=E(pidx); end end if (Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff>Eout(Ei)/invEff % charging EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+((Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff- Eout(Ei)/invEff)*battCharEff; else % discharging EbatSeasonal{li, pidx, Ei}(i+1)=EbatSeasonal{li, pidx, Ei}(i)+(Ein{li, pidx}(i)+genOnSeason*genOutput)*charConEff-Eout(Ei)/invEff; end if EbatSeasonal{li, pidx, Ei}(i+1)>E(pidx) % Cannot be more than 100% charged.

78

EbatSeasonal{li, pidx, Ei}(i+1)=E(pidx); end if EbatSeasonal{li, pidx, Ei}(i+1)<=(E(pidx) * DoDlimit) % Detect when minimum is reached genOnSeason = true; % Turn on generator end if genOnSeason == true genRunSeasonal{li, pidx, Ei}(i+1)=true; % Record that the generator is on if EbatSeasonal{li, pidx, Ei}(i+1)>= genShutOffLevel*E(pidx) % We have reached shutoff level genOnSeason = false; % Stop Generator end end else % Not operating runningSeason{li, pidx, Ei}(i) = false; end end end

genHoursPerYear{li}(Ei, pidx)=(sum(genRun{li, pidx, Ei})/21); genHours{li}(Ei, pidx) = sum(genRun{li, pidx, Ei});

genHoursPerYearSeasonal{li}(Ei, pidx)=(sum(genRunSeasonal{li, pidx, Ei})/21); genHoursSeasonal{li}(Ei, pidx) = sum(genRunSeasonal{li, pidx, Ei}); end end end tableheaderYear = ["Generator Operating Hours over 21 Years (All Year)", "Batt Size (kWh) [E]", "panel size (kw)", Eout+" kW Load"]'; tableheaderSeason = [["Generator Operating Hours over 21 Years (Seasonal: Month "+StartMonth+" to "+EndMonth+")"], "Batt Size (kWh) [E]", "panel size (kw)", Eout+" kW Load"]';

79

tableSeasonStart = size(tableheaderYear,1)+2;

for li = 1:size(location,2) filename = 'GenTables.xlsx'; % All year data writematrix(tableheaderYear,filename,'Sheet', location(li),'Range','A1','WriteMode','overwritesheet'); % Erases the worksheet first writematrix(E,filename,'Sheet', location(li),'Range','B2'); writematrix(panel,filename,'Sheet', location(li),'Range','B3'); writematrix(genHours{li},filename,'Sheet', location(li),'Range','B4');

writematrix(tableheaderSeason,filename,'Sheet', location(li),'Range',["A"+tableSeasonStart]); writematrix(E,filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+1)]); writematrix(panel,filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+2)]); writematrix(genHoursSeasonal{li},filename,'Sheet', location(li),'Range',["B"+(tableSeasonStart+3)]); end

% The outputs will be in four tables. The tables are cell arrays with each % corresponding to a location. They are in order of the location list. % For each location, the rows are loads (Eout), and the columns are the % panel/battery system sizes (panel+E). The first column will be for the % first value in both the panel and battery size (E) array. The second will % be for the second set of values in each array (E(2) and panel(2), and so % on

80

REFERENCES

Abedi, S., Alimardani, A., Gharehpetian, G. B., Riahy, G. H., & Hosseinian, S. H. (2012). A comprehensive method for optimal power management and design of hybrid RES-based autonomous energy systems. Renewable and Sustainable Energy Reviews, 16(3), 1577-1587. doi:10.1016/j.rser.2011.11.030 Andrews, J. W., & Matsuda, Y. (1975). The Influence of Various Culture Conditions on the Oxygen Consumption of Channel Catfish. Transactions of the American Fisheries Society, 104(2), 322-327. doi:10.1577/1548- 8659(1975)104%3C322:tiovcc%3E2.0.co;2 Andrews, J. W., Murai, T., & Gibbons, G. (1973). The Influence of Dissolved Oxygen on the Growth of Channel Catfish. Transactions of the American Fisheries Society, 102(4), 835-838. doi:10.1577/1548- 8659(1973)102%3C835:TIODOO%3E2.0.CO;2 Ang, K. H., Chong, G., & Li, Y. (2005). PID Control System Analysis, Design, and Technology. IEEE Transactions on Control Systems Technology, 13(4), 559-576. doi:10.1109/TCST.2005.847331 Bajpai, P., & Dash, V. (2012). Hybrid renewable energy systems for power generation in stand-alone applications: A review. Renewable and Sustainable Energy Reviews, 16(5), 2926-2939. doi:10.1016/j.rser.2012.02.009 Baurzhan, S., & Jenkins, G. P. (2016). Off-grid Solar PV: Is it an affordable or appropriate solution for rural electrification in Sub-Saharan Africa? Renewable and Sustainable Energy Reviews, 60, 1405-1418. doi:10.1016/j.rser.2016.03.016 Bearon Manufacturing. (2021). Vari-Speed: An Axial Flow System for Oysters. Retrieved from Bearon Aquatics: https://thepowerhouseinc.com/vari-speed-22/ Bernardez, R. G. (1995). Evaluation of an in-pond raceway system and its economic feasibility for fish production. Auburn, AL: Auburn University. Blair, N., DiOrio, N., Freeman, J., Gilman, P., Janzou, S., Neises, T., & Wagner, M. (2018). System Advisor Model General Description (Version 2017.9.5). Golden, CO: National Renewable Energy Laboratory. Retrieved from https://www.nrel.gov/docs/fy18osti/70414.pdf Boyd, C. E. (2015, November 5). Efficiency of Mechanical Aeration. Retrieved from Global Aquaculture Alliance Web Site: https://www.aquaculturealliance.org/advocate/efficiency-of-mechanical- aeration/?headlessPrint=AAAAAPIA9c8r7gs82oWZBA Boyd, C. E., & Tucker, C. S. (1998). Pond Aquaculture Water Quality Management. Kluwer Academic Publishers: Boston, MA. Boyd, C. E., Romaire, R. P., & Johnston, E. (1978). Predicting Early Morning Dissolved Oxygen Concentrations in Channel Catfish Ponds. Transactions of the American Fisheries Society, 107(3), 484-492. doi:10.1577/1548- 8659(1978)107%3C484:PEMDOC%3E2.0.CO;2

81

Boyd, C. E., Torrans, E. L., & Tucker, C. S. (2018). Dissolved Oxygen and Aeration in Ictalurid Catfish Aquaculture. Journal of the World Aquaculture Society, 49(1), 7- 70. doi:10.1111/jwas.12469 Brown, T. W., & Tucker, C. S. (2013). Pumping Performance of a Slow-Rotating Paddlewheel for Split-Pond Aquaculture Systems. North American Journal of Aquaculture, 75(2), 153-158. doi:10.1080/15222055.2012.743935 Brown, T. W., & Tucker, C. S. (2014). Pumping Performance of a Modified Commercial Paddlewheel Aerator for Split-Pond Aquaculture Systems. North American Journal of Aquaculture, 76, 72-78. doi:10.1080/15222055.2013.860067 Brown, T. W., Boyd, C. E., & Chappell, J. A. (2015). Organic carbon and dissolved oxygen budgets for a commercial-size, in-pond raceway. Journal of the World Aquaculture Society, 46(5), 448-539. doi:10.1111/jwas.12208 Brown, T. W., Chappell, J. A., & Boyd, C. E. (2011). A commercial-scale, in-pond raceway system for Ictalurid catfish production. , 44, 72-79. doi:10.1016/j.aquaeng.2011.03.003 Brown, T. W., Chappell, J. A., & Hanson, T. R. (2010). In-Pond Raceway Demonstrates Economic Benefits for Catfish Production. Global Aquaculture Advocate, 18-21. Brown, T. W., Hanson, T. R., Chappell, J. A., Boyd, C. E., & Wilson, Jr., D. S. (2014). Economic Feasibility of an In-Pond Raceway System for Commercial Catfish Production in West Alabama. North American Journal of Aquaculture, 76(1), 79- 89. doi:10.1080/15222055.2013.862195 Brown, T. W., Tucker, C. S., & Rutland, B. L. (2016). Performance evaluation of four different methods for circulating water in commercial-scale, split-pond aquaculture systems. Aquacultural Engineering, 70, 33-41. doi:10.1016/j.aquaeng.2015.12.002 Brune, D. E., Kirk, K., & Eversole, A. G. (2004). Autotrophic intensification of pond aquaculture; shrimp production in a partitioned aquaculture system. Proceedings of the Fifth International Conference on Recirculating Aquaculture, (pp. 201-210). Roanoke, VA. Brune, D. E., Schwartz, G., Eversole, A. G., Collier, J. A., & Schwedler, T. E. (2003). Intensification of pond aquaculture and high rate photosynthetic systems. Aquacultural Engineering, 28, 65-86. doi:10.1016/S0144-8609(03)00025-6 Brune, D. E., Tucker, C., Massingill, M., & Chappell, J. (2012). Partitioned Aquaculture Systems. In J. Tidwell, Aquaculture Production Systems (pp. 308-342). Hoboken, NJ: John Wiley & Sons. Carbon Disclosure Project. (2009). Global Electric Utilities Building business resilience to inevitable climate change. Oxford: Carbon Disclosure Project. Celik, A. N., Muneer, T., & Clarke, P. (2008). Optimal Sizing and Life Cycle Assessment of Residential Photovoltaic Energy Systems with Battery Storage. Progress in Photovoltaics: Research and Applications, 16, 69-85. doi:10.1002/pip.774 Chang, C.-k. (2019). Factors Affecting Capacity Design of Lithium-Ion Stationary Batteries. Batteries, 5(3), 58. doi:10.3390/batteries5030058

82

Chen, J. H., Sung, W. T., & Lin, G. Y. (2015). Automated monitoring system for the fish farm aquaculture environment. IEE International Conference on Systems, Man, and Cybernetics (pp. 1161-1166). Piscataway, N.J.: IEEE. doi:10.1109/SMC.2015.208 CLICK PLCs. (2021). CLICK Programming Help File v. 3.00. Cumming, GA: AutomationDirect.com. Comello, S., & Reichelstein, S. (2019). The emergence of cost-effective battery storage. Nature Communications, 10, 2038-9. doi:10.1038/s41467-019-09988-z Corbett, D. M. (1943). Stream Gauging Procedure: A Manual Describing Methods and Practices of the Geologic Survey. Washington, D.C.: USGS. Cripps, S. J., & Bergheim, A. (2000). Solids management and removal for intensive land- based aquaculture systems. Aquacultural Engineering, 22(1-2), 33-56. doi:10.1016/S0144-8609(00)00031-5 de Blas, I., Muniesa, A., Vallejo, A., & Ruiz-Zarzuela, I. (2020). Assessment of Sample Size Calculations Used in Aquaculture by Simulation Techniques. Frontiers in Veterinary Science, 7(253), 9. doi:10.3389/fvets.2020.00253 Dobos, A. P. (2014). PVWatts Version 5 Manual. Golden, CO: National Renewable Energy Laboratory. Dong, J., Zhu, L., Su, Y., Ma, Y., Wang, F., Tolbert, L. M., . . . Bruce, L. (2018). Battery and backup sizing for a resilient microgrid under stochastic extreme events. IET Generation, Transmission & Distribution, 12(20), 4443-4450. doi:10.1049/iet- gtd.2018.5883 Dvorak, J. S., Shockley, J., Mason, R., & McNeill, S. (2018). Modelling the costs associated with high-moisture grain for mobile apps. Computers and Electronics in Agriculture, 153, 313-317. doi:https://doi.org/10.1016/j.compag.2018.06.038 Ellabban, O., Abu-Rub, H., & Blaabjerg, F. (2014). Renewable energy resources: Current status, future prospects and their enabling technology. Renewable and Sustainable Energy Reviews, 39, 748-764. doi:10.1016/j.rser.2014.07.113 EnerSys. (2016). Odyssey Battery Technical Manual - Ninth Edition. Reading, PA: Odyssey Battery. Ericson, S., & Olis, D. (2019). A Comparison of Fuel Choice for Backup Generators. Washington, D.C.: National Renewable Energy Laboritory. FAO. (2020). State of World Fisheries and Aquaculture 2020. Sustainability in action. Rome. doi:10.4060/ca9229en Fathima, H., & Palanisamy, K. (2015). Optimized Sizing, Selection, and Economic Analysis of Battery Energy Storage for Grid-Connected Wind-PV Hybrid System. Modelling and Simulation in Engineering, 2015, 16. doi:10.1155/2015/713530 Feldman, D., Barbose, G., Margolis, R., Wiser, R., Darghouth, N., & Goodrich, A. (2012). Photovoltaic (PV) Pricing Trends: Historical, Recent, and Near-Term Projections. Washington, D.C.: NREL.

83

Ferjani, A., Mann, S., & Zimmermann, A. (2018). An evaluation of Swiss agriculture's contribution to food security with decision support system for food security strategy. British Food Journal, 120(9), 2116-2128. doi:10.1108/BFJ-12-2017-0709 Fernandez-Infantes, A., Contreras, J., & Bernal-Agustin, J. L. (2006). Design of grid connected PV systems considering electrical, economical, and environmental aspects: A practical case. Renewable Energy, 31, 2041-2061. doi:10.1016/j.renene.2005.09.028 Filgueira, R., Grant, J., & Strand, Ø. (2014). Implementation of marine spatial planning in shellfish aquaculture management: modeling studies in a Norwegian fjord. Ecological Applications, 24(4), 832-843. doi:10.1890/13-0479.1 Fizer, E., Gray, D. D., & Semmens, K. J. (2013). Effect of screens on the turbulent flow in the quiescent zone of a rectangular aquaculture raceway. Aquacultural Engineering, 57, 48-53. doi:10.1016/j.aquaeng.2013.07.001 Fornshell, G., Hinshaw, J., & Tidwell, J. H. (2012). Flow-Through Raceways. In J. Tidwell, Aquaculture Production Systems (pp. 173-190). Hoboken, NJ: John Wiley & Sons. Fukase, E., & Martin, W. (2017). Economic Growth, Convergence, and World Food Demand and Supply. Washington, D.C.: World Bank Group. Fullerton, G. C. (2016). Economic Viability of Floating In-Pond Raceway Systems for Commercial Hybrid Catfish Production. Auburn, AL: Auburn University. Füllner, G., Gottschalk, T., & Pfeifer, M. (2007). Experiments for the production of hybrid striped bass in in-pond circulation systems. Aquaculture International, 15, 241- 248. doi:10.1007/s10499-007-9093-2 Gao, G. (2012). World Food Demand. American Journal of Agricultural Economics, 94(1), 25-51. doi:10.1093/ajae/aar133 Genconnex. (2019). Propane & Natural Gas Modified EU2200i. Hingham, MA: New England Gen-Connect, LLC. Retrieved from Genconnex Direct: https://genconnexdirect.net/honda_propane_modified_generators.htm Gordon, N. D., McMahon, T. A., Finlayson, B. L., Gippel, C. J., & Nathan, R. J. (2004). Stream Hydrology: An Introduction for Ecologists (2nd ed.). Chichester, West Sussex: John Wiley & Sons. Han, E., Ines, A. V., & Baethgen, W. E. (2017). Climate-Agriculture-Modeling and Decision Tool (CAMDT): A software framework for climate risk management in agriculture. Environmental Modelling & Software, 95, 102-114. doi:10.1016/j.envsoft.2017.06.024 Hartleb, C. (2004). Floating raceways used to raise yellow perch at cranberry farms. Aquaculture Magazine, 30(1), 1-4. Hendarti, R., & Septiafani, L. (2020). Economic analysis on the application of solar panels on an aquaculture. IOP Conf. Series: Earth and Environmental Science, 426, 012073. doi:10.1088/1755-1315/426/1/012073

84

Hendarti, R., Wangidjaja, W., & Septiafani, L. (2018). A study of solar energy for an aquaculture in Jakarta. IOP Conf. Series: Earth and Environmental Science, 195, 012096. doi:10.1088/1755-1315/195/1/012096 Hesse, H. C., Martins, R., Musilek, P., Naumann, M., Truong, C. N., & Jossen, A. (2017). Economic Optimization of Component Sizing for Residential Battery Storage Systems. Energies, 10, 835. doi:10.3390/en10070835 Honda Motor Company. (2019). Owner's Manual Generator EU2200i. Samut Prakan, TH: Honda Motor Company. Hossain, M., Mekhelif, S., & Olatomiwa, L. (2017). Performance evaluation of a stand- alone PV-Wind-Diesel-Battery hybrid system feasible for a large resort center in South China Sea, Malaysia. Sustainable Cities Society, 28, 358-366. Huang, J., Qu, X., Wan, M., Ying, J., Li, Y., Zhu, F., . . . Li, W. (2015). Investigation on the performance of raceway ponds with internal structures by means of CFD simulations and experiments. Algal Research, 10, 64-71. doi:10.1016/j.algal.201504.012 Huang, J., Yang, Q., Chen, J., Wan, M., Ying, J., Fan, F., . . . Li, Y. (2016). Design and optimization of a novel airlift-driven sloping raceway pond with numerical and practical experiments. Algal Research, 20, 164-171. doi:10.1016/j.algal.2016.09.023 Irfan, M., Zhao, Z.-y., Ahmad, M., & Rehman, A. (2019). A techno-economic analysis of off-grid solar PV system: A case study for Punjab Province in Pakistan. Processes, 7(10), 708. doi:10.3390/pr7100708 ISO. (2007). Hydrometry-Measurement of Liquid Flow in Open Channels Using Current- Meters or Floats. Brussels: European Committee for Standardization. Kasco Marine. (2021). Surface Aerators. Retrieved from Kasco Marine: https://kascomarine.com/product/surface-aerators/ Kentucky Geological Survey. (2020). Water Fact Sheet. Lexington, KY: University of Kentucky. Kramer, D. L. (1987). Dissolved oxygen and fish behavior. Environmental Biology of Fishes, 18(2), 81-92. Krishna, K. S., & Kumar, K. S. (2015). A review on hybrid renewable energy systems. Renewable and Sustainable Energy Reviews, 52, 907-916. doi:10.1016/j.rser.2015.07.187 Kumar, G., Engle, C. R., Hanson, T. R., Tucker, C. S., Brown, T. W., Bott, L. B., . . . Torrans, E. L. (2018). Economics of alternative catfish production technologies. Journal of the World Aquaculture Society, 49, 1039-1057. doi:10.1111/jwas.12555 Lai-Fa, Z., & Boyd, C. E. (1988). Nightly Aeration to Increase the Efficiency of Channel Catfish Production. The Progressive Fish-Culturist, 50(4), 237-242. doi:10.1577/1548-8640(1988)050%3C0237:NATITE%3E2.3.CO;2

85

Larson, J., & Roney, J. M. (2013, June 12). Farmed Fish Production Overtakes Beef. Retrieved from Earth Policy Institute: http://www.earth- policy.org/plan_b_updates/2013/update114 Leavitt, D., Abbate, D., Brice, J.-P., DeBellis, M., Phothisane, S., White, C., . . . Riley, L. (2008). Can we upgrade the flupsy for improved performance with reduced operating costs? Journal of Shellfish Research, 27(4), 1023-1023. Lorenz, P., Pinner, D., & Seitz, T. (2008). The economics of solar power. Houston: McKinsey & Company. Madziga, M., Rahil, A., & Mansoor, R. (2018). Comparison between Three Off-Grid Hybrid Systems (Solar Photovoltaic, Diesel Generator and Battery Storage System) for Electrification for Gwakwani Village, South Africa. Environments, 5, 57. doi:10.3390/environments5050057 Mandelli, S., Merlo, M., & Colombo, E. (2016). Novel procedure to formulate load profiles for off-grid rural areas. Energy for Sustainable Development, 31, 130-142. doi:10.1016/j.esd.2016.01.005 Masser, M. P. (2012). In-Pond Raceways. In J. Tidwell, Aquaculture Production Systems (pp. 387-394). Hoboken, N.J.: John Wiley & Sons. Mbugwa, G. w., Prager, S. D., & Krall, J. M. (2015). Utilization of spatial decision support systems decision-making in dryland agriculture: A Tifton burclover case study. Computers and Electronics in Agriculture, 118, 215-224. doi:10.1016/j.compag.2015.09.008 Missouri State University OEWRI. (2007). Standard Operating Procedure for: Operation of the SonTek/YSI FlowTracker Handheld Acoustic Doppler Velocimeter (ADV). Springfield: Missouri State University. Mohanty, P., Sharma, K. R., Gujar, M., Kolhe, M., & Azmi, A. N. (2016). PV System Design for Off-Grid Applications. In P. Mohanty, T. Muneer, & M. Kolhe, Solar Photovoltaic System Applications (pp. 49-83). Springer. doi:10.1007/978-3-319- 14663-8_3 Muceka, R., Kukeera, T., Alokore, Y., Noara, K., & Groh, S. (2018). Integrating a Solar PV System with a Household Based Backup Generator for Hybrid Swarm Electrification: A Case Study of Nigeria. Africa-EU Renewable Energy Research and Innovation Symposium 2018 (RERIS 2018) (pp. 43-58). Chame: Springer. doi:10.1007/978-3-319-93438-9_4 Mukherjee, S., Nateghi, R., & Hastak, M. (2018). A multi-hazard approach to assess severe weather-induced major power outage risks in the U.S. Reliability Engineering & System Safety, 175, 283-305. doi:10.1016/j.ress.2018.03.015 Narayan, N., Qin, Z., Popovic-Gerber, J., Diehl, J.-C., Bauer, P., & Zeman, M. (2020). Stochastic load profile construction for the multi-tier framework for household electricity access using off-grid DC appliances. Energy Efficiency, 13, 197-215. doi:10.1007/s12053-018-9725-6

86

National Marine Fisheries Service. (2020). Fisheries of the United States, 2018. Silver , MD: U.S. Department of Commerce. Retrieved from https://www.fisheries.noaa.gov/national/commercial-fishing/fisheries-united- states-2018 Navarro-Hellín, H., Martínez-del-Rincon, J., Domingo-Miguel , R., Soto-Valles, F., & Torres-Sánchez, R. (2016). A decision support system for managing irrigation in agriculture. Computers and Electronics in Agriculture, 124, 121-131. doi:10.1016/j.compag.2016.04.003 Nguyen, T., Sines, S., Smith, B., & Whitlock, H. (2019, April 30). Final Report. Retrieved from InPondRaceway: https://sites.google.com/view/inpondraceway/final- report?authuser=0 Nobre, A. M., Musango, J. K., de Wit, M. P., & Ferreira, J. G. (2009). A dynamic ecological-economic modeling approach for aquaculture management. Ecological Economics, 68, 3007-3017. doi:10.1016/j.ecolecon.2009.06.019 NREL. (2019, December 5). What is the NSRDB? Retrieved from NSRDB: National Solar Radiation Database: https://nsrdb.nrel.gov/about/what-is-the-nsrdb.html Ouyang, Y., Nkedi-Kizza, P., Wu, Q. T., Shinde, D., & Huang, C. H. (2006). Assessment of Seasonal Variation in Surface Water Quality. Water Resources, 40(20), 3800- 3810. doi:10.1016/j.watres.2006.08.030 Parker, N. C., & Suttle, M. A. (1987). Design of Airlift Pumps for Water Circulation and Aeration in Aquaculture. Aquacultural Engineering, 6, 97-110. doi:10.1016/0144- 8609(87)90008-2 Rahil, A., Gammon, R., & Brown, N. (2018). Techno-economic assessment of dispatchable hydrogen production by multiple electrolysers in Libya. Journal of Energy Storage, 16, 46-60. Rantz, S. E. (1982). Measurement and Computation of Streamflow: Volume 1. Measurement of Stage and Discharge. Washington, D.C.: USGS. Rehmel, M. (2007). Application of Acoustic Doppler Velocimeters for Streamflow Measurements. Journal of Hyrdraulic Engineering, 133(12), 1433-1438. doi:10.1061/(ASCE)0733-9429(2007)133:12(1433) Rivara, G., Tetrault, K., & Patricio, R. M. (2002). A Low Cost Floating Upweller Shellfish Nursery System Construction and Operations Guide. Ithaca: Cornell Cooperative Extension Marine Program. Round, F. E. (1981). The Ecology of Algae. Cambridge: Cambridge University Press. Ru, Y., Kleissel, J., & Martinez, S. (2013). Storage Size Determination for Grid-Connected Photovoltaic Systems. IEEE Transactions on Sustainable Energy, 4(1), 68-81. doi:10.1109/TSTE.2012.2199339 Rupnik, R., Kukar, M., Vračar, P., Košir, D., Pevec, D., & Bosnić, Z. (2019). AgroDSS: A decision support system for agriculture and farming. Computers and Electronics in Agriculture, 161, 260-271. doi:10.1016/j.compag.2018.04.001

87

Shang, Y. C., & Costa-Pierce, B. A. (1983). Integrated Aquaculture-Agriculture Farming Systems: Some Economic Aspects. Journal of the World Society, 14(1- 4), 523-530. doi:10.1111/j.1749-7345.1983.tb00104.x Solaun, K., & Cerdá, E. (2019). Climate change impacts on renewable energy generation. A review of quantitative projections. Renewable and Sustainable Energy Reviews, 116, 109415. doi:10.1016/j.rser.2019.109415 SonTek. (2019). FlowTracker 2 User's Manual. San Diego: SonTek. Stevens, J. W., & Corey, G. P. (1996). A Study of Lead-Acid Battery Efficiency Near Top- of-Charge and the Impact on PV System Design. Photovoltaics Specialists Conference (pp. 1485-1488). Washington, D.C: IEEE. doi:0-7803-3166- 4/96/$5.00 Superior Raceway Systems. (2018). Size & Shape. Retrieved from Superior Raceway Systems: https://www.superiorraceways.com/sizeandshape/ Thornhill, N. F., & Hägglund, T. (1997). Detection and diagnosis of oscillation in control loops. Control Engineering Practice, 5(10), 1343-1354. doi:10.1016/S0967- 0661(97)00131-7 Torrans, E. L. (2005). Effect of Oxygen Management on Culture Performance of Channel Catfish in Earthen Ponds. North American Journal of Aquaculture, 67, 275-288. doi:10.1577/A04-071.1 Trojan Battery Company. (2019). User's Guide. Santa Fe Springs, CA: Trojan Battery Company. Tucker, C., & Hargreaves, J. (2012). Ponds. In J. Tidwell, Aquaculture Production Systems (pp. 191-244). Hoboken, NJ: John Wiley & Sons. U.S. EIA. (17, June 2020). How much carbon dioxide is produced when different fuels are burned? Retrieved from U.S. Energy Information Administration Web Site: https://www.eia.gov/tools/faqs/faq.php?id=73&t=11 United States Energy Information Administration. (2020). Levelized Cost and Levelized Avoided Cost of New Generation Resources in the Annual Energy Outlook 2020. Washington, D.C. Retrieved from https://www.eia.gov/outlooks/aeo/pdf/electricity_generation.pdf Viadero, Jr., R. C., Rumberg, A., Gray, D. D., Tierney, A. E., & Semmens, K. J. (2006). Acoustic doppler velocimetry in aquaculture research: Raceway and quiescent zone hydrodynamics. Aquacultural Engineering, 34, 16-25. doi:10.1016/j.aquaeng.2005.03.002 Vivas, F. J., De las Heras, A., Segura, F., & Andújar, J. M. (2018). A review of energy management strategies for renewable energy systems with hydrogen backup. Renewable and Sustainable Energy Reviews, 82, 126-155. doi:10.1016/j.rser.2017.09.014 WMO. (2010). Manual of Field Work: Volume 1. Stream Gauging. Geneva: WMO.

88

Wurts, W. A., & Durborow, R. M. (1992). Interactions of pH, Carbon Dioxide, Alkalinity and Hardness in Fish Ponds. Stoneville, MS: Southern Regional Aquaculture Center. Retrieved from https://appliedecology.cals.ncsu.edu/wp- content/uploads/SRAC-0464.pdf Yang, B., Li, W., Zhao, Y., & He, X. (2010). Design and Analysis of a Grid-Connected Photovoltaic Power System. IEEE Transactions on Power Electronics, 25(4), 992- 1000. doi:10.1109/TPEL.2009.2036432 Zhai, Z., Martínez, J. F., Beltran, V., & Martínez, N. L. (2020). Decision support systems for agriculture 4.0: Survey and Challenges. Computers and Electronics in Agriculture, 170, 105256. doi:https://doi.org/10.1016/j.compag.2020.105256

89

VITA

Bo H. Smith was raised in Nortonville, KY. He has a B.S. from the University of Kentucky

in Biosystems and Agricultural Engineering, and he worked as both a Graduate Research

Assistant at the University of Kentucky in the Department of Biosystems and Agricultural

Engineering and a Teaching Assistant in the same department. He was very involved in

professional development societies such as the Biosystems and Agricultural Engineering

Graduate League of Student and the Alpha Epsilon Honors Society for Agricultural and

Biological Engineers, serving on the officer board of both. During his graduate tenure, he gave an oral presentation virtually at the American Society of Agricultural and Biological

Engineers’ (ASABE) 2020 Annual International Meeting.

Bo H. Smith

90