Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Results”

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 451 records · Page 25

Developing an experimental-computational framework to investigate the deformation mechanisms and mechanical properties of Al-8Ce-10Mg alloys at micro and macroscales

There is a promising future for the use of aluminum-cerium-magnesium alloys in a broad range of applications, including devices that operate at high temperatures. With cerium currently considered a waste product of rare earth mining validation studies of possible applications are essential to reduce the environmental waste. Here, a computational-experimental framework is developed to investigate the role both the intermetallic and matrix have on the mechanical properties of these alloys. A set of experiments, including SEM/EBSD imaging, nanoindentation, in-situ SEM tensile testing, and in-situ SEM-DIC tests are performed to characterize the microstructure and mechanical properties of these alloys. Furthermore, the elastic, plastic, and failure deformation mechanisms of the microstructure, and their correlation with the bulk scale mechanical properties are investigated. Experimental results are also used to calibrate parameters for a crystal plasticity finite element model, by performing a computational framework that minimizes the error between the computational and experimental results. This model is then utilized to investigate how the area percentage of intermetallics and the crystallographic texture control the mechanical properties of the alloy. Simulation results show that an increase in the percentage of intermetallics increase the strength but decrease the ductility of the alloy. Also, a change in material texture improves strength and reduces damage that leads to material failure. The development of the crystal plasticity model, as discussed in this work, opens opportunities for future investigations of similar aluminum-cerium-magnesium alloys.

36 MATERIALS SCIENCE↗

Rigid-flexible coupling multi-body dynamics modeling of a semi-submersible floating offshore wind turbine

A 14-degree-of-freedom (14-DOF) flexible multibody dynamics model is developed and verified for a semi-submersible floating offshore wind turbine (FOWT). The model considers the coupled dynamics of the platform, tower, nacelle, blades, and mooring subjected to external wind and wave loads. The platform is simplified as a rigid body connected to the seabed by the mooring system. The tower and blade are modeled as flexible cantilever beams. The energy method is used to derive the governing equations of motion, where the kinetic energy, potential energy, and work done by external wind-wave forces are all deduced in a global coordinate system. The 5-MW baseline semi-submersible FOWT is used to verify the derived model against the results simulated from FAST developed by the US National Renewable Energy Laboratory (NREL) at two scenarios: free decay state and different wind-wave load cases. Results show that the established FOWT model can well reflect the vibration characteristics of FAST model. The application of the model to the control of platform pitch with a tuned mass damper is studied. Furthermore, the simplified model could provide a low-order method for the structural dynamics analysis and advanced vibration control design for the multi-body components of the semi-submersible FOWT in the future.

17 WIND ENERGY↗

Accelerated Zymonic Acid Formation from Pyruvic Acid at the Interface of Aqueous Nanodroplets

To explore the role of the liquid interface in mediating reactivity in small compartments, the formation kinetics of zymonic acid (ZA) is measured in submicron aerosols (average radius = 240 nm) using mass spectrometry. The formation of ZA, from a condensation reaction of two pyruvic acid (PA) molecules, proceeds over days in bulk solutions, while in submicron aerosols, it occurs in minutes. The experimental results are replicated in a kinetic model using an apparent interfacial reaction rate coefficient of krxn = (0.9 ± 0.2) × 10-3 M-1 s-1. The simulation reveals that surface activity of PA coupled with an enhanced interfacial reaction rate drives accelerated ZA formation in aerosols. Experimental and simulated results provide compelling evidence that the condensation reaction of PA occurs exclusively at the aerosol interface with a reaction rate coefficient that is enhanced by 4 orders of magnitude (∼104) relative to what is estimated for macroscale solutions.

Kim, Pyeongeun↗

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗

Mesh Failure Prediction Using Deep Learning Techniques

The goal of this joint effort between SRI and LLNL has been to investigate the use of state-of-the-art deep learning architectures to effectively explore the correlations between a set of mesh zone attributes and manifestations of failure in a given mesh simulation and draw useful conclusions towards a practical failure detection system. Each zone attribute is a scalar feature that is recorded separately for each zone during simulation, resulting in a multivariate time series signal for each zone. Currently, these features are used to manually decide where to apply manual relaxations so that simulations run without any failures. Our long-term vision is to build a machine learning system that is capable of reliably detecting upcoming failures and is able to apply finer-grained targeted relaxations automatically to relevant zones to avoid an expected failure.

97 MATHEMATICS AND COMPUTING↗

Testing of Helium-Cooled Metal Molybdenum Disk Target

NorthStar Medical Radioisotopes LLC is planning to produce an important medical radioisotope, molybdenum-99 (Mo-99), through photonuclear reaction on molybdenum-100 (Mo-100). In this approach, molybdenum metal will be bombarded with a 40-MeV electron beam. Because enriched Mo-100 is expensive, it is desired to use as much beam power as possible to achieve maximum production yield and minimize target mass. This objective leads to very high beam power density (heat deposition in the target), which sets challenging requirements for cooling. Together with scientists at Los Alamos National Laboratory, a team at Argonne National Laboratory has developed and demonstrated a cooling approach using pressurized helium, which allows for efficient heat removal. One of the main challenges in this approach is the management of the heat load on the target window. The target window separates the high-pressure helium inside the target from the vacuum in the beamline, so it is constantly under stress from differential pressure. Also, the window is cooled only by the helium gas flowing on one side, making the window cooling even more challenging. High heat deposition in the target disks also imposes a strict requirement on performance of the helium cooling system and thickness of the target disks. The target disks are produced from metal powder via a press-and-sinter process. The resulting disks do not possess tensile strength as high as solid molybdenum and might not survive the vibration from the high-velocity helium coolant and high thermal stress from beam heating. An Argonne team of scientists performed a series of tests at Argonne’s Low Energy Accelerator Facility (LEAF) [3-7]. This report describes two series of tests for scale down production target designs that utilize full-scale 29 millimeters diameter, 0.75 mm thick press-and-sintered disks. We performed two thermal tests with different beam parameters and configurations of the disk laminations. We compared the results of the window temperature measurements and cooling system parameters obtained in the experiments with those predicted by analytical calculations and Computation Flow Dynamic (CFD) simulations. Results of the experiments and calculations are presented below.

36 MATERIALS SCIENCE↗

DC-Link Capacitor Design for a Neutral-Point-Less Three-Level Dual-Phase Inverter for Traction Application

Conventional multilevel inverters, such as neutral￾point-clamped and T-type inverters, have gained popularity in electric vehicle applications due to their advantages, including high voltage range, high power capability, low switching losses, low total harmonic distortion, and low electromagnetic interference. However, these traditional multilevel inverters require a neutral point connection to generate a zero-voltage vector. The neutral point current oscillates at three times the fundamental frequency, leading to voltage imbalance and overvoltage stress on the power modules. Additionally, the use of two stacked DC-link capacitors increases the volume required for the same overall capacitance and complicates packaging due to ripple current and heat dissipation from separate components. This is a significant concern for traction drive units, where space is limited. In this paper, a neutral-point-less multilevel dual three￾phase inverter topology is investigated for traction inverter applications. Simulation results demonstrate that the proposed topology effectively retains the benefits of multilevel operation while utilizing a single DC-link capacitor. The inverter model was simulated in conjunction with an industry-standard battery model to evaluate the potential for capacitor size reduction compared to a conventional three-level inverter.

33 ADVANCED PROPULSION SYSTEMS↗

Understanding the Charging Flexibility of Shared Automated Electric Vehicle Fleets

The combined anticipated trends of vehicle sharing (ride-hailing), automated control, and powertrain electrification are poised to disrupt the current paradigm of predominately owner-driven gasoline vehicles with low levels of utilization. Shared, automated, electric vehicle (SAEV) fleets offer the potential for lower cost and emissions and have garnered significant interest among the research community. While promising, unmanaged operation of these fleets may lead to unintended negative consequences. One potentially unintended consequence is a high quantity of SAEVs charging during peak demand hours on the electric grid, potentially increasing the required generation capacity. This research explores the flexibility associated with charging loads demanded by SAEV fleets in response to servicing personal mobility travel demands. Travel demand is synthesized in four major United States metropolitan areas: Detroit, MI; Austin, TX; Washington, DC; and Miami, FL. In each of these four cities, SAEV simulations are performed using local projected electricity prices from the Regional Energy Deployment System (ReEDS) for a handful of supply side scenarios. A coordinated charging model is introduced that seeks to reduce fleet charging costs in response to time-varying electricity prices and increasing shares of personal trip demand served (between 1% and 25% of all metro trips served by the SAEV fleet). Simulation results are presented across various scenarios assuming fleetwide coordination to minimize charging energy costs while constrained by offering adequate mobility service to fleet customers. The results indicate that the SAEV charging load is highly flexible; energy costs were shown to reduce between 13% and 46% across a wide range of simulated scenarios. In addition, these savings were realized without detrimentally impacting the fleet’s ability to service trips.

autonomous vehicle↗

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗

Quasi-brittle porous material: Simulated effect of stochastic air void structure on compressive strength

