Research article

Chebyshev pseudospectral approximation of two dimensional fractional Schrodinger equation on a convex and rectangular domain

  • In this article, the authors report the Chebyshev pseudospectral method for solving twodimensional nonlinear Schrodinger equation with fractional order derivative in time and space both. The modified Riemann-Liouville fractional derivatives are used to define the new fractional derivatives matrix at CGL points. Using the Chebyshev fractional derivatives matrices, the given problem is reduced to a diagonally block system of nonlinear algebraic equations, which will be solved using Newton's Raphson method. The proposed methods have shown error analysis without any dependency on time and space step restrictions. Some model examples of the equations, defined on a convex and rectangular domain, have tested with various values of fractional order α and β. Moreover, numerical solutions are demonstrated to justify the theoretical results.

    Citation: A. K. Mittal, L. K. Balyan. Chebyshev pseudospectral approximation of two dimensional fractional Schrodinger equation on a convex and rectangular domain[J]. AIMS Mathematics, 2020, 5(3): 1642-1662. doi: 10.3934/math.2020111

    Related Papers:

    [1] Haifa Bin Jebreen, Beatriz Hernández-Jiménez . An approach based on the pseudospectral method for fractional telegraph equations. AIMS Mathematics, 2023, 8(12): 29221-29238. doi: 10.3934/math.20231496
    [2] Haifa Bin Jebreen, Beatriz Hernández-Jiménez . Pseudospectral method for fourth-order fractional Sturm-Liouville problems. AIMS Mathematics, 2024, 9(9): 26077-26091. doi: 10.3934/math.20241274
    [3] M. Mossa Al-Sawalha, Azzh Saad Alshehry, Kamsing Nonlaopon, Rasool Shah, Osama Y. Ababneh . Fractional view analysis of delay differential equations via numerical method. AIMS Mathematics, 2022, 7(12): 20510-20523. doi: 10.3934/math.20221123
    [4] Iman Ben Othmane, Lamine Nisse, Thabet Abdeljawad . On Cauchy-type problems with weighted R-L fractional derivatives of a function with respect to another function and comparison theorems. AIMS Mathematics, 2024, 9(6): 14106-14129. doi: 10.3934/math.2024686
    [5] Zahra Pirouzeh, Mohammad Hadi Noori Skandari, Kamele Nassiri Pirbazari, Stanford Shateyi . A pseudo-spectral approach for optimal control problems of variable-order fractional integro-differential equations. AIMS Mathematics, 2024, 9(9): 23692-23710. doi: 10.3934/math.20241151
    [6] Dina Abdelhamid, Wedad Albalawi, Kottakkaran Sooppy Nisar, A. Abdel-Aty, Suliman Alsaeed, M. Abdelhakem . Mixed Chebyshev and Legendre polynomials differentiation matrices for solving initial-boundary value problems. AIMS Mathematics, 2023, 8(10): 24609-24631. doi: 10.3934/math.20231255
    [7] Alessandra Jannelli, Maria Paola Speciale . On the numerical solutions of coupled nonlinear time-fractional reaction-diffusion equations. AIMS Mathematics, 2021, 6(8): 9109-9125. doi: 10.3934/math.2021529
    [8] Md. Asaduzzaman, Md. Zulfikar Ali . Existence of positive solution to the boundary value problems for coupled system of nonlinear fractional differential equations. AIMS Mathematics, 2019, 4(3): 880-895. doi: 10.3934/math.2019.3.880
    [9] Lakhlifa Sadek, Tania A Lazǎr . On Hilfer cotangent fractional derivative and a particular class of fractional problems. AIMS Mathematics, 2023, 8(12): 28334-28352. doi: 10.3934/math.20231450
    [10] Adisorn Kittisopaporn, Pattrawut Chansangiam . Approximate solutions of the $ 2 $D space-time fractional diffusion equation via a gradient-descent iterative algorithm with Grünwald-Letnikov approximation. AIMS Mathematics, 2022, 7(5): 8471-8490. doi: 10.3934/math.2022472
  • In this article, the authors report the Chebyshev pseudospectral method for solving twodimensional nonlinear Schrodinger equation with fractional order derivative in time and space both. The modified Riemann-Liouville fractional derivatives are used to define the new fractional derivatives matrix at CGL points. Using the Chebyshev fractional derivatives matrices, the given problem is reduced to a diagonally block system of nonlinear algebraic equations, which will be solved using Newton's Raphson method. The proposed methods have shown error analysis without any dependency on time and space step restrictions. Some model examples of the equations, defined on a convex and rectangular domain, have tested with various values of fractional order α and β. Moreover, numerical solutions are demonstrated to justify the theoretical results.


    1. Introduction

    The development of new materials is often the beginning of innovative and future-oriented technologies. This is necessary to meet tomorrow’s engineering applications, for example, to increase efficiency, reliability and productivity or to reduce energy consumption and with this to raise standard of living [1,2].

    It is also possible to understand the term“new material”as a combination of materials. Joining metals and ceramics for example, allows the combination of the positive properties of both. Metallic materials offer magnetic properties, a high ductility or an electrical conductivity, while ceramic materials have a high strength even at elevated temperatures, a high corrosion resistance and, in the case of zirconia, an oxygen-ion conductivity.

    These combinations are necessary in solid oxide fuel cells, gas separation membranes, thermoelectric generators, bipolar surgical tools or in mechanical engineering applications. Major obstacle for using these composites is, that existing manufacturing technologies cannot fulfil requirements in terms of costs, complex geometries or operation conditions.

    Conventional frictional and form-locking joining techniques, like screwing or shrinking, offer a simple possibility to combine both materials. However, these joining techniques have many disadvantages regarding wall thicknesses and geometrical design.

    Conventional cohesive joining techniques, like gluing or brazing, allow chemical bonding between both materials. Especially gluing allows a simple joining, while for brazing a special development of brazing material is necessary [3]. However, within both processes a further material is inserted which can limit functionality, operation temperature or strength. Nevertheless, brazing is already well established for joining power electronics or stacks of solid oxide fuel cells.

    It is generally known that in addition to conventional joining techniques several powder technological approaches exist. One of these approaches is thermal spraying. This is used for applying a metallic or ceramic layer on top of ceramic or metallic substrate [4]. This allows the production of electrical circuits on insulating substrates, wear-resistant layers or thermal barrier coatings. Layer thicknesses thereby range from 30 µm to some millimetres. The high layer thickness is a disadvantage as well as the high porosity and high surface roughness. Furthermore, the coating of complex geometries is quite difficult.

    One promising approach to simplify the production of composites and achieving a high degree of freedom regarding design is a co-shaping, continued with co-sintering. Especially the co-sintering step is the challenge of this manufacturing route. Differences in coefficient of thermal expansion (CTE) or in densifying behaviour (shrinkage) can lead to stresses during heating or cooling [5]. These stresses can cause flaws as cracks, delaminations, pores, changes in microstructure or in the case of asymmetric composites a warpage. Cai et al. has theoretically analysed occurring stresses during sintering of composites [5]. He considered that the materials exhibit during sintering pure viscous behaviour and developed, based on this assumption, equations for maximum stresses σK in bi-layered asymmetric (equation 1) and three-layered symmetric (equation 2) composites.

    σK=[m4+mnn2+2mn(2m2+3m+2+m(4))]ηΔ˙ε (1)
    σK=11+mnηΔ˙ε (2)

    Especially his equation for the rate of warping (equation 3), which was derived from equation 1, is commonly used for a comparison of theory and practice [6,7,8].

    ˙k=6(m+1)2mnm4n2+2mn(2m2+3m+2)+1Δ˙ε (3)

    It can be derived from these equations that besides of the difference in CTE and sintering behaviour, which is represented as difference in sintering strain rate Δ˙ε̇, the ratio of layer thickness m and the material viscosity η respectively ratio of the material viscosities n influence stresses and warpage.

    Furthermore, it is known that non-or low-shrinking substrates can harm and slow down densifying of thin ceramic layers during co-sintering, which result in lower densities in comparison to free shrinking materials [9,10,11]. In this context Yamaguchi et al. and Muecke et al. showed that high densities in co-sintered ceramic layers can only be reached, if differences in strain rate are small enough or if the negative effect is compensated by an increased isothermal dwell time [11,12].

    The first step to avoid flaws during co-sintering is the choice of materials with well adapted CTE’s. The second step is the reduction of mismatches in sintering strain rate. The main focus of investigations for the second step is the adaptation of green densities, use of specific particle sizes, changing heating rates, using doping elements or integrating a material gradient [13,14,15,16,17]. In this work, sintering behaviour was influenced and adjusted by treating metallic powder with a high energy milling step [18,19].

    The aim of this work was to achieve a multi-layered composite consisting of a thin ceramic layer on top or between two thicker porous metallic layers. These composites are suitable as semi-finished metal supported fuel cells, filter elements or gas separation membranes. Thereby, investigations about the co-sintering behaviour of adapted and un-adapted composites were done. A combination of zirconia and stainless steel as materials was used. The composites were produced by using a tape casting process, which is the standard process in ceramic industry for the production of large-area, thin and flat ceramic substrates [20].


    2. Materials and Method

    For the investigations 3 mol% yttrium stabilized zirconia TZ-3YS-E (Tosoh Inc., Japan) was selected as ceramic material and the gas atomized high temperature corrosion resistant ferritic iron-chromium alloy Crofer®22APU (H. C. Starck GmbH, Germany) as metallic one. The metallic powder was treated with a high energy milling step, to influence its sintering behaviour. Further information about this can be found in earlier publications [18,19,21]. Overall, four different metallic powders (metal powder 1-4) were taken, which differentiate in terms of their sintering behaviour. Crofer®22APU was chosen due to its melting temperature of 1510°C which is above the sintering temperature of 1450°C of the zirconia powder.

    Different metal-ceramic (asymmetric) and metal-ceramic-metal (symmetric) composites were produced by tape casting. For this a slurry consisting of water, polyvinylalcohol as binder and glycerine as plasticizer was prepared. Polyvinylalcohol and glycerine were taken due to their good decomposition behaviour under hydrogen atmosphere, which is important to avoid a harmful increase of carbon content in steel after sintering [22]. Besides binder and plasticizer, the defoamer Foammaster F111 (BASF SE, Germany) and the surfactant Glycol N 109 (Zschimmer & Schwarz GmbH % Co. KG, Germany) were used. The composition of the individual slurries is shown in table 1. More details concerning the production process are published in [18,22,23].

    Table 1.Composition of ceramic and metallic slurries for tape casting in wt%.
    powderwaterbinderplasticizerde-foamersurfactant
    TZ-3YS-E44.545.14.35.70.20.2
    metal powder 143.335.74.36.50.10.1
    metal powder 280.515.61.620.10.2
    metal powder 3 65.328.133.30.10.2
    metal powder 480.515.61.620.10.2
     | Show Table
    DownLoad: CSV

    Tape casting allows the production of asymmetric and symmetric composites. Asymmetric composites were produced with the combination TZ-3Y-SE and metal powder 1 (composite 1a) respectively powder 2 (composite 2). Symmetric composites were produced with metal powder 3 (composite 3), powder 4 (composite 4) and again with powder 1 (composite 1s). The ceramic green layer thickness was in all cases 11.6±1 µm while the metallic layer thickness varied from 200 to 400 µm. Hence, the metallic layer thickness is 17 to 34 times thicker than the ceramic one. It has to be noted that the specific metallic layer thickness has no influence on the shown results, for which reason it is not specified in detail. Figure 1 exhibits an asymmetric metal-ceramic multi-layered green tape out of composite 2. The ceramic layer thickness was measured before and after sintering using field emission electron scanning microscopy (FESEM) (NVision 40, Carl Zeiss SMT GmbH, Germany).

    Figure 1. Metal-ceramic asymmetric green tape with the thin ceramic on top of the thicker metallic layer.

    Round specimens with a diameter of 20 mm were punched out of the tapes for co-sintering experiments. These tapes were debinded and sintered on top and between two porous, 1 mm thick Keralpor 99 alumina sintering supports (Kerafol Keramische Folien GmbH, Germany) respectively. Covering the tapes with a sintering support was necessary to avoid warpage during sintering. Furthermore, two smaller specimens with a rectangular shape and dimensions of 5×7 mm out of composite 1a and 2 were sintered without a coverage.

    The specimens were debinded under hydrogen atmosphere with a heating rate of 1 K/min up to 600°C. After a dwell time of 2 h the furnace (RRO280/300-900V, MUT advanced heating, Germany) was cooled down to room temperature with a rate of 3 K/min. Sintering took place in a separated furnace (HTBL20W22-2G, Carbolite Gero GmbH & Co. KG, Germany) with a heating rate of 3 K/min up to the sintering temperature between 1200 and 1400°C and a cooling rate of 4 K/min. The dwell time was 2 h and the sintering atmosphere was a mixture of argon (80%) and hydrogen (20%).

    The sintering behaviour of the materials was characterized with high temperature heating microscope (EHM 201-17K, Hesse Instruments, Germany). Thereby a camera records pictures of the specimen during sintering in a tube furnace. A following image analyses allows a calculation of linear shrinkage as well as sintering strain rate. During this measurement the tube furnace was flushed with gas (95% argon, 5% hydrogen). Heating rate was 1 K/min up to 600°C and then 3 K/min up to 1400°C. The used specimens were 1 mm in thickness and had a diameter of 5 mm. Using a heating microscope instead of pushing rod dilatometer, had the advantage to get free shrinkage curves of the materials. When using a pushing rod dilatometer, the pressure of the pushing rod can especially falsify the measurement of the metallic material due to its low viscosity at high temperatures.

    The pushing rod dilatometer Dil 402 D (Netsch Gerätebau GmbH, Germany) was used for measuring the CTEs. The CTE was measured by means of sintered laminates (density > 99%) of each material with dimensions of 4×4×20 mm3.

    The sintered composites were cut in the middle and were prepared by grinding and polishing with abrasive papers for microstructural investigations with the field emission scanning microscope NVision 40 (Carl Zeiss SMT GmbH, Germany), . The green tapes, however, were prepared using broad ion beam method. 10-15 pictures with a magnification of 7000 and a backscattered electron detector, which shows the material contrast, were taken to analyse porosity. The porosity was afterwards determined by image analysing with the software ImageJ V1.6.0 (National Institute of Health, USA).

    The viscosity of the ceramic and metallic materials during sintering was measured using the bending beam method. This method is usually used to measure the viscosity of glasses [24] but was adapted by Lee et al. [25] and Lame et al. [26] to measure the viscosity of oxide ceramics and steels during sintering. Thereby, a beam is sintered on top of a special support which has only two contact points at the left and right side of the tape. It is then possible to calculate the viscosity η from the rate of deflection ˙δ using equation 4 [27].

    η=5 ρgL432 ˙δh2 (4)

    The span length L during the measurements amounted 10 mm. After preliminary tests, the thickness h of the tapes was set to 100 µm for the ceramic and 500 µm for the metallic. These tapes were pre-sintered up to 1000°C to avoid a bending during de-binding and to give the tapes some strength.

    Before and after co-sintering, the ceramic layer of the asymmetric composites was investigated in terms of its phase composition by using x-ray diffraction (Bruker D8 ID 3003TT, Bruker Corporation, USA).

    The tensile adhesive strength between the ceramic and metallic layers was measured using the symmetric structured composite 1s and, furthermore, composites with dense metallic layers. The powder treatment and the composition for the dense sintering metallic tapes is published in [19]. For testing round specimens with a diameter of approximately 25 mm (after sintering) were prepared and sintered at 1270, 1330 and 1400°C. The test was carried out in accordance to DIN EN 582 with a velocity of 0.5 mm/min [28]. It is known from Dourandish et al. [29] that boron can increase tensile strength of composites, due to the formation of a Fe-B eutectic at 1170°C, which enables an improved wetting of the ceramic material. To investigate this 0.03 and 0.63 wt% amorphous boron was added to the ceramic layer. The test was done with a 100 kN universal testing machine (Zwick GmbH, Germany) (figure 2). For this, the metallic layers were glued to the testing stamps with FM®1000 gluing pads (Cytec Industries Corp., USA).

    Figure 2. Universal testing machine (left) and scheme of the experimental build up for testing the adhesive tensile strength.

    3. Results and Discussion


    3.1. CTE and sintering behaviour

    Figure 3 shows the CTE of the used materials up to a temperature of 1200°C. CTEs match well to each other until 700°C, with a maximum mismatch below 1.0×10−6 1/K. With increasing temperature especially the CTE of the metallic material increases, which leads to a mismatch of 2.4×10−6 1/K at 1200°C [18]. However, it has to be assumed that the gap between the CTEs at higher temperatures isn’t as bad as at lower temperatures because the metallic material has a low elasticity/viscosity at higher temperatures. Hence, the metallic layer can absorb some occurring stresses by deformation. Furthermore, other metallic materials and especially the austenitic steels 17-4PH or 316L, which are often used for the combination with zirconia, have much higher CTE’s and with this, higher mismatches [30].

    Figure 3. Coefficient of thermal expansion of the zirconia and the Crofer22APU.

    Figure 4 compares the sintering behaviour of ceramic and the four different metallic powders as a function of temperature and dwell time. The TZ-3YS-E powder starts sintering at 1150°C. Its linear shrinkage at the end of the dwell time amounts 23.5%.

    Figure 4. Sintering behaviour of the ceramic and metallic materials up to 1365°C with a heating rate of 3 K/min and a dwell time of 2 h under argon (95%)/hydrogen (5%) atmosphere.

    Metal powder 1 starts sintering in the same temperature region (1150°C) as the TZ-3YS-E. Hereby, also the course of shrinkage exhibits a good agreement until the beginning of the dwell time. This leads to a good accordance of the sintering strain rates (figure 5) and, furthermore, only a small maximum mismatch in strain rate of±0.25×10−3 1/min (figure 6) until the beginning of the dwell time. The mismatch in strain rate increases up to a maximum value of −0.75×10−3 1/min due to the lower strain rate of metal powder 1 in the first 30 minutes of the dwell time.

    Figure 5. Sintering strain rate of the TZ-3YS-E and the 4 metallic powders up to 1365°C with a heating rate of 3 K/min and a dwell time of 2 h under argon (95%)/hydrogen (5%) atmosphere.
    Figure 6. Differences in sintering strain rate of the 4 metallic powders minus TZ-3YS-E up to 1365°C with a heating rate of 3 K/min and a dwell time of 2 h under argon (95%)/hydrogen (5%) atmosphere.

    In comparison to powder 1, all other metallic powders start sintering at higher temperatures. Furthermore, their strain rates are significantly lower which results in lower total shrinkage values of 16.5% for powder 2, 10% for powder 3 and 7.5% for powder 4, respectively.

    The later onset of sintering, as well as the lower maximum sintering strain rates, result in higher differences in strain rate. Composite 3 has a maximum mismatch at 1280°C of −2×10−3 1/min, while composite 2 and 4 have the highest mismatch values at −3×10−3 1/min (figure 6).


    3.2. Sintering results

    Figure 7 presents FESEM images of the microstructure of pure 1000 µm thick zirconia tapes after sintering for 2 h at 1350 and 1400°C. Analysis exhibited that the porosity at 1350°C amounts 4.2%, while increasing sintering temperature to 1400°C decreases porosity below 1%.

    Figure 7. Electron scanning microscopy pictures of the ceramic microstructure sintered at 1350°C and 1400°C with a heating rate of 3 K/min and a dwell time of 2 h under argon (95%)/hydrogen (5%) atmosphere.

    The asymmetric composites 1a and 2 were sintered with and without a coverage. While the mismatches in strain rate led to a warpage of the uncovered composites (figure 8), the coverage completely avoided a warpage. Thereby the coverage led to no macroscopic cracks on the surface of the composites. Considering figure 8, it is visible that the higher mismatch in strain rate of composite 2 in comparison to composite 1a led to a more significant warpage, despite its higher metallic layer thickness.

    Figure 8. SEM images of asymmetric composites 1 (left) and 2 (right) after sintering for 2 h at 1350°C without a coverage.

    It is reducible to differences in viscosity that the much thinner ceramic layer (green thickness 11.6 µm) can bend the whole composite with a metallic layer thickness between 200 µm (composite 1a) and 400 µm (composite 2). Figure 9 shows the temperature-and density-dependent viscosity of the ceramic layer and metal powder 1 during sintering. The ceramic material starts with a viscosity of 1100 GPa*s at 1150°C. At lower temperatures, the viscosity was too high so that no measurements were possible. With increasing sintering temperature, viscosity decreased to 20 GPa*s. During dwell time, further sintering and densification led to an increase of the viscosity up to 250 GPa*s. The same response exhibits the metallic material, whereby its viscosity is the whole time 1 to 2 magnitudes lower, which is the reason for that the thin ceramic layer can bend the whole composite. The measured viscosity values, thereby match well with already published values from literature (24, 25).

    Figure 9. Viscosity of TZ-3YS-E and metal powder 1 during sintering.

    Figure 10 shows the porosity of the free sintered ceramic layer (TZ-3YS-E) as well as the porosity of the co-sintered ceramic layers of composite 1a and 2 at different temperatures with a dwell time of 2h, covered and uncovered. Comparing porosities of the uncovered composites at 1350°C with the free sintered zirconia layer, it is obvious that the sintering behaviour of the zirconia layer is influenced by the metallic substrate. The porosity of composite 2 amounts 7.5%, which is 3% higher than the porosity of the free sintered zirconia layer, while the porosity of composite 1a is 1.5% lower than the porosity of the free sintered layer. This is caused by compression stresses (composite 1a) and tensile forces (composite 2), respectively, which improve respectively reduce sinter activity of the zirconia layer.

    Figure 10. Porosity after sintering of a single zirconia tape (TZ-3YS-E) and the zirconia layer sintered as composite at different temperatures, a heating rate of 3 K/min and a dwell time of 2 h (left) and sintered at 1350°C (right).

    The negative influence is increased when the composite is covered during sintering. In this case, tensions cannot be prevented or relieved by a warping of the composite. Regarding the covered composite 2, it shows already at 1200°C an 8% higher porosity than the free sintered layer. With increasing sintering temperature the mismatch is slightly increased up to 9%. Comparisons of the covered and uncovered composites show that the porosity is 4% higher when it is covered during sintering. In contrast to these results, the coverage shows no negative influence in the case of composite 1. Here, both ceramic layers have porosities around 1%.

    Further investigations, using symmetric structured composites out of the combination 1s, 3 and 4, confirm these results (figure 11). After sintering at 1350°C, the porosity of the zirconia layer of composite 3 and 4 amount 11%, respectively 15%, while composite 1s has a porosity of 2%. The porosity of composite4 is higher than the porosity of composite 3, due to its lower lateral shrinkage.

    Figure 11. Porosity of zirconia layer of composite 1, 3 and 4 after sintering for 2 h at 1350°C respectively 1400°C with a heating rate of 3 K/min under argon (95%)/hydrogen (5%) atmosphere.

    Nevertheless, an increase of the sintering temperature up to 1400°C minimizes the negative effect of the metallic layers on densification of zirconia layers. At this temperature, the porosity is in all cases below 1%. This is also visible in the scanning microscopy images of the composites in figure 12 on the right hand side.

    Figure 12. Scanning electron microscopy image of composite 1, 3 and 4 sintered for 2 h at 1400°C, a heating rate of 3 K/min under Argo (95%)/Hydrogen (5%) atmosphere.

    The lateral shrinkage of the ceramic layers in the case of composite 3 and 4 is constraint due to the metallic substrate. Accordingly, densification of the ceramic layer at 1400°C is only possible with an increased shrinkage in thickness direction. This is evidenced by measurements of the thickness before and after sintering (table 2). In green state the ceramic layer thickness amounts in all cases 11.6 µm. At 1400°C the ceramic layer thickness is between 9.5 and 8.2 µm, depending on the composite composition. Composite 1s has the biggest layer thickness after sintering, because the metallic substrate has nearly the same lateral shrinkage as the pure zirconia layer. This leads to a calculated shrinkage in thickness-direction of 18.1%. In the case of composite 3 and 4 the metallic layers have lower lateral shrinkage values, which leads to an increased shrinkage in z-direction of 27% and 29% respectively.

    Table 2.Ceramic layer thickness of symmetric composites 1, 3 and 4 before and after co-sintering at 1400°C.
    compositegreen1400 °C
    layer thickness/µmlayer thickness/µmshrinkage/%
    19.5 ± 1.418.1
    311.6 ± 1.08.5 ± 1.426.7
    48.2 ± 1.229.3
     | Show Table
    DownLoad: CSV

    The results show that it is possible to densify the zirconia layer, even in combination with a high mismatch in strain rate. Despite these results, an adaptation of the sintering behaviour is indispensable. This is visible in the FESEM-images with lower magnifications (left hand side and middle in figure 12). While composite 1s has a completely defect free zirconia layer which is necessary for the different applications, sintering of both other composites is possible but they exhibit (micro-) cracks in the ceramic layer. Hence, a good adaption of sintering behaviour is necessary to achieve defect free composites.


    3.3. X-ray diffraction measurements

    Figure 13 shows x-ray diffraction measurements of the zirconia layer of composite 1a before and after sintering. As the measurements show, before sintering the powder has some monoclinic phases which are non-existent after sintering. This result is independent of the investigated composite (composite 1a or 2) and demonstrates that a mismatch in CTE above 700°C, does not lead to a stress induced transformation of the tetragonal in the monoclinic phase.

    Figure 13. XRD-analysis of the zirconia layer of the asymmetric composite 1a after sintering at 1350°C.

    3.4. Interphase and adhesive tensile strength

    The formation of mixed oxide phases [18] as well as a mechanical grouting is the bases of bonding between both materials.

    Figure 14 left shows a metal-ceramic-metal composite consisting of a dense metallic layer (porosity < 3%) in combination with the zirconia layer. As the higher magnification at the right side of the figures shows, there is a 0.5 µm thick mixed oxide phase in the interface of both materials. The EDX measurement of this phase shows, that in addition to the main elements of both materials there exist an accumulation of different alloying elements as Al, La, Si. These elements can form phases as La2O3, AlLaO3 or for example ZrSiO4.

    Figure 14. Scanning electron microscopy image of a metal-ceramic-metal compound with dense metallic layers (left) and the EDX-measurement of the interphase phase (right).

    Measurements of the adhesive tensile strength of composite 1s (porous metallic layer) exhibit a high dependency on sintering temperature (figure 15). Starting with a tensile strength between 1 and 2 MPa at 1270°C it is increased up to 4-5 MPa at 1330°C. Thereby, also a slight increase of the strength with increased boron content is visible. An increase of the sintering temperature up to 1400°C doubled tensile strength to 8 MPa without boron. Hereby, the strongest influence of the boron content is visible with an increase up to 9.6 MPa with 0.03% boron and up to 10.8 MPa with 0.63%. Regarding the use of boron, it has to be assumed that boron can lead to the formation of glass phases in the zirconia material which can negatively influence its properties. So, according to the specific application, the use of boron has to be reconsidered.

    Figure 15. Tensile strength of the symmetric composite 1 in terms of sintering temperature and boron content.

    The experiments exhibited that the composite breaks preferably at one metal ceramic interface, while the ceramic layer stays intact. This is visible in figure 16, which shows a top view on the ceramic layer after tensile tests. Some metallic particles are on top of the ceramic layer, which were pulled out of the metallic layer.

    Figure 16. FESEM picture of the surface of ceramic layer of the symmetric composite 1 after tensile strength test.

    Against the expectation tests of composites with dense metallic layers, sintered for 2 h at 1400°C, showed, with a maximum tensile strength up to 3 MPa (without boron), much lower strengths. Measurements at lower sintering temperatures were not feasible, because they delaminated during clamping in the testing machine. Actually, a higher tensile strength was expected, due to the larger contact area between both materials.

    A possible reason for this result can be assumed comparing the interface of the composites with dense metallic layers and porous metallic layers. Both show mixed oxide phase, while the only visible difference is, that the composites with the dense metallic layer have an interface with a lower roughness and less micro undercuts. Hence, the authors assume that especially the micro undercuts give strength and not the phases. So regarding design, undercuts or a materials gradient should be included in the composite. However, these results shall be proven in further investigations.


    4. Conclusion

    The presented work shows co-sintering results of metal-ceramic symmetric and asymmetric composites. The investigations exhibited that a mismatch in sintering behaviour due to a less shrinking metallic layer, led to a reduced density of the zirconia layer in comparison to a free sintered zirconia layer after sintering. Furthermore, mismatches in the strain rate cause cracks in the ceramic layer.

    On the other hand, adapting of sintering behaviour allows manufacturing of defect-free metal-ceramic composites with dense ceramic layers. Largest dimensions of manufactured composites were 10 cm×10 cm. X-ray diffraction measurements exhibited that there was no stress induced transformation from the tetragonal to the monoclinic phase due to mismatches in the CTE.

    Adhesive tensile strength measurements showed a strong dependence of strength in terms of sintering temperature and boron content while maximum measured strength was 10.8 MPa.


    Acknowledgement

    This project was funded by the German Ministry of Economic Affairs and Energy (BMWi) through the German Federation of Industrial Research Associations (“Arbeitsgemeinschaft industrieller Forschungsvereinigung-AiF”) under the IGF-project number 18520 BR.


    Conflict of Interest

    The authors declare no conflicts ofinterest regarding this paper.




    [1] E. A. B. Abdel-Salam, E. A. Yousif, M. A. El-Aasser, Analytical solution of the space-time fractional nonlinear schrödinger equation, Rep. Math. Phys., 77 (2016), 19-34. doi: 10.1016/S0034-4877(16)30002-7
    [2] Z. Asgari, S. Hosseini, Efficient numerical schemes for the solution of generalized time fractional burgers type equations, Numer. Algorithms, 77 (2018), 763-792. doi: 10.1007/s11075-017-0339-4
    [3] A. Bhrawy, M. Abdelkawy, A fully spectral collocation approximation for multi-dimensional fractional schrödinger equations, J. Comput. Phys., 294 (2015), 462-483. doi: 10.1016/j.jcp.2015.03.063
    [4] A. Bhrawy, M. A. Zaky, Highly accurate numerical schemes for multi-dimensional space variable-order fractional schrödinger equations, Comput. Math. Appl., 73 (2017), 1100-1117. doi: 10.1016/j.camwa.2016.11.019
    [5] J. P. Boyd. Chebyshev and Fourier Spectral Methods, Courier Corporation, 2001.
    [6] A. Chechkin, R. Gorenflo, I. Sokolov, Retarding subdiffusion and accelerating superdiffusion governed by distributed-order fractional diffusion equations, Phys. Rev. E, 66 (2002), 046129.
    [7] J. B. Chen, M. Z. Qin, Y. F. Tang, Symplectic and multi-symplectic methods for the nonlinear schrödinger equation, Comput. Math. Appl., 43 (2002), 1095-1106. doi: 10.1016/S0898-1221(02)80015-3
    [8] X. Cheng, J. Duan, D. Li, A novel compact adi scheme for two-dimensional riesz space fractional nonlinear reaction-diffusion equations, Appl. Math. Comput., 346 (2019), 452-464.
    [9] M. Dehghan, A. Taleei, A compact split-step finite difference method for solving the nonlinear schrödinger equations with constant and variable coefficients, Comput. Phys. Commun., 181 (2010), 43-51. doi: 10.1016/j.cpc.2009.08.015
    [10] K. Diethelm, The Analysis of Fractional Differential Equations: An Application-oriented Exposition Using Differential Operators of Caputo Type, Springer Science and Business Media, 2010.
    [11] J. Dong, Scattering problems in the fractional quantum mechanics governed by the 2d spacefractional schrödinger equation, J. Math. Phys., 55 (2014), 032102.
    [12] J. Dong, M. Xu, Some solutions to the space fractional schrödinger equation using momentum representation method, J. Math. Phys., 48 (2007), 072105.
    [13] W. Fan, F. Liu, A numerical method for solving the two-dimensional distributed order spacefractional diffusion equation on an irregular convex domain, Appl. Math. Lett., 77 (2018), 114-121. doi: 10.1016/j.aml.2017.10.005
    [14] W. Fan, H. Qi, An efficient finite element method for the two-dimensional nonlinear time-space fractional schrödinger equation on an irregular convex domain, Appl. Math. Lett., 86 (2018), 103-110. doi: 10.1016/j.aml.2018.06.028
    [15] B. Jin, B. Li, Z. Zhou, Subdiffusion with a time-dependent coefficient: Analysis and numerical solution, Math. Comput., 88 (2019), 2157-2186. doi: 10.1090/mcom/3413
    [16] G. Jumarie, An approach to differential geometry of fractional order via modified riemannliouville derivative, Acta Math. Sin., 28 (2012), 1741-1768. doi: 10.1007/s10114-012-0507-3
    [17] D. Li, J. Wang, J. Zhang, Unconditionally convergent l1-galerkin fems for nonlinear timefractional schrodinger equations, SIAM J. Sci. Comput., 39 (2017), A3067-A3088.
    [18] D. Li, C. Wu, Z. Zhang, Linearized galerkin fems for nonlinear time fractional parabolic problems with non-smooth solutions in time direction, J. Sci Comput., 80 (2019), 403-419. doi: 10.1007/s10915-019-00943-0
    [19] L. Li, D. Li, Exact solutions and numerical study of time fractional burgers' equations, Appl. Math. Lett., 100 (2020), 106011.
    [20] M. Li, A high-order split-step finite difference method for the system of the space fractional cnls, Eur. Phys. J. Plus, 134 (2019), 244.
    [21] M. Li, X. M. Gu, C. Huang, et al. A fast linearized conservative finite element method for the strongly coupled nonlinear fractional schrödinger equations, J. Comput. Phys., 358 (2018), 256-282. doi: 10.1016/j.jcp.2017.12.044
    [22] M. Li, C. Huang, W. Ming, A relaxation-type galerkin fem for nonlinear fractional schrödinger equations. Numer. Algorithms, 83 (2019), 99-124.
    [23] M. Li, C. Huang, P. Wang, Galerkin finite element method for nonlinear fractional schrödinger equations, Numer. Algorithms, 74 (2017), 499-525. doi: 10.1007/s11075-016-0160-5
    [24] M. Li, C. Huang, Z. Zhang, Unconditional error analysis of galerkin fems for nonlinear fractional schrödinger equation, Appl. Anal., 97 (2018), 295-315. doi: 10.1080/00036811.2016.1262947
    [25] M. Li, C. Huang, Y. Zhao, Fast conservative numerical algorithm for the coupled fractional kleingordon-schrödinger equation, Numer. Algorithms, 82 (2019), 1-39. doi: 10.1007/s11075-018-0591-2
    [26] M. Li, D. Shi, J. Wang, et al. Unconditional superconvergence analysis of the conservative linearized galerkin fems for nonlinear klein-gordon-schrödinger equation, Appl. Numer. Math., 142 (2019), 47-63. doi: 10.1016/j.apnum.2019.02.004
    [27] Y. Lin, C. Xu, Finite difference/spectral approximations for the time-fractional diffusion equation, J. Comput. Phys., 225 (2007), 1533-1552. doi: 10.1016/j.jcp.2007.02.001
    [28] Y. F. Luchko, M. Rivero, J. J. Trujillo, et al. Fractional models, non-locality, and complex systems, Comput. Math. Appl., 59 (2010), 1048-1056. doi: 10.1016/j.camwa.2009.05.018
    [29] J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci., 16 (2011), 1140-1153. doi: 10.1016/j.cnsns.2010.05.027
    [30] Z. Mao, G. E. Karniadakis, Fractional burgers equation with nonlinear non-locality: Spectral vanishing viscosity and local discontinuous galerkin methods, J. Comput. Phys., 336 (2017), 143-163. doi: 10.1016/j.jcp.2017.01.048
    [31] A. K. Mittal, A stable time-space Jacobi pseudospectral method for two-dimensional sine-Gordon equation, J. Appl. Math. Comput., 63 (2020), 1-26.
    [32] A. K. Mittal, L. K. Balyan, A highly accurate time-space pseudospectral approximation and stability analysis of two dimensional brusselator model for chemical systems, Int. J. Appl. Comput. Math., 5 (2019), 140.
    [33] A. Mohebbi, Analysis of a numerical method for the solution of time fractional burgers equation, B. Iran. Math. Soc., 44 (2018), 457-480. doi: 10.1007/s41980-018-0031-z
    [34] A. Mohebbi, M. Dehghan, The use of compact boundary value method for the solution of twodimensional schrödinger equation, J. Comput. Appl. Math., 225 (2009), 124-134. doi: 10.1016/j.cam.2008.07.008
    [35] K. Oldham, J. Spanier, The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order, Elsevier, 1974.
    [36] B. Ross, The development of fractional calculus, Hist. Math., 4 (1977), 75-89. doi: 10.1016/0315-0860(77)90039-8
    [37] H. Rudolf, Applications of Fractional Calculus in Physics, World Scientific, 2000.
    [38] L. N. Trefethen, Finite Difference and Spectral Methods for Ordinary and Partial Differential Equations, 1996.
    [39] L. Wei, Y. He, X. Zhang, et al. Analysis of an implicit fully discrete local discontinuous galerkin method for the time-fractional schrödinger equation, Finite Elem. Anal. Des., 59 (2012), 28-34. doi: 10.1016/j.finel.2012.03.008
    [40] E. Yousif, E. B. Abdel-Salam, M. El-Aasser, On the solution of the space-time fractional cubic nonlinear schrödinger equation, Results Phys., 8 (2018),702-708. doi: 10.1016/j.rinp.2017.12.065
    [41] F. Zeng, C. Li, F. Liu, et al. The use of finite difference/element approaches for solving the timefractional subdiffusion equation, SIAM J. Sci. Comput., 35 (2013), A2976-A3000.
    [42] G. Zhang, C. Huang, M. Li, A mass-energy preserving galerkin fem for the coupled nonlinear fractional schrödinger equations, Eur. Phys. J. Plus, 133 (2018), 155.
    [43] H. Zhang, X. Jiang, C. Wang, et al. Galerkin-legendre spectral schemes for nonlinear space fractional schrödinger equation, Numer. Algorithms, 79 (2018), 337-356. doi: 10.1007/s11075-017-0439-1
    [44] X. Zhao, Z. Z. Sun, Z. P. Hao, A fourth-order compact adi scheme for two-dimensional nonlinear space fractional schrodinger equation, SIAM J. Sci. Comput., 36 (2014), A2865-A2886.
  • This article has been cited by:

    1. Avinash K. Mittal, Lokendra K. Balyan, Numerical solutions and stability analysis for solitary waves of complex modified Korteweg–de Vries equation using Chebyshev pseudospectral methods, 2020, 36, 0749-159X, 1662, 10.1002/num.22497
    2. Avinash K. Mittal, Parnika Shrivastava, Manoj K. Panda, Time–space Jacobi pseudospectral simulation of multidimensional Schrödinger equation, 2021, 37, 0749-159X, 1725, 10.1002/num.22605
    3. Avinash K. Mittal, Lokendra K. Balyan, Manoj K. Panda, Parnika Shrivastava, Harvindra Singh, Pseudospectral analysis and approximation of two‐dimensional fractional cable equation, 2021, 0170-4214, 10.1002/mma.7362
    4. M. T. Darvishi, Mohammad Najafi, Lanre Akinyemi, Hadi Rezazadeh, Gaussons of some new nonlinear logarithmic equations, 2023, 32, 0218-8635, 10.1142/S0218863523500133
    5. A. K. Mittal, A spectrally accurate time–space pseudospectral method for reaction–diffusion Malaria infection model, 2022, 41, 2238-3603, 10.1007/s40314-022-02094-9
    6. Waleed Mohamed Abd-Elhameed, Amr Kamel Amin, New Formulas and Connections Involving Euler Polynomials, 2022, 11, 2075-1680, 743, 10.3390/axioms11120743
    7. A. K. Mittal, L. K. Balyan, K. K. Sharma, A spectrally accurate time ‐ space pseudospectral method for viscous Burgers' equation, 2023, 0749-159X, 10.1002/num.23011
    8. A. K. Mittal, Pseudospectral quadrature simulation technique for spatio-temporally parabolic multi-dimensional nonlinear fractional evolution equation, 2024, 32, 0971-3611, 2647, 10.1007/s41478-024-00750-3
    9. A. K. Mittal, Two-dimensional Jacobi pseudospectral quadrature solutions of two-dimensional fractional Volterra integral equations, 2023, 60, 0008-0624, 10.1007/s10092-023-00545-1
    10. A.K. Mittal, A space-time pseudospectral method for solving multi-dimensional quasi-linear parabolic partial differential (Burgers') equations, 2024, 195, 01689274, 39, 10.1016/j.apnum.2023.09.005
    11. Fahir Talay Akyildiz, Kuppalapalle Vajravelu, Cemil Tunç, John Abraham, Chebyshev delta shaped and Chebyshev pseudo-spectral methods for solutions of differential equations, 2025, 236, 03784754, 52, 10.1016/j.matcom.2025.03.034
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4363) PDF downloads(509) Cited by(11)

Article outline

Figures and Tables

Figures(3)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog