Geology, gravimetric and numerical modeling of the Nova Collinas impact structure, Parnaiba Basin, Brazil
Explore the source record for details and available documents.
Avona, P. D.↗
SEARCH · Engineering Papers
Results for “numerical modelling”
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.
Influence of Numerical Modeling Approaches on Damped Behavior of Flexible Beams: Preprint
Composites structures are widely used in aerospace and wind energy applications for their excellent stiffness and strength-to-weight properties. In these structures, structural damping is critical to predict vibration amplitudes, performance, and reliability. Structural damping is of particular interest for slender wings, rotorcraft blades, and wind turbine blades that can exhibit complex vibration phenomena and are frequently modeled with geometrically exact beam theory (GEBT). Standard approaches of stiffness proportional or modal damping merely assign user defined values and cannot predict damping behavior. This work compares stiffness proportional damping to two more advanced damping approaches: modal strain energy and Prony series. The modal strain energy approach uses a sectional analysis tool to calculate the beam stiffness and postprocess internal stresses from GEBT simulations. The internal stresses are then used to calculate modal damping factors. The Prony series is implemented within GEBT to directly model viscoelastic behavior of the composites. These approaches are compared by modeling the evolution of the damping factors of a realistic flexible wind turbine blade with varying rotational speed. Discrepancies between the approaches suggest areas for future modeling development, but differences in nonlinear damping values are less than current uncertainties about the magnitude of structural damping.
17 WIND ENERGY↗
A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage
The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.
14 SOLAR ENERGY↗
Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas
The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.
54 ENVIRONMENTAL SCIENCES↗
Development of Cr, Se, U, Sb, and Te isotopes as indicators of redox reactions, contaminant fate, and contaminant transport in aqueous systems: A review
Cr, Se, U, Sb, and Te are toxic, redox-active elements that are more mobile and environmentally problematic in their oxidized forms, and less mobile and bioavailable in their reduced forms. This chapter reviews the development of Cr, Se, U, Sb, and Te isotope ratio measurements as new indicators of redox reactions and contaminant migration. Reliable analytical methods exist, but are still evolving. Understanding of isotopic fractionation induced by various (bio)geochemical processes has been explored in dozens of publications, yet is far from complete: Reduction reactions, the major driver of isotopic variation, have been relatively well studied. However, the magnitude of fractionation is variable and the systematics of that variation are still being explored. Isotopic fractionation induced by oxidation reactions is not well understood. Non-redox reactions, which involve smaller changes in bonding of these elements, tend to induce less isotopic fractionation, but can nonetheless cause significant isotopic shifts. Field applications of Cr, Se, U isotope ratios have demonstrated that they are useful as indicators of reduction in natural systems. A few studies suggest they are also useful as indicators of oxidation and contaminant sources. The physical and chemical complexity of groundwater systems hinders accurate quantitative interpretation of Cr, Se, U isotope data using simple models. Numerical models have been developed that capture the behavior of complex, coupled systems and enable the most effective extraction of information from field data sets.
58 GEOSCIENCES↗
Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models
Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.
54 ENVIRONMENTAL SCIENCES↗
Comparative investigations of multi-fidelity modeling on performance of electrostatically-actuated cracked micro-beams
Silicon is a commonly used material for the fabrication of beams for use in micro-electrical-mechanical systems (MEMS). Although silicon is a brittle material, it has been shown to accumulate fatigue damage at the micro-scale. Understanding the effect this has on the overall device performance is critical to the design of reliable devices. Analytical methods for modeling damage provide expedient results but are limited by broad modeling assumptions. Numerical models account for more detailed physical phenomena but can be computationally intensive. In this work, two different crack scenarios are modeled using both analytical techniques and 3D computational simulations. First, the effects of a single surface crack on the static deflection and natural frequency of an electrostatically actuated micro-beam are formulated and compared. Then, a new method for approximating damage associated with realistic distributed crack networks is formulated for use in an analytical model and numerical simulations. A method for utilizing experimentally derived crack statistics to inform the analytical and numerical distributed crack models is developed. Good agreement between the analytical and numerical models is obtained for both crack scenarios. Altogether, these models can be used to effectively simulate a variety of damage and fatigue behaviors in silicon-based MEMS devices.
42 ENGINEERING↗
CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications
Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.
Murdoch, Larry↗
Code-to-code comparison between FLASH and HYDRA in gas-puff Z-pinch modeling
The numerical modeling of gas-puff Z pinches involves the nonlinear coupling of a broad range of complex, multi-physics phenomena that makes such simulations challenging. The challenge is further compounded by nonlinear instabilities that can impact the dynamics of imploding gas-puff Z pinches, such as the magneto Rayleigh–Taylor instability (MRTI). If the growth rate and amplitude of the latter is comparable to the relevant timescales and properties of the imploding plasma, the MRTI can dramatically alter implosion dynamics, dictate pinch stability, and govern the plasma properties achievable in pulsed-power-driven laboratory experiments. National Laboratories and academic teams have developed numerical tools that can accurately model Z-pinch configurations and provide reliable design capabilities that can guide experimental choices and assist in interpreting experimental results. Most such tools, however, are not broadly available. Here, we present newly developed Z-pinch simulation capabilities of the publicly available FLASH code, applied in the study of MRTI growth and dynamical effects in gas-puff implosions. To verify the new implementations, we perform a comparison of FLASH gas-puff implosion simulations with previously published calculations with the HYDRA code from Lawrence Livermore National Laboratory, which have been validated with experimental data from the CESZAR pulsed-power driver at the University of California, San Diego. The experiments involved double- and triple-nozzle configurations, in an experimental attempt to stabilize the pinch to the MRTI. The code-to-code comparison shows similar results between the FLASH and HYDRA simulations, supporting the use of FLASH in the modeling of future gas-puff Z-pinch experiments at CESZAR.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Optimizing Prediction Error for Time-dependent Solar Radiation Modeling
Numerical weather forecasting models and statistical methods have found wide use to help power companies estimate renewable output, but better methods are needed, particularly for extended forecasts. Machine learning approaches have been used here as well, but so far a major limitation is the ability to also predict the corresponding uncertainty in a forecast. Here we show that both can be done and demonstrate this using a long-term-short memory neural network where the difference between predicted and ground truth data are used to train a model for the corresponding forecast uncertainties.
97 MATHEMATICS AND COMPUTING↗
System and method for wave prediction
A method and system for prediction of wave properties include collecting time-series data streams from one or more wave measurement devices and processing the data to identify data parameters to establish boundary conditions of a numerical model. The numerical model may be used to compute a predicted wave field of time-series data for a variety of wave properties at a target location.
Previsic, Mirko↗
Proppant Transport and Coverage in Rock Fractures –A Computational Modeling Approach
Numerical Modeling of Proppant Transport and Coverage in Rock Fractures
Rousta, Farid↗
Machine learning based algorithms for uncertainty quantification in numerical weather prediction models
Complex numerical weather prediction models incorporate a variety of physical processes, each described by multiple alternative physical schemes with specific parameters. The selection of the physical schemes and the choice of the corresponding physical parameters during model configuration can significantly impact the accuracy of model forecasts. There is no combination of physical schemes that works best for all times, at all locations, and under all conditions. It is therefore of considerable interest to understand the interplay between the choice of physics and the accuracy of the resulting forecasts under different conditions. This paper demonstrates the use of machine learning techniques to study the uncertainty in numerical weather prediction models due to the interaction of multiple physical processes. The first problem addressed herein is the estimation of systematic model errors in output quantities of interest at future times, and the use of this information to improve the model forecasts. The second problem considered is the identification of those specific physical processes that contribute most to the forecast uncertainty in the quantity of interest under specified meteorological conditions. In order to address these questions we employ two machine learning approaches, random forests and artificial neural networks. The discrepancies between model results and observations at past times are used to learn the relationships between the choice of physical processes and the resulting forecast errors. Numerical experiments are carried out with the Weather Research and Forecasting (WRF) model. The output quantity of interest is the model precipitation, a variable that is both extremely important and very challenging to forecast. The physical processes under consideration include various micro-physics schemes, cumulus parameterizations, short wave, and long wave radiation schemes. The experiments demonstrate the strong potential of machine learning approaches to aid the study of model errors.
97 MATHEMATICS AND COMPUTING↗
Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)
In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.
03 NATURAL GAS↗
Dynamic responses of submerged components in advanced reactors: experimental and numerical studies
The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.
Engineering↗
A numerical study of process complexity in permafrost dominated regions
Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.
Calibration↗
Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity
Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .
79 ASTRONOMY AND ASTROPHYSICS↗
Wave Energy Converter Power Take-Off Modeling and Validation From Experimental Bench Tests
This article describes the implementation of a new numerical model of the power take-off system installed in the Monterey Bay Aquarium Research Institute wave energy converter, a device developed to provide power to various oceanic research missions. The simultaneous presence of hydraulic, pneumatic, and electrical subsystems in the power take-off system represents a significant challenge in forging an accurate model able to replicate the main dynamic characteristics of the system. The validation of the new numerical model is addressed by comparing simulations with the measurements obtained during a series of bench tests. Data from the bench tests show good agreement with the numerical model. The validated model provides deeper insights into the complex nonlinear dynamics of the power take-off system and will support further performance improvements in the future.
16 TIDAL AND WAVE POWER↗