Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear programming”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Verification of GENE and GYRO with L-mode and I-mode plasmas in Alcator C-Mod

Here, verification comparisons are carried out for L-mode and I-mode plasma conditions in Alcator C-Mod. We compare linear and nonlinear ion-scale calculations by the gyrokinetic codes GENE and GYRO to each other and to the experimental power balance analysis. The two gyrokinetic codes' linear growth rates and real frequencies are in good agreement throughout all the ion temperature gradient mode branches and most of the trapped electron mode branches of the kyρs spectra at r/a = 0.65, 0.7, and 0.8. The shapes of the toroidal mode spectra of heat fluxes in nonlinear simulations are very similar for k y ρ s ≤ 0.5, but in most cases GENE has a relatively higher heat flux than GYRO at higher mode numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Identifying key factors of the seawater intrusion model of Dagu river basin, Jiaozhou Bay

Highlights: • The key factors of seawater intrusion model is identified by mutual entropy analysis. • Precipitation rate and groundwater pumping are the key factors for seawater intrusion. • The hydraulic conductivity of zone 1 has significant influence on seawater intrusion. Seawater intrusion is a complex groundwater - seawater interaction process, and it is influenced by many factors from ground surface to underground, from groundwater to seawater. Generally, for seawater intrusion model, some model parameters and boundary conditions are always specified by model users’ personal experiences or literature's reference value. The defective model would damage the groundwater management for controlling and preventing seawater intrusion when making decisions are based on this model. In order to improve the reliability of seawater intrusion model, the influences of model inputs on output should be identified prior at optimizing model inputs. Dagu river basin, Jiaozhou Bay is one of the most serious areas of seawater intrusion in China, and it is chosen as the study area in this study. The seawater intrusion model of Dagu river basin is built based on a general program SEAWAT4. The key influence factors of model output are analyzed by two sensitivity analysis methods, i.e., stepwise regression and mutual entropy. The results demonstrated that the most important influence factors which have largest sensitivities to groundwater Cl{sup -} concentration are the precipitation rate and groundwater pumping in agriculture area. In addition, the hydraulic conductivity of zone 1 has a non-negligible influence on seawater intrusion process. Stepwise regression analysis is capable of identifying most important influence factor, and it can’t handle complicated nonlinear input-output relationship. Mutual entropy analysis is reliable for identifying the influence factors for complex seawater intrusion model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Optimization of Superconducting Linac for Proton Improvement Plan-II (PIP-II)

PIP-II is an essential upgrade of the Fermilab complex that will enable the worlds most intense high-energy beam of neutrinos for the international Deep Underground Neutrino Experiment at LBNF and support a broad physics program at Fermilab. Ultimately, the PIP-II superconducting linac will be capable of accelerating the $H-$ CW beam to 800 MeV with an average power of 1.6 MW. To operate the linac with such high power, beam losses and beam emittance growth must be tightly controlled. In this paper, we present the results of global optimization of the Linac options towards a robust and efficient physics design for the superconducting section of the PIP-II linac. We also investigate the impact of the nonlinear field of the dipole correctors on the beam quality and derive the requirement on the field quality using statistical analysis. Finally, we assess the need to correct the quadrupole focusing produced by Half Wave, and Single Spoke accelerating cavities. We assess the feasibility of controlling the beam coupling in the machine by changing the polarity of the field of linac focusing solenoids

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

The 2020 Nonlinear Mechanics and Dynamics Research Institute

The 2020 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 15 to July 30, 2020. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers – especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 11 students participated in the seven-week long program held virtually due to the COVID-19 health pandemic. The students collaborated on one of four research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and other academic and research institutions. In addition to the research activities, the students attended weekly technical seminars, various virtual tours, and socialized at virtual gatherings. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD 2020 are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2022

Simultaneous Observations of Surface Layer Profiles of Humidity, Temperature, and Wind Using Scanning Lidar Instruments

