Incident Wave Angle Research Articles

Advanced terahertz refractive index sensor in all-dielectric metasurface utilizing second order toroidal quasi-BIC modes for biochemical environments

We present a novel design for an all-dielectric metasurface ultra-sensitive refractive index (RI) sensor, characterized by a high-quality factor (Qf) and figure of merit (FOM). This design incorporates two high-order toroidal modes, including electric toroidal quadrupole (\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{e}$$\\end{document}​) and magnetic toroidal quadrupole (\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{m}$$\\end{document}​) based on quasi‑bound states in the continuum (q-BIC) resonances. Our findings demonstrate the feasibility of switching between \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{e}$$\\end{document} and \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{m}$$\\end{document}​ through various symmetry-breaking mechanisms. To excite these high-order toroidal modes, we explored several symmetry-breaking strategies, including complex structural designs, symmetry disruption, and variations in the incident wave angle. Symmetry breaking in the structure induces new modes in the transmission spectrum, which are highly advantageous for sensing applications due to the presence of dark modes. The designed metasurface exhibits the capability to sense a broad range of RI in diverse environments, particularly in biochemical fields. Sensitivity (S) is significantly enhanced by the excitation of new resonance modes and adjustments in the incidence angle, increasing from 217 GHz/RIU in a symmetrical structure to 512.3 GHz/RIU for \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{e}$$\\end{document}​. The FOM improves from 197 RIU 1 to 8538 RIU− 1 for \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{e}$$\\end{document}​ and 152,395 RIU− 1 for \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{m}$$\\end{document}​. Additionally, the Qf rises from 872 to 17,983 and 921,351 for \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{e}$$\\end{document}​ and \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\:{Q}_{T}^{m}$$\\end{document}​, respectively.

Second harmonic generation based on graphene hyperstructure for higher resolution and performance on top of achieving fundamental wave detection in theory

In this paper, an innovative one-dimensional graphene hyperstructure (GHS) is proposed, allowing for the concurrent detection of multiple physical parameters in both the fundamental and second harmonic generation. The sensing characteristics of GHS pertaining to magnetic field strength (B), incident electromagnetic wave angle (θ), and graphene thickness (dgt) are systematically investigated. Moreover, through the incorporation of second harmonic generation alongside fundamental detection, higher resolution and performance are achieved. The findings indicate an expansion of the measurement range for B, θ, and dgt, from 0.3∼0.5 T, 35∼55°, and 1∼6 layers to 0.3∼1 T, 35∼65°, and 1∼10 layers, providing increased flexibility and adjustability. Additionally, by leveraging nonlinear effects and widening the Fabry-Perot cavity width, this structure effectively enhances the quality factor (Q) from 2.94 × 102 to 1.95 × 105, resulting in a substantial improvement in sensing performance. This development holds tremendous promise in surpassing the diffraction limit and addressing high-Q value sensing requirements. In comparison to conventional detectors, the GHS not only enhances detection efficiency but also harbors the potential for multiple physical quantities detection. This forward-looking research is pivotal in its successful resolution of detector performance limitations, ushering in novel possibilities across diverse domains.

Characterization of sound absorption by a foliage hedge at oblique incidence

Global warming requires public authorities to introduce legislation to protect urban air quality. Among existing environmental solutions, green walls meet this objective by reducing micro-particles, while improving biodiversity, sustainable development, thermal aspects and noise reduction. Numerous studies have demonstrated the noise-reducing benefits of vegetation. Despite the many published works on street modelling, most of them do not consider the effect of the incidence angle of rays on green walls. However, this can have very significant effects on wave absorption. The impact of the incidence angle of a sound wave passing through a foliage media on acoustic absorption is proposed in this paper by considering the Horoshenkov-Miki model. The effective parameters (density and compressibility) of a foliage sample for several leaves arrangements are first compared with those obtained experimentally by the transfer matrix method. These parameters are then used to map the absorption coefficient of this media as a function of the incidence angle.

Propagation of Free Infragravity Waves Generated at Distant Beaches

AbstractInfragravity (IG) waves are important drivers of extreme run‐up, morphological changes, and seiches. While locally forced IG waves have been extensively investigated, recent studies have revealed the significance of free IG waves generated on distant beaches. This study focuses on free IG waves generated at a coast facing southeast in Japan during the passage of a typhoon. The relationship between incident short waves and free IG waves as well as their contributions to nearshore IG waves, in particular to seiches in a small port, are analyzed using unique measurement data and numerical experiments. During the typhoon passage, more than 75% of the observed IG wave energy originates from free IG waves at an observatory located at a depth of 23 m, and their peak direction is alongshore. Six‐year measurement data demonstrate that peak directions of free IG waves strongly depend on the incident wave angles of short waves and that swells from the south generate alongshore propagating free IG waves. A numerical model can reproduce the alongshore propagating free IG waves accurately when using a large computational domain. Moreover, numerical experiments performed using the model demonstrate that the alongshore propagating free IG waves are IG waves reflected from distant beaches. The free IG waves from distant beaches are small outside the port, but seiches in the port are amplified by more than 10%. The relationship between seiches and incoming free IG waves is further discussed based on the numerical experiments.