Highlights: • A simulation procedure to study the effect of porosity on the compressive strength of a quasi-brittle material is introduced • Based on the simulated results, the relationship between compressive strength and air void fraction (porosity) is presented. • A linear relationship was achieved between the cubic root of porosity and the simulated compressive strength. The effect of porosity comprised of spherical air voids on the compressive strength of quasi-brittle material was studied via simulations. The simulated porous structures were based on pore size distributions of two mortar samples measured by X-ray microtomography. While the simulation method set practical limits on the size of sample, the base of the statistics was established by simulating 128 small structures generated by sampling from pore structures of two mortars. By studying the application of the classical strength-porosity formulas to the simulated data, a new simple model was formed. A linear relationship was achieved between the cubic root of air void fraction (porosity) and the simulated compressive strength. The reasons for scattering of simulated strength around fitted trend remained unresolved in this study; no clear dependence on pore number or other distribution properties was observed. With the presented simulation approach, the dependence of compressive strength on porosity is achieved independently of disturbances that occur in experimental studies creating understanding of compressional behavior of low porosity materials.

36 MATERIALS SCIENCE↗

The role of focusing geometry in MeV x-ray production from petawatt laser–solid interaction

Relativistic laser–plasma interactions provide a compact and flexible route to generating bright, ultrashort pulses of MeV x rays, with applications in high-energy density science, nuclear physics, and radiography. Despite extensive study of intensity scaling in laser-driven electron acceleration, the role of focusing geometry and focal-volume effects in MeV radiation production remains insufficiently understood. Here, we present experimental and kinetic simulation results from the Texas Petawatt Laser (120 J, 140 fs) in which the focusing geometry (f/3 or f/1.5) is varied, while the laser energy and pulse duration are fixed. Experimentally, the f/3 geometry produces approximately three times more MeV radiation than the f/1.5 geometry, despite its twofold lower nominal vacuum intensity. This result is based on the time-integrated radiation per unit solid angle, as measured along the diagnostic line of sight. Three-dimensional particle-in-cell simulations reproduce this trend when a modest (10s of μm) effective focal plane shift is introduced, demonstrating that relativistic laser–plasma coupling is highly sensitive to focal geometry in the presence of a preplasma. This behavior is consistent with differences in interaction length and effective intensity at the critical surface between the two configurations. The sensitivity of the tightly focused f/1.5 configuration ($z_R ≈$⁠ 5 μ m) reflects the combined influence of thermal lensing and an extended preplasma, while the f/3 geometry (⁠$z_R ≈$ 80 μ m) remains comparatively robust. These results demonstrate a breakdown of conventional intensity scaling and identify focusing geometry as a critical control parameter for MeV electron and x-ray generation at petawatt powers. Implications for laser–plasma-based radiographic facilities are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

This data package contains input files for TOUGHREACT for a modeling study examining the effects of managed aquifer recharge on agricultural lands on nitrate cycling and transport in the Central Valley of CA near Modesto. The files contain all the geochemical species, reactions, and hydrological parameters for the model. The files are text files used for the TOUGH family of code created by LBNL. To use the files a license is required. https://tough.lbl.gov/licensing-download/toughreact-licensing-download/Accompanying Paper Abstract: Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO3-) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying water all-at-once rather than in increments for a fixed volume of recharge transports higher concentrations of NO3- deeper into the profile, which has the potential to exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as riparian corridors, floodplains, wetlands, and other managed landscapes like dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs

This paper presents a mooring system design featuring polymer springs for the VolturnUS-S 15-MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modeling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15-MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, whereas the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads - on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗

Design and Analysis of a Floating-Wind Shallow-Water Mooring System Featuring Polymer Springs: Preprint

In this paper, a mooring system featuring polymer springs is designed for the VolturnUS-S 15 MW reference floating wind turbine in site conditions for the New York Bight at a 50-m water depth. Polymer springs have a nonlinear stress-strain curve that allows a stiffer response at low loads and a more flexible response at higher loads, potentially reducing peak mooring line tensions. The mooring dynamics model MoorDyn has been extended to model springs with nonlinear tension-strain curves from a user-inputted look-up table. This MoorDyn modelling advancement is verified against OrcaFlex simulation results. Using MoorDyn's updated capabilities, a spring-equipped catenary mooring system is designed for the 15 MW floating system, along with a baseline catenary mooring system that does not use polymer springs. The floating wind turbine simulator OpenFAST is used to simulate the mooring systems in design-driving load cases to show the effect of polymer springs on key dynamic behaviours. The results show that the spring-equipped design reduces peak tensions by up to 60%, while the turbine offsets stay within a maximum of 7.2 m, which is still a reasonable offset limit for cable considerations. The reduction in peak tensions results in a significant decrease in damage equivalent loads, on the order of 50% for upwind lines in fully loaded conditions. These results show that polymer springs can effectively reduce peak tensions and fatigue loads in mooring systems at shallow water depths.

floating offshore wind↗