Next Article in Journal
Optimization of Distillation Sequences with Nonsharp Separation Columns
Next Article in Special Issue
Ultra-Pure Hydrogen via Co-Valorization of Olive Mill Wastewater and Bioethanol in Pd-Membrane Reactors
Previous Article in Journal
Reactive Power Optimization of Large-Scale Power Systems: A Transfer Bees Optimizer Application
Previous Article in Special Issue
Long-Term Stability of Thin-Film Pd-Based Supported Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microstructure Control of Tubular Micro-Channelled Supports Fabricated by the Phase Inversion Casting Method

by
Yuliang Liu
1,
Arash Rahimalimamaghani
1,
Martin van Sint Annaland
2 and
Fausto Gallucci
1,*
1
Inorganic Membranes and Membrane Reactors, Sustainable Process Engineering, Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, De Rondom 70, 5612 AZ Eindhoven, The Netherlands
2
Chemical Process Intensification, Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, De Rondom 70, 5612 AZ Eindhoven, The Netherlands
*
Author to whom correspondence should be addressed.
Processes 2019, 7(6), 322; https://doi.org/10.3390/pr7060322
Submission received: 6 May 2019 / Revised: 19 May 2019 / Accepted: 28 May 2019 / Published: 31 May 2019
(This article belongs to the Special Issue Catalysis in Membrane Reactors)

Abstract

:
Thin-film membrane layers coated onto porous supports is widely considered as an efficient way to obtain high-performance oxygen transport membranes with both good permeability and high mechanical strength. However, conventional preparation methods of membrane supports usually result in highly tortuous channels with high mass transfer resistance. Tubular porous MgO and MgO/CGO supports were fabricated with a simple phase inversion casting method. Long finger-like channels were obtained inside the dual-phase supports by adjusting the ceramic loading, polymer concentration and particle surface area, as well as by introducing ethanol inside the casting slurries. Slurries that exhibit lower viscosity in the zero-shear viscosity region resulted in more pronounced channel growth. These supports were used to produce thin supported CGO membranes for possible application in O2 separation. Similar shrinkage speeds for the different layers during the sintering process are crucial for obtaining dense asymmetric membranes. The shrinkage of the support tube at a high temperature was greatly affected by the polymer/ceramic ratio and compatible shrinkage behaviours of the two layers were realized with polymer/ceramic weight ratios between 0.175 and 0.225.

1. Introduction

Ceramic-based membrane separation is widely accepted as a clean and cost-effective way to produce high purity oxygen from air [1,2,3,4,5]. Additionally, these membranes can be integrated in reactors for partial oxidations, making these membrane reactors more efficient in terms of operational expenditures (OPEX) and capital expenditures (CAPEX). Since the first introduction of the mixed conducting ceramic oxides concept as a base for oxygen transport membranes by Cales and Baumard [6,7], considerable efforts have been invested to increase the oxygen conductivity and long-term stability of these membranes. The bulk-diffusion-based oxygen permeability through the membranes could be described by Wagner Equation [8]:
j O 2 = R T 4 2 F 2 1 L ln P O 2 ln P O 2 σ el σ ion σ el + σ ion d ln p O 2
where jO2 is the oxygen permeation flux, F is the Faraday constant, R is the gas constant, L is the membrane thickness, P’’O2 and PO2 are the low and high oxygen partial pressure and σel and σion are partial electronic and ionic conductivity, respectively. Accordingly, the oxygen permeability could be improved by increasing the partial pressure difference across the membrane, improving the operating temperature, choosing material with high conductivity or reducing the membrane thickness. As a simple way to achieve high oxygen fluxes, thin-film membranes are preferred for separation processes. Since a reduced membrane thickness is accompanied with lower mechanical stability, the thin-film membrane layers were usually coated onto porous supports resulting in asymmetric membrane structures [9,10,11,12,13,14,15]. The support is a key component for the asymmetric membranes, which should provide sufficient mechanical stability, smooth surface for the deposition of thin dense membrane layers and good porosity to circumvent additional mass transfer resistances for oxygen transport.
Selection of proper materials for the support and selective layers is of great importance for the formation of asymmetric membranes with high performance. Permeability of the oxygen selective oxides, material cost, chemical and mechanical stabilities, as well as the compatibility of different layers are the main factors that should be taken into consideration during the material selection. Ce0.9Gd0.1O1.95-δ (CGO) is a fluorite material with high oxygen ionic conductivity and it shows good long term stability under both reducing and oxidizing atmospheres, which makes it an attractive membrane layer material for operation under harsh conditions [16]. To form gas-tight membranes the porous supports and the selective layers should exhibit identical shrinkage behaviour during the sintering process; for this reason most of the reports on multilayer membranes are dedicated to the fabrication of different layers with similar materials [9,10,12,13]. However, the membrane preparation could be more cost-effective if low-price alternatives are found for the supports of the membranes. MgO is a candidate because of its good mechanical and chemical stability and lower cost compared with most of the oxygen selective oxides used for the fabrication of the functional layers. Moreover, due to the similar thermal expansion coefficient between MgO (13.9·10−6 °C-1) and Ce0.9Gd0.1O1.95δ (CGO, 12.7·10−6 °C−1) [17,18], MgO is preferred for the porous substrate when CGO based compounds are employed as the membrane layer material. Ramachandran et al. [16] prepared a tubular MgO support by thermoplastic extrusion and CGO-based layers were deposited onto the support with a repeated dip-coating method. Ovtar et al. [19] also applied the combined extrusion and dip-coating process for fabricating a tubular MgO supported CGO-La0.6Sr0.4FeO3δ (LSF) dual phase membrane.
Thermoplastic extrusion is frequently used to produce tubular substrates of asymmetric membranes. Pores in the support tubes are formed by burning of the sacrificial phases, which are added to the paste prepared before the extrusion process. The resulting support tubes usually show a sponge-like microstructure. Although compared with the single-layer membranes, improved oxygen permeability has been achieved by coating thin-film membrane layers on to porous supports, the highly tortuous channels in the sponge-like support layer (Figure 1a) still result in mass transfer resistances that reduce the oxygen flux. For further decreasing the gas transfer resistance, novel fabrication methods need to be developed. Phase inversion is a promising and low-cost method for preparation of membrane substrates showing low tortuosity. In this process, no pore formers need to be introduced into the ceramic paste. After casted into desired shapes, the ceramic pastes should be submerged into water to perform the phase inversion process. Instead of connected pores, finger like channels are formed because of the solvent replacement with water. The channels provide comparatively short path for gas transfer through the membrane support, which is shown in Figure 1b. He et al. [20] fabricated a Zr0.84Y0.16O1.92 (YSZ)-La0.8Sr0.2MnO3δ (LSM) asymmetric membrane with phase inversion process, which shows improved oxygen permeability than the one with a sponge-like substrate prepared by a sequential two-step tape casting method. The phase inversion method is mostly applied for fabricating asymmetric membranes with hollow fibre [3,4,21,22,23] or planar disc configurations [20,24,25,26], while its application to the manufacture of tubular membranes is quite limited. A porous alumina support tube was successfully prepared by Zhu et al. [27] with a novel phase inversion casting process but the resulting support tubes mainly showed sponge-like microstructure. Length of the channels in the tubular supports still need further improvement.
The formation of uniformly distributed channels of enough length in tubular membrane supports is more challenging than inside the hollow fibres or planar discs, because long micro-channels usually favour low viscosity of the slurries, while a too low slurry viscosity would bring difficulties in shaping the slurry into the tubular form. Optimization of the processing conditions in the phase-inversion-based preparation process to increase the microstructure of the support tubes is thus an area of great interest. In this work, the low cost phase inversion casting process designed by Zhu et al. [27] is applied to fabricate MgO and MgO-CGO dual-phase support tubes. CGO layers were coated onto the substrates with dip-coating method. The effect of ceramic and polymer concentrations, particle size of the ceramic powders and the introduction of ethanol inside the casting slurries on the micro-channel size and tube shrinkage speed during sintering process were investigated to obtain MgO-CGO supports with a low mass transfer resistance and to achieve compatible shrinkage behaviour of the membrane layers and supports.