Experimental Study of the Random Wave-Induced Hydrodynamics and Soil Response in a Porous Seabed Around Double Piles

The evaluation of the wave-induced pore pressures around the offshore piles has attracted great attentions among coastal engineers, because they have been commonly used as foundations of numerous marine infrastructures. This paper presents comparative studies of the random wave-induced transient seabed response around single and double piles in a sandy seabed through a series of wave flume experiments. The influences of relative spacing ratios, wave incidence angles, and front pile diameters under different random wave parameters on oscillatory pore pressures in the vicinity of double piles are examined. In addition, variations in wave profiles and dynamic wave pressures surrounding single and double piles are quantitatively analyzed. Based on the experimental results, the following conclusions can be drawn: (1) under the influence of random waves, the wave profiles around the double piles exhibit obvious irregularity and nonlinearity; (2) the shielding effect existing in the tandem piles results in lower dynamic wave pressures around the rear pile compared to the front pile; (3) the pore pressures on the front surface of the double piles decrease with increasing soil depth, with a decreasing attenuation rate at each layer; (4) when the relative spacing ratio G/D2=3, the group-pile effect weakens, leading to an increase in the pore pressures around the rear pile, approaching the results of a single pile under conditions of lower significant wave heights or periods; (5) the intense disturbance effect caused by large wave incidence angles exacerbates the pore pressure response around the double piles; (6) when the diameter of the front pile in the tandem piles increases, it enhances the shielding effect, thus suppressing the seabed response around the rear pile. In contrast, it causes an increase in the wave surface around the double piles, exacerbating the pore pressure response in the seabed. The latter effect becomes more pronounced when the significant wave height is larger.

A frequency domain analysis method for the dynamic response of a submerged floating tunnel under wave action

To investigate the dynamic response of a submerged floating tunnel (SFT) under the action of waves, a modal superposition method in the frequency domain was established, and the amplitudes of the modals were determined using the generalized wave loads and hydrodynamic coefficients of the structural modals. In this model, the SFT was simplified as a constant-cross-section Euler–Bernoulli beam with multiple anchor cables, and its structural modals were resolved using the finite element method (FEM). The generalized wave loads and generalized hydrodynamic coefficients were obtained using the Fourier expansion strategy of the modal function, by which the hydrodynamic problem can be converted into 2D problems with the modified Helmholtz controlling equation and solved using the matched higher-order boundary element method (HOBEM) and eigenfunction expansion. The effects of the diameter of the anchor cables and length of the SFT on the stiffness of the structure and the dynamic response of an SFT with anchor cables under wave action were investigated. Furthermore, the effects of different added mass calculation methods on the structural response were studied. In addition, the wave exciting forces were calculated using potential flow theory or linearized Morison equations according to the wave scale. Finally, the effects of the incident wave frequencies and incident angles were investigated in detail to reveal the combined resonance mechanism.

GWECS: An indicator to describe the energy absorption characteristics of wave energy converter arrays in real sea states

Grouping wave energy converters (WECs) into an array has become a trend to improve power production and achieve economies of scale. Therefore, the performance of a WEC array in real sea states is of vital importance. To intuitively reflect an array's frequency-based energy absorption characteristics instead of only the total extractable power in given irregular wave conditions, a new indicator of group wave energy capture spectrum (GWECS) is proposed by extending the concept of wave energy capture spectrum (WECS) to an array. A practical method to measure the indicator by estimating the absorbed wave surface from the array's power signal is also given. Simulations in two-dimensional conditions with a 2-cuboid WEC array and in three-dimensional conditions with a 2-cylinder WEC array are conducted to study the concept. The results examine the consistency in capture efficiency between the newly proposed indicator and conventional ones, and present the effect of various parameters, including linear power take-off (PTO) damping, buoy permutation, separation distance, and incident wave angle on the GWECS.

Enhancement of broadband terahertz trace absorption spectrum based on angular multiplexing of trapezoidal dielectric grating