We demonstrate the combination of three Doppler lidars (DLs), a water vapor differential absorption lidar and a temperature rotational Raman lidar for the investigation of the interactions between the land-surface and the atmospheric boundary layer. This combination of scanning lidars was operated for the first time during the Land-Atmosphere Feedback Experiment at the Atmospheric Radiation Measurement program's Southern Great Plain site, Oklahoma, USA, in August 2017, and provided simultaneous surface layer profiles of horizontal wind, humidity and temperature. The horizontal wind profiles were determined using the dual-Doppler method with two DLs. The scans were performed above four towers providing atmospheric variables and fluxes. These combined lidar data allowed for the estimation of the friction velocity as well as the surface latent and sensible heat fluxes. For this purpose, profiles calculated with Monin-Obukhov similarity theory were fitted with the Levenberg-Marquardt nonlinear least squares curve-fitting method to the measured surface-layer lidar profiles. We present case studies of three 50-min periods to illustrate the new method. The measurements resulted in friction velocities of 0.31, 0.29, and 0.38 ms –1 , sensible surface heat fluxes of 312.1, 234.1, and 183.3 Wm –2 and latent surface heat fluxes of 251.8, 227.3, and 274.4 Wm –2 . These values were compared with the in-situ measurements of the towers. Considering the sampling differences and the error analyses, the results agree within 25.1%, 42.5%, and 28.8%, respectively, demonstrating that this synergy of scanning active remote sensing systems can be used to derive surface fluxes with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Microwave Imaging Reflectometry for the study of Edge Harmonic Oscillations on DIII-D [Microwave Imaging Reflectometry (MIR) for the study of Edge Harmonic Oscillations (EHOs) on DIII-D]

Quiescent H-mode (QH) is an ELM free mode of operation in which edge-localized harmonic oscillations (EHOs) are believed to enhance particle transport, thereby stabilizing ELMs and preventing damage to the divertor and plasma facing components. Microwave Imaging Reflectometer (MIR) enabling direct comparison between the measured and simulated 2D images of density fluctuations near the edge can determine the 2D structure of density oscillation which can help to explain the physics behind EHO modes. MIR data sometimes indicates a counter-propagation between higher (n>1) and dominant (n=1) harmonics of coherent EHOs in the steep gradient regions of the pedestal. To preclude diagnostic artifacts, we have performed forward modeling that includes possible optical misalignments to show that offsets between transmitting and receiving antennas do not account for this feature. We have also simulated the non-uniform rotation of the EHO structure, which induces multiple harmonics that are properly characterized in the synthetic diagnostic. Excluding these possible explanations for the data, the counter-propagation observed in MIR data, which is not corroborated by external Mirnov coil array measurements, may be due to subtleties of the eigenmode structure, such as an inversion radius consistent with a magnetic island. Similar effects are observed in analysis of internal ECE-Imaging and BES data. Furthermore, the identification of a non-ideal structure motivates further exploration of nonlinear models of this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Application of Variance-Based Sensitivity Analysis to a Large System Dynamics Model

Variance-based sensitivity methods can provide insights into large computational models. We present a novel application of sensitivity analysis to the Biomass Scenario Model (BSM) a large and complex system dynamics model of the developing biofuels industry in the United States. We apply a two-stage sensitivity approach consisting of an initial sensitivity screening, followed by a variance decomposition approach. Identifying key system levers and quantifying their strength is not straightforward in complex system dynamics models that have numerous feedbacks and nonlinear results. Variance-based sensitivity analysis (VBSA) offers a systematic, global approach to assessing system dynamics models because it addresses nonlinear responses and interactive effects. Especially when a large model's size makes manual exploration of the input space difficult and time-consuming, the approach can help to provide a comprehensive understanding of interactions that drive model behaviors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The 2019 Nonlinear Mechanics and Dynamics Research Institute

The 2019 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 17 to August 1, 2019. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 20 students came to Albuquerque, New Mexico to participate in the seven-week long program held at the Mechanical Engineering building on the University of New Mexico campus. The students collaborated on one of seven research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and academic institutions. In addition to the research activities, the students attended weekly technical seminars, various tours, and socialized at various off-hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗
DOE OSTI2019

The 2018 Nonlinear Mechanics and Dynamics Research Institute

The 2018 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 18 to August 2, 2018. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers -- especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 17 students came to Albuquerque, New Mexico to participate in the seven-week long program held at the Mechanical Engineering building on the University of New Mexico campus. The students collaborated on one of six research projects that were developed by various mentors from Sandia National Laboratories, University of New Mexico, and academic institutions. In addition to the research activities, the students attended weekly technical seminars, various tours, and socialized at various off-hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗
DOE OSTI2015

RELAP-7 Development Updates

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory, and is the next generation tool in the RELAP reactor safety/systems analysis application series. RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). During Fiscal Year (FY) 2015, the RELAP-7 code has been further improved with expanded capability to support boiling water reactor (BWR) and pressurized water reactor NPPs analysis. The accumulator model has been developed. The code has also been coupled with other MOOSE-based applications such as neutronics code RattleSnake and fuel performance code BISON to perform multiphysics analysis. A major design requirement for the implicit algorithm in RELAP-7 is that it is capable of second-order discretization accuracy in both space and time, which eliminates the traditional first-order approximation errors. The second-order temporal is achieved by a second-order backward temporal difference, and the one-dimensional second-order accurate spatial discretization is achieved with the Galerkin approximation of Lagrange finite elements. During FY-2015, we have done numerical verification work to verify that the RELAP-7 code indeed achieves 2nd-order accuracy in both time and space for single phase models at the system level.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

