Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “behavior modeling”

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.

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

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic Modeling of Latent Heat Thermal Battery

A dynamic thermodynamic model is constructed in Modelica for the thermal battery portion of a heat pipe integrated thermal battery and is demonstrated using a shakedown test. The model is also leveraged to inform aspects of an experiment configuration. The theory section discusses modeling assumptions of the thermal battery system, establishes how melting and fusing is calculated, and how internal convection is modeled within the thermal battery. The shakedown test uses three charge-standby-discharge-standby cycles in which the charging cycles are 10 hours in length and the discharge cycles are 8 hours. System behavior including heat loss, temperature distributions and divergences, and internal convection behavior are presented. The model is then leveraged to show how experiment design can be informed through testing the model for insulation thickness impact on heat loss, round trip efficiency calculation, and heat tracing control.

25 ENERGY STORAGE↗
DOE OSTI2022

Modeling and Analysis of a Micro-reactor using the DireWolf Code Suite

Micro-reactors, specifically those cooled with heat pipes, are tightly coupled systems which require multi-physics tools to accurately model both steady-state and transient behavior. Idaho National Laboratory’s DireWolf code suite is tailor-built to model heat-pipe reactors. DireWolf’s ability to model the coupled steady-state behavior of a heat-pipe micro-reactor are demonstrated herein. The DireWolf modeling methodology is described, and a micro-reactor design is introduced. A full-core model is developed, and the steady-state coupled neutronics-thermal analysis is performed. Results of the analyses, including flux, power, and temperature distributions are presented and discussed. This is a Westinghouse Electric Company (WEC) led publication.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Sensitivity analysis for characterizing the impact of HNGD model on the prediction of hydrogen redistribution in Zircaloy cladding using BISON code

Hydrogen in zirconium cladding is able to precipitate into zirconium hydrides which impacts cladding integrity. The Hydride Nucleation-Growth-Dissolution (HNGD) model in the BISON code accounts for the precipitation and dissolution kinetics of hydride in Zircaloy material. This paper presents global sensitivity analyses of the HNGD model aiming to enhance our understanding of the hydride precipitation phenomena by quantifying the variance that key parameters have on the prediction of hydrogen behavior under various environmental conditions. Here, model predictions are compared to experimental data obtained under two different conditions: 1) with uniformly precharged specimens subjected to a linear thermal gradient, and 2) specimens precharged with a cathodically applied hydride rim at one end of the sample and subjected to an asymmetric thermal gradient. The Sobol sensitivity analysis identifies the key parameters in the HNGD model for both types of specimens. For linear temperature cases, the heat of transport dominates the accuracy of predictions when no precipitation occurs at the cold end, while Terminal Solid Solubility for Dissolution (TSSD) is the most important parameter when precipitation occurs. A large variation in the predicted hydrogen concentration profiles is found in the range of high TSSD due to the occurrence of precipitation. For asymmetric temperature cases, the solubility coefficient gives the largest impact on the predicted hydrogen distribution, as it determines the amount of solute hydrogen dissolved from the initially applied hydride rim. A large discrepancy in hydrogen distribution between simulations and experiments exists with the asymmetric specimens because BISON simulations fail to predict the precipitation of hydride at the cooler end. Comparative studies using former and updated models verifies the significant impact of the hydride growth mechanism on predicted hydrogen concentration profiles. In particular, when hydride initially exists, changes in TSSD generate a large variation in the predicted amount of precipitation by hydride growth, giving large uncertainty in predicting the hydrogen distribution over the sample length. The outputs characterize the significant impact of the hydride growth mechanism in the HNGD model on predicting hydrogen behavior, and improve the understanding of the precipitation of hydride in Zircaloy cladding within a range of expected environmental conditions. The analyses indicate work is still needed to improve the hydride solvus models in the BISON code to accurately predict experimentally observed hydride concentrations and distributions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Utilizing CellBox to Describe Varied Dynamical Systems

CellBox is a hybrid modeling framework that integrates machine-learning approaches with explicit mathematical modeling to determine the behavior of dynamical systems, intended for determining the complex interactions of molecule and phenotype within a cellular environment using perturbation/response-based data. The original work utilized analysis of a melanoma cell line, SK-Mel-133, which we sought to expand to and evaluate effectiveness in describing other nonlinear systems both within and outside of biology. We evaluated with two models. One model is based upon the free-fall of an object under a velocity-dependent drag force; this produced inherent perturbation/response data. The other model is based upon a form of mass-action molecular kinetics; for this, we ran the simulation to steady state in each species for each perturbation and select the steady state as input for CellBox. We determined that CellBox, outside of its original use case, can still strongly predict even non-biological dynamical systems in the small-perturbation regime, with r = 0.99977 to 0.85123, p < 0.05 for freefall. However, CellBox cannot accurately predict behavior within large-perturbation regimes, p >> 0.05 for freefall. In addition, we performed a naïve comparison of CellBox network structure to known network structure in the mass-action model to assess accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Simplified Numerical Approach to Characterize the Thermal Response of a Moving Bed Solar Reactor

Abstract Concentrated solar thermochemical storage in the form of a zero-emission fuel is a promising option to produce long-duration energy storage. Solar fuel is produced using a cavity reactor that captures concentrated solar radiation from a solar field of heliostats. In this paper, heat transfer model of a tubular plug-flow reactor designed and manufactured for a solar fuel production is presented. Experimental data collected from a fixed bed tubular reactor testing are used for model comparison. The system consists of an externally heated tube with counter-current flowing gas and moving solid particles as the heated media. The proposed model simulates the dynamic behavior of temperature profiles of the tube wall, gas, and particles under various gas flowrates and residence times. The heat transfer between gas–wall, solid particle–wall, and gas–solid particle is numerically studied. The model results are compared with the results of experiments done using a 4 kW furnace with a 150 mm heating zone surrounding a horizontal alumina tube (reactor) with 50.8 mm outer diameter and thickness of 3.175 mm. Solid fixed particles of magnesium manganese oxide (MgMn2O4) with the size of 1 mm are packed within the length of 250 mm at the center of the tube length. Simulation results are assessed with respect to fixed bed experimental data for four different gas flowrates, namely, 5, 10, 15, and 20 standard liters per minute of air, and furnace temperatures in the range of 200–1200 °C. The simulation results showed good agreement with maximum steady state error that is less than 6% of those obtained from the experiments for all runs. The proposed model can be implemented as a low-order physical model for the control of temperature inside plug-flow reactors for thermochemical energy storage applications.

Engineering↗
DOE OSTI2020

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗
DOE OSTI2024

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Data‐Driven Safety Risk Prediction of Lithium‐Ion Battery

Abstract Inevitable safety issues have pushed battery engineers to become more conservative in battery system design; however, battery‐involved accidents still frequently are reported in headlines. Identifying, understanding, and predicting safety risks have become priorities to further accelerate technology and industry development. However, diverse loading scenarios, significantly varied stress‐induced short circuit mechanisms, and highly coupled mechanical–electrochemical safety behaviors have remained grand challenges. Herein, the safety risk is termed as the probability of the mechanical triggering of an internal short circuit, to reflect the safety related behaviors of lithium‐ion batteries. Based on a mechanical model and experimental results, a sufficient dataset is generated consisting of strain states and their corresponding safety risks, covering both cylindrical and pouch cells, various states of charges, and loading conditions. Machine‐learning tools combined with the established finite element mechanical model are applied to predict the safety risks of the cells. The results achieve a high level of accuracy on the test data (the relative error of the average short circuit prediction deviation is less than 6.2%.). This work underpins the safety risk concept and highlights the promise of physics combined with data‐driven modeling methodology to predict the safety behaviors of energy storage systems.

Jia, Yikai↗
DOE OSTI2026

Microcanonical Kinetics of Water-Mediated Proton Transfer in 4ABAH + ·(H 2 O) n = 4–6 Clusters (ABA = Aminobenzoic Acid): A Model System for Size-Dependent Relaxation to Ergodic Behavior

Here, we leverage the unique properties of the 4ABAH + · (H 2 O) n clusters (ABA = 4-aminobenzoic acid, n = 4−6) to quantitatively address how a finite, isolated system evolves into an ergodic condition starting from localized arrangements in configuration space. This system adopts two distinct structural isomers in which water molecules cluster around the cationic centers of its two protomers with widely separated positive charge centers. These isomers arise from excess proton attachment to either the acid (O) or amino (N) group on opposite sides of the benzene ring. Both forms are captured and kinetically trapped using cryogenic ion methods and then selectively vibrationally excited through their mutually exclusive IR bands involving NH and OH stretching fundamentals. Because the IR excitation lies below the water binding energy, the system can evolve to explore slow, rare events that lead to the interconversion between the two isomers. The rates of these intracluster reactions are determined by using a pump−probe scheme involving ∼5 ns IR pump and UV probe lasers. The rates occur on the microsecond time scale, leading to steady state populations of the isomers, thus revealing the cluster size-dependent fractionation between the two species at microcanonical equilibrium. The steady state distributions are correlated with the expected trend in the cluster size-dependent reaction energetics, which are in turn consistent with changes in the relative densities of states of the two species. These results thus provide an unusually clear example in which complex, protic-solvent-mediated chemical transformations are captured within a finite system at a precisely determined internal energy.