Terahertz spectroscopy provide a fast, label-free, and non-destructive method for detecting bio-medical molecules, with wide applications in national security, pharmacy, and healthcare. However, detecting terahertz absorption spectra of trace analytes still faces significant challenges. In recent years, metamaterial parameter multiplexing technology has been used to detect the terahertz fingerprint spectra of trace substances. This paper proposes a new structure based on the angle-multiplexing enhancement principle of a trapezoidal dielectric grating to enhance the terahertz absorption spectra of thin-film analytes. By scanning the incident wave angle, a series of resonance absorption peaks can be obtained. The enhanced absorption spectrum formed by the envelope of absorption spectrum peaks increases by 141 times compared to directly measuring the absorption spectrum of the same film. This structure provides a new and effective means for terahertz trace detection using metamaterial parameter multiplexing technology.

All-Dielectric Metasurface-Based Terahertz Molecular Fingerprint Sensor for Trace Cinnamoylglycine Detection.

Terahertz (THZ) spectroscopy has emerged as a superior label-free sensing technology in the detection, identification, and quantification of biomolecules in various biological samples. However, the limitations in identification and discrimination sensitivity of current methods impede the wider adoption of this technology. In this article, a meticulously designed metasurface is proposed for molecular fingerprint enhancement, consisting of a periodic array of lithium tantalate triangular prism tetramers arranged in a square quartz lattice. The physical mechanism is explained by the finite-difference time-domain (FDTD) method. The metasurface achieves a high quality factor (Q-factor) of 231 and demonstrates excellent THz sensing capabilities with a figure of merit (FoM) of 609. By varying the incident angle of the THz wave, the molecular fingerprint signal is strengthened, enabling the highly sensitive detection of trace amounts of analyte. Consequently, cinnamoylglycine can be detected with a sensitivity limit as low as 1.23 μg·cm-2. This study offers critical insights into the advanced application of THz waves in biomedicine, particularly for the detection of urinary biomarkers in various diseases, including gestational diabetes mellitus (GDM).

Modeling obliquely incident seismic wave propagation in rocks via the FDM-DEM coupling scheme: Algorithms and applications

The incidence angle of seismic waves plays a pivotal role in the stability analysis of structures in near-field zones. Despite its importance, the application of obliquely incident seismic waves within models of discontinuous media has been limited. This study introduces a novel approach that integrates the oblique incidence of seismic waves into a coupled model employing both the Finite Difference Method (FDM) and the Discrete Element Method (DEM). We validate the accuracy of seismic wave propagation under oblique incidence in two FDM-DEM coupling models—overlapping coupling and interface coupling—through comparisons with theoretical solutions. For vertically incident seismic P-waves, both models yield similar results. However, under oblique incidence, the overlapping coupling model demonstrates superior accuracy, with a maximum error of approximately 4.217 %, compared to about 10 % in the interface coupling model. Our findings also show that while the relative sizes of particle diameter and overlap length minimally affect accuracy, the arrangement of particles markedly influences the results. Furthermore, we apply this method to a slope model to investigate the failure process, revealing that the nature of the sliding body shifts from tensile sliding to tensile ejection as the incidence angle increases.

Flow characteristics of the rip current system adjacent to a coastal vertical structure for irregular waves

Rip currents are common hazards found on many beaches around the world. The present research examines the rip current flow features adjacent to a coastal vertical structure with total reflection function to the incident waves. The laboratory experiments were conducted over the planar beach for both regular and irregular waves to explore the specific phenomenon for the irregular waves, and a higher-order Boussinesq equation model was adopted in the simulations to analyze the formation mechanism of the flow features. The study results show that, the standing wave pattern for the irregular wave is much weaker than that for the regular wave with inducing only one apparent node rip at the first node line in the reflection wave area. Moreover, the specific “single-vortex” flow pattern was formed in the crossing wave field for the irregular wave tests. This specific flow pattern is found essentially induced by the larger wave reflection rate of the vertical structure, which leads to the strong lateral flow against the downstream longshore current. This specific flow pattern prefers to present for the smaller wave incident angles and larger wave heights and periods, and leads to the feeding flow volume rates discharged more in the form of lateral flow. The C-shaped lateral flow is found mainly induced by the longshore gradient of the wave set-up instead of the non-uniform wave radiation stress nearshore.

Numerical investigation of rip currents near a barred beach induced by irregular waves

Rip currents pose a common natural hazard at coastal tourist beaches, presenting a significant threat to the safety of beachgoers and beach management. To better understand the characteristics of rip currents, particularly those induced by irregular waves on a barred beach, this study utilized a Boussinesq-type, phase-resolving hydrodynamic model for numerical simulations. We investigated various factors, including significant wave height, peak period, incident angle, and bottom friction, to assess their impact on rip currents formation. The results of our research reveal some key insights. Firstly, an increase in significant wave height and peak period fosters the development of rip currents. However, beyond a certain threshold, the intensity of the rip currents decreases. Conversely, a larger incident angle tends to reduce the intensity of rip currents. Additionally, higher bottom friction leads to increased energy dissipation during wave propagation, consequently resulting in decreased rip currents intensity. In an effort to create a more precise representation of real sea conditions, multi-directional waves were employed as incident waves in this study. An analysis of different standard deviations of directional distributions revealed that as the directional spread of incident waves widens, the waves become more susceptible to breaking, ultimately leading to a reduction in the intensity of rip currents.

Microwave Hyperthermia Technology Based on Near‐Field Focused Metasurfaces: Design and Implementation

AbstractThis paper presents a near‐field focused metasurface design method for microwave hyperthermia. Utilizing full‐wave electromagnetic simulation, a novel dual‐layer metasurface unit with the following characteristics is developed: transmission efficiency higher than −1 dB, adjustable transmission phase from 0° to 360°, and insensitivity to the polarization and incident angle of the incoming wave. Leveraging phase compensation theory, the units at different positions in the array to achieve near‐field focusing functionality is meticulously designed. To further optimize transmission efficiency, the transmission phase of the central unit of the metasurface to 76° is tuned. When applying this near‐field focused metasurface to microwave hyperthermia, the electromagnetic and thermal effects in biological tissues is conducted simulations and experiments to analyzed. The results demonstrate that the proposed metasurface can effectively focus microwave energy, achieving efficient microwave transmission in pork tissue and generating localized high temperatures in the target area. This research introduces new design concepts and implementation strategies for the development of microwave hyperthermia technology.

Study on the interaction mechanism of laser-generated Rayleigh waves and subsurface inclined cracks

Abstract In this paper, the finite element method was used to study the reflected and transmitted waves of laser-generated Rayleigh waves from subsurface inclined cracks, the propagation paths and mode conversion mechanisms of different characteristic waves are determined. The Rayleigh wave will interact with the crack top tip and propagate back and forth along the crack surface, and be converted to shear waves at the crack top tip. The shear waves will mode-convert to Rayleigh waves at the free surface when the incidence angle of the shear wave is larger than 60°. Moreover, for the Rayleigh wave interacting with the crack bottom tip, when the crack inclined angle is less than 60°, some Rayleigh waves will travel along the crack surface to the crack top tip. When the crack inclination angle is greater than 60°, in addition to the Rayleigh waves propagating upwards along the crack surface, some Rayleigh waves convert to shear waves at the crack bottom tip and then incident on the free surface of the workpiece. Experiments were carried out to validate some of the Rayleigh wave propagation paths. The experimental results matched the theoretical arrival time well, thus verifying the reliability of the analytical wave path. The results are helpful for the quantitative detection of subsurface inclined cracks using laser ultrasonic techniques.

Broadband THz metasurface bandpass filter/antireflection coating based on metalized Si cylindrical rings

Abstract The operation of the metasurface based on silicon cylindrical rings coated by a gold as a terahertz (THz) bandpass filter/antireflection structure is studied. The decrease in the reflectance is conditioned by the destructive interference of electromagnetic waves reflected from structural layers of the metasurface. An efficient antireflection band with the reflectance below 10% is formed in the frequency spectrum of 0.71–0.92 THz with a relative bandwidth of 26%. It is shown that the operating spectrum of the suggested metasurface can be varied by changing the total radius of cylindrical rings, whereas the filter’s performance is rather insensitive to the variations in cylinder height and inner radius. The dependence of the antireflection band on the polarization and incidence angle of the THz waves is also analyzed. The antireflection band is sensitive to changes in the surrounding medium, hence enabling control of the transmittance band by exploiting refractive-index-changing materials.

Computational study on complex wave hydrodynamics of multidirectional extreme waves at fringing reef