2. Materials and Methods

2.1. Materials

MgO powders purchased from VWR were used as the support tube material. Polyvinylpyrrolidone (PVP, Sigma-Aldrich, Missouri, USA), polyethersulfone (PES, Veradel 3000P, Solvay Advanced Polymers, Brussels, Belgium) and 1-methyl-2-pyrrolidone (NMP, 99.5%, Sigma-Aldrich, Missouri, USA) were dispersant, polymer and solvent, respectively, for the preparation of the tube-casting slurries. CGO powders supplied from Ningbo Sofcman Energy were used as the second phase in the support tube, after milling for 1 h in a planetary ball mill (Pulverisette 6, Fritsch, Idar-Oberstein, Germany). The influence of particle surface area on the channel growth was studied by applying two other kinds of CGO powders provided by Solvay and Marion Technologies, respectively, which were also ball milled for 1 h before adding them into the slurries.
The milled CGO powders from Ningbo Sofcman Energy were used for the membrane layer fabrication. PVP and PVB purchased from Sigma-Aldrich were used as the dispersant and binder in the suspensions for the dip-coating process, while ethanol (absolute, VWR, Pennsylvania, USA) was employed as the solvent.

2.2. Phase Inversion Casting for Support Tube Preparation

Phase inversion casting is a slurry-based method. Polymer and ceramic concentrations in the slurries are key factors that may influence the microstructure and performance of the support tubes. For preparing the slurries, first, PES and PVP were fully dissolved into NMP by mixing for 20 min using a planetary centrifugal mixer (ARE-250, Thinky, Tokyo, Japan) to form polymer solutions. The rotational speed of the mixer was set at 1850 rpm. MgO and CGO powders (or only MgO powders) were added successively into the solutions and mixed again for 30 min, which resulted in homogeneous slurries. Bubbles inside the slurries may influence mechanical strength and surface smoothness of the support tubes, so de-foaming procedures were applied using a rotational speed of 1850 rpm for 15 min after mixing the slurries in the same planetary centrifugal mixer. The PVP concentration (weight ratio of PVP over PVP+NMP) was fixed at 3.2 wt%, while the ceramic loadings (the weight ratio of CGO+MgO to CGO+MgO+NMP) were varied from 43 wt% to 60 wt% and a range of polymer concentrations (weight ratio of PES to PES+NMP) from 11 wt% to 25 wt% was selected to investigate their influence on the microstructure and shrinkage of the support tubes.
The phase-inversion casting process employed for the preparation of the micro-channelled support tubes is schematically depicted in Figure 2. The simple casting device used for this process is a cylindrical mould with one side connected to a screw cap. An appropriate amount of the de-foamed slurry was used in the mould through the bell mouth and pushed inside the device with a piston for driving air out. The cap of the device was then taken off before placing it on top of water. Then the slurry was pushed out into water with the piston at a controlled speed. The phase inversion process started at the moment the slurry entered water. The slurry was kept inside water for 12 h at a temperature of 20 °C. Solidified green tubular supports were thus obtained in the water bath. Afterward, the green tube was taken off from the piston and dried at room temperature for 2 days before deposition of the membrane layer, because wet surfaces of the porous supports would result in non-uniform membrane layer formation. Outer and inner diameters of the obtained green support tubes were ~12 mm and ~8 mm, respectively.

2.3. Fabrication of Membrane Layer

CGO membrane layers were coated onto the green tubular supports by a dip-coating method. Because the micro-channels grew from outside surface of the porous supports in the water bath, CGO membrane layers were coated inside, while the external surface of the tubes was covered with plastic films during the dip-coating process. Before the dip-coating process, an ethanol-based suspension with a concentration of 23.9 wt% CGO, 2.4 wt% PVP and 1.9 wt% PVB was prepared with the planetary centrifugal mixer. Similar as for the preparation of the casting slurries, PVP and PVB were firstly added into ethanol to make a solution and subsequently the CGO powders were mixed with the solution at a speed of 1850 rpm for 30 min. An automatic dip-coater was employed to dip-coat membrane layers onto the green support tubes. The thickness of the membranes could be precisely controlled by programming the dip-coating speed and times with the dip-coater. The tubes were coated three times to achieve desirable gas-tightness by increasing the thickness of the coating.
The sintering process is of great importance for the fabrication of asymmetric membranes. After drying at room temperature, the membrane layers were sintered together with the green support tubes at 1350 °C for 2 h to form ceramic membranes with good mechanical strength. A heating rate of 30 °C/h was employed to increase the oven temperature to 1350 °C. The low heating rate is required because of two reasons. Firstly, uneven heating could break the ceramic supports and lower heating rate provides enough time for each part of the tubes to reach a uniform temperature. Besides, faster heating rate could result in a larger difference in strain rate between the membrane layers and the support tubes, which would cause cracks on the membrane layers. After sintering, the samples were cooled down inside the oven with a speed of 60 °C/h.

