
In this paper, we study an adaptive edge finite element method for time-harmonic Maxwell's equations in metamaterials. A-posteriori error estimators based on the recovery type and residual type are proposed, respectively. Based on our a-posteriori error estimators, the adaptive edge finite element method is designed and applied to simulate the backward wave propagation, electromagnetic splitter, rotator, concentrator and cloak devices. Numerical examples are presented to illustrate the reliability and efficiency of the proposed a-posteriori error estimations for the adaptive method.
Citation: Hao Wang, Wei Yang, Yunqing Huang. An adaptive edge finite element method for the Maxwell's equations in metamaterials[J]. Electronic Research Archive, 2020, 28(2): 961-976. doi: 10.3934/era.2020051
[1] | Gunduz Caginalp . Surface tension, higher order phase field equations, dimensional analysis and Clairaut’s equation. AIMS Mathematics, 2017, 2(2): 207-214. doi: 10.3934/Math.2017.2.207 |
[2] | Yellamma, N. Manjunatha, Umair Khan, Samia Elattar, Sayed M. Eldin, Jasgurpreet Singh Chohan, R. Sumithra, K. Sarada . Onset of triple-diffusive convective stability in the presence of a heat source and temperature gradients: An exact method. AIMS Mathematics, 2023, 8(6): 13432-13453. doi: 10.3934/math.2023681 |
[3] | Khalil Ur Rehman, Wasfi Shatanawi, Zeeshan Asghar, Haitham M. S. Bahaidarah . Neural networking analysis for MHD mixed convection Casson flow past a multiple surfaces: A numerical solution. AIMS Mathematics, 2023, 8(7): 15805-15823. doi: 10.3934/math.2023807 |
[4] | Jonathan D. Evans, Morgan L. Evans . Stress boundary layers for the Giesekus fluid at the static contact line in extrudate swell. AIMS Mathematics, 2024, 9(11): 32921-32944. doi: 10.3934/math.20241575 |
[5] | Jean-Paul Pascal, Serge D’Alessio, Syeda Rubaida Zafar . The instability of liquid films with temperature-dependent properties flowing down a heated incline. AIMS Mathematics, 2019, 4(6): 1700-1720. doi: 10.3934/math.2019.6.1700 |
[6] | Khalil Ur Rehman, Wasfi Shatanawi, Zead Mustafa . Artificial intelligence (AI) based neural networks for a magnetized surface subject to tangent hyperbolic fluid flow with multiple slip boundary conditions. AIMS Mathematics, 2024, 9(2): 4707-4728. doi: 10.3934/math.2024227 |
[7] | Melek Erdoğdu, Ayșe Yavuz . On differential analysis of spacelike flows on normal congruence of surfaces. AIMS Mathematics, 2022, 7(8): 13664-13680. doi: 10.3934/math.2022753 |
[8] | Nevin Gürbüz, Dae Won Yoon . Geometry of curve flows in isotropic spaces. AIMS Mathematics, 2020, 5(4): 3434-3445. doi: 10.3934/math.2020222 |
[9] | T. Mahesh Kumar, Nehad Ali Shah, V. Nagendramma, P. Durgaprasad, Narsu Sivakumar, B. Madhusudhana Rao, C. S. K. Raju, Se-Jin Yook . Linear regression of triple diffusive and dual slip flow using Lie Group transformation with and without hydro-magnetic flow. AIMS Mathematics, 2023, 8(3): 5950-5979. doi: 10.3934/math.2023300 |
[10] | Haifaa Alrihieli, Musaad S. Aldhabani, Ghadeer M. Surrati . Enhancing the characteristics of MHD squeezed Maxwell nanofluids via viscous dissipation impact. AIMS Mathematics, 2023, 8(8): 18948-18963. doi: 10.3934/math.2023965 |
In this paper, we study an adaptive edge finite element method for time-harmonic Maxwell's equations in metamaterials. A-posteriori error estimators based on the recovery type and residual type are proposed, respectively. Based on our a-posteriori error estimators, the adaptive edge finite element method is designed and applied to simulate the backward wave propagation, electromagnetic splitter, rotator, concentrator and cloak devices. Numerical examples are presented to illustrate the reliability and efficiency of the proposed a-posteriori error estimations for the adaptive method.
Free-surface flow over a submerged obstacles is one of the well-known classical problem in fluid mechanics. Many researchers have investigated free surface flows over an obstacle for different bottoms topography, for examples, Forbes and Schwartz [13] used the boundary integral method to find fully non-linear solutions of subcritical and supercritical flows over a semi-circular obstacle. Supercritical and critical flows over a submerged triangular obstacle were investigated by Dias and Vanden-Broeck [11]. They employed a series truncation methods to calculate the solutions. Abd-el-Malek and Hanna [1] studied flow over a triangular obstacle by using the Hilbert method with gravity effect. When the fluid is subjected to the interaction of gravity and surface tension, in this case, the problem is generally difficult to solve. Forbes [12] was among the first to propose numerical solutions of non-linear flows over a semi-circular obstruction under the effect of gravity and surface tension. Later, many authors have studied this problem, for example, Grandison [14], Vanden-Broeck [24]. In the case of flows over two obstacles, Pratt [23] investigated this problem experimentally and theoretically using weakly non-linear analysis. Later, Belward [3] computed numerical solutions of a critical flow for which the hydraulic fall occurred at the leftmost obstacle with downstream supercritical flow. Recently, Binder, Vanden-Broeck and Dias [9] showed that there exist two types of solution in subcritical flow regime, and one type in supercritical flow regime. This paper is concerned with the numerical calculations of flow of finite depth over a successive obstacles. The purpose of this work is to examine further flows with many obstacles in order to classify the possible solutions. The effects of surface tension and gravity are included in the boundary conditions, the problem is solved numerically by using the boundary integral equation methods. These methods are based on a reformulation of the problem as a system of non-linear integro-differential equations for the unknown quantities on the free surface. These equations are then discretized and the resultant non-linear algebraic equations is solved by iteration. Such boundary integral equation methods have been used extensively by many researchers [3,4,6,7,15,17] and others. It is assumed that there is uniform flow far upstream where the flow approaches a uniform stream with constant velocity and depth (see Figure 1). The problem is characterized by the two parameters the Froude number defined by
(1.1) |
and the inverse Weber number where
(1.2) |
Here is the surface tension, is the gravity and is the fluid density. When , the flow is called subcritical and for the flow far upstream is called supercritical. In this work, we calculate waveless solutions for both supercritical and subcritical flows by introducing the effects of surface tension.
Formulation of the problem and numerical procedure are given in section and section respectively. In section we discuss the numerical results of free surface flows over a successive triangular obstacles with different angles , and for various values of the two parameters and . Solution diagrams for all flow regimes are presented.
We consider steady two-dimensional potential free surface flows past a submerged obstacles at the bottom of a channel (see Figure 1). The flow is assumed to be inviscid and incompressible. Fluid domain is bounded below by a horizontal rigid wall and the successive obstructions forming the angles , with the horizontal, where , and above by the free surface . Let us introduce Cartesian coordinates with the axis along the bottom and the axis directed vertically upwards, gravity is acting in the negative direction. Let's introduce the velocity potential and the stream function by defining the complex potential function as
(2.1) |
The complex velocity can be written as
(2.2) |
Here and are velocity components in the and directions, and . For convenience, we define dimensionless variables by taking as the reference length and as the reference velocity. Without loss of generality, we choose on the free surface . By the choice of our dimensionless variables, we have on the bottom and at the point (see Figure 2).
The problem is formulated in terms of the velocity potential . This function satisfies Laplace's equation
The Bernoulli's equation on the free surface can be written
(2.3) |
Here is curvature of the free surface, and are defined in (1.1) and (1.2) respectively.
The kinematic boundary conditions in -plane are given by
(2.4) |
Now we reformulate the problem as an integral equation. We define the function by
(2.5) |
and we map the flow domain onto the upper half of the plane by the transformation
(2.6) |
The flow in the plane is shown in Figure 3. The curvature of a streamline, in terms of , is given by
(2.7) |
Substituting (2.7) into (2.3), Bernoulli's equation becomes
(2.8) |
We apply the Cauchy's integral formula to the function in the complex -plane with a contour consisting of the axis and the semicircle of arbitrary large radius in the upper half plane. After taking the real part and , we obtain
(2.9) |
Where and denote the value of and on the free surface. The integral in is a Cauchy principal value type.
The kinematic boundary conditions (2.4) in -plane become
(2.10) |
By using (2.10), Eq (2.9) becomes :
(2.11) |
Rewriting this equation in terms of by substituting , ,
this gives
(2.12) |
Here and . The Eq (2.8) is now rewritten in terms of and as
(2.13) |
Integrating the identity
(2.14) |
We obtain the following parametric representation of the free surface
(2.15) |
(2.16) |
By substituting (2.16) into (2.13), an integro-differential equation is created and it is solved numerically.
In this section, we describe numerical approach for the nonlinear problem derived in previous section. This numerical procedure has been successfully used by B. J. Binder [9], P. Guayjarernpanishk [15] and others for solving nonlinear integral equations. Firstly, the free surface must be truncated at and for the corresponding far upstream , and far downstream , respectively. The truncated free surface is then discretized into equally segments with
(3.1) |
and the unknown variables on the free surface are
Here is the mesh spacing. We evaluate the values at the midpoints
(3.2) |
by applying the trapezoidal rule to the integral in (2.12) with summations over such that is the midpoints. We evaluate and by applying the Euler's method and by using (2.14). This yields
and
Here We now satisfy (2.13) at the midpoints (3.2). This yields non-linear algebraic equations for the unknowns . The derivative, , at the mesh points , is approximated by a finite difference, whereby
The system of equations with unknowns is solved by Newton's method.
The numerical procedure of section is used to compute solutions for free surface flows over a successive triangular obstacles. For simplicity, we assume that the triangles are isosceles forming the angles with the horizontal (see Figure 1). Also, we choose . Most of the calculations in this paper are obtained with and . For a given values of at the points , we compute waveless solutions for various values of the angles , Froude number and the inverse Weber number . We denote by which represents the length of the sides of the triangles. The problem is essentially characterized by four parameters; The Froude number , the inverse of Weber , and L. In supercritical () or subcritical flow () and a fixed values of and , the effect of surface tension on the shape of free surface, is shown in Figures 4 and 5. It should be noted that the free surface elevation increases when the inverse Weber number decreases. The Figures 6 and 7 show the effect of the Froude number for fixed values of , and . It can be seen that the elevation of the free surfaces increases as Froude number increases. When the surface tension is neglected and (without gravity) and for an arbitrary values of and ; the problem has an exact solution that can be computed via the streamline method due to Kirchhoff [2], in this case, the elevation of free surfaces reaches its maximum. The effect of varying the length , whilst , and are fixed is shown in Figure 8. Figure 9 illustrates the effect of varying the angles where and are fixed.
In this paper, the problem of irrotational, two-dimensional free-surface flow over a successive obstacles has been presented. The fluid is assumed to be incompressible and inviscid. The fully non-linear problem is formulated by using a boundary integral equation technique. The numerical solutions are obtained, in the presence of surface tension and gravity. For supercritical flow , there is a three-parameters family of solutions depending on the height of obstacle, the inverse Weber number and the Froude number which is similar to the case of subcritical flow . We have seen the effect of surface tension on free surface profiles for supercritical and subcritical flows. It noted that when the inverse Weber number decreases or the Froude number increases, the free surface elevation increases. The same observation is made when or decreases the elevation of the free surface decreases and vice versa. The maximum free-surface elevation is obtained in the absence of the surface tension and gravity, in this case, the exact solution can be found via the hodograph transform due to Kirchhoff [2].
The authors declare no conflicts of interest in this paper.
[1] |
Plasmon resonance with finite frequencies: A validation of the quasi-static approximation for diametrically small inclusions. SIAM J. Appl. Math. (2016) 76: 731-749. ![]() |
[2] |
Nearly cloaking the full Maxwell equations: Cloaking active contents with general conducting layers. J. Math. Pures. Appl. (9) (2014) 101: 716-733. ![]() |
[3] |
A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. (1994) 114: 185-200. ![]() |
[4] |
Localization and geometrization in plasmon resonances and geometric structures of Neumann-Poincaré eigenfunctions. ESAIM Math. Model. Numer. Anal. (2020) 54: 957-976. ![]() |
[5] | Analysis of a Cartesian PML approximation to the three dimensional electromagnetic wave scattering problem. Int. J. Numer. Anal. Model. (2012) 9: 543-561. |
[6] |
An adaptive finite element method for two-dimensional Maxwell's equations. J. Sci. Comput. (2013) 55: 738-754. ![]() |
[7] |
An adaptive finite element method for two-dimensional transverse magnetic time harmonic Maxwell's equations with general material properties and general boundary conditions. J. Sci. Comput. (2016) 68: 848-863. ![]() |
[8] |
A recovery-based a posteriori error estimator for (curl) interface problems. Comput. Methods. Appl. Mech. Engrg. (2015) 296: 169-195. ![]() |
[9] |
Recovery-based error estimators for interface problems: Mixed and nonconforming finite elements. SIAM J. Numer. Anal. (2010) 48: 30-52. ![]() |
[10] |
Flux recovery and a posteriori error estimators: Conforming elements for scalar elliptic equations. SIAM J. Numer. Anal. (2010) 48: 578-602. ![]() |
[11] |
Multigrid methods for two-dimensional Maxwell's equations on graded meshes. J. Comput. Appl. Math. (2014) 255: 231-247. ![]() |
[12] |
Full and partial cloaking in electromagnetic scattering. Arch. Ration. Mech. Anal. (2017) 223: 265-299. ![]() |
[13] |
A convergent adaptive algorithm for Poisson's equation. SIAM J. Numer. Anal. (1996) 33: 1106-1124. ![]() |
[14] | Y. Hao and R. Mittra, FDTD modeling of metamaterials: Theory and applications, Artech. House., (2008). |
[15] |
B. He, W. Yang and H. Wang, Convergence analysis of adaptive edge finite element method for variable coefficient time-harmonic Maxwell's equations, J. Comput. Appl. Math., 376 (2020), 16pp. doi: 10.1016/j.cam.2020.112860
![]() |
[16] |
The averaging technique for superconvergence: Verification and application of 2D edge elements to Maxwell's equations in metamaterials. Comput. Methods Appl. Mech. Engrg. (2013) 255: 121-132. ![]() |
[17] |
Interior penalty DG methods for Maxwell's equations in dispersive media. J. Comput. Phys. (2011) 230: 4559-4570. ![]() |
[18] |
Y. Huang, J. Li and W. Yang, Modeling backward wave propagation in metamaterials by the finite element time-domain method, SIAM J. Sci. Comput., 35 (2013), B248–B274. doi: 10.1137/120869869
![]() |
[19] |
The superconvergent cluster recovery method. J. Sci. Comput. (2010) 44: 301-322. ![]() |
[20] |
On quasi-static cloaking due to anomalous localized resonance in . SIAM J. Appl. Math. (2015) 75: 1245-1260. ![]() |
[21] |
Analysis of electromagnetic scattering from plasmonic inclusions beyond the quasi-static approximation and applications. ESAIM Math. Model. Numer. Anal. (2019) 53: 1351-1371. ![]() |
[22] |
Analysis and application of the nodal discontinuous Galerkin method for wave propagation in metamaterials. J. Comput. Phys. (2014) 258: 915-930. ![]() |
[23] |
An adaptive edge finite element method for electromagnetic cloaking simulation. J. Comput. Phys. (2013) 249: 216-232. ![]() |
[24] |
H. Liu, Virtual reshaping and invisibility in obstacle scattering, Inverse Problems, 25 (2009), 16pp. doi: 10.1088/0266-5611/25/4/045006
![]() |
[25] |
(2003) Finite Element Methods for Maxwell's Equations. New York: Oxford University Press. ![]() |
[26] |
A posteriori error estimates based on the polynomial preserving recovery. SIAM J. Numer. Anal. (2004) 42: 1780-1800. ![]() |
[27] |
Hybridizable discontinuous Galerkin methods for the time-harmonic Maxwell's equations. J. Comput. Phys. (2011) 230: 7151-7175. ![]() |
[28] |
Controlling electromagnetic fields. Science (2006) 312: 1780-1782. ![]() |
[29] |
Adjoint recovery of superconvergent functionals from PDE approximations. SIAM. Rev. (2000) 42: 247-264. ![]() |
[30] | A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, Artech House, Inc., Boston, MA, 2000. |
[31] |
Adaptive finite element method for the sound wave problems in two kinds of media. Comput. Math. Appl. (2020) 79: 789-801. ![]() |
[32] |
D. H. Werner and D.-H. Kwon, Transformation Electromagnetics and Metamaterials. Fundamental Principles and Applications, Springer-Verlag, London, 2014. doi: 10.1007/978-1-4471-4996-5
![]() |
[33] |
Developing a time-domain finite element method for the Lorentz metamaterial model and applications. J. Sci. Comput. (2016) 68: 438-463. ![]() |
[34] |
Modeling and analysis of the optical black hole in metamaterials by the finite element time-domain method. Comput. Methods Appl. Mech. Engrg. (2016) 304: 501-520. ![]() |
[35] |
Mathematical analysis and finite element time domain simulation of arbitrary star-shaped electromagnetic cloaks. SIAM J. Numer. Anal. (2018) 56: 136-159. ![]() |
[36] |
Developing finite element methods for simulating transformation optics devices with metamaterials. Commun. Comput. Phys. (2019) 25: 135-154. ![]() |
[37] |
The superconvergence patch recovery (SPR) and adaptive finite element refinement. Comput. Methods. Appl. Mech. Engrg. (1992) 101: 207-224. ![]() |
1. | Andrew Nichols, Matteo Rubinato, Yun-Hang Cho, Jiayi Wu, Optimal Use of Titanium Dioxide Colourant to Enable Water Surfaces to Be Measured by Kinect Sensors, 2020, 20, 1424-8220, 3507, 10.3390/s20123507 | |
2. | Abdelkader Laiadi, 2020, Numerical solution of flows over a trapezoidal obstacles, 978-1-7281-2580-0, 121, 10.1109/ICMIT47780.2020.9046993 | |
3. | Abdelkader Laiadi, The influence of surface tension and gravity on cavitating flow past an inclined plate in a channel, 2022, 0033-569X, 10.1090/qam/1617 | |
4. | Abdelkader Laiadi, Free surface flow over a trapezoidal cavity with surface tension effect, 2025, 3, 2972-4805, 2311, 10.59400/jam2311 |