Research article Special Issues

The instability of liquid films with temperature-dependent properties flowing down a heated incline

  • We investigate the stability of free-surface flow on a heated incline. We develop a complete mathematical model for the flow which captures the Marangoni effect and also accounts for changes in the properties of the fluid with temperature. We apply a linear stability analysis to determine the stability of the steady and uniform flow. The associated eigenvalue problem is solved numerically by means of a spectral collocation method. Comparisons with nonlinear simulations are also made and the agreement is good.

    Citation: Jean-Paul Pascal, Serge D’Alessio, Syeda Rubaida Zafar. The instability of liquid films with temperature-dependent properties flowing down a heated incline[J]. AIMS Mathematics, 2019, 4(6): 1700-1720. doi: 10.3934/math.2019.6.1700

    Related Papers:

    [1] Zeid Al Qaisi, Qais Alsafasfeh, Ahmad Harb . Stability impact of integrated small scale hybrid (PV/Wind) system with electric distribution network. AIMS Energy, 2018, 6(5): 832-845. doi: 10.3934/energy.2018.5.832
    [2] K. M. S. Y. Konara, M. L. Kolhe, Arvind Sharma . Power dispatching techniques as a finite state machine for a standalone photovoltaic system with a hybrid energy storage. AIMS Energy, 2020, 8(2): 214-230. doi: 10.3934/energy.2020.2.214
    [3] Amer Braik, Asaad Makhalfih, Ag Sufiyan Abd Hamid, Kamaruzzaman Sopian, Adnan Ibrahim . Impact of photovoltaic grid-tied systems on national grid power factor in Palestine. AIMS Energy, 2022, 10(2): 236-253. doi: 10.3934/energy.2022013
    [4] Katherin Indriawati, Bambang L. Widjiantoro, Ali Musyafa . Design of sensor and actuator fault tolerant control system on wind turbine benchmark for Region II. AIMS Energy, 2019, 7(2): 111-126. doi: 10.3934/energy.2019.2.111
    [5] Mohamed G Moh Almihat . An overview of AC and DC microgrid energy management systems. AIMS Energy, 2023, 11(6): 1031-1069. doi: 10.3934/energy.2023049
    [6] Rashid Al Badwawi, Mohammad Abusara, Tapas Mallick . Speed control of synchronous machine by changing duty cycle of DC/DC buck converter. AIMS Energy, 2015, 3(4): 728-739. doi: 10.3934/energy.2015.4.728
    [7] Carla Viveiros, Rui Melicio, Victor Mendes, Jose Igreja . Adaptive and predictive controllers applied to onshore wind energy conversion system. AIMS Energy, 2018, 6(4): 615-631. doi: 10.3934/energy.2018.4.615
    [8] Hari Charan Nannam, Atanu Banerjee . A novel control technique for a single-phase grid-tied inverter to extract peak power from PV-Based home energy systems. AIMS Energy, 2021, 9(3): 414-445. doi: 10.3934/energy.2021021
    [9] Nagaraj C, K Manjunatha Sharma . Fuzzy PI controller for bidirectional power flow applications with harmonic current mitigation under unbalanced scenario. AIMS Energy, 2018, 6(5): 695-709. doi: 10.3934/energy.2018.5.695
    [10] Arqum Shahid, Ahsan Ali, Faheem Ashiq, Umer Amir Khan, Muhammad Ilyas Menhas, Sajjad Manzoor . Comparison of resistive and inductive superconductor fault current limiters in AC and DC micro-grids. AIMS Energy, 2020, 8(6): 1199-1211. doi: 10.3934/energy.2020.6.1199
  • We investigate the stability of free-surface flow on a heated incline. We develop a complete mathematical model for the flow which captures the Marangoni effect and also accounts for changes in the properties of the fluid with temperature. We apply a linear stability analysis to determine the stability of the steady and uniform flow. The associated eigenvalue problem is solved numerically by means of a spectral collocation method. Comparisons with nonlinear simulations are also made and the agreement is good.



    Due to fossil fuel environmental impacts, renewable energy sources (RESs) have been increasingly used across the globe [1]. Solar energy, especially photovoltaic (PV) energy, is considered one of the most attractive renewable energy technologies due to its abundant energy source and availability [2,3]. According to the International Energy Agency Renewable 2021 report, an annual addition of 305 GW per year of renewable energy capacity is expected to be utilized between 2021 and 2026 [4]. The solar PV annual market is expected to grow by 17% to approximately 160 GW in 2021, reaching 200 GW in 2026. Even though renewable energy sources provide better alternatives to conventional power plants, the variation of renewable energy outputs poses major technical issues in integrating these sources into the power grid [5,6]. Moreover, the output nature of these sources, such as DC, AC, and frequency level, makes these sources unsuitable to be interconnected with conventional power generation plants. Distributed generator systems (DGSs) enabled the idea of integrating RESs into the existing power grids. In general, a DGS usually consists of a power generator such as a solar PV or wind generator connected to the utility grid via a power electronic device. Depending on the nature of the systems' output (i.e., DC or AC), the power electronic devices that are utilized for the integration are chosen. These devices provide the opportunity to perform different control tasks, such as flexible power flow, voltage and frequency control, or maximum power point tracking (MPPT) [7,8].

    Different topologies have been used for integrating solar PV systems into existing power grids [10,11,12,13]. One of the most commonly used topologies is the two-stages transformerless [9]. Such topology consists of a DC-DC converter and a DC-AC three-phase two-level inverter. The converter extracts the PV maximum power and boosts the output voltage of the PV generator, while the inverter transfers the DC voltage to AC voltage that is compatible with the power grid voltage level. Other topologies that reduce the leakage current of a single-phase grid connected PV system were introduced in [14,15,16]. These works implemented different power electronics devices such as H5, H6, and HERIC inverters, which show a remarkable reduction of PV solar system leakage current. An H7 three-phase inverter for integrating PV with leakage current reduction was proposed in [17]. Besides leakage current, many issues (e.g., nonlinear load and grid fault) can affect the power grid stability. In the case of grid fault, it requires the system to be capable of restoring the system operation to the normal conditions.

    Various control techniques have been implemented to control the PV systems' outputs such as voltage and power. The classical proportional integral (PI) method has been widely used for integrating PV systems. A fuzzy logic MPPT combined with a PI controller for charging batteries was presented in [18]. A small-signal modeling and designing of a multi-variable PI controller for a voltage source inverter in a high-voltage DC system was introduced in [19]. In [20], an adaptive PI control method to maintain the DC link capacitor voltage balance in a solar-PV-fed grid-connected neutral point clamped inverter was proposed. PI control combined with artificial neural networks for controlling a PV system supplying a direct current water pump was presented in [21]. The authors in [22] proposed an optimal PI tuning method using a genetic algorithm for controlling a single-phase grid-connected PV system. The research work in [23] introduced a PI controller-based active filter for a grid-connected PV system's current harmonics compensation. The previous works show that PI controllers have a fast dynamic response as well as satisfactory steady-state behavior. However, the PI controller requires external signals modulation generation and excessive time and effort for tuning parameters. Moreover, severe faults may deteriorate the system's stability.

    Sliding mode controllers (SMC) offer robustness to uncertainties, upper bounded disturbances and parameter variations, and the derivation of such controllers relies on the availability of system models [24]. These controllers have been frequently reported in renewable energy applications, such as super twisting SMC for grid connected PV systems [25], non integer SMC for DC micro-grid stability problems [26], adaptive non integer SMC for wind energy [27], industrial grade robust controllers [28] and fractional order control for wind energy [29]. Although SMC offer robust performance, however, such controllers compromise robustness when practically implemented using digital processors.

    Another control technique method that has been implemented for controlling grid-connected PV systems is model predictive control (MPC). Compared to other types of controllers, MPC provides some advantages, including ease of implementation, the ability to control different variables at the same time, and the ability to add limitations in the control objective. All these features and others make MPC an attractive control method not only for grid-connected PV systems but also for different applications in power systems. In [30], an MPPT and leakage current reduction technique using MPC for the H5 inverter of the PV system was proposed. Voltage regulation and an MPPT method based on MPC for a two-stage grid-connected PV system were introduced in [31]. The work in [32] proposed a current sensor-less model predictive control of an MPPT algorithm. A new MPPT technique of a PV module based on MPC was presented in [33]. The authors in [34] proposed an MPC strategy and a high-efficiency two-stage topology for a single-phase PV system.

    This paper investigates the performance of the current model predictive control (CMPC) for a two-stage three-phase PV system under grid fault conditions. To verify the CMPC capability, the results of the CMPC under grid fault were compared to the PI controller. Moreover, the results were compared with our previously proposed sliding mode control (SMC) [35]. The rest of the paper is organized as follows: Section 2 describes the grid connected PV system, section 3 derives the converter controller, and section 4 discusses the simulation results. Finally, the conclusion is drawn based on the presented data and results.

    The PV system considered in this work is shown in Figure 1. A two-stage transformerless topology is implemented to integrate the PV system into the main grid. The system consists of a PV panel connected to a DC-DC converter and then is connected to a DC-AC three-phase two-level inverter which delivers the power to the grid through an RL filter. The DC-DC converter is controlled via an MPPT algorithm to ensure maximum output power during the PV panel operation. Then, the DC link voltage is controlled via a PI controller in the d-q reference frames to generate the reference currents for the CMPC controller. After transferring the current reference values from the d-q frames to α-β, the reference values are then sent to the CMPC, which controls the output current of the three-phase two-level inverter. Figure 1 shows the schematic diagram of the PV system control process implemented in this work.

    Figure 1.  PV system topology and control strategy.

    Maximum power point tracking (MPPT) controllers are widely utilized for harvesting peak power from stochastic renewable energy sources using high-frequency power converters interface. In the existing literature, several MPPT methods have been discussed, such as perturb and observed, incremental conductance, ripple correlation, and fuzzy logic based maximum power harvesting methods. However, the incremental conductance algorithm is the most widely utilized technique and is easy to implement. In this work, the incremental conductance algorithm is used to extract maximum power from the energy source. The detailed diagram of the incremental conductance algorithm is shown in Figure 2.

    Figure 2.  Flow chart of incremental conductance algorithm.

    In this subsection, the DC link voltage model is realized from the system block diagram of Figure 1, and then a proportional integral (PI) is utilized to stabilize the DC link voltage to a reference value. As given in Figure 1, a DC link capacitor stores input power from the PV panel in the form of energy and then delivers the stored energy to the inverter. So, the dynamics across the DC link capacitor are expressed as follows:

    ˙Vc=1CVc(Pin32Vgdifd) (3.1)

    where Vc represents DC link voltage, Pin is the PV power stored in the DC link capacitors, Vgd represents grid voltage, and ifd is the grid current.

    In order to control the voltage dynamics of Eq 3.1, a PI controller is utilized and expressed as follows:

    idref=kpeVc+kieVc (3.2)

    where eVc is the voltage error, and kp and ki are the PI proportional and integral gains, respectively.

    The model predictive control strategy for controlling power electronic devices is based on the fact that there is only a finite number of switching states related to the voltage vectors of a specific power electronic device [36]. Through an optimization process, the system model is used to predict the future behavior of the control variables for each voltage vector. Then, the optimal voltage vector that results in the least error between the reference and controlled value is selected. The optimal voltage vector selection is done via a cost function that evaluates all the possible switching states and selects the optimal solution. To implement the MPC strategy, the following steps must be identified:

    ● Continuous-time system model,

    ● Discrete-time system model,

    ● Cost function to select the optimal voltage vector,

    ● Generation of switching signals related to the selected optimal voltage vector.

    The block diagram in Figure 3 illustrates the MPC process. The current values that have been measured x(k) are used in the system model to forecast future values x(k+1). Then, by minimizing the error between the future and reference values, a cost function assesses all future values. Finally, the switching signals are generated by the MPC controller.

    Figure 3.  Model predictive control process.

    The three-phase inverter topology is shown in Figure 4. It is made up of three legs, each having two switches. Since the upper switches are complementary with the lower ones, each leg results in two control actions, which means that the three-phase two-level inverter produces eight different switching state combinations [37]. The three-phase two-level inverter's voltage state-space vector can be written as

    Vx=23Vc(S1+aS2+a2S3) (3.3)

    where Vx is the output voltage of the inverter, and a equals ej2π3. The state-space vectors of the three-phase two-level inverter are illustrated in Table 1.

    Figure 4.  Three-phase two-level inverter.
    Table 1.  Three-phase two-level inverter voltage vectors.
    Mode S1 S2 S3 Vc
    1 0 0 0 0
    2 0 0 1 23Vdc0
    3 0 1 1 23Vdc60
    4 0 1 0 23Vdc120
    5 1 1 0 23Vdc180
    6 1 0 0 23Vdc240
    7 1 0 1 23Vdc300
    8 1 1 1 0

     | Show Table
    DownLoad: CSV

    The mathematical model of the system AC side in Figure 1 can be derived as

    Vx=VR+VL+Vg (3.4)

    where VL is the voltage across the filter inductor, VR is the voltage across the equivalent series resistance (ESR) of the inductor, and Vg is the main grid voltage. The current load dynamics can be expressed as

    Vx=Rif+Ldifdt+Vg (3.5)

    where the current flow from the inverter to the main grid is denoted by if. Eq 3.5 can be rearranged in the following way:

    Ldifdt=Vx(Rif+Vg). (3.6)

    To predict the future values of inverter output current, the system discrete-time must be optioned [38]. Using the approximation derivative equation

    dxdtx(k+1)x(k)Ts (3.7)

    where x(k+1) is the future value, x(k) is the present value, and Ts is the sampling time, the discrete-time of Eq 3.6 is

    Lif(k+1)if(k)Ts=Vx(k)(Rif(k)+Vg(k)). (3.8)

    After rearrangement of Eq 3.8, the future value of the inverter output current is predicted as

    if(k+1)=if(k)+TsL(Vx(k)(Rif(k)+Vg(k))). (3.9)

    Since there are eight voltage vectors of the three-phase two-level inverter Vs (Table 1), the CMPC predicts eight values of the current future value in each sampling time Ts.

    The square error cost function (3.10) is used to evaluate each predicted value and select the optimal solution. In other words, the cost function minimizes the error for each future value of the eight predicted values in Eq 3.9. After that, the best voltage vector is chosen, and the CMPC generates the best switching states to send to the inverter. Note that the CMPC does this procedure every Ts sampling period. Figure 5 illustrates the CMPC strategy algorithm.

    g=(i(k+1)i(k+1))2 (3.10)
    Figure 5.  CMPC algorithm.

    In order to investigate the performance of the CMPC under grid fault, the system in Figure 1 was simulated using a MATLAB/Simulink environment. A comparative analysis between CMPC and PI control was conducted. For further investigation, the results are also compared to our previously proposed SMC. The details of implementing the SMC strategy for the PV two-stage transformerless topology were presented in [29]. The parameters of the system in Figure 1, MPC and PI controller are shown in Table 2. Figure 6 illustrates the input sun irradiance waveform to the PV power generator. A short circuit fault between grids' phase A and ground is applied at t = 1 s. The fault duration is set to one cycle, which approximately equals 40 ms. The fault resistance is set to 1.31 Ω, while the ground resistance is 0.01 Ω.

    Table 2.  System parameters.
    Parameter Symbol Value
    PV Rated Power VPV 5 kW
    DC-link Reference Voltage VDC 800 V
    Filter Inductance L 2.8 mH
    ESR of L R 0.02 Ω
    Utility Grid Voltage Vg 380 VLL,RMS
    Grid Frequency fgrid 50 Hz
    Sampling Time Ts 10 e6
    Proportional Gain Kp 1.6
    Integral Gain Ki 20

     | Show Table
    DownLoad: CSV
    Figure 6.  Irradiance waveform.

    Figure 7 shows the MPPT reference value and the actual PV output power using CMPC, PI, and SMC controllers. As can be seen, using CMPC, the PV output power follows the MPPT reference value smoothly even under grid fault. The CMPC does not produce any overshoot or undershooting of PV output power signals. Even though the PV output powers using PI and SMC controllers follow the reference value, they both suffer from undershooting, which almost equals 15 kW in both cases of PI and SMC controlling strategies. Therefore, Figure 7 proves that CMPC has a satisfactory performance in stabilizing PV output power under grid fault compared to PI and SMC controllers.

    Figure 7.  MPPT reference value and PV output powers.

    Similar to the PV output power, the PV output current tracks the radiation value smoothly using the three controllers. Figure 8 shows the PV output current of the CMPC, PI, and SMC. Under the grid fault, CMPC maintains the PV output current value at the same level without any oscillation, overshooting, or undershooting. In contrast, a jump of 45 A is introduced using PI and SMC during the grid fault.

    Figure 8.  PV output currents.

    The DC link voltage is shown in Figure 9. As is illustrated in the figure, the recovery time for the DC link voltage is approximately 150 ms using all three controllers. However, CMPC results in less overshoot and undershoot of the DC link voltage under grid fault. A zoomed version of Figure 9 is presented in Figure 10. As can be seen, CMPC introduces a maximum of 75 V undershooting, while PI and SMC result in 180 V and 185 V, respectively. Moreover, CMPC produces less overshooting, with 20 V; however, PI and SMC result in almost 65 V overshooting. As a result, Figures 9 and 10 show the capability of the CMPC in stabilizing the DC link voltage compared to both PI and SMC during grid fault.

    Figure 9.  Zoomed version of Figure 8.
    Figure 10.  DC link voltages.

    Figures 11 and 12 present the output active power of the inverter. It is evident that the three controllers are capable of tracking the reference values smoothly. Even though CMPC has a bit more undershoot than the PI and SMC controllers, its recovery time is considered the lowest, at 50 ms. However, the recovery times of SMC and PI controllers are 60 ms and 80 ms, respectively. The inverter's reactive powers are shown in Figures 13 and 14. Since the q axis reference current is set to zero, the three-phase two-level inverter's output reactive power should be zero. It is evident from the figures that both CMPC and SMC are capable of maintaining the reactive power at zero even under the grid fault, while the PI controllers introduce some noticeable oscillations in the reactive power. That is because of the parameters decoupling issues that the PI controller imposes.

    Figure 11.  Active powers of the three-phase two-level inverter.
    Figure 12.  Zoomed version of Figure 10.
    Figure 13.  Reactive powers of the three-phase two-level inverter.
    Figure 14.  Zoomed version of Figure 12.

    The three-phase two-level inverters' output currents using the CMPC, PI, and SMC are shown in Figure 15. The comparative analysis shows that the recovery time of the output current under grid fault is 50 ms, while the recovery times for PI and SMC are 80 ms and 60 ms, respectively. The total harmonic distortion of the three-phase two-level inverter output current using the three control strategies is shown in Figure 16. As can be seen, CMPC results in a lower total harmonic distortion percentage and also a lower settling time compared to PI and SMC.

    Figure 15.  Three-phase two-level output currents.
    Figure 16.  Total harmonic distortion.

    This research work analyzes and investigates the CMPC strategy performance for a PV two-stages transformerless topology under grid fault conditions. A comparative analysis between CMPC, PI, and SMC is conducted under grid fault. The results illustrate that CMPC has a satisfactory performance compared to PI and SMC under faulty conditions. First, CMPC does not result in any overshoot or undershoot for the PV output current and power. In contrast, PI and SMC result in undershoot, with almost 15 kW and 45 A for the PV output power and current, respectively. Second, the settling time of the three-phase inverter's output active power and current is 50 ms using CMPC. However, the settling times of the PI and SMC are 80 ms and 60 ms, respectively. Finally, CMPC results in less voltage dip (75 V) of the DC link voltage compared to a voltage dip of 180 V for both PI and SMC.

    This work was supported by the Research Groups Program funded by the Deanship of Scientific Research, Taif University, Ministry of Education, Saudi Arabia, under Grant 1-441-106.

    The authors declare no conflict of interest.



    [1] P. L. Kapitza, Wave flow of thin layers of a viscous fluid: II. Fluid flow in the presence of continuous gas flow and heat transfer, Zh. Eksp. Teor. Fiz., 18 (1948), 19-28.
    [2] P. L. Kapitza and S. P. Kapitza, Wave flow of thin layers of a viscous fluid: III. Experimental study of undulatory flow conditions, Zh. Eksp. Teor. Fiz., 19 (1949), 105-120.
    [3] T. B. Benjamin, Wave formation in laminar flow down an inclined plane, J. Fluid Mech., 2 (1957), 554-574. doi: 10.1017/S0022112057000373
    [4] C.-S. Yih, Stability of liquid flow down an inclined plane, Phys. Fluids, 6 (1963), 321-334. doi: 10.1063/1.1706737
    [5] D. J. Benney, Long waves on liquid films, J. Math. and Phys., 45 (1966), 150-155.
    [6] B. Gjevik, Occurrence of Finite-Amplitude Surface Waves on Falling Liquid Films, Phys. Fluids, 13 (1970), 1918-1925.
    [7] S. Kalliadasis, E. A. Demekhin, C. Ruyer-Quil, et al. Falling liquid films, Springer-Verlag, London, 2012.
    [8] S. Kalliadasis, E. A. Demekhin, C. Ruyer-Quil, et al. Thermocapillary instability and wave formation on a film flowing down a uniformly heated plane, J. Fluid Mech., 492 (2003), 303-338. doi: 10.1017/S0022112003005809
    [9] C. Ruyer-Quil, B. Scheid, S. Kalliadasis, et al. Thermocapillary long waves in a liquid film flow. part 1. low-dimensional formulation, J. Fluid Mech., 538 (2005), 199-222.
    [10] B. Scheid, C. Ruyer-Quil, S. Kalliadasis, et al. Thermocapillary long waves in a liquid film flow. part 2. linear stability and nonlinear waves, J. Fluid Mech., 538 (2005), 223-244.
    [11] P. Trevelyan, B. Scheid, C. Ruyer-Quil, et al. Heated falling films, J. Fluid Mech., 592 (2007), 295-334.
    [12] S. Lin, Stability of a liquid flow down a heated inclined plane, Lett. Heat Mass Transfer, 2 (1975), 361-369. doi: 10.1016/0094-4548(75)90002-8
    [13] Y. O. Kabova and V. V. Kuznetsov, Downward flow of a nonisothermal thin liquid film with variable viscosity, J. Appl. Mech. Tech. Phys., 43 (2002), 895-901.
    [14] D. A. Goussis and R. A. Kelly, Effects of viscosity on the stability of film flow down heated or cooled inclined surfaces: Long-wavelength analysis, Phys. Fluids, 28 (1985), 3207-3214.
    [15] C.-C. Hwang and C.-I. Weng, Non-linear stability analysis of film flow down a heated or cooled inclined plane with viscosity variation, Int. J. Heat Mass Transfer, 31 (1988), 1775-1784.
    [16] K. A. Ogden, S. J. D. D'Alessio, J.P. Pascal, Gravity-driven flow over heated, porous, wavy surfaces, Phys. Fluids, 23 (2011), 122102.
    [17] J. P. Pascal, N. Gonputh and S. J. D. D'Alessio, Long-wave instability of flow with temperature dependent fluid properties down a heated incline, Intl J. Eng. Sci., 70 (2013), 73-90. doi: 10.1016/j.ijengsci.2013.05.003
    [18] S. J. D. D'Alessio, C. J. M. P. Seth and J. P. Pascal, The effects of variable fluid properties on thin film stability, Phys. Fluids, 26 (2014), 122105.
    [19] M. Chhay, D. Dutykhb, M. Gisclonb, et al. New asymptotic heat transfer model in thin liquid films, Appl. Math. Model., 48 (2017), 844-859.
    [20] J. H. Spurk and N. Aksel, Fluid Mechanics, 2nd edition, Springer-Verlag, Berlin, 2008.
    [21] C. Ruyer-Quil and P. Manneville, Improved modeling of flows down inclined planes, Eur. Phys. J. B, 15 (2000), 357-369.
    [22] C. Ruyer-Quil and P. Manneville, Further accuracy and convergence results on the modeling of flows down inclined planes by weighted-residual approximations, Phys. Fluids, 14 (2002), 170-183.
    [23] N. Gonputh, J. P. Pascal and S. J. D. D'Alessio, Inclined flow of a heated fluid film with temperature dependent fluid properties, In: Proceedings of the 2012 International Conference on Scientific Computing, pp. 247-253, CSREA Press, 2012.
  • This article has been cited by:

    1. Waleed M. Zapar, Khalaf Gaeid , Hazli Bin Mokhlis Mokhlis, Takialddin A. Al Smadi , Review of the Most Recent Articles in Fault Tolerant Control of Power Plants 2018 – 2022, 2023, 30, 2312-7589, 103, 10.25130/tjes.30.2.11
    2. Mebarek LAHBIB, Mohammed BOUSSAID, Houcine MOUNGAR, Ahmed TAHRI, Analytical assessment of the filed results on the PV system connectors performance under operating temperature, 2023, 11, 2333-8334, 555, 10.3934/energy.2023029
  • Reader Comments
  • © 2019 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(5318) PDF downloads(492) Cited by(3)

Figures and Tables

Figures(9)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog