SEARCH · Search NASA

Results for “STRESS AND LOAD”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides

Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗

DOE OSTI2026

Electrospun Ti–Zr Oxide Heterostructures Enable Strongly Anchored Ultralow-Ir Anodes for Durable Acidic Oxygen Evolution

Proton-exchange-membrane water electrolysis (PEMWE) requires acidic oxygen-evolution-reaction (OER) anodes that combine high activity, high durability, and low Ir loading. Here, we report a Ti-Zr composite electrospun oxide (ESO) nanorod support that enables ultralow-Ir anodes for high-performance PEMWE. Zr-containing Ti oxide heterostructures stabilize anatase-rich TiO2, tune the local oxygen-coordination environment, and strengthen interfacial anchoring of IrOx under acidic anodic conditions. The electrospun nanorod network further creates an open, mechanically coherent catalyst layer that improves Ir utilization, ionomer penetration, and mass transport. At an anode loading of 0.2 mgIr cm-2, the optimized Ir/TiZr20-ESO anode delivers a PEMWE mass activity of 0.99 A mgIr-1 at 1.45 V, 28.3 and 43.0 times higher than commercial Ir black and commercial IrO2/TiO2, respectively. The same anode reaches 3.0 and 4.0 A cm-2 at 1.75 and 1.83 V, respectively, and sustains 2000 h operation at 2.0 A cm-2. Also, accelerated stress tests up to 525 hours over 31,500 cycles confirm promising long-term durability, with an insignificant performance decay of 0.4 μV per cycle. Density functional theory indicates that the Ti-Zr oxide heterostructure suppresses Ti demetallation and strengthens IrO2 interfacial binding, rationalizing the improved high-current-density stability.

25 ENERGY STORAGE↗

DOE OSTI2025

Compression–tension cell with sample manipulator for in situ X‐ray nanotomography experiments

In situ X-ray nanotomography experiments where tensile or compressive force is applied on the sample require specialized equipment. A compression-tension device with fluid flow-through capability has been designed for X-ray nanotomography beamlines. The compression-tension cell is equipped with a triaxial stage for sample alignment and a high sensitivity loadcell for measurement of applied force. To handle the <100 µm samples used for X-ray nanotomography imaging and for loading samples on the compression-tension cell a sample manipulator has been built. The sample manipulator is capable of selecting a single <100 µm particle for nanotomography scanning while viewing multiple samples under an optical microscope. To test the functionality of these two devices an initial compression experiment involving two glass beads was performed. To demonstrate instrument stability two spherical glass beads were compressed from a no load condition until one of the beads fractured. Nanotomography data were collected at each step of increasing compressive force. The experimentally observed contact area of the spherical glass beads was compared with the theoretical estimate using the Hertz analysis. To demonstrate the fluid flow capability, two calcite grains were compressed against each other under a calcite saturated solution. Surface topological changes were observed for the stressed grain contact area.

X-ray tomography↗

DOE OSTI2026

Modeling Framework for Data Center

This chapter highlights the critical need for advanced modeling of data centers due to their rapidly increasing energy consumption and impact on grid reliability. Driven by the demand for AI applications, data centers are projected to consume a significant portion of US energy by 2028, putting stress on an already challenged power grid. The chapter emphasizes the importance of "fast" time-scale models to understand the dynamic interactions between data centers and the grid, especially given the rapid power fluctuations of AI workloads. It outlines a modeling framework that includes both offline and real-time EMT domain simulations, detailing the necessary representations for various components like utility interfaces, transformers, IT loads, UPS, cooling loads, Battery Energy Storage Systems (BESS), generators, protection systems, and higher-level control systems. While standard simulation tools like PSCAD offer basic models, custom development is often required to accurately capture the unique and fast-changing behaviors of modern data centers. The chapter also discusses key metrics and test cases for validating these models, focusing on transient load responses, protection relay coordination, and demand flexibility. Finally, it addresses the challenges of modeling large-scale data centers, such as computational complexity and the trade-off between model fidelity and practicality, suggesting hybrid modeling approaches as a solution. The overarching goal is to create a robust framework that helps assess data center impacts on grid stability, identify vulnerabilities, and inform the development of standards for reliable integration of these large loads into the bulk power system.

25 ENERGY STORAGE↗

DOE OSTI2024

Friction stir processing on a strontium modified, thin-wall, vacuum-assisted high-pressure die-cast Aural-5 alloy to improve tensile and fatigue performance

Here, this study explores the application of friction stir processing (FSP) to enhance the material properties of Sr-modified Aural-5 alloy, with a focus on improved tensile and fatigue properties. Aural-5 is a well-known vacuum-assisted high-pressure die-cast (HPDC) Al-Si7-Mg alloy used in the automotive industry to reduce vehicle weight, enhance fuel efficiency, and lower carbon emissions. This alloy modifies its material chemistry with Sr for fine fibrous networks of eutectic silicon and manganese (Mn) to reduce die soldering. It has significantly less iron (Fe) content resulting in the elimination of detrimental needle-shaped Fe-bearing ß-phase intermetallic and improving ductility. The initial microstructure of as-received HPDC Aural-5 exhibits shrinkage porosity in the middle section, a dendritic microstructure with fibrous Al-Si eutectic colonies, a shear-band structure beneath the die-wall, large dendritic externally solidified crystals (ESCs), needle-shaped Mg 2 Si phase and significant second-phase particulates. Some of those microstructural features, such as porosity, ESCs, needle-shaped Mg 2 Si phase, and large second-phase particles, serve as initiation sites for cracks under mechanical loading, resulting in adverse effects on tensile properties, particularly ductility. FSP effectively transforms the microstructure into a wrought configuration with uniform particle distribution by eliminating porosity and disintegrating dendrites, eutectic colonies, ESCs, second-phase particles, and shear-band structures. FSP-driven microstructure modification enhances yield strength and tensile ductility by ~30% and ~35%, respectively. The fatigue life of the material in a bending mode configuration (stress ratio R = 0.1) after FSP exhibits enhancements ranging from 2.0 to 3.9 times that of the original HPDC Aural-5 alloy, depending on the applied stress level.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Potential of Data Center Controls in Grid Services

The rapid proliferation of large data centers brings both challenges and opportunities for grid reliability. The data center resources and their potential flexibility have the potential to contribute resources to grid operations. Through capabilities like energy shifting and resource coordination, data centers can help reduce their net demand on the transmission network, as well as provide additional grid services to support reliable operation on the grid. While transient and long-term grid planning and operations are the scenarios that draw most attention, the quasi-steady state timeseries (QSTS) operation of data centers and grid bring interesting scenarios that can help evaluate the data center controls to aid grid services. This work is focused on modeling data centers for QSTS applications – incorporating the AI data center load profiles and building on the PNNL digital twin model for the thermal management loads to enable simulation studies to reveal the impact of data center controls on grid performance. This includes the integration of a QSTS battery and natural gas generator model to incorporate local resource impacts to the system. The simulation study is performed with a modified IEEE 24-Bus transmission system. Scenarios are focused on evaluating the data center load impacts on the transmission system and leveraging both data center and local generation controls to mitigate those impacts and provide additional grid services. The data center controls revealed the ability to contribute to two main kinds of grid services: preventing congestion on a weak grid by coordinating the data center resources with the collocated BESS and onsite generation; and the ability to help the grid operations during stressed times of operation like during a contingency. Leveraging these and other capabilities has the potential to help data centers become grid responsive assets, aiding in both their integration into the power system and grid reliability.

power grid simulation↗

DOE OSTI2025

Determining anisotropic slip system rate sensitivities of Ti-6Al-4V using high-energy X-ray diffraction microscopy

Variation of strain rate sensitivity among different families of slip systems in the hexagonal close-packed (α) phase of titanium (Ti) alloys has the potential to alter microscale load redistribution during both creep and dwell fatigue loading. However, existing literature contains conflicting reports regarding the degree of anisotropy present in the strain rate sensitives between α slip system families across Ti alloys. Here, we quantify the strain rate sensitivity of α slip system families in Ti-6Al-4V using high-energy X-ray diffraction microscopy (HEDM). We present a novel procedure in which we utilize the HEDM-measured grain-scale stress states during stress relaxation to determine the strain rate sensitivities of the basal $\langle a \rangle$ , prismatic $\langle a \rangle$, and first-order pyramidal $\langle c + a \rangle$ slip system families. In addition, we measure the strain rate sensitivity exponent of the macroscopic response and the α phase for comparison. We find rate sensitivities for the different slip system families to range from 0.02 to 0.04, which—while varied—are relatively isotropic in comparison with some values presented in the literature. We also demonstrate the effects of the measured anisotropic rate sensitivities using crystal plasticity finite element simulations.

Peterson, Kenneth M. [Pennsylvania State Univ., Un↗

DOE OSTI2024

Performance and Durability of Heavy-Duty Fuel Cell Systems with an Advanced Ordered Intermetallic ORR Alloy Catalyst and Novel Support

Ordered PtCo intermetallic (OIM) catalyst (L1 0 -PtCo/C) is a promising candidate as the oxygen reduction reaction (ORR) catalyst in hybrid fuel cell systems (FCS) for class-8 heavy duty (HD) trucks. Compared to a baseline annealed Pt on high surface area carbon (a-Pt/HSC) catalyst, its mass activity (MA) is 71% higher initially and 144% higher after 90,000 potential cycles in an accelerated stress test (AST). Analysis of the AST data indicates that the ORR kinetic constants do not change with aging and the degradation in the OIM catalyst activity is linearly proportional to the loss in the electrochemically active surface area (ECSA). Several operational strategies are investigated to mitigate catalyst degradation and achieve 25,000-h electrode lifetime and 2.5 kW g −1 Pt utilization on a HD truck duty cycle including load sharing with the hybrid battery, regulating the radiator fan power to maintain the coolant temperature close to 60 °C, clipping the maximum cell voltage below 850 mV, limiting the ECSA loss to 55%, and oversizing the active area of the membrane electrode assemblies by 20%. Drive cycle simulations indicate that the lifetime average voltage degradation rate is about 1.8 μV h −1 and the integrated stack and FCS drive cycle efficiencies decrease by 3.5 to 3.9%.

25 ENERGY STORAGE↗

DOE OSTI2024

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be time-consuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the tower-base damage equivalent loading and better tracking of the rated generator speed and power.

axial flow↗

DOE OSTI2025

FY25 Status Update on A709 Code Case Testing at ANL

This report provides an update on the status of the creep rupture testing on the precipitation-treated (PT) Alloy 709 samples fabricated from the first, the second and the third commercial heats, in support of the first, 100,000-hour Alloy 709 Code Case development under the American Society of Mechanical Engineers (ASME) Section III, Division 5 for high-temperature reactor construction. In Fiscal Year (FY) 25, no new tests were initiated, while 5 previously loaded tests were completed. This report presents an overview of the creep-rupture data gathered at Argonne National Laboratory (ANL) since FY21 that totaled 52 completed tests. Metallographic images showing crack morphologies in some representative samples ruptured in FY25 are also included. Key conclusions include: 1) the creep behavior of the three heats of Alloy 709 can be well described by the Larson-Miller relationship; 2) higher temperatures and larger stresses lead to shorter creep life and higher minimum creep rate; 3) Alloy 709 is overall very ductile under creep; 4) the crack morphologies near fracture surfaces are dependent on material’s grain size and test conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling

Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

AGC-4 Specimen Post-Irradiation Examination Data Interim Report

This interim report documents preliminary results of the post-irradiation examination material property testing from the fourth advanced graphite creep (AGC), AGC 4, capsule specimens. This is the fourth of a series of six irradiation test trains planned as part of the AGC experiment to fully characterize the neutron irradiation effects and radiation creep behavior of current nuclear graphite grades to moderate dose levels (=7 dpa). The AGC 4 capsule was irradiated in the Idaho National Laboratory Advanced Test Reactor at a nominal temperature of 800°C and to a peak dose of 8 dpa. Half of the AGC-4 specimens were subjected to compressive stresses to induce irradiation creep. Post-irradiation testing and measurement results are reported with the exception of thermal testing, which is still in progress, and irradiation mechanical strength testing. Additionally, some specimens initially deemed too hot to be examined in the ART Graphite laboratory may still be measured. The data reported includes specimen dimensions for both stressed and unstressed specimens to establish the irradiation creep rates, mass and dimensional data necessary to derive density, elastic constants (Young?s modulus, shear modulus, and Poisson?s ratio) from ultrasonic time of flight velocity measurements, Young?s modulus from the fundamental frequency of vibration, and electrical resistivity. A more complete evaluation of trends in the material property changes, as well as irradiation-induced creep due to the irradiation environment and applied load on the specimens, will be discussed later in AGC 4 post-irradiation examination analysis reports.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

AGC-4 Specimen Post-Irradiation Examination Data Interim Report

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗

DOE OSTI2024

Glass production rate in an electric melter: Melting rate correlation and primary foam stability

Waste loading and melting rate are major factors determining the lifecycle of nuclear waste disposal by vitrification in electric melters in which the heat to the cold cap is delivered from the melt pool. One of the crucial parameters that affect the rate of melting is the temperature at which the primary foam collapses at the cold cap bottom. Further, apart from the transient glass-forming melt viscosity, primary foam stability is affected by the shear stress imposed by the flow in the melt pool and the presence of residual solid particles.

36 MATERIALS SCIENCE↗

DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Community structure and function during periods of high performance and system upset in a full-scale mixed microalgal wastewater resource recovery facility

Microalgae have the potential to exceed current nutrient recovery limits from wastewater, enabling water resource recovery facilities (WRRFs) to achieve increasingly stringent effluent permits. The use of photobioreactors (PBRs) and the separation of hydraulic retention and solids residence time (HRT/SRT) further enables increased biomass in a reduced physical footprint while allowing operational parameters (e.g., SRT) to select for desired functional communities. However, as algal technology transitions to full-scale, there is a need to understand the effect of operational and environmental parameters on complex microbial dynamics among mixotrophic microalgae, bacterial groups, and pests (i.e., grazers and pathogens) and to implement robust process controls for stable long-term performance. Here, we examine a full-scale, intensive WRRF utilizing mixed microalgae for tertiary treatment in the US (EcoRecover, Clearas Water Recovery Inc.) during a nine-month monitoring campaign. We investigated the temporal variations in microbial community structure (18S and 16S rRNA genes), which revealed that stable system performance of the EcoRecover system was marked by a low-diversity microalgal community (D INVSIMPSON = 2.01) dominated by Scenedesmus sp. (MRA = 55 %-80 %) that achieved strict nutrient removal (effluent TP < 0.04 mg·L -1 ) and steady biomass concentration (TSS monthly avg . = 400–700 mg·L −1 ). Operational variables including pH, alkalinity, and influent ammonium (NH 4 + ), correlated positively (p < 0.05, method = Spearman) with algal community during stable performance. Further, the use of these parameters as operational controls along with N/P loading and SRT allowed for system recovery following upset events. Importantly, the presence or absence of bacterial nitrification did not directly impact algal system performance and overall nutrient recovery, but partial nitrification (potentially resulting from NO 2 − accumulation) inhibited algal growth and should be considered during long-term operation. The microalgal communities were also adversely affected by zooplankton grazers (ciliates, rotifers) and fungal parasites (Aphelidium), particularly during periods of upset when algal cultures were experiencing culture turnover or stress conditions (e.g., nitrogen limitation, elevated temperature). Altogether, the active management of system operation in order to maintain healthy algal cultures and high biomass productivity can result in significant periods (>4 months) of stable system performance that achieve robust nutrient recovery, even in winter months in northern latitudes (WI, USA).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