2.4. Analysis of the Membranes

The microstructure of the tubular asymmetric membranes was observed and compared with scanning electron microscopy (SEM, Quanta 3D FEG, Oregon, USA) operating in the secondary-electron mode at an excitation voltage of 5 kV. Before the SEM study, membrane samples were attached to the metallic SEM stages with graphite tapes and gold was sputtered on the membranes using a high-resolution sputter coater (Emitech K575X Peltier-cooled, Quorum Technologies, East Sussex, England), at the constant sputtering current of 65 mA for 30 s, with the purpose of preventing beam damage and allowing charge dissipation on the samples during the SEM study. The element distribution of the CGO-MgO dual-phase supported CGO membrane was obtained using the EDX mapping of a desktop SEM (Phenom ProX, Eindhoven, The Netherlands) at the voltage of 15 kV. A pixel time of 100ms was selected while the map resolution was set at 128×128 pixels. The specific surface area of the ceramic powders was measured with BET nitrogen adsorption isotherm (TriStar II, Micromeritics, Georgia, USA). The slurry viscosity was determined with a rheometer (Physica MCR 301, Anton Paar, Graz, Austria). Nitrogen permeation measurements were conducted to test the mass transfer resistance of the support tubes. The setup for the permeability test mainly consisted of a chamber, a gas supply unit and a film flow meter (STEC SF, Horiba, Kyoto, Japan). The membranes were introduced into the chamber after being sealed with some metal tubes. Nitrogen was fed into the membranes through the tubes using a mass flow controller and the nitrogen permeated through the porous supports was collected at the shell side and sent to the flow meter for the flow rate measurement. In this work, we are only interested in the gas tightness of the membranes, thus no perm-selectivities and oxygen fluxes have been measured. The chemical stability of the dual-phase supported membrane was tested in pure CO2 atmosphere at 850 °C for 500 min. The flowrate of CO2 was set at 480 mL/min. For better contact of the membrane material with the feeding gas, the membrane was crashed into small pieces. The XRD patterns of the samples before and after the stability test were obtained using an X-ray diffractometer (Rigaku MiniFlex, Tokyo, Japan).

3. Results and Discussion

SEM images of the prepared porous MgO-supported CGO membrane are presented in Figure 3a. Partial delamination of the sintered CGO layer from the support tube can be observed from surface and cross section of the asymmetric membrane. This result was reproducible for several membranes and is mainly caused by differences in shrinkage behaviour of the membrane support and the selective layer. To increase attachment of the membrane layer onto the substrate and to improve gas tightness of the supported membrane, a mixture of 60 wt% MgO and 40 wt% CGO powders was used instead of the pure MgO powders for the preparation of the casting slurries for phase inversion. A crack-free CGO membrane layer with a thickness of ~20 µm was obtained on the dual-phase support and strong attachment was achieved between the two layers (Figure 3b). This result was again quite reproducible for several membranes. However, the figure also shows a very homogeneous membrane support, which would result in pronounced mass transfer resistances. Thus, the next step is to create a porous structure in the supports before coating the selective layers. The composition of the casting slurry was varied to investigate the effect on the microstructure and performance of the dual-phase ceramic support.
The element distribution on the cross section of the CGO-MgO dual-phase supported membrane is presented in Figure 4. It can be observed that in the membrane support both the element of Ce and Mg were detected, while only Ce was found in the membrane layer and the distribution of Ce in the membrane layer is denser than that in the support, suggesting that the CGO concentration is higher in the membrane layer than inside the support.

3.1. Influence of Ceramic Loading

The polymer concentration in the solution is the key factor that can influence the channel length in the polymer membranes by affecting the resistance for water transfer into the solution. In addition to the polymer concentration, also the ceramic loading in the slurry is an important parameter that can influence the length of the micro-channels in the ceramic tubes. Moreover, the ceramic concentration also influences the shrinkage behaviour of the support tubes during the sintering process and thereby the adherence of the membrane layers on the supports. When the ceramic loading was lower than 43 wt%, the slurry exhibited inadequate plasticity to maintain the tubular shape during the casting process, while a ceramic loading above 60 wt% caused cracks on the sintered support tube because the chosen PES concentration in the solution was insufficient to handle this ceramic loading. According to this, the ceramic loading was varied from 43 wt% to 60 wt% to investigate its influence on quality of CGO coatings and growth of the micro-channels, while the polymer concentration was fixed at 18 wt%. For each 20 g of NMP, 15 g, 17.5 g, 20 g, 25 g and 30 g PES were added, respectively, in the polymer solution, resulting in ceramic loadings of 43 wt%, 47 wt%, 50 wt% 56 wt% and 60 wt%. Bigger step was applied for the change of ceramic loading from 50 wt% to 60 wt%, because the ceramic content in the slurry shows smaller influence on the slurry viscosity and microstructure of the support at high ceramic loading than at the range from 43 wt% to 50 wt%.
Figure 5 shows SEM images of porous supports prepared at different ceramic loadings and the membrane layers coated on top. It can be observed from cross sections of the membrane supports, that all the support tubes have two layers: one layer with finger-like pores (channels) which provides low mass transfer resistance and the other one with sponge-like microstructure. It is also clear that the ceramic loading in the preparation recipe has a large effect on the thickness of the layer with channels. Sufficiently long and very homogeneous channels were obtained in the support tube, when the ceramic loading was 43 wt%, while the “pore-size” decreased when adding ceramic powders to the casting slurries, which could be attributed to the improved resistance for the channel growth because of the increased viscosity of the slurries (Figure 6). It should be noted that the slurry with a ceramic concentration of 43 wt% showed a rather high fluidity, which means a high pushing speed was required to keep a tubular shape of the slurry when introducing the piston into water.
Large shrinkages of the porous support tubes were observed after the sintering, because of burning of the polymers and sintering of the ceramic particles. The ceramic loading had a great influence on the tube shrinkage. In this work the shrinkage of a support tube was quantified by measuring the tube length before and after the sintering process. Table 1 shows the effect of the ceramic loading in the casting slurry on the shrinkage of the porous supports. The results can be explained by the decrease in the weight ratio of polymers to ceramic powders, where the smaller shrinkage of the support is attributed to reduced polymer burning during the sintering process. Dense membrane layers can be achieved only when the porous supports and membrane layers exhibit uniform shrinkage behaviours. From SEM images of the surfaces of the membrane layers (Figure 5), it is clear that with a ceramic loading of 60 wt%, several cracks were formed on the membrane layer which indicates that the shrinkage of the membrane was faster than that of the support tube. However, when the ceramic loading was below 47 wt%, the MgO-CGO supports shrank faster than the layers, which caused the membrane layers to be corrugated and thus also fragile. Accordingly, in terms of the shrinkage behaviour compatibility between the membranes and supports, ceramic loadings between 50 wt% and 56 wt% (polymer/ceramic ratios between 0.225 and 0.180) were the best choice for preparing the casting slurries, which resulted in selective layers with no connected pores or cracks but some closed pores and surface defects with a pore size of 0.2-1µm. No obvious effect of the closed pores and surface defects was observed on the gas-tightness of the selective layers.
Nitrogen permeation measurements were conducted for the porous supports fabricated with different ceramic loadings. Pressure differences (ΔP) from 10 kPa to 50 kPa were applied across the support tubes and the gas permeation through the supports are plotted as a function of ΔP in Figure 7. A small increase in the nitrogen permeability with ΔP was observed for all the test membranes, because of viscous flow through the small pores [28]. Figure 7 also reveals that the effect of reducing ceramic loading inside the slurry is an increase in gas permeability through the support tubes. Especially when the ceramic loading was decreased from 47 wt% to 43 wt%, a very significant increase in the gas permeation flux was observed, which could be attributed to decrease of the tortuosity of the micro-channels.