RELAP-7 Pressurizer Component Development Updates

RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). As part of the efforts to expand the capability for PWR simulation, an equilibrium single-region pressurizer model has been implemented in RELAP-7. The pressurizer component can simulate pressure and water level change through insurge, spray, and heating processes. Two simple tests – one for insurge process and another for outsurge process – have been reported to demonstrate and verify the functions of the pressurizer model. The typical single-phase PWR system model presented in the first RELAP-7 milestone report has been updated, as part of system level test for the new pressurizer model. The updated PWR system model with the pressurizer component can be used for more realistic transient simulations. The addition of the equilibrium single-region pressurizer model represents the first step of developing a suite of pressurizer models with different levels of fidelities. Non-equilibrium two-region two-zone model and non-equilibrium five-region three zone model will be implemented in the near future. Cross comparison and validation will be pursued when higher fidelity pressurizer models are completed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Design Considerations for Proposed Fermilab Integrable RCS

Integrable optics is an innovation in particle accelerator design that provides strong nonlinear focusing while avoiding parametric resonances. One promising application of integrable optics is to overcome the traditional limits on accelerator intensity imposed by betatron tune-spread and collective instabilities. The efficacy of high-intensity integrable accelerators will be undergo comprehensive testing over the next several years at the Fermilab Integrable Optics Test Accelerator (IOTA) and the University of Maryland Electron Ring (UMER). We propose an integrable Rapid-Cycling Synchrotron (iRCS) as a replacement for the Fermilab Booster to achieve multi-MW beam power for the Fermilab high-energy neutrino program. We provide a overview of the machine parameters and discuss an approach to lattice optimization. Integrable optics requires arcs with integer-pi phase advance followed by drifts with matched beta functions. We provide an example integrable lattice with features of a modern RCS - long dispersion-free drifts, low momentum compaction, superperiodicity, chromaticity correction, separate-function magnets, and bounded beta functions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

NID Copper Sample Analysis

The current focal point of the nuclear physics program at PNNL is the MAJORANA DEMONSTRATOR, and the follow-on Tonne-Scale experiment, a large array of ultra-low background high-purity germanium detectors, enriched in 76Ge, designed to search for zero-neutrino double-beta decay (0νββ). This experiment requires the use of germanium isotopically enriched in 76Ge. The MAJORANA DEMONSTRATOR is a DOE and NSF funded project with a major science impact. The DEMONSTRATOR will utilize 76Ge from Russia, but for the Tonne-Scale experiment it is hoped that an alternate technology, possibly one under development at Nonlinear Ion Dynamics (NID), will be a viable, US-based, lower-cost source of separated material. Samples of separated material from NID require analysis to determine the isotopic distribution and impurities. DOE is funding NID through an SBIR grant for development of their separation technology for application to the Tonne-Scale experiment. The Environmental Molecular Sciences facility (EMSL), a DOE user facility at PNNL, has the required mass spectroscopy instruments for making isotopic measurements that are essential to the quality assurance for the MAJORANA DEMONSTRATOR and for the development of the future separation technology required for the Tonne-Scale experiment. A sample of isotopically separated copper was provided by NID to PNNL in January 2011 for isotopic analysis as a test of the NID technology. The results of that analysis are reported here. A second sample of isotopically separated copper was provided by NID to PNNL in August 2011 for isotopic analysis as a test of the NID technology. The results of that analysis are also reported here.

copper isotopics↗
DOE OSTI2018

The 2017 Nonlinear Mechanics and Dynamics Research Institute.

The 2017 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 19 to July 28, 2017. NOMAD seeks to bring together participants with diverse technical backgrounds to work in small teams to utilize an interactive approach to cultivate new ideas and approaches in engineering . NOMAD provides an opportunity for researchers - especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 17 students from around the world came to Albuquerque, New Mexico to participate in the six - week long program held at the University of New Mexico campus. The students collaborated on one of six research projects that were developed by various mentors from Sandia National Laboratories, academia, and other government laboratories. In addition to the research activities, the students attended weekly technical seminars, toured the National Museum of Nuclear Science & History, and socialized at various off - hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation o n their research findings that was broadcast via Skype. Many of the research discoveries made at NOMAD are published as proceedings at technical conference s and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