Floating solar renewable energy is of enormous potential in Indonesia. This paper presents a comprehensive study of the design of Floating Photovoltaic (FPV) systems with Battery Energy Storage Systems (BESS) for three islands in Indonesia. These islands represent three typical scenarios in Indonesia (a) using a national grid powered by fossil fuel generators, (b) using a local grid powered by diesel generators, and (c) no grid at all. In-person surveys were conducted at these islands to collect data, and then FPV and BESS were designed to meet the demands of each island. Subsequently, the systems’ energy simulations were conducted using the System Advisor Model, demonstrating daily energy demand and supply in hour variation. Based on the results, a series of sustainability analyses were created from the aspects of economics, society, and the environment. The economic analysis demonstrated cost savings by using FPV to replace contemporary energy methods. The social analysis provides valuable insights into the local community, forming a demographic profile and obtaining perceptions and opinions regarding the new energy approach. The environmental analysis quantifies the potential CO₂ emissions. Overall, the work provides valuable insights into the roadmap for implementing floating solar technologies in Indonesia which can also inform global ocean-based solar energy developments.

The problem of wave interaction with multiple adjacent floating solar panels with arbitrary types and numbers of constraints is considered. All the solar panels are assumed to be homogeneous, with the same physical properties, as well as modeled by using the Kirchhoff-Love plate theory. The motion of the fluid is described by the linear velocity potential theory. The domain decomposition method is employed to obtain the solutions. In particular, the entire fluid domain is divided into two types, the one below the free surface, and the other below elastic plates. The velocity potential in the free surface domain is expressed into a series of eigenfunctions. By contrast, the boundary integral equation and the Green function are employed to construct the velocity potential of fluid beneath the entire elastic cover, with unknowns distributed along two interfaces and jumps of physical parameters of the plates. All these unknowns are solved from the system of linear equations, which is established from the matching conditions of velocity potentials and edge conditions. This approach is confirmed with much higher computational efficiency compared with the one only involving eigenfunction expansion for the fluid beneath each plate. Extensive results and discussions are provided for the reflection and transmission coefficients of water waves, maximum deflection, and principal strain of the elastic plates; especially, the influence of different types and numbers of edge constraints are investigated in detail.