3.2. Influence of Polymer Concentration

A high water-solvent exchange rate promotes the growth of micro-channels during the phase inversion. For polymer membranes, the number and size of finger-like channels formed in the water bath tend to decrease when the PES concentration in the solution was increased, which indicates that a high polymer concentration would frustrate the conversion between water and solvent [29]. Based on the research of polymer membranes, in this work the effect of the PES concentration on the channel growth in the ceramic membrane supports was investigated. As with excessive ceramic loadings inside the slurries, when the polymer concentrations were below 11 wt%, the supports generally broke during the sintering process, because the polymer was insufficient to handle the ceramic powders. As a result, a range of polymer concentrations from 11 wt% to 25 wt% was selected to investigate its influence on the microstructure of the supports and membrane layers while the ceramic loading was kept at 50 wt%. For each 20 g of NMP, 2.5 g, 3.5 g, 4.5 g, 5.5 g and 6.5 g PES was added, respectively, in the polymer solution.
As shown in Figure 8, for the chosen ceramic loading, the channel length increases when increasing the polymer concentration in the recipe, which is opposite to the result obtained for polymer membranes. Adding ceramic powders to the polymer solutions caused a large increase in the resistance for water exchange with the solvent, while PES in the suspension performed as a surfactant for the ceramic powders by adsorbing on the surface of the particles [30]. Aggregation of the ceramic particles could be avoided due to the presence of the polymers, so increased PES contents resulted in a reduced resistance for channel growth. In addition, the finger-like pores were found to disappear gradually with increasing the phase inversion time [31] but a higher content of PES, which also performed as binder of the particles, would result in stronger walls of the channels. The rheological behaviour of the slurries with different PES concentrations at a fixed ceramic loading of 50 wt% is shown in Figure 9a. It is evident that all the slurries exhibited shear thinning behaviour (viscosities decreased with increased shear rate) and this behaviour became more pronounced when the polymer concentration was decreased. As a result, at low shear rate the viscosities of the ceramic slurries decreased with increased PES concentrations, while at high shear rate the slurries with low PES concentrations were less viscous than the ones with more PES. The rheological behaviours in the zero-shear viscosity plateau region were consistent with the decrease of water transfer resistances when increasing the PES concentrations.
Table 2 presents how the shrinkages of the porous supports after the sintering process are affected by the polymer concentrations in the slurries. It is clear that increased polymer concentrations result in larger length shrinkages of the support tubes. During sintering, shrinkage of the ceramic substrate can be divided into two stages: the initial stage starts at relatively low temperature, which is caused by the release of the organic polymers, while the second and more significant shrinkage occurs at a much higher temperature due to sintering of the ceramic particles and the reduction of the residual porosity formed by the burning of the polymers. Increase of the polymer content in the green substrates caused larger initial shrinkages during release of the polymers and larger pore volumes, which also resulted in increased shrinkages during the second stage. Microstructures of membrane layers formed on substrates with different polymer concentrations are shown in Figure 8. Good membrane layers were obtained with polymer concentrations between 15 wt% and 18 wt% (polymer/ceramic ratios between 0.175 and 0.225), which is because of the compatible shrinkage behaviour of the substrates and the inside films. No connected pores or cracks but some surface defects and closed pores with a pore size of 0.2-1µm were observed in these selective layers. However, these closed pores and surface defects did not show obvious effect on the gas-tightness of the selective layers. Detachment and cracks of the membrane were formed with a polymer concentration of 25 wt% in the casting slurry, because of over-shrinkage of the support tube.
Figure 10a shows the gas permeation as a function of ΔP for the support tubes prepared with different polymer concentrations, when ceramic loadings fixed at 50 wt%. The slurry with 25 wt% PES resulted in a porous support with the highest permeability and a declined permeation flux was observed when the polymer concentration was reduced to 15 wt%, which is attributed to the reduced channel length and the decreased porosity formed by burning of the polymers. However, decreasing the polymer concentration to 11 wt% did not further reduce the gas permeability, because some micro-cracks were formed on the support tubes (Figure 10b), which could be expected from the insufficient polymer concentration for handling the ceramic loading.
To further examine the role of PES in the phase inversion process for ceramic support fabrication, another group of slurries with a ceramic loading of 43 wt% was prepared and the polymer concentration was increased from 11 wt% to 25 wt% to investigate the growth of micro-channels inside the tubular supports. From the SEM images of the support tube cross sections (Figure 11), it is observed that the longest and most homogeneous channels were achieved when the polymer concentration was 18 wt%, which means that this slurry provided the lowest resistance for water-solvent exchange. Below the polymer concentration of 18 wt%, the channel length increased with the PES content, which could be attributed to the dispersing effect of PES. However, a reduced channel size was observed when more PES was added into the slurry with polymer concentration of 18 wt% and this can be understood by considering that, when the surfaces of the particles were fully covered with polymer, excess polymer increases the resistance for phase inversion, as for the polymer membranes. Consistent results were also observed for the rheological behaviour of the slurries at the zero-shear viscosity plateau region (Figure 9b), which means that the lowest viscosity of the slurries (lowest channel growth resistance) was obtained at the polymer concentration of 18 wt%. Dong et al. [30] reported similar rheological behaviour for NiO/YSZ ceramic slurries by adjusting the dispersant concentration in the suspension. By comparing the results of these two groups of experiments with the ceramic loadings fixed at 50 wt% and 43 wt%, it is evident that an increased ceramic loading requires more polymer as surfactant to achieve sufficiently long channels in the ceramic supports because of the improved specific surface areas of the powders.