Due to the special topography of coral reefs, most of the energy of incoming waves can be reduced through the wave breakings at the reef edge, thus playing a vital role in protecting the coastal regions behind coral reefs. In severe weather conditions, wind-generated waves tend to propagate in multiple directions in real ocean environment. The abrupt multidirectional wave-focusing mechanisms may result in the formation of huge waves in both deep and shallow waters, threatening the safety of coastal infrastructures and well-being of coastal communities at the reef islands. Therefore, it is crucial to thoroughly examine the complex hydrodynamics of multidirectional extreme waves at the fringing reefs. Previous studies on reef wave hydrodynamics predominantly focus on unidirectional waves. Only a small number of studies have considered the wave hydrodynamics of fringing reefs under multidirectional waves, let alone multidirectional extreme waves. To address the knowledge gap in previous studies, this study analyzes the complex hydrodynamic processes of multidirectional extreme waves at the fringing reef by applying a nonhydrostatic numerical wave model (NHWAVE). Influences of several key factors are analyzed, such as significant wave height, peak wave period, water depth, incident wave angle, and reef topography. It is expected that the research findings of this study will contribute to a deeper understanding of the hydrodynamics of multidirectional extreme waves at the fringing reef.

Influence of ground motion parameters on seismic response of a large-longitudinal-slope and small-radius curved girder bridge.

The mechanical performance of curved bridges under the action of an earthquake is complex. To obtain the influence of seismic parameters on the seismic response of curved girder bridges, this paper relies on a large slope small-radius curved steel box girder bridge (LSCGB) and selects seismic wave incidence angle, vertical component of ground motion, and site category as seismic parameters to carry out nonlinear time history analysis. Based on the analysis results of the case bridge, it is shown that the torsional vibration of the first 10 modes of LSCGB is significant, the modes are dispersed, and the contribution of high-order modes of vibration cannot be ignored. The most unfavorable seismic wave incidence angle is in the direction of 45°∼60° counterclockwise Angle from the central connection line of Pier No. 1 and Pier No. 4 of the bridge. The seismic response of the curved bridge components increases with the vertical seismic intensity, and the influence on displacement responses is more significant. The basic vibration period of curved girder bridges built on soft soil sites is extended by approximately 18.23%, and the seismic response of key components increases with the softening of the site soil. Therefore, when analyzing the seismic response of LSCGBs, the influence of vertical component of ground motion and site category should not be ignored.

A bidirectional broadband multifunctional terahertz device based on vanadium dioxide

We propose a switchable and bidirectional metamaterial terahertz absorber/reflector based on vanadium dioxide. Simulation results show that the device can achieve broadband absorption and reflection by adjusting the phase state of vanadium dioxide. When the vanadium dioxide is in the insulating phase, the device can be a perfect bidirectional broadband absorber with more than 90% absorptivity under normal forward/backward TE/TM-polarized incident wave in the frequency range from 1.9 THz to 10 THz. This excellent absorption characteristic is insensitive to the polarization angle of the incident wave. When the vanadium dioxide is in the metal phase, the device achieves more than 90% reflectivity in the low-frequency range under forward/backward TE-polarized incident wave. In the high-frequency range, the device has both reflection and absorption effects on the TE-polarized wave. When vanadium dioxide is in two different phase states, the performance of device is not sensitive to the electromagnetic wave incidence angle. The change in vanadium dioxide conductivity does not affect the stable performance of the device at high (low) frequencies when used as an absorber (reflector). The proposed device has potential applications in terahertz energy harvesting, detection, sensing, optical switching, and shielding.

Propagation characteristics of obliquely incident terahertz waves in inhomogeneous microplasma

The transmission characteristics of terahertz (THz) waves in a non-uniform microplasma are investigated by using the scattering matrix method. The electron density distribution in microplasma is simulated by Epstein and parabolic models. The effects of physical parameters, such as the incidence angle of THz waves, microplasma size, electron density, and collision frequency, on the propagation of THz waves are numerically analyzed. The results show that lower frequency THz waves are difficult to penetrate the microplasma with high electron density and high collision frequency. The microplasma density distribution, especially the gradient variation of the density in the first layer, has a large effect on the reflection of THz waves. Thus, THz waves can be used to diagnose the physical parameters of microplasmas.

Tunable terahertz filter based on graphene photonic crystals with defective layers

In this paper, we design a high-precision tunable terahertz filter by using transfer matrix method. The filter structure mainly consists of graphene embedded photonic crystals (GPCs). The front part of the GPCs contains artificial synthetic material and air layer, the back part of the GPCs is composed by and periodic stack of isotropic dielectric slabs (MgF2) embedded with graphene sheets, where air defect layer is located in the middle of the GPC as a central layer. Our simulation reveals that graphene layer and air defective layer strongly affect the filter performance. And we can get a relatively pure transmission peak in a wide frequency region. Additionally, the influence of incidence angle of terahertz wave, thickness of air layer, the unit number of front periodic structure and chemical potential of the graphene sheets can also modulate the function of the filter. And the filter has strong stability when the temperature changes from 150 K to 350 K.The results indicate that single channel, dual and multiple channels filter in a narrow frequency can be obtained by optimizing the structure parameter.