Rana, Abhijit [Yale Univ., New Haven, CT (United S↗
DOE OSTI2020

Linking the Dynamic Response of the Carbon Dioxide-Concentrating Mechanism to Carbon Assimilation Behavior in Fremyella diplosiphon

ABSTRACT Cyanobacteria use a carbon dioxide (CO 2 )-concentrating mechanism (CCM) that enhances their carbon fixation efficiency and is regulated by many environmental factors that impact photosynthesis, including carbon availability, light levels, and nutrient access. Efforts to connect the regulation of the CCM by these factors to functional effects on carbon assimilation rates have been complicated by the aqueous nature of cyanobacteria. Here, we describe the use of cyanobacteria in a semiwet state on glass fiber filtration discs—cyanobacterial discs—to establish dynamic carbon assimilation behavior using gas exchange analysis. In combination with quantitative PCR (qPCR) and transmission electron microscopy (TEM) analyses, we linked the regulation of CCM components to corresponding carbon assimilation behavior in the freshwater, filamentous cyanobacterium Fremyella diplosiphon . Inorganic carbon (C i ) levels, light quantity, and light quality have all been shown to influence carbon assimilation behavior in F. diplosiphon . Our results suggest a biphasic model of cyanobacterial carbon fixation. While behavior at low levels of CO 2 is driven mainly by the C i uptake ability of the cyanobacterium, at higher CO 2 levels, carbon assimilation behavior is multifaceted and depends on C i availability, carboxysome morphology, linear electron flow, and cell shape. Carbon response curves (CRCs) generated via gas exchange analysis enable rapid examination of CO 2 assimilation behavior in cyanobacteria and can be used for cells grown under distinct conditions to provide insight into how CO 2 assimilation correlates with the regulation of critical cellular functions, such as the environmental control of the CCM and downstream photosynthetic capacity. IMPORTANCE Environmental regulation of photosynthesis in cyanobacteria enhances organismal fitness, light capture, and associated carbon fixation under dynamic conditions. Concentration of carbon dioxide (CO 2 ) near the carbon-fixing enzyme RubisCO occurs via the CO 2 -concentrating mechanism (CCM). The CCM is also tuned in response to carbon availability, light quality or levels, or nutrient access—cues that also impact photosynthesis. We adapted dynamic gas exchange methods generally used with plants to investigate environmental regulation of the CCM and carbon fixation capacity using glass fiber-filtered cells of the cyanobacterium Fremyella diplosiphon . We describe a breakthrough in measuring real-time carbon uptake and associated assimilation capacity for cells grown in distinct conditions (i.e., light quality, light quantity, or carbon status). These measurements demonstrate that the CCM modulates carbon uptake and assimilation under low-C i conditions and that light-dependent regulation of pigmentation, cell shape, and downstream stages of carbon fixation are critical for tuning carbon uptake and assimilation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗
DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Adaptive Stiffness Power Take-Off System for a Hydraulic Hose-Pump Point Absorber Wave Energy Converter

Two crucial challenges faced in the design of ocean wave energy converters (WECs) are maintaining high capacity factors and remaining robust in everchanging, harsh ocean environments. This study introduces a novel mechanism for achieving reactive WEC control with a soft hydraulic power take-off, thereby addressing both challenges. The device of interest involves a hydraulic pump, constructed from a fluidic flexible matrix composite (F2MC), which is cyclically stretched by driving waves and pressurizes working fluid through a turbine-generator. This class of fiber-reinforced stretch hoses is established for ocean use as robust moorings and shock absorbers. The dynamic properties of the F2MC pump depend on internal pressure, which can be controlled by a continuously variable transmission between the turbine and generator, allowing the system to be adapted for better performance over a wide range of wave conditions. This work presents a dynamic model for operation and adaptive behavior of this power take-off system, experimentally validates components of the model, and performs dynamic simulations of the system that show increases in power production and capacity factor compared to non-adaptive operation for a variety of driving wave conditions. This work develops the basis for future large-scale experimental testing and controller design of this system.

Wave energy conversion, fluidic flexible matrix co↗