3.3. Influence of Particle Surface Area

Three kinds of ceramic powders with BET specific surface areas (SSA) of 12.53 m2/g, 24.11 m2/g and 26.33 m2/g, respectively, were applied to the preparation of ceramic casting slurries for phase inversion. The PES concentration was fixed at 25 wt%, while the ceramic loading was 50 wt%. The microstructures of the support tubes fabricated with these three slurries are shown in Figure 12. The channel length decreased significantly when the SSA of the ceramic particles was increased from 12.53 m2/g to 24.11 m2/g, leading to an abrupt decrease in the gas permeability through the porous support (Figure 13). This result further confirms the surfactant role of PES by coating on the surface of the ceramic particles. Larger SSA of the powders requires more surfactant to reduce the resistance for the water-solvent conversion. Accordingly, for a specific polymer concentration, the slurry with the lower SSA might have a lower resistance for channel growth.
Figure 14 reports the effect of the ceramic powder SSA on the slurry viscosities. At low shear rate, the CGO powders with an SSA of 12.53 m2/g resulted in a much lower slurry viscosity than the ones with an SSA of 24.11 m2/g and 26.33 m2/g.

3.4. Addition of Ethanol to the Ceramic Slurry

Smolders et al. [32] proposed a mechanism for the formation of macrovoids (channels), by attributing the occurrence of the channels inside polymer membranes to the instantaneous liquid-liquid demixing and the formation of membranes without channels to the delayed demixing. They also concluded that the low non-solvent (water) concentration in the diffusion front would cause delayed demixing, which would hinder growth of the channels, while adding a certain amount of water to the polymer solution could change delayed demixing to instantaneous demixing and increase the channel length.
In this work, also the effect of adding a non-solvent to the casting slurries on the microstructures of the ceramic membrane supports was also investigated. Ethanol was applied instead of water as the non-solvent addictive because of the low solubility of PES in water, which would cause precipitation of the PES before the phase inversion casting process when adding a small amount of water into the polymer solutions. Excess ethanol in the mixture would also cause precipitation of the polymers, which means that the amount of ethanol added to the solution should be precisely controlled. Various NMP-ethanol-PES-PVP solutions with ethanol contents (weight ratios of ethanol to ethanol+NMP) ranging from 0 wt% to 27 wt% were prepared. Precipitation of PES occurred at the surface of the polymer solution when ethanol was dripped into the mixture because of the lower density of ethanol compared to NMP. However, the solution tended to become pellucid after mixing for 5 min in the centrifugal mixer, when ethanol mixed thoroughly with the solution. Then, the ceramic powders were added into the solution to obtain the ceramic casting slurries. The PES concentration (weight ratio of PES to PES+ethanol+NMP) was fixed at 18 wt%, while the ceramic loading (weight ratio of ceramic powders to CGO+MgO+ethanol+NMP) was 50 wt% in this group of experiments.
From the cross sections of the obtained ceramic supports (Figure 15), it can be observed that the length of the channels tended to increase when the ethanol contents were increased from 0 wt% to 27 wt%. Jin et al. [31] employed the solubility difference parameter (|∆δ|) for the evaluation of the driving force to form finger-like channels. By applying liquids with different solubilities as the coagulants to adjust the driving force, they verified that the figure-like channels are easier to form when |∆δ| between NMP and the coagulant is improved. The solubility parameters for NMP, water and ethanol are 22.9, 47.9 and 26.0 respectively. In the channel growth front, the water concentration is lower (NMP concentration is higher) than at the initial interface between the slurry and water bath, which may result in a declined solubility of the coagulant and the resulting |∆δ| is insufficient for the formation of long enough channels. Since the solubility parameter of ethanol is higher than that of NMP, adding a certain amount of ethanol to the casting slurry can help increasing the solubility in the diffusion front and improve the driving force for channel growth. Moreover, in Figure 16 it is shown that the slurry viscosity at the zero-shear viscosity plateau region declined when the ethanol concentration was increased. According to the former results, a decreased low shear rate viscosity would promote channel growth inside the ceramic membrane supports.

3.5. Chemical Stability of the Membrane

The XRD patterns of the dual-phase supported CGO membrane pieces before and after the chemical stability test in CO2 atmosphere, together with the XRD patterns of MgO and CGO powders are presented in Figure 17. It indicated that CGO and MgO remain unchanged in the CGO-MgO dual-phase material after the heat treatment during the membrane preparation process. Besides, it can also be observed that after the stability test, the membrane pieces completely kept their dual phases (CGO and MgO). No carbonate was formed. The chemical stability of the membrane in CO2 atmosphere was confirmed to be good during the testing period (500 min).

4. Conclusions

Micro-channelled MgO and MgO-CGO dual-phase support tubes were fabricated with a low-cost phase inversion casting method. Gas-tight CGO membrane layers were obtained on the dual-phase supports, while the layers tended to peel off from the MgO tubes after sintering. The composition of the casting slurries was a key factor that influenced channel length and shrinkage of the MgO-CGO supports. Increased ceramic loading frustrated the channel growth, causing decreased gas permeability through the support tubes. PES performed as surfactant of the ceramic powders, which means that addition of a certain amount of PES to the casting slurry or a decrease in the specific surface area of the ceramic particles can help reducing the channel-growth resistances. However, saturated adsorption of polymers onto the ceramic powders resulted in a lower water-solvent conversion speed and reduced channel length. Moreover, addition of ethanol to the casting slurry was found to be an effective way to increase the channel length in the casted support tubes. Slurries that provided lower channel-growth resistances exhibited lower viscosities at the zero-shear plateau region. An increase in the PES/ceramic ratio caused larger shrinkage of the tubes during the sintering process and compatible shrinkage behaviours between the dense membranes and substrates were achieved when PES/ceramic ratios were between 0.175 and 0.225. The chemical stability of the dual-phase supported membrane in CO2 atmosphere was confirmed to be good during a testing period of 500 min.

Author Contributions

Design of the experiment, Y.L. and F.G.; performing the experiment, Y.L. and A.R.; data analysis, Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, F.G.; supervision, M.v.S.A. and F.G.

Funding

This research received no external funding.

Acknowledgments

We would like to acknowledge China Scholarship Council for supporting this work by providing PhD scholarship to Y.L.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dyer, P.N.; Richards, R.E.; Russek, S.L.; Taylor, D.M. Ion transport membrane technology for oxygen separation and syngas production. Solid State Ion. 2000, 134, 21–33. [Google Scholar] [CrossRef]
  2. Liu, Y.; Tan, X.; Li, K. Mixed conducting ceramics for catalytic membrane processing. Catal. Rev. 2006, 48, 145–198. [Google Scholar] [CrossRef]
  3. Liu, H.; Tan, X.; Pang, Z.; Diniz da Costa, J.C.; Lu, G.Q.; Liu, S. Novel dual structured mixed conducting ceramic hollow fibre membranes. Sep. Purif. Technol. 2008, 63, 243–247. [Google Scholar] [CrossRef]
  4. Bi, X.; Meng, X.; Liu, P.; Yang, N.; Zhu, Z.; Ran, R.; Liu, S. A novel CO2-resistant ceramic dual-phase hollow fiber membrane for oxygen separation. J. Memb. Sci. 2017, 522, 91–99. [Google Scholar] [CrossRef]
  5. Badwal, S.P.S.; Ciacchi, F.T. Ceramic membrane technologies for oxygen separation. Adv. Mater. 2001, 13, 993–996. [Google Scholar] [CrossRef]
  6. Cales, B.; Baumard, J.F. Mixed conduction and defect structure of ZrO2-CeO2-Y2O3 solid solutions. J. Electrochem. Soc. 1984, 131, 2407–2413. [Google Scholar] [CrossRef]
  7. Cales, B.; Baumard, J.F. Oxygen semipermeability and electronic conductivity in calcia-stabilized zirconia. J. Mater. Sci. 1982, 17, 3243–3248. [Google Scholar] [CrossRef]
  8. Bouwmeester, H.J.M.; Burggraaf, A.J. Dense Ceramic Membranes for Oxygen Separation. J. Mater. Sci. 1996, 4, 435–528. [Google Scholar]
  9. Souza, D.F.; Nunes, E.H.M.; Vasconcelos, W.L. Preparation of Ba0.5Sr0.5Co0.8Fe0.2O3-δ asymmetric structures by freeze-casting and dip-coating. Ceram. Int. 2018, 44, 1002–1006. [Google Scholar] [CrossRef]
  10. Rachadel, P.L.; Souza, D.F.; Nunes, E.H.M.; da Costa, J.C.D.; Vasconcelos, W.L.; Hotza, D. A novel route for manufacturing asymmetric BSCF-based perovskite structures by a combined tape and freeze casting method. J. Eur. Ceram. Soc. 2017, 37, 5249–5257. [Google Scholar] [CrossRef]
  11. Xing, W.; Fontaine, M.L.; Li, Z.; Polfus, J.M.; Larring, Y.; Denonville, C.; Nonnet, E.; Stevenson, A.; Henriksen, P.P.; Bredesen, R. Asymmetric tubular CaTi0.6Fe0.15Mn0.25O3-δ membranes: Membrane architecture and long-term stability. J. Memb. Sci. 2018, 548, 372–379. [Google Scholar] [CrossRef]
  12. Liu, T.; Chen, Y.; Fang, S.; Lei, L.; Wang, Y.; Ren, C.; Chen, F. A dual-phase bilayer oxygen permeable membrane with hierarchically porous structure fabricated by freeze-drying tape-casting method. J. Memb. Sci. 2016, 520, 354–363. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, Z.; Zhang, G.; Dong, X.; Jiang, W.; Jin, W.; Xu, N. Fabrication of asymmetric tubular mixed-conducting dense membranes by a combined spin-spraying and co-sintering process. J. Memb. Sci. 2012, 415–416, 313–319. [Google Scholar] [CrossRef]
  14. Yin, X.; Hong, L.; Liu, Z. Oxygen permeation through the LSCO-80/CeO2 asymmetric tubular membrane reactor. J. Memb. Sci. 2006, 268, 2–12. [Google Scholar] [CrossRef]
  15. Teraoka Yasutake, N.Y.; Fukuda, T.; Miura, N. Development of oxygen semipermeable membrane using mixed conductive perovskite-type oxides (part 2): Preparation of dense film of perovskite-type oxide on porous substrate. J. Ceram. Soc. Jpn. 1989, 97, 533–538. [Google Scholar] [CrossRef]
  16. Ramachandran, D.K.; Søgaard, M.; Clemens, F.; Gurauskis, J.; Kaiser, A. Fabrication and performance of a tubular ceria based oxygen transport membrane on a low cost MgO support. Sep. Purif. Technol. 2015, 147, 422–430. [Google Scholar] [CrossRef] [Green Version]
  17. Tietz, F. Thermal expansion of SOFC materials. Ionics 1999, 5, 129–139. [Google Scholar] [CrossRef]
  18. Shiratori, Y.; Tietz, F.; Buchkremer, H.P.; Stöver, D. YSZ-MgO composite electrolyte with adjusted thermal expansion coefficient to other SOFC components. Solid State Ion. 2003, 164, 27–33. [Google Scholar] [CrossRef]
  19. Ovtar, S.; Gurauskis, J.; Haugen, A.B.; Chatzichristodoulou, C.; Kaiser, A.; Hendriksen, P.V. Oxygen transport properties of tubular Ce0.9Gd0.1O1.95-La0.6Sr0.4FeO3−δ composite asymmetric oxygen permeation membranes supported on magnesium oxide. J. Memb. Sci. 2017, 523, 576–587. [Google Scholar] [CrossRef]
  20. He, W.; Huang, H.; Gao, J.; Winnubst, L.; Chen, C. Phase-inversion tape casting and oxygen permeation properties of supported ceramic membranes. J. Memb. Sci. 2014, 452, 294–299. [Google Scholar] [CrossRef]
  21. Tan, X.; Liu, Y.; Li, K. Preparation of LSCF ceramic hollow-fiber membranes for oxygen production by a phase-inversion/sintering technique. Ind. Eng. Chem. Res. 2005, 44, 61–66. [Google Scholar] [CrossRef]
  22. Li, K.; Tan, X.; Liu, Y. Single-step fabrication of ceramic hollow fibers for oxygen permeation. J. Memb. Sci. 2006, 272, 1–5. [Google Scholar] [CrossRef]
  23. Wu, Z.; Hidayati Othman, N.; Zhang, G.; Liu, Z.; Jin, W.; Li, K. Effects of fabrication processes on oxygen permeation of Nb2O5-doped SrCo0.8Fe0.2O3-δ micro-tubular membranes. J. Memb. Sci. 2013, 442, 1–7. [Google Scholar] [CrossRef]
  24. Meng, Y.; He, W.; Li, X.; Gao, J.; Zhan, Z.; Yi, J.; Chen, C.; Bouwmeester, H.J.M. Asymmetric La0.6Sr0.4Co0.2Fe0.8O3-δ membrane with reduced concentration polarization prepared by phase-inversion tape casting and warm pressing. J. Memb. Sci. 2017, 533, 11–18. [Google Scholar] [CrossRef]
  25. Huang, H.; Cheng, S.; Gao, J.; Chen, C.; Yi, J. Phase-inversion tape-casting preparation and significant performance enhancement of Ce0.9Gd0.1O1.95-La0.6Sr0.4Co0.2Fe0.8O3-δ dual-phase asymmetric membrane for oxygen separation. Mater. Lett. 2014, 137, 45–248. [Google Scholar] [CrossRef]
  26. Shao, X.; Dong, D.; Parkinson, G.; Li, C.Z. Thin ceramic membrane with dendritic microchanneled sub structure and high oxygen permeation rate. J. Memb. Sci. 2017, 541, 653–660. [Google Scholar] [CrossRef]
  27. Zhu, Z.; Xiao, J.; He, W.; Wang, T.; Wei, Z.; Dong, Y. A phase-inversion casting process for preparation of tubular porous alumina ceramic membranes. J. Eur. Ceram. Soc. 2015, 35, 3187–3194. [Google Scholar] [CrossRef]
  28. Ramachandran, D.K.; Søgaard, M.; Clemens, F.; Sudireddy, B.R.; Kaiser, A. Low cost porous MgO substrates for oxygen transport membranes. Mater. Lett. 2016, 169, 254–256. [Google Scholar] [CrossRef] [Green Version]
  29. Barzin, J.; Madaeni, S.S.; Mirzadeh, H.; Mehrabzadeh, M. Effect of polyvinylpyrrolidone on morphology and performance of hemodialysis membranes prepared from polyether sulfone. J. Appl. Polym. Sci. 2004, 92, 3804–3813. [Google Scholar] [CrossRef]
  30. Dong, D.; Gao, J.; Liu, X.; Meng, G. Fabrication of tubular NiO/YSZ anode-support of solid oxide fuel cell by gelcasting. J. Power Sources 2007, 165, 217–223. [Google Scholar] [CrossRef]
  31. Jin, C.; Yang, C.; Chen, F. Effects on microstructure of NiO-YSZ anode support fabricated by phase-inversion method. J. Memb. Sci. 2010, 363, 250–255. [Google Scholar] [CrossRef]
  32. Smolders, C.A.; Reuvers, A.J.; Boom, R.M.; Wienk, I.M. Microstructures in phase-inversion membranes. Part 1. Formation of macrovoids. J. Memb. Sci. 1992, 73, 259–275. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic illustration of gas transport through porous supports with (a) sponge-like microstructure and (b) micro channels.
Figure 1. Schematic illustration of gas transport through porous supports with (a) sponge-like microstructure and (b) micro channels.
Processes 07 00322 g001
Figure 2. Schematic representation of the support tube preparation process by the phase-inversion casting method.
Figure 2. Schematic representation of the support tube preparation process by the phase-inversion casting method.
Processes 07 00322 g002
Figure 3. Scanning electron microscopy (SEM) images of CGO layers fabricated on (a) MgO and (b) MgO-CGO dual-phase support tubes.
Figure 3. Scanning electron microscopy (SEM) images of CGO layers fabricated on (a) MgO and (b) MgO-CGO dual-phase support tubes.
Processes 07 00322 g003
Figure 4. Element distribution on the cross section of the CGO-MgO dual-phase supported membrane: (a) Combined mapping and SEM image; (b) Combined element mapping; (c) Mg mapping; (d) Ce mapping.
Figure 4. Element distribution on the cross section of the CGO-MgO dual-phase supported membrane: (a) Combined mapping and SEM image; (b) Combined element mapping; (c) Mg mapping; (d) Ce mapping.
Processes 07 00322 g004
Figure 5. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 60 wt%, (b) 56 wt%, (c) 50 wt%, (d) 47 wt% and (e) 43 wt% ceramic loadings together with the CGO layers coated inside.
Figure 5. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 60 wt%, (b) 56 wt%, (c) 50 wt%, (d) 47 wt% and (e) 43 wt% ceramic loadings together with the CGO layers coated inside.
Processes 07 00322 g005aProcesses 07 00322 g005b
Figure 6. Rheological behaviour of the slurries with different ceramic loadings.
Figure 6. Rheological behaviour of the slurries with different ceramic loadings.
Processes 07 00322 g006
Figure 7. Nitrogen permeability through MgO-CGO supports fabricated from slurries with different ceramic loadings.
Figure 7. Nitrogen permeability through MgO-CGO supports fabricated from slurries with different ceramic loadings.
Processes 07 00322 g007
Figure 8. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 11 wt%, (b) 15 wt%, (c) 18 wt%, (d) 22 wt% and (e) 25 wt% polymer concentrations together with the CGO layers coated inside.
Figure 8. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 11 wt%, (b) 15 wt%, (c) 18 wt%, (d) 22 wt% and (e) 25 wt% polymer concentrations together with the CGO layers coated inside.
Processes 07 00322 g008aProcesses 07 00322 g008b
Figure 9. Rheological behaviour of the slurries with different polymer concentrations when the ceramic loadings fixed at (a) 50 wt% and (b) 43 wt%.
Figure 9. Rheological behaviour of the slurries with different polymer concentrations when the ceramic loadings fixed at (a) 50 wt% and (b) 43 wt%.
Processes 07 00322 g009
Figure 10. (a) Nitrogen permeability through MgO-CGO supports fabricated from slurries with different polymer concentrations; (b) SEM image of cracks on the support tube made from a slurry with 11 wt% polymer concentration.
Figure 10. (a) Nitrogen permeability through MgO-CGO supports fabricated from slurries with different polymer concentrations; (b) SEM image of cracks on the support tube made from a slurry with 11 wt% polymer concentration.
Processes 07 00322 g010
Figure 11. SEM images of MgO-CGO supports fabricated from slurries with (a) 11 wt%, (b) 15 wt%, (c) 18 wt%, (d) 22 wt% and (e) 25 wt% polymer concentration with a fixed ceramic loading of 43 wt%.
Figure 11. SEM images of MgO-CGO supports fabricated from slurries with (a) 11 wt%, (b) 15 wt%, (c) 18 wt%, (d) 22 wt% and (e) 25 wt% polymer concentration with a fixed ceramic loading of 43 wt%.
Processes 07 00322 g011
Figure 12. SEM images of MgO-CGO porous supports fabricated using ceramic powders with specific surface areas of (a) 12.53 m2/g, (b) 24.11 m2/g and (c) 26.33 m2/g.
Figure 12. SEM images of MgO-CGO porous supports fabricated using ceramic powders with specific surface areas of (a) 12.53 m2/g, (b) 24.11 m2/g and (c) 26.33 m2/g.
Processes 07 00322 g012
Figure 13. Nitrogen permeability through MgO-CGO supports fabricated with ceramic powders of different specific surface areas.
Figure 13. Nitrogen permeability through MgO-CGO supports fabricated with ceramic powders of different specific surface areas.
Processes 07 00322 g013
Figure 14. Rheological behaviour of the slurries prepared with ceramic powders of different specific surface areas.
Figure 14. Rheological behaviour of the slurries prepared with ceramic powders of different specific surface areas.
Processes 07 00322 g014
Figure 15. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 0 wt%, (b) 7 wt%, (c) 13 wt%, (d) 20 wt% and (e) 27 wt% ethanol concentrations.
Figure 15. SEM images of MgO-CGO porous supports fabricated from slurries with (a) 0 wt%, (b) 7 wt%, (c) 13 wt%, (d) 20 wt% and (e) 27 wt% ethanol concentrations.
Processes 07 00322 g015aProcesses 07 00322 g015b
Figure 16. Rheological behaviour of the slurries prepared with different ethanol concentrations.
Figure 16. Rheological behaviour of the slurries prepared with different ethanol concentrations.
Processes 07 00322 g016
Figure 17. X-ray diffraction (XRD) patterns of (a) CGO powders, (b) MgO powders, (c) CGO-MgO dual-phase supported CGO membrane before the chemical stability test, (d) CGO-MgO dual-phase supported CGO membrane after the chemical stability test.
Figure 17. X-ray diffraction (XRD) patterns of (a) CGO powders, (b) MgO powders, (c) CGO-MgO dual-phase supported CGO membrane before the chemical stability test, (d) CGO-MgO dual-phase supported CGO membrane after the chemical stability test.
Processes 07 00322 g017
Table 1. Shrinkage comparison of the coated support tubes fabricated from slurries with different ceramic loadings.
Table 1. Shrinkage comparison of the coated support tubes fabricated from slurries with different ceramic loadings.
Ceramic Loading (wt%)Polymer/Ceramic RatioShrinkage of Support Tube (%)Microstructure of CGO Layer
600.15033.6Cracked
560.18035.6Smooth
500.22536.7Smooth
470.25737.1Corrugated
430.30038.6Corrugated
Table 2. Shrinkage comparison of the support tubes fabricated from slurries with different polymer concentrations.
Table 2. Shrinkage comparison of the support tubes fabricated from slurries with different polymer concentrations.
Polymer Concentration (wt%)Polymer/Ceramic RatioShrinkage of Support Tube (%)Microstructure of CGO Layer
110.12531.8Cracked
150.17535.3Smooth
180.22536.7Smooth
220.27538.3Corrugated
250.32538.7Corrugated

Share and Cite

MDPI and ACS Style

Liu, Y.; Rahimalimamaghani, A.; van Sint Annaland, M.; Gallucci, F. Microstructure Control of Tubular Micro-Channelled Supports Fabricated by the Phase Inversion Casting Method. Processes 2019, 7, 322. https://doi.org/10.3390/pr7060322

AMA Style

Liu Y, Rahimalimamaghani A, van Sint Annaland M, Gallucci F. Microstructure Control of Tubular Micro-Channelled Supports Fabricated by the Phase Inversion Casting Method. Processes. 2019; 7(6):322. https://doi.org/10.3390/pr7060322

Chicago/Turabian Style

Liu, Yuliang, Arash Rahimalimamaghani, Martin van Sint Annaland, and Fausto Gallucci. 2019. "Microstructure Control of Tubular Micro-Channelled Supports Fabricated by the Phase Inversion Casting Method" Processes 7, no. 6: 322. https://doi.org/10.3390/pr7060322

APA Style

Liu, Y., Rahimalimamaghani, A., van Sint Annaland, M., & Gallucci, F. (2019). Microstructure Control of Tubular Micro-Channelled Supports Fabricated by the Phase Inversion Casting Method. Processes, 7(6), 322. https://doi.org/10.3390/pr7060322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop