Array Of Cylinders Research Articles

A Simulation of Calibration and Map-making Errors of the Tianlai Cylinder Pathfinder Array

The Tianlai cylinder array is a pathfinder for developing and testing 21 cm intensity mapping techniques. In this paper, we use numerical simulation to assess how its measurement is affected by thermal noise and the errors in calibration and map-making process, and the error in the sky map reconstructed from a drift scan survey. Here we consider only the single frequency, unpolarized case. The beam is modeled by fitting to the electromagnetic simulation of the antenna, and the variations of the complex gains of the array elements are modeled by Gaussian processes. Mock visibility data are generated and run through our data processing pipeline. We find that the accuracy of the current calibration is limited primarily by the absolute calibration, where the error comes mainly from the approximation of a single dominating point source. We then studied the m-mode map-making with the help of Moore–Penrose inverse. We find that discarding modes with singular values smaller than a threshold could generate visible artifacts in the map. The impacts of the residue variation of the complex gain and thermal noise are also investigated. The thermal noise in the map varies with latitude, being minimum at the latitude passing through the zenith of the telescope. The angular power spectrum of the reconstructed map show that the current Tianlai cylinder pathfinder, which has a shorter maximum baseline length in the North–South direction, can measure modes up to l ≲ 2π b NS/λ ∼ 200 very well, but would lose a significant fraction of higher angular modes when noise is present. These results help us to identify the main limiting factors in our current array configuration and data analysis procedure, and suggest that the performance can be improved by reconfiguration of the array feed positions.

Impacts of Biomimetic Mangrove Structures on Wave Attenuation in Kuala Sepetang, Perak

Biomimicry is introduced as an approach to develop an ecologically-informed design solution. The study proposes a regenerative design solution through hydrodynamic effects of the complex mangrove roots system that is applied to a biomimetic structural column design, with consideration on interactions between root porosity, water flows, and sediment transportation. Sangga Besar River, Perak is chosen as the key parameters of the flow. The roots are modelled with a circular array of cylinders (patch) with different porosities and varying spacing ratios. 15 patches are derived from 3 main sets where each set consisted of five patches with varying porosities. The simulation of each case to understand the impacts of different patch porosity on the flow structure is carried out using Ansys Fluent Student Version 2021. The study found that turbulence is impacted, but not directly proportional to the patch porosity where there could be an optimal porosity of the patch that generates maximum yield of turbulent kinetic energy dissipation rate with minimum sediment erosion for the onset of sediment deposition. The patch with φ=40% has the highest turbulence dissipation rate which relatively coherent with the range of porosity that would likely be encountered in natural mangrove ecosystem.

Simulation of laminar to transitional wakes past cylinders with a discontinuous Galerkin inviscid shallow water model

Laminar to transitional wakes occur in slow, quasi-steady flows past cylinders at low cylinder Reynolds numbers (R ed ≤ 250). Inviscid numerical solvers of the depth-averaged shallow water equations (SWE) introduce numerical dissipation that, depending on Red , may imitate the mechanisms of viscous turbulent models. However, the numerical dissipation rate in a second-order finite volume (FV2) SWE solver is so large at a practical resolution that this can instead hide these mechanisms. The extra numerical complexity of the second-order discontinuous Galerkin (DG2) SWE solver results in a lower dissipation rate, making it a potential alternative to the FV2 solver to reproduce cylinder wakes. This paper compares the DG2 and FV2 solvers, initially for wake formation behind one cylinder. The findings confirm that DG2 can reproduce the expected wake formations, which FV2 fails to capture, even at a 10-fold finer resolution. It is further demonstrated that DG2 is capable of reproducing key features of the flow fields observed in a laboratory random cylinder array.

Direct numerical simulations of flows around an array of rough cylinders partially coated with porous media at Re = 3900

In recent years, interest has been growing in the stability control of cylindrical structures in a flow field using full/partial coatings of porous media. However, the partial application of porous coatings to arrays of cylinders with rough surfaces has rarely been studied. This paper numerically investigates the passive flow control of an array of rough cylinders partially coated with porous media on their leeward side in the subcritical flow regime (Re = 3900). This is achieved through two-dimensional direct numerical simulations and the Cartesian cut-cell method, which facilitate flow simulations around complex porous media and rough surfaces in an accurate and flexible manner. In contrast to macroscopic models, this study employs a quasi-microscopic model to simplify the porous structure for each cylinder, providing an accurate and intuitive means of describing the inflow details of porous media near the central cylinder. The porosities and pores per inch of the porous media are defined in terms of the alignment angles and radii. Similarly, the surface roughness of a cylinder is simulated by the attachment of small circular appendages. The accuracy and capacity of the presented numerical approach are demonstrated via a numerical case study with a single smooth cylinder. The influence of the porous media on the flows around a single rough cylinder and an array of rough cylinders is thoroughly investigated and discussed.

Near trapping phenomena in arrays of porous vertical cylinders

The near trapping phenomena are associated with the wave trapping in arrays of large number of solid bodies, which experience increased amplitudes of free-surface elevation between adjacent bodies, and corresponding large wave loads on each array element. In relevant situations of infinite arrays, the wave energy is trapped within the array allowing only a small amount of energy to radiate to the far field, forming the so-called trapped mode phenomenon. The question which this paper tries to answer is whether trapped waves are simulated by arrays of porous vertical cylinders. To tackle this task, a solution method that solves the corresponding diffraction problem by the application of the multiple scattering approach, for the simulation of the hydrodynamic interactions between the waves and the elements of the array, is proposed. Numerical results for several geometrical configurations are presented and discussed, highlighting the effect of porous cylinders on the exerted hydrodynamic loads on the array.

A Bottom-up Transient Theory for Concentration Polarization with Applicability to over-Limiting Current Conditions

Renewable power sources have drawn interest due to scarcity of fossil fuels and associated pollution during their combustion-based power production. Renewables are projected to produce approximately 60% of US electricity generation by 2050, with an expected contribution of 2000 TW-h from wind power alone. To overcome the mismatch between electricity demand and intermittent supply, such power plants require large-scale energy-storage device integrated with them. In contrast to other technologies, redox flow batteries (RFBs) offer independently scalable power and energy capacity that is demanded for such applications. Despite significant achievement in their designs through experimentation and modeling, limited understanding has been established to date concerning their transient response while operating under irregular charge/discharge scenarios in renewable-powered grids and micro-grids. Under such conditions pore-scale mass transport physics and different mechanisms of concentration polarization can occur, but the conventional film law of mass transfer (FLoMT) using a constant Sherwood number cannot capture time-varying diffusion layer thickness near electrode surfaces. Further, the mean solute advection rate through porous RFB electrodes can deviate from the value expected from product of superficial velocity and average concentration inside the electrode if pore-scale concentration polarization is substantial. Consideration of these effects in the present bottom-up theory makes it applicable to over-limiting conditions, while predicting concentration polarization in agreement with experimental observations. We accomplish this by creating theory with similarity to the theory of Ralph White and co-workers for solid-state, particle-scale diffusion in Li-ion batteries [J. Power Sources, 2012, Vols. 214 and 218], while our theory is instead applied to an arbitrary representative volume element of flowing liquid electrolyte inside electrode pores.Here, we present a theory to incorporate these transient effects into a bottom-up transient model (B-UTM) by using frequency-dependent transfer functions (TFs) that are derived from the Fourier transformed pore-scale mass conservation equations (MCEs). These TFs preserve pore-scale mass transfer physics during their embedding into an up-scaled model. One of these TFs is the spectral Sherwood number that extends the film law of mass transfer to transient conditions and that correlates reactive flux with concentration polarization through an output-input relationship in the frequency domain. Another TF named the advective-flux transfer function captures the acceleration/suppression of solute advection rate that arises from the inhomogeneity in pore-scale concentration distributions. We first present validation of this theory using transient RFB experiments where the introduced B-UTM predicts the time evolution of concentration polarization in good agreement. Meanwhile the conventional FLoMT model that uses time-independent Sherwood number and neglects the solute acceleration/suppression effect shows systematically larger polarization (up to 400%), which restricts such models from simulating charging/discharging with applied current near and above the limiting current.The utility of the introduced bottom-up theory is demonstrated by incorporating it into a multi-scale RFB model in which redox-active electrolyte flows through a porous electrode comprised of an array of circular cylinders in crossflow that mimics commonly used carbon felt. The frequency dependence of the introduced TFs obtained from numerical solution of pore-scale MCEs are analyzed for different electrode porosity and Péclet number, where both the TFs reach a non-zero finite value in the limit of vanishing frequency. The non-zero value of the spectral Sherwood number in this limit is identical to the Sherwood number obtained in the pseudo-steady limit (PSL) where transient effects are negligible. With increasing frequency both the TFs show a transition from lagless response to semi-infinite Warburg response with orders higher values compared to that of the PSL. Finally, these TFs are embedded into an up-scaled model to obtain the time-domain response of RFBs operating under a step impulse applied current with different step durations. Numerically obtained concentration polarization and the average reactant concentration inside electrodes is used to construct non-dimensional regime maps showing a regime of over-limiting current. The results show that fast current fluctuations are sustained even when current exceeds the limiting current, which provides a means to increase the utilization of charge capacity of existing flow-based electrochemical devices that are operated under transient conditions. The theory introduced here can find applicability in other flow-based electrochemical devices including electrochemical separations and CO2 capture. Figure 1

On the flow conditions requiring detailed geometric modeling for multiscale evaluation of coastal forests

The multiscale evaluation method is applied to assess the influence of detailed geometric modeling of trees on their macroscopic attenuation effect against tsunami-like flow. Specifically, we conduct a series of numerical flow tests (NFTs), i.e., 3D flow simulations in a local test domain (LTD), under various inflow conditions to evaluate the macroscopic flow characteristics in the LTD accommodating an array of either simple cylinder or detailed tree models that mimic a coastal forest. After the procedure of NFTs in the multiscale evaluation method is briefly summarized and the corresponding governing equations and analysis conditions are presented, we introduce two indices for evaluating the macroscopic flow characteristics within the framework of multiscale modeling. Based on the NFT results, we discuss how the modeling scheme for trees influences the macroscopic flow characteristics in terms of these indices and clarify the microscopic mechanisms that influence the macroscopic attenuation property. Additionally, the NFT results are compared with the experimental results to justify these discussions, and key factors are explored in terms of reproducing real phenomena. In addition, to apply the multiscale evaluation method to assess the disaster mitigation performance of various types of ecosystems in a realistic situation, we discuss under what circumstances and with what level of detail vegetation should be modeled.

Numerical study of heat transfer across in-line circular cylinders for low to moderate Reynolds numbers

A computational fluid dynamics (CFD) model is generated to study the fluid flow and heat transfer from a series of six cylinders arranged in a horizontal row. The momentum and energy equations were solved using the finite volume method, with convective terms discretized using a second-order upwind method and diffusion terms discretized using central differencing. A second-order implicit technique was used for time integration. Simulations were run with four different values of the free stream Reynolds numbers 100, 200, 300, 400 and three different values of the longitudinal pitch (SL ) 1.5D, 2D, 2.5D to analyze the heat transfer behavior with fluid flow over a series of horizontal cylinders. First the grid independence test and CFD model validation are performed with the exact relation available in the literature. Next, the numerical results obtained from the developed CFD model for the row of cylinders were compared with the previously published analytical model, and correlation available in the literature, and were found to be in good agreement. The average Nusselt number calculated from the analytical and numerical models falls in the range of 8 to 10 for free stream Reynolds number equal to 100 for all values of dimensionless longitudinal pitches. Moreover, for the free stream Reynolds number equal to 400 with longitudinal pitch of 2.5D, the values of the average Nusselt number reach up to double the above-mentioned range. The maximum percentage error between the average Nusselt number obtained from the numerical and analytical solutions was less than 20% for larger free stream Reynolds number. The study found that the average Nusselt number for the array of cylinders increases with both the free stream Reynolds number and the dimensionless longitudinal pitch ratio increments. This information can be used to design more efficient heat exchangers or other fluid systems involving arrays of cylinders.

Generalized Lagrangian heterogeneous multiscale modelling of complex fluids

We introduce a fully Lagrangian heterogeneous multiscale method (LHMM) to model complex fluids with microscopic features that can extend over large spatio/temporal scales, such as polymeric solutions and multiphasic systems. The proposed approach discretizes the fluctuating Navier–Stokes equations in a particle-based setting using smoothed dissipative particle dynamics (SDPD). This multiscale method uses microscopic information derived on-the-fly to provide the stress tensor of the momentum balance in a macroscale problem, therefore bypassing the need for approximate constitutive relations for the stress. We exploit the intrinsic multiscale features of SDPD to account for thermal fluctuations as the characteristic size of the discretizing particles decreases. We validate the LHMM using different flow configurations (reverse Poiseuille flow, flow passing a cylinder array and flow around a square cavity) and fluid (Newtonian and non-Newtonian). We show the framework's flexibility to model complex fluids at the microscale using multiphase and polymeric systems. We also show that stresses are adequately captured and passed from micro to macro scales, leading to richer fluid response at the continuum. In general, the proposed methodology provides a natural link between variations at a macroscale, whereas accounting for memory effects of microscales.

Experimental study on fabrication of large-scale micro cylinders array by mask jet electrochemical deposition

Experimental study on fabrication of large-scale micro cylinders array by mask jet electrochemical deposition

Effect of the surface peak-valley features on droplet splash dynamics

In this study, we fabricated two types of functional surfaces with the same roughness (Sa) but entirely opposite surface morphological features, which are defined as the positively skewed surface filled with protruding cylinder array (Ssk > 0) and the negatively skewed surface filled with circular pit array (Ssk < 0). The effect of surface morphology peak-valley features on droplet splash is analyzed, and the formation mechanism of the prompt splash and corona splash is also indicated based on the Kelvin–Helmholtz instability. Our results demonstrate that, under the same roughness conditions of Sa, the interaction between the liquid lamellae and the thin air layer is much stronger on the negatively skewed surface, which would inhibit droplet spreading and promote the generation of droplet splash. Increasing the depth of microstructures, resulting in more pronounced peak-valley features, has been found to facilitate both prompt and corona splash phenomena to some extent. Additionally, it is found that the ease of splash formation on each surface is related to the initial spreading speed variation, with the degree of reduction in the initial spreading speed indirectly reflecting the instability of the liquid lamellae. The findings from our study contribute to the development of advanced surface engineering strategies for controlling droplet splash and enhancing the performance of various industrial applications.

Developing an Active Microfluidic Pump and Mixer Driven by AC Field-Effect-Mediated Induced-Charge Electro-Osmosis of Metal–Dielectric Janus Micropillars: Physical Perspective and Simulation Analysis

We propose herein a novel microfluidic approach for the simultaneous active pumping and mixing of analytes in a straight microchannel via the AC field-effect control of induced-charge electro-osmosis (ICEO) around metal–dielectric solid Janus cylinders of inherent inhomogeneous electrical polarizability immersed in an electrolyte solution. We coin the term “Janus AC flow field-effect transistor (Janus AC-FFET)” to describe this interesting physical phenomenon. The proposed technique utilizes a simple device geometry, in which one or a series of Janus microcylinders are arranged in parallel along the centerline of the channel’s bottom surface, embedding a pair of 3D sidewall driving electrodes. By combining symmetry breaking in both surface polarizability and the AC powering scheme, it is possible, on demand, to adjust the degree of asymmetry of the ICEO flow profile in two orthogonal directions, which includes the horizontal pump and transversal rotating motion. A comprehensive mathematical model was developed under the Debye–Hückel limit to elucidate the physical mechanism underlying the field-effect-reconfigurable diffuse-charge dynamics on both the dielectric and the metal-phase surfaces of the Janus micropillar. For innovation in applied science, an advanced microdevice design integrating an array of discrete Janus cylinders subjected to two oppositely polarized gate terminals is recommended for constructing an active microfluidic pump and mixer, even without external moving parts. Supported by a simulation analysis, our physical demonstration of Janus AC-FFET provides a brand-new approach to muti-directional electro-convective manipulation in modern microfluidic systems.

Time-dependent water wave scattering by a marine structure consisting of an array of compound porous cylinders

Under the assumption of linear wave theory, a semi-analytical model is developed to address the time-dependent water wave scattering problem involving a marine structure consisting of several circular rigid vertical cylinders, each of which is surrounded by a thin cylindrical porous wall with a water region between the inner cylinder and the wall. The problem is tackled by applying the eigenfunction expansion approach. The energy dissipation relation is also derived for the system of compound cylinders. The principal focus of this work lies in locating the optimum geometrical configurations for which the wave forces acting on the structure are minimal along a given wave direction in the frequency domain. Subsequently, the time-dependent response of the structure under the impact of a Gaussian wave pulse is examined. The wavenumber at which the force acting on the inner cylinder is at its maximum is obtained for a given value of the porous-effect parameter. Moreover, the value of that corresponding wavenumber decreases as the magnitude of the porous-effect parameter increases. The model is also carefully validated numerically against results available in the existing literature.

Large-eddy simulation of non-isothermal flow over an array of cylinders placed inside a converging channel

This work evaluates the flow and thermal characteristics of a gas-turbine blade trailing-edge. The study is carried out using large-eddy simulation (LES) for Reynolds number, Re=5900. In view of the realistic conditions, a converging channel with a staggered arrangement of fins with varying aspect ratio is considered for the higher pressure of the coolant and the larger temperature difference between the blade and the coolant. The present results are compared with the parallel channel case to identify the effects of converging flow passage of the channel on the flow and thermal characteristics. The computational setup and methodology are thoroughly validated using reference data from experiments as well as LES available in the literature. The results are analyzed in terms of Q-criterion, streamlines, vorticity and turbulent heat flux contours, and the Nusselt number profiles at several locations in the flow domain. The present results shows, unlike in the parallel channel case, a dominant spanwise component of velocity that helps in the formation of a strong vortex behind the first row itself, and the Nusselt number consistently increases from the first to the last row. The maximum gain and diffusion of turbulent heat flux occur within a recirculation zone and along the free-shear layer, respectively. The heat flux from the end-wall is extracted mainly by the spanwise velocity component, which is also responsible for its advection to the central region. Heat transfer is observed to be higher from fins compared to the end-walls, which underscores the importance of the active participation of fins in the analysis of heat transfer at the trailing-edge.

Influence of submergence ratio on flow and drag forces generated by a long rectangular array of rigid cylinders at the sidewall of an open channel

This paper discusses how the submergence ratio, defined as the ratio between the flow depth, D, and the height, h, of the solid rigid cylinders forming the array affects flow and turbulence structure inside and around a rectangular array of cylinders placed adjacent to one of the channel sidewalls. As the array becomes submerged, a vertical shear layer develops in between the top face of the array and the free surface, which strongly increases flow three-dimensionality and modifies how the momentum exchange between the array and the surrounding open water regions occurs with respect to the case of an emerged array where only a horizontal shear layer forms as part of the incoming flow approaching the array is deflected laterally. Eddy-resolving simulations are conducted for several values of the solid volume fraction, ϕ, and of the submergence ratio, 1.0 ≤ D/h ≤ 4.0. Similar to the limiting case of an emerged array (D/h = 1.0), the width- and depth-averaged streamwise velocity inside the array reaches a constant value after an initial adjustment region in the submerged-array cases. For D/h ≥ 1.33, the mean normal velocities through the top and side faces of the array do not become equal to zero downstream of the initial adjustment region. The flow inside the array reaches an equilibrium regime where the local flux of fluid leaving the array through its side face is balanced by the local flux of fluid entering the array through its top face. This regime is observed until close to the end of the array. For D/h ≥ 2.0, the horizontal shear layer vortices do not generate successive regions of high and low streamwise velocity and bed friction velocity inside the array, as is observed in the emerged cases. With respect to the emerged case, the size of the shear layer vortices and the shear layer width increase for low submergence ratios before decreasing rapidly for D/h ≥ 1.33. Significant three-dimensional effects are present inside the array and the horizontal shear layer for cases with both emerged and submerged arrays. In the ϕ = 0.08 cases, strong upwelling and downwelling motions are observed inside the array for D/h ≥ 1.33, while the circulation of the streamwise-oriented cell of secondary flow forming close to the lateral face of the array peaks when the submergence ratio is close to 1.33. For constant ϕ, the total streamwise drag force normalized with the height and width of the array increases with increasing submergence ratio. As the array submergence increases, the cylinders near the front of the array contribute less to the total force acting on the cylinders forming the array. For constant D/h, the total streamwise drag force acting on the array increases with ϕ for the submerged and emerged cases.

Flow-induced forces and vortex transportation characteristics of three circular cylinders at subcritical Reynolds number

Three-dimensional numerical simulation for the three circular cylinders at Re = 3900 is carried out by using the large eddy simulation (LES) in this paper. The present numerical model and parameters are verified by solving single cylinder and two cylinders cases. Then, the parameters consisted of spacing ratio (L/D) and incidence angle (α) are checked in three circular cylinders array. The fluid force and near-wake behavior of the two- and three-cylinder systems are analyzed with emphasis. Under a tight arrangement, there are three flow patterns: reattachment flow pattern, large-incidence-angle flow pattern, and biased flow pattern, corresponding to cases with incidence angles of 0, 45, and 90 degrees, respectively. The Strouhal number (St) of the three cylinders is essentially overlapping at α= 0° and α= 90°, the distinctive secondary peaks can be observed at the small spacing ratio in power spectral density (PSD). The mean pressure coefficient of midstream cylinder (MC) influenced by the number of cylinders is negative in the tandem arrangement at the stagnation point, accompanied by reattachment phenomenon the upstream surface of MC with L/D= 1.5. Placing the third cylinder between the two cylinders in a tandem arrangement significantly suppresses their lift oscillations. The phase difference (φum) between the fluctuating lift of the upstream two cylinders is highly sensitive to the incidence angle of the system.

Flow structure and dynamics behind cylinder arrays at Reynolds number ∼100

The flow behind nine different arrays of cylinders is experimentally investigated via Particle Image Velocimetry (PIV) at a Reynolds number of Re ∼100 based on the diameter of the cylinders. Each array consists of a column of four cylinders in front and three in the rear. The horizontal distance between the two columns and the vertical distance between the cylinders within each column are varied for H/D=[2,4,8] and V/D=[2,4,6], resulting in nine different arrays denoted as mVnH, where m corresponds to V/D and n stands for H/D. The PIV measurements are conducted for 15 s at 200 Hz frequency, corresponding to 39 to 360 vortex shedding events for the wakes in this study. Then, proper orthogonal decomposition is applied to the velocity fields to analyze the flow dynamics. All arrays show unsteady flow, and based on their flow structures, they are classified in to three main categories of single bluff body (SBB), transitional (TR), and co-shedding (CS) flow. SBB characteristics can be seen for 2V2H and 2V4H arrays, but the latter has more steady vortex shedding as the H/D increases from 2 to 4. Then, 2V8H and 4V2H have an asymmetric flow with several vortex streets and act as an intermediary stage in the shift from SBB to CS flow structure when the distances are increased. The highest total kinetic energy values and widest probability density functions of the velocity components are observed for this group. The five remaining arrays in the CS group have symmetric flow, with three or five vortex streets present behind. However, based on the distances, the frequency and phase synchronization of the vortex streets change considerably, which might have an important effect on, for example, the heat transfer or the structural load of the cylinders.

TiN-based metasurface absorber for efficient solar energy harvesting

Solar absorber, which is widely used in concentrated solar thermal systems and solar photovoltaics, has attracted great attention in recent years. Concentrated solar irradiation not only saves production costs, but also effectively improves the photothermal conversion efficiency of the system. However, the huge radiative heat flux will cause the system temperature to rise rapidly, which poses significant challenges for designing solar absorbers. Here, a metasurface absorber based on cylinder array structure is proposed for efficient solar energy harvesting, using the high-temperature resistant materials TiN and SiO2. The total solar absorption of the absorber is up to 0.94 at wavelengths of 300–2500 nm. The high absorption performance of the metasurface absorber can be explained as the coupling effect of surface plasmon resonance and Fabry-Pérot resonance, which is confirmed by the electric field distribution. Moreover, the effect of geometric parameters on absorption performance is analyzed. Finally, we discuss the influence of incident angle on the solar absorber. We believe this work could deepen the understanding of coupling resonance mode and guide the design of high-temperature solar absorbers.

Waves in frequency-independent single-mode 2-D artificial media with forced refraction

Investigations of the waves in planar structures with anisotropic gratings proved that structures of this type demonstrated properties of 2D media with forced refraction. Parasitic modes and limited frequency range were typical for them. This study investigated the propagation of guided waves along metal grating lying between two metal screens inside parallel plate waveguide (PPW). Dominant T-mode was the focus of this study. This mode was studied theoretically by the method of integral equations (IE) and numerically with the help of an electromagnetic Ansys HFSS simulation. Both methods demonstrated that the investigated structure had the properties of a 2D medium with forced refraction independent of frequency. The modified structure contained metal grating and an additional 2D array of metal cylinders was also studied numerically. The numerical simulation demonstrated that an array of metal cylinders suppressed parasitic PPW mode and thus modified structure demonstrated properties of a single-mode frequency independent artificial media with forced refraction. The results have the potential for application in microwave and millimeter-wave planar quasi-optical lens antennas.

Achieving Complete Band Gaps of Water Waves with Shallow-Draft Cylinder Arrays

Conventionally, band gaps of water waves can be obtained with periodic arrays of bottom-mounted obstacles or vertical cylinders. However, such structures need to possess a large ratio of volume in water, and may not be easy to build and move in practical ocean engineering. Here, we theoretically propose and experimentally demonstrate that when the water surface is pierced by a fixed, rigid cylinder array, a complete band gap of water waves can be achieved, in which the propagation of water waves is forbidden in all the directions. Such a band gap can exist even when the cylinders have a draft much smaller than the water depth ($d\ensuremath{\ll}h$). Based on the complete band gap, a semicircular shape of the fixed-surface disk array is further constructed, which can block omnidirectional water waves generated by a point source.