Next Article in Journal
A Novel Application on the Drive Elements of Using Electrical Contact Resistance and Friction Coefficient for Evaluating Induction Heat Treatment
Previous Article in Journal
Nano-Scale Modifications of Amniotic Membrane Induced by UV and Antibiotic Treatment: Histological, AFM and FTIR Spectroscopy Evidence
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 14
Issue 4
10.3390/ma14040864
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Facile and Scalable Approach in the Fabrication of Tailored 3D Graphene Foam via Freeze Drying

by
Tony Thomas
and
Arvind Agarwal
*
Plasma Forming Laboratory, Mechanical and Materials Engineering, Florida International University; Miami, FL 33174, USA
*
Author to whom correspondence should be addressed.
Materials 2021, 14(4), 864; https://doi.org/10.3390/ma14040864
Submission received: 1 January 2021 / Revised: 5 February 2021 / Accepted: 8 February 2021 / Published: 11 February 2021
(This article belongs to the Special Issue Graphene Foam Based Composites)
Browse Figures
Versions Notes

Abstract

:
One of the challenges in the processing of advanced composite materials with 2D reinforcement is their extensive agglomeration in the matrix. 3D architecture of 2D graphene sheets into a Graphene Foam (GrF) assembly has emerged as an effective way to overcome agglomeration. The highly reticulated network of branches and nodes of GrF offers a seamless pathway for photon and electron conduction in the matrix along with improved mechanical properties. 3D GrF nano-filler is often fabricated by chemical vapor deposition (CVD) technique, which demands high energy, slow deposition rate, and restricting production to small scale. This work highlights freeze-drying (FD) technique to produce 3D graphene nanoplatelets (GNP) foam with a similar hierarchical structure to the CVD GrF. The FD technique using water as the main chemical in 3D GNP foam production is an added advantage. The flexibility of the FD in producing GNP foams of various pore size and morphology is elucidated. The simplicity with which one can engineer thermodynamic conditions to tailor the pore shape and morphology is presented here by altering the GNP solid loading and mold geometry. The FD 3D GNP foam is mechanically superior to CVD GrF as it exhibited 1280 times higher elastic modulus. However, thermal diffusivity of the FD GNP foam is almost 0.5 times the thermal diffusivity of the CVD GrF due to the defects in GNP particles and pore architecture. The versatility in GNP foam scalability and compatibility to form foam of other 1D and 2D material systems (e.g., carbon nanotubes, boron nitride nanotubes, and boron nitride nanoplatelets) brings a unique dimensionality to FD as an advanced engineering foam development process.

Graphical Abstract

1. Introduction

Over the past five decades, extensive progress has been achieved in understanding the effect of reinforcing materials in enhancing the matrix properties [1,2,3,4]. Reinforced composite materials are broadly adapted in all engineering applications due to a unique combination of material properties they can manifest, for example lightweight yet strong, thermally insulating yet electrically conducting, and tailored fracture mechanics. Reinforcement materials are categorized into fiber, flakes, particulates, and fillers [5]. The choice of reinforcement materials is based on the desired intrinsic property enhancement in the material matrix. For example, particulates are not effective in enhancing the toughness or fracture resistance in a material, as they randomly orient themselves in the material matrix. Typically, flakes or two-dimensional (2D) platelet material forms of reinforcement improve the composites’ physical and mechanical properties. 2D platelets are packed parallel to one another, resulting in higher density packing than fiber reinforcements. Hence, 2D reinforcements can produce composite materials with boosted material properties. Carbonaceous nanomaterials have shown promising potential as a reinforcement material in developing advanced engineering composites. Graphene (Gr), a 2D material, which is a single layer of sp2 hybridized carbon atom, has inspired material researchers and industry owing to its impressive arsenal of properties, as summarized in Figure 1 [6,7,8].
Gr offers ~3 times the specific surface area than the carbon nanotubes (CNT), and hence its addition has proved to improve the overall material properties such as strength, stiffness, and thermal and electrical conductivity [9]. Such augmentation in the matrix material properties due to Gr addition has led to its potential application in energy devices, semiconductors, heat sinks, protective coating, etc. Figure 2 summarizes potential applications of Gr based materials [6,7,10].
2D Gr sheets have a higher tendency to agglomerate due to intermolecular π–π bonding [11]. These agglomerates form clusters that can restack to form graphitic structure and act as a stress concentrator leading to an adverse effect on the mechanical properties. Moreover, lack of homogeneous distribution of 2D Gr in the non-conducting (both thermal and electrical) material matrix, as shown in Figure 3a [low temperature co-fired dense ceramic reinforced with graphene nano-platelets (GNP)] produced by spark plasma sintering), causes a discontinuity in the contact between Gr flakes, as represented in Figure 3b, affecting the phonon and electron conduction path, thus limiting the thermal and electrical conductivity in the composite [12,13]. Although various physical and chemical dispersion techniques such as sonication, surface functionalization, and ball milling have been employed to overcome the agglomeration, Gr tends to align randomly. Other strategies such as crumbling the graphene sheets and adding spacers have also been explored to overcome the agglomeration challenge [14,15].
Graphene foam with a highly reticulated hierarchical structure, as shown in Figure 3c, has proven to be an effective solution against agglomeration and restacking of 2D Gr flakes. This macroporous, free-standing, 3D graphene foam (GrF) has emerged as a prevalent nanofiller material [16,17,18,19,20]. 3D GrF has a homogenous microstructure and can be readily introduced into a low viscous material matrix without complicated dispersion techniques. The reticulated architecture comprises an interconnected network of branches and nodes, as shown in Figure 3d, which facilitates an uninterrupted pathway for phonon and electron conduction, unlike the 2D Gr flakes (refer to Figure 3b). The pristine free-standing GrF shown in Figure 3c offers filler properties such as ultra-low density (<4 mg/cm3), enhanced surface area (~850 m2/g), and improved electron and phonon conduction due to reduced inter Gr sheet contact resistance [21,22,23,24]. As a result, 3D GrF has been used in producing several composite materials for applications such as scaffolds [21], strain sensors [16], vibration dampeners [12,25], supercapacitors [26], fracture-resistant materials [27], and thermal interfacing [28]. The pristine 3D GrF is often fabricated via a template-directed chemical vapor disposition (CVD) technique [24,25,26,27,28]. Since its advent in 2011, from the trend analysis of scientific publications on CVD GrF (Figure 4, data from the web of science), it is evident that there is a tremendous interest in tapping into intrinsic properties of graphene arranged in a 3D hierarchical structure.
The procedure of producing 3D GrF by CVD technique is shown in the schematic in Figure 5 [29,30,31]. A porous metal foam with a reticulated structure, typically nickel with the desired pore shape and size, is chosen as the graphene deposition template. A hydrocarbon gas, usually methane (CH4), is used as the carbon source. At a high temperature of 1000 °C, carbon released from the methane decomposition is deposited on the porous metal template. Due to a catalytic reaction, thin layers of Gr films are precipitated on the nickel template. Once the required layers of Gr film are deposited on the template, a poly methyl methacrylate (PMMA) is deposited on the precipitated Gr film. PMMA prevents the Gr network’s collapse upon removing the metal template by chemical etching, as shown in Figure 5. Hydrochloric or ferric chloride acid is the typical etchant used. Finally, the etched GrF foam structure undergoes an acetone bath to dissolve the PMMA and obtain a monolithic 3D GrF with a highly reticulated system [24].
Although the CVD technique produces high-quality GrF with controlled layer numbers and crystallinity [24], it has some limitations. One of the main challenges in producing 3D GrF with CVD is scalability, as the size of the foam is limited to the reaction chamber capacity. Having the metal substrate with the desired architecture brings further challenges. The additional challenge includes obtaining a completely uniform graphene layer on the substrate, as the kinetics of gas diffusion change within the reaction chamber’s volume [32,33]. This affects the chemical reactions on the metal substrate. The current CVD technique to produce free-standing GrF requires specialized equipment and is a highly energy driven process as it demands elevated temperature for chemical reactions [24]. The byproducts of the process, such as reaction gases, are extremely combustible. The metal etchant acid used to separate the substrate from Gr film is corrosive and toxic. Special reformers are required to neutralize these toxic byproducts. As a result, the pristine GrF produced is expensive, and not easily scalable. To address the limitations mentioned above and the disadvantages of the CVD technique in the fabrication of GrF, an eco-friendly process known as Freeze Drying (FD) is explored in this research. The procedure involves freezing an aqueous suspension, followed by sublimation of the frozen aqueous solvent at low pressure to obtain porous architecture of suspended materials. A variety of materials can be subjected to FD, which suggests that the process relies on the physical interaction of the slurry materials rather than their chemical interactions, unlike CVD [34,35,36,37].
Over the past two decades, FD has emerged as a promising technique to fabricate porous materials with high surface area for biomedical, structural, and energy applications [34,35,36,37,38]. FD was first developed as a single-step forming process to produce dense ceramic and only later it was realized that the thermodynamic control of the freezing vehicle could yield porous material with hierarchical structure [39]. In FD, the particles in the slurry are rejected from the propagating crystal front, collecting the particles between the growing crystal front. Thus, by controlling the heat transfer rate and direction, isotropic or anisotropic solidification can be introduced in the slurry system to get directional pores with uniform size or random pores of varying size [40]. The solidification time influences the pore size and the wall thickness. Solidification time increases linearly with the thickness of the sample [34]. Hence, optimal temperature control and control of formulation is paramount to obtain a 3D foam structure with isotropic properties. Many graphene and graphene based composite 3D foam structures of various chemical binders and dispersants have been produced by FD for applications such as flexible supercapacitors, electromagnetic shielding, sensors and catalyst support [40,41,42,43,44,45].
Hence, herein the FD technique’s simplicity is demonstrated in artificially stacking the graphene nanoplatelets (GNP) into a 3D hierarchical structure similar to the CVD produced GrF. The FD process’s simplicity in designing 3D GNP foam with controlled pore size and morphology is explored as this is the main challenge with the CVD technique. Here, two methodologies were used successfully to produce GNP foam with different pore size and morphology: (i) regulating the heat transfer rate by changing the amount of GNP (solid loading) in the aqueous slurry; and (ii) controlling the heat transfer direction by changing the mold geometry. The 3D GNP foam produced by FD is compared with the CVD GrF to understand the effect of processing on the mechanical and thermal properties. This study also demonstrates that 3D foams of other nanomaterials can be produced by the simple, eco-friendly, and cost-effective FD technique.

2. Experimental

2.1. Slurry and 3D GNP Foam Preparation

Graphene nanoplatelets used in this study were purchased from XG Science, Lancing, MI, USA (Avg. particle size 15 µm and thickness of 6–14 nm, bulk density 0.0–0.1 g cm−3). Deionized water was the primary solvent used. Carboxymethyl cellulose (CMC, ~1.1 wt.%, Mw ~700,000, Millipore Sigma, Burlington, MA, USA) was used as the dispersant. Styrene-butadiene rubber (SBR, ~2.1 wt.%, MTI Corporation, Richmond, CA, USA) was the binder used to hold the GNP particles together once the solvent is sublimed. Two slurry compositions were prepared to regulate the pore size by controlling the length of the solidifying waterfront. Two mold shapes, circular and rectangular, were used to regulate the pore morphology by controlling the heat transfer direction and rate. Both the molds were made from Aluminum (Al). Initially, CMC and SBR were vortex mixed in DI water for 30 min to obtain the solvent. GNP slurry batches of 8 wt.% and 9.5 wt.% was vortex mixed in the solvent (~89 wt.% and 87 wt.% respectively de-ionized water) for 120 min. The slurry was frozen in the Al molds for 4 h at −56 °C and dried (sublimation of water) for 24 h in 1 Pa ambient pressure in a freeze drier (Pro-Freeze Dryer PLT300, Pro-Lab Inc, Fort Worth, TX, USA).

2.2. Microstructure Analysis

A field emission scanning electron microscope (FE-SEM JEOL JSM-6330F, JEOL Ltd. Tokyo, Japan) was used to analyze the FD GNP foam’s pore size and morphology. An acceleration voltage of 20 kV was used throughout the secondary electron imaging. The pore size was evaluated from the SEM images using ImageJ 1.52a an open source image processing tool.

2.3. Compression Test

Displacement controlled compression tests were performed on FD GNP foam and commercially available CVD GrF (Graphene Supermarket, Calverton, NY, USA). All samples were strained to 55% of their original length at a displacement rate of 0.016 mm/second using Electroforce 3200 mechanical tester (Bose corporation, Eden Prairie, MN, USA) equipped with a 25 N load cell. Both FD (8 wt.% GNP, cast in Al mold) and commercial CVD GrF had dimension of 10 mm in diameter and a thickness of ~1.2–1.5 mm. They were sliced into required dimensions using a razor blade.

2.4. Thermal Analysis

Thermal diffusivity (α) of the FD GNP foam and commercial CVD GrF (Graphene Supermarket, Calverton, NY, USA, USA) across the thickness (in-plane) was measured by laser flash technique (NETZSH LFA 467 HT HyperFlash, Selb, Germany) at temperatures between 25 and 75 °C. Both FD (8 wt.% GNP, cast in Al mold) and commercial CVD GrF had dimension of 10 mm in diameter and a thickness of ~1.2–1.5 mm.

2.5. Raman Analysis

Micro-Raman spectroscopy analysis was performed using Spectra-Physics 3900S (Newport Corporation, Irvine, CA, USA) equipped with Ti-sapphire crystal (514 nm) as the target. The Spectra-Physics also consists of a detector with 4 cm−1 spectral resolution from Kaiser Optical Systems, Inc. (Ann Arbor, MI, USA), a laser power (18 mW), with a spot size of 2 µm.

3. Results and Discussion

3.1. Controlling the Foam’s Pore Size and Morphology

3.1.1. Solid Loading

To comprehend the FD technique’s theoretical understanding, one should consider the system’s thermodynamic free energy as a critical parameter [38]. Since the process involves solidifying the colloidal slurry (GNP mixed in DI water, CMC, and SBR solvent), the freezing solvent front must segregate the GNP particles as they propagate. For the GNP particles to be pushed away from the propagating ice crystal, the free energy of the system ∆σ must be positive. If ∆σ is negative, the GNP particles will be trapped inside the freezing ice crystal, inducing the formed foamy structures’ low structural integrity. Hence, a slurry for FD should satisfy the following thermodynamic condition, which is given by Equation (1) [36]:
Δ σ = σ s p − ( σ l p + σ s l ) > 0
where σsp, σlp, and σsl are the interfacial free energies related with the solid(ice-crystal)–particle (GNP), liquid (DI water) –particle (GNP), and solid(ice-crystal)–liquid (DI-water) interfaces, respectively. If Equation (1) is satisfied, then the solidification front rejects the GNP particles in the slurry system. If the system’s overall free energy decreases (become negative), then, based on Equation (1), the particles in the slurry are engulfed by the propagating ice crystal, which is not ideal for developing porous 3D foam with a hierarchical structure. The solid loading is one such variable that drastically affects the ∆σ of the system. ∆σ is directly influenced by the thermal gradient generated in the slurry system, which dictates the freezing ice crystals’ growth rate or kinetics. The temperature gradient induced in a slurry system can be represented as Equation (2) [36]:
( T 1 − T 2 ) L
where T1 and T2 are the time-dependent temperatures of the atmosphere outside the mold and inside the mold (GNP slurry), respectively, as shown in the schematic in Figure 6a,c. L is the length of the mold along which the majority of the heat transfer takes place. The practical demonstration of Equations (1) and (2) is shown in Figure 6. Figure 6b is the 3D GNP foam produced with 8 wt.% GNP and Figure 6d is the 3D GNP foam made with 9.5 wt.% GNP. Solid loading was selected based on a general guideline that slurry with solid loading <5 vol.% and >45 vol.% often results in 3D porous structures with random pore size and morphologies [38]. As mentioned above, the solid loading affects the ∆σ of the slurry system, which alters the propagating ice crystals’ kinetics, influencing the size of the pores formed.
3D GNP foam with 8 wt.% GNP (Figure 6b) has an average pore size of ~110 µm with oblong pore morphology, while 3D GNP foam with 9.5 wt.% GNP (Figure 6d) has an average pore size of ~50 µm with near-circular pore morphology. In the slurry with 8 wt.% GNP, the magnitude of ∆σ is higher than slurry with 9.5 wt.% GNP. With the increased GNP content in the slurry, the thermal conductivity of the slurry increases. As a result, T1 equilibrates with T2 at a higher rate (low ∆σ), forming ice crystals with a smaller radius and length (refer to Figure 6c) that manifests into pores of smaller size and thick wall, as shown in Figure 6d and its respective inset. Increased pore walls or strut thickness is the result of mass conservation [46,47]. In the slurry with 8 wt.% GNP, its thermal conductivity is comparatively lower than 9.5 wt.% slurries. In this case, the slurry system experiences a thermal gradient for a prolonged period (high ∆σ), resulting in a longer ice crystal with a higher radius (refer to Figure 6a). After sublimation, the resultant foam will have increased pore size with a thin wall, as shown in Figure 6b and its respective inset. Researchers have also looked into other complex methods to better control the pore size and morphology by regulating the viscous stresses and pressure drops experienced by the particles as it interacts with the propagating solidification front [48].

3.1.2. Effect of Mold Geometry on Pores

Another easily adaptable technique in tailoring desired pore size and morphology is by using molds of different geometry. The mold into which the slurry is poured is the critical source to impact the heat transfer between the freezing atmosphere and the slurry during solidification. This heat transfer is mainly convective type and takes the form, Equation (3):
Q = h c A ( T 2 − T 1 )
where Q is the heat transferred per unit time, A is the mold surface area for the heat transfer, hc is the heat transfer thin film co-efficient of the heat transfer medium, T2 is the initial mold surface temperature, and T1 is the temperature of the freezing environment.
The solidification conditions control pore morphology and size as described in the previous section. The ice crystals grow parallel but in the opposite direction to the thermal gradient (red arrows in Figure 7 represent the heat transfer direction). In the cylindrical mold (Figure 7a), most of the heat transfers along the radial surface (radial surface area ~640 mm2).
Hence, the ice crystals grow from the circumference of the mold inwards. Therefore, in the SEM in Figure 7a, pores can be observed aligned along the rim (represented by the red dashed curve in Figure 7a and radially grown inwards. In addition, there is heat transfer from the bottom of the mold (surface area ~80 mm2). However, heat transfer along the length of the mold is not as dominant as across the circumference. Here, the thermal gradient can be considered as homogenous, resulting in oriented hierarchical porous architecture. Due to higher heat transfer along the mold’s circumference, pore channels run through the entire sample forming a reticulated porous structure. In the rectangular mold (drawing in Figure 7b), there is a heterogeneous thermal gradient as the entire mold surface area is of significant dimension. Along the length, the surface area is ~500 mm2, the surface area is ~200 mm2 along the width, and the bottom mold surface area is ~1020 mm2. As seen in the SEM in Figure 7b, the difference in a thermal gradient across the rectangular mold surface resulted in ice crystals’ growth in a random direction. Although long-range lamellar porous structures can be observed, they do not form a reticulated system with a continuous pore channel. Localized small range lamellar structures can also be observed in the SEM in Figure 7b. Unlike in cylindrical mold, the pores are not aligned along the edges of the rectangular mold (red dashed line in Figure 7b), but, instead, they are at an angle due to the difference in a thermal gradient along the length, width, and height of the mold. Therefore, the temperature field must be carefully controlled to obtain homogenous pore size and morphology in FD samples. In CVD, the pore size and morphology depend on the metal template used to deposit the graphene film. Hence, using CVD processing, tailoring the pore morphology and size in the 3D reticulated GrF is not straightforward as in the FD technique. The new porous metal substrate with the desired pore shape and morphology needs to be fabricated every time, whenever customization is required, which is later etched away, leading to wastage of metal and generation of toxic acid byproducts.
Another advantage of producing 3D GNP foam by FD is obtaining a foam with aligned GNP along its basal plane. Since the direction of the solidification crystal front’s growth can be controlled, the GNP’s basal plane aligns along the solidification direction, as shown in Figure 8a. With such a reinforcement nanofoam, it is possible to produce advanced engineering materials with directional properties, such as a material system that conducts electrons along the thickness, preferably across the material’s bulk. In the case of CVD produced GrF, although one can obtain high purity GrF with a preferred number of graphene films, as the basal plane rests on the deposited surface, they are at different angles concerning each other. Hence, producing 3D GrF with the aligned basal plane is challenging via the CVD technique. Figure 8c shows a pristine 3D GrF produced by the CVD technique, and the respective inset shows the wall or the strut of CVD GrF. It can be seen that the CVD technique results in foam with a hollow strut due to the removal of the metal substrate by acid etching, whereas 3D GNP foam produced via FD technique has foam with solid strut with highly oriented GNP, as shown in Figure 8b, of similar strut architecture to CVD GrF. FD GNP foam packs higher GNP content in a given strut as compared to the CVD produced GrF.

3.2. Mechanical and Thermal Properties of Graphene Foams

Since GrF is an emerging nanofiller used to augment the mechanical and phonon conduction properties of the matrix material, the FD GNP foam’s mechanical and thermal properties are analyzed and compared with CVD GrF. The FD GNP used in the property characterization has 8 wt.% of GNP and was cast in a cylindrical mold, as shown in Figure 7a.
Figure 9a is the compressive stress vs. strain plot of FD GNP foam and CVD GrF. The inset in Figure 9a shows the test sample (within the red oval) placed between the flat plates of the mechanical testing machine. The sample merely has a large surface area compared to the surface area of the flat plate used in the compression test. Approximately 0.15 mm2 of the sample area is not in contact with the flat plate. FD GNP foam records a compressive strength of ~120 kPa with an elastic modulus of 283 ± 33 kPa, whereas CVD exhibited a compressive strength of ~10 kPa and an elastic modulus of 0.22 ± 0.024. FD GNP foam exhibits ~12 times the compressive strength and ~1280 times more elastic modulus than CVD GrF. It should be noted that the GNP foams in this study were not subjected to any type of heat treatment process after FD to increase its rigidity. The increase in the mechanical properties can be attributed to the presence of a solid wall architecture (refer Figure 8b) exhibited by FD GNP foam unlike CVD GrF, where the wall is hollow (refer Figure 8c). Trace amount of polymeric rubber in FD GNP foam also contributes to the increase in mechanical properties. The substantial increase in the slope of the FD GNP foam with increased stress can be due to the efficient stress transfer along the node–branch architecture of the foam and bending of the solid struts of the foam [25,28]. Other reasons for strengthening in graphene-based foams can be the flattening of intrinsic corrugations, inter layer van der Waals spring-like actions, kink band formation, and membrane vibrations [25,28]. It is well known that pristine CVD GrF is challenging to handle, as even a minuscule force, as small as the soft touch of a human, disintegrates the foam. As a result, integrating CVD GrF into a highly viscous material matrix without compromising its hierarchical structure is challenging. The use of large CVD GrF (the foam buckles due to its weight and disintegrates) is still a challenge. Thus, commercially available pristine CVD GrF are small with a typical size of 2 inches × 2 inches and a thickness of 1.2 mm. The presence of a small amount of elastomeric binder in FD GNP foam produces a rigid foam that can be easily cut into any shape and thickness. The production of FD GNP foam of various shapes, sizes, and thicknesses with controlled pore morphology and directionality, as shown in Figure 10, is possible with ease, and no special tooling or equipment is needed.
Figure 9b is the thermal diffusivity exhibited by FD GNP foam and CVD GrF in the temperature range 25–75 °C. CVD GrF has 2.25 times higher thermal diffusivity than FD GNP foam at 25 °C, and it increases to three times at 75 °C. The elastomeric binder that increased FD GNP foam’s mechanical property is the one of the reasons for the decreased thermal diffusivity. SBR has a negative glass transition temperature (−60 °C), and, above this temperature, the thermal conductivity of SBR decreases asymptotically [49]. Hence, the thermal diffusivity of the FD GNP foam decreases with increased temperature. Another reason for decreased thermal diffusivity is due to the lack of solid-phase continuity across the strut wall of the FD sample (refer to Figure 8a,b). Raman spectroscopy was employed to evaluate the graphene defects in both foams, as defects can influence the intrinsic phonon conduction property [24]. Table 1 tabulates the intensity of D, G, and 2D peaks deduced from the Raman spectra shown in Figure 9c.
The ratio of ID/IG shines a light on the extent of defect in graphene. Table 1 shows that GNP used in FD GNP foam has three times more defects than graphene forming the CVD GrF. This resulted in the distortion of the phonon conduction band in FD GNP foam, as confirmed by the I2D/IG ratio, which is 1.1 times higher than that of CVD GrF. The higher number of defects and distortion of the phonon conduction band in the GNP used in FD GNP foam can also explain the reduced thermal diffusivity in FD GNP foam.
Due to its chemical purity, CVD GrF is advantageous for high thermal conductivity and small-scale non-structural applications. However, when the nano-filler 3D graphene foam’s structural strength is essential, FD GNP foam is a great alternative. Although using FD GNP foam, the thermal conductivity may be compromised to a certain extent due to the polymeric binder, the benefits, such as particle alignment and tailorable pore size, morphology, and scalability, make FD GNP foam versatile in this class of nano-fillers. This work provides an outlook on the FD process in producing novel reinforcement material with a 3D hierarchical structure. One of FD’s primary advantages over the CVD technique in producing 3D foam structures is that the FD process is independent of material chemistry, unlike CVD. Hence, any material with proper slurry formulation can be artificially stacked into hierarchical reticulated structures (Figure 10).
Figure 11 demonstrates that any material system, whether 1D, 2D, or 3D, can be cast into highly reticulated structures easily without the need for any expensive equipment. Figure 11a shows GNP/CNT foam’s SEM image with the GNP decorated by CNT (shown in the inset). Figure 11b is the SEM image of 2D hexagonal boron nitride platelet (BNNP) foam with a thin leaf-like wall made up of highly aligned BNNP particles (refer to inset in Figure 11b). Finally, Figure 11c shows 1D hexagonal boron nitride nanotubes (BNNT) formed into porous structures of untangled BNNT (refer to inset in Figure 11c). All the 3D structures shown in Figure 11 were produced using similar slurry formation chemicals described in Section 2.1 (more details can be obtained upon request). Using the energy-efficient FD technique, the material science community can produce various novel engineering materials that can open multiple engineering applications such as 3D scaffolds for tissue engineering, high-density energy storage devices, high strength lightweight composite structural materials, and many more.

4. Conclusion

FD’s versatility in producing mechanically rigid yet highly reticulated GNP foam is presented. The FD process is advantageous when controlling the pore size and morphology of GrF. Using FD, without any need for specialized equipment, the intrinsic thermodynamics in the freezing solvent’s solidification can be engineered to achieve control over the reticulated architectural GNP foam design. In this work, the critical thermodynamic parameter that dictates the freezing ice crystal’s length and radius, free energy, was regulated by formulating slurry with different GNP solid loading and varying geometry molds. It was demonstrated that both these approaches influence the extent of thermal gradient and heat transfer direction to produce GNP foam with controlled pore size and morphology. A unique feature in the GNP foam made by FD is the alignment of GNP’s basal plane along the direction of the propagating ice crystal. This can induce intrinsic directional properties of GNP in the reinforced material matrix. In addition, unlike CVD produced GrF, the FD produced GNP foam walls are solid and not hollow, providing high graphene density for similar wall thickness. Compression test of FD GNP foam exhibited ~1280 times higher elastic modulus than CVD GrF, reassuring that FD GNP foam is a rigid and robust reinforcing nano-filler that would not disintegrate or collapse when infiltrated with high viscous slurry for new material development. Although FD GNP foam’s thermal property is 0.5 times the CVD GrF, the ease of process scalability and ability to artificially stack any material into a highly reticulated 3D structure makes FD a unique foam producing technology.

Author Contributions

Conceptualization, A.A. and T.T.; methodology, A.A. and T.T.; software, T.T.; validation, A.A. and T.T.; formal analysis, T.T.; investigation, A.A. and T.T.; resources A.A. and T.T.; data curation, T.T.; writing—original draft preparation, T.T.; writing—review and editing, A.A. and T.T.; visualization, T.T.; supervision, A.A.; project administration, A.A.; funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded in part by NASA under Award # 80NSSC19M0201.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data reported in this study can be obtained by requesting the authors.

Acknowledgments

The authors acknowledge the Advanced Materials Engineering Research Institute (AMERI) support at Florida International University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Voigt, W. Ueber die Beziehung zwischen den beiden Elasticitätsconstanten isotroper Körper. Ann. der Phys. 1889, 274, 573–587. [Google Scholar] [CrossRef] [Green Version]
  2. Hill, R. Elastic properties of reinforced solids: Some theoretical principles. J. Mech. Phys. Solids 1963, 11, 357–372. [Google Scholar] [CrossRef]
  3. Qiu, Y.; Weng, G. On the application of Mori-Tanaka’s theory involving transversely isotropic spheroidal inclusions. Int. J. Eng. Sci. 1990, 28, 1121–1137. [Google Scholar] [CrossRef]
  4. Sheng, Y.Q.; Xu, T.; Hui, Y. A stepping scheme for predicting effective properties of the multi-inclusion composites. Int. J. Eng. Sci. 2007, 45, 997. [Google Scholar]
  5. Qin, Q. Introduction to the composite and its toughening mechanisms. In Toughening Mechanisms in Composite Materials; Woodhead Publishing: Cambridge, UK, 2015; pp. 1–32. [Google Scholar]
  6. Geim, A.K.; Novoselov, K.S. The rise of graphene. Nat. Mater. 2007, 6, 183–191. [Google Scholar] [CrossRef] [PubMed]
  7. Lee, C.; Wei, X.; Kysar, J.W.; Hone, J. Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 2008, 321, 385–388. [Google Scholar] [CrossRef]
  8. Balandin, A.A.; Ghosh, S.; Bao, W.; Calizo, I.; Teweldebrhan, D.; Miao, F.; Lau, C.N. Superior Thermal Conductivity of Single-Layer Graphene. Nano Lett. 2008, 8, 902–907. [Google Scholar] [CrossRef] [PubMed]
  9. Kinloch, I.A.; Suhr, J.; Lou, J.; Young, R.J.; Ajayan, P.M. Composites with carbon nanotubes and graphene: An outlook. Science 2018, 362, 547–553. [Google Scholar] [CrossRef] [Green Version]
  10. Song, S.H.; Park, K.H.; Kim, B.H.; Choi, Y.W.; Jun, G.H.; Lee, D.J.; Kong, B.-S.; Paik, K.-W.; Jeon, S. Enhanced Thermal Conductivity of Epoxy-Graphene Composites by Using Non-Oxidized Graphene Flakes with Non-Covalent Functionalization. Adv. Mater. 2013, 25, 732–737. [Google Scholar] [CrossRef] [PubMed]
  11. Li, D.; Müller, M.B.; Gilje, S.; Kaner, R.B.; Wallace, G.G. Processable aqueous dispersions of graphene nanosheets. Nat. Nanotechnol. 2008, 3, 101–105. [Google Scholar] [CrossRef] [PubMed]
  12. Nieto, A.; Boesl, B.; Agarwal, A. Multi-scale intrinsic deformation mechanisms of 3D graphene foam. Carbon 2015, 85, 299–308. [Google Scholar] [CrossRef]
  13. Nautiyal, P.; Mujawar, M.; Boesl, B.; Agarwal, A. In-situ Mechanics of 3D Graphene Foam Based Ultra-stiff and Flexible Metallic Metamaterial. Carbon 2018, 137, 502. [Google Scholar] [CrossRef]
  14. Yao, X.; Zhao, Y. Three-Dimensional Porous Graphene Networks and Hybrids for Lithium-Ion Batteries and Supercapacitors. Chem 2017, 2, 171–200. [Google Scholar] [CrossRef]
  15. Li, T.; Li, N.; Liu, J.; Cai, K.; Foda, M.F.; Lei, X.; Han, H. Synthesis of functionalized 3D porous graphene using both ionic liquid and SiO2 spheres as “spacers” for high-performance application in supercapacitors. Nanoscale 2014, 7, 659–669. [Google Scholar] [CrossRef] [PubMed]
  16. Idowu, A.; Boesl, B.; Agarwal, A. 3D graphene foam-reinforced polymer composites–A review. Carbon 2018, 135, 52–71. [Google Scholar] [CrossRef]
  17. Zhu, C.; Han, T.Y.-J.; Duoss, E.B.; Golobic, A.M.; Kuntz, J.D.; Spadaccini, C.M.; Worsley, M.A. Highly compressible 3D periodic graphene aerogel microlattices. Nat. Commun. 2015, 6, 6962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Meng, Y.; Wang, K.; Zhang, Y.; Wei, Z. Hierarchical Porous Graphene/Polyaniline Composite Film with Superior Rate Performance for Flexible Supercapacitors. Adv. Mater. 2013, 25, 6985–6990. [Google Scholar] [CrossRef] [PubMed]
  19. Qin, Z.; Jung, G.S.; Kang, M.J.; Buehler, M.J. The mechanics and design of a lightweight three-dimensional graphene assembly. Sci. Adv. 2017, 3, e1601536. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, M.; Duan, X.; Xu, Y.; Duan, X. Functional Three-Dimensional Graphene/Polymer Composites. ACS Nano. 2016, 10, 7231–7247. [Google Scholar] [CrossRef]
  21. Bustillos, J.; Zhang, C.; Boesl, B.; Agarwal, A. Three-Dimensional Graphene Foam–Polymer Composite with Superior Deicing Efficiency and Strength. ACS Appl. Mater. Interfaces 2018, 10, 5022–5029. [Google Scholar] [CrossRef] [PubMed]
  22. Nieto, A.; Dua, R.; Zhang, C.; Boesl, B.; Ramaswamy, S.; Agarwal, A. Three Dimensional Graphene Foam/Polymer Hybrid as a High Strength Biocompatible Scaffold. Adv. Funct. Mater. 2015, 25, 3916–3924. [Google Scholar] [CrossRef]
  23. Wang, C.; Zhang, C.; Chen, S. The microscopic deformation mechanism of 3D graphene foam materials under uniaxial compression. Carbon 2016, 109, 666–672. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, Z.; Ren, W.; Gao, L.; Liu, B.; Pei, S.; Cheng, H.-M. Three-dimensional flexible and conductive interconnected graphene networks grown by chemical vapour deposition. Nat. Mater. 2011, 10, 424–428. [Google Scholar] [CrossRef]
  25. Nautiyal, P.; Boesl, B.; Agarwal, A. The mechanics of energy dissipation in a three-dimensional graphene foam with macroporous architecture. Carbon 2018, 132, 59–64. [Google Scholar] [CrossRef]
  26. Jia, J.; Sun, X.; Lin, X.; Shen, X.; Mai, Y.; Kim, J. Exceptional Electrical Conductivity and Fracture Resistance of 3D Interconnected Graphene Foam/Epoxy Composites. ACS Nano. 2014, 8, 5774–5783. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, Z.; Xu, C.; Ma, C.; Ren, W.; Cheng, H.-M. Lightweight and Flexible Graphene Foam Composites for High-Performance Electromagnetic Interference Shielding. Adv. Mater. 2013, 25, 1296–1300. [Google Scholar] [CrossRef]
  28. Embrey, L.; Nautiyal, P.; Loganathan, A.; Idowu, A.; Boesl, B.; Agarwal, A. Three-Dimensional Graphene Foam Induces Multifunctionality in Epoxy Nanocomposites by Simultaneous Improvement in Mechanical, Thermal, and Electrical Properties. ACS Appl. Mater. Interfaces 2017, 9, 39717. [Google Scholar] [CrossRef]
  29. Jiang, W.; Xin, H.; Li, W. Microcellular 3D graphene foam via chemical vapor deposition of electroless plated nickel foam templates. Mater. Lett. 2016, 162, 105–109. [Google Scholar] [CrossRef]
  30. Huang, M.; Wang, C.; Quan, L.; Nguyen, T.H.-Y.; Zhang, H.; Jiang, Y.; Byun, G.; Ruoff, R.S. CVD Growth of Porous Graphene Foam in Film Form. Matter 2020, 3, 487–497. [Google Scholar] [CrossRef]
  31. Lee, J.-S.; Ahn, H.-J.; Yoon, J.-C.; Jang, J.-H. Three-dimensional nano-foam of few-layer graphene grown by CVD for DSSC. Phys. Chem. Chem. Phys. 2012, 14, 7938–7943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Manjakkal, L.; Garcia, C.G.; Dang, W.; Dahiya, R. Flexible self-charging supercapacitor based on gra-phene-Ag-3D graphene foam electrodes. Nano Energ. 2018, 51, 604. [Google Scholar] [CrossRef]
  33. Krueger, E.; Chang, A.N.; Brown, D.; Eixenberger, J.; Brown, R.; Rastegar, S.; Yocham, K.M.; Cantley, K.D.; Estrada, D. Graphene Foam as a Three-Dimensional Platform for Myotube Growth. ACS Biomater. Sci. Eng. 2016, 2, 1234–1241. [Google Scholar] [CrossRef] [PubMed]
  34. Deville, S. Freeze-Casting of Porous Biomaterials: Structure, Properties and Opportunities. Materials 2010, 3, 1913–1927. [Google Scholar] [CrossRef]
  35. Gaudillere, C.; Serra, J.M. Freeze-casting: Fabrication of highly porous and hierarchical ceramic supports for energy applications. Boletín Soc. Española Cerámica Vidr. 2016, 55, 45–54. [Google Scholar] [CrossRef] [Green Version]
  36. Deville, S.; Saiz, E.; Tomsia, A.P. Ice-templated porous alumina structures. Acta Mater. 2007, 55, 1965–1974. [Google Scholar] [CrossRef] [Green Version]
  37. Deville, S. Freeze-Casting of Porous Ceramics: A Review of Current Achievements and Issues. Adv. Eng. Mater. 2008, 10, 155–169. [Google Scholar] [CrossRef]
  38. Nelson, I.; Naleway, S.E. Intrinsic and extrinsic control of freeze casting. J. Mater. Res. Technol. 2019, 8, 2372–2385. [Google Scholar] [CrossRef]
  39. Sofie, S.W.; Dogan, F. Freeze Casting of Aqueous Alumina Slurries with Glycerol. J. Am. Ceram. Soc. 2004, 84, 1459–1464. [Google Scholar] [CrossRef]
  40. Fukasawa, T.; Ando, M.; Ohji, T.; Kanzaki, S. Synthesis of Porous Ceramics with Complex Pore Structure by Freeze-Dry Processing. J. Am. Ceram. Soc. 2001, 84, 230–232. [Google Scholar] [CrossRef]
  41. Li, W.; Li, Y.; Meng, Q.; Xu, Y.; Yang, B.; Li, T. Facile Preparation of a Ni/Mn-LDH Nanosheet @Graphene Foam Hierarchical Structure for High-performance Flexible Supercapacitors. Mater. Res. Bull. 2020, 129, 110889. [Google Scholar] [CrossRef]
  42. Zhu, S.; Cheng, Q.; Yu, C.; Pan, X.; Zuo, X.; Liu, J.; Chen, M.; Li, W.; Li, Q.; Liu, L. Flexible Fe3O4/graphene foam/poly dimethylsiloxane composite for high-performance electromagnetic interference shielding. Compos. Sci. Technol. 2020, 189, 108012. [Google Scholar] [CrossRef]
  43. Geng, Y.; Li, Z.; Chen, M.; Zhao, H.; Zhang, L. The preparation of graphene foam by one-step reduction and air-drying for oil–water separation. J. Sol. Gel Sci. Technol. 2020, 94, 375–383. [Google Scholar] [CrossRef]
  44. Zhao, Z.; Shen, T.; Liu, Z.; Zhong, Q.; Qin, Y. Facile fabrication of binder-free reduced graphene oxide/MnO2/Ni foam hybrid electrode for high-performance supercapacitors. J. Alloys Compd. 2020, 815, 152124. [Google Scholar] [CrossRef]
  45. Wang, Z.; Zhao, Z.; Baucom, J.; Wang, D.; Dai, L.; Chen, J.-F. Nitrogen-Doped Graphene Foam as a Metal-Free Catalyst for Reduction Reactions under a High Gravity Field. Engineering 2020, 6, 680–687. [Google Scholar] [CrossRef]
  46. Ghazanfari, S.M.H.; Zamanian, A. Phase transformation, microstructural and mechanical properties of hydroxyapatite/alumina nanocomposite scaffolds produced by freeze casting. Ceram. Int. 2013, 39, 9835–9844. [Google Scholar] [CrossRef]
  47. Munch, E.; Saiz, E.; Tomsia, A.P.; Deville, S. Architectural Control of Freeze-Cast Ceramics Through Additives and Templating. J. Am. Ceram. Soc. 2009, 92, 1534–1539. [Google Scholar] [CrossRef] [Green Version]
  48. Saint-Michel, B.; Georgelin, M.; Deville, S.; Pocheau, A. Interaction of Multiple Particles with a Solidification Front: From Compacted Particle Layer to Particle Trapping. Langmuir 2017, 33, 5617–5627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Saxena, N.; Pradeep, P.; Mathew, G.; Thomas, S.; Gustafsson, M.; Gustafsson, S. Thermal conductivity of styrene butadiene rubber compounds with natural rubber prophylactics waste as filler. Eur. Polym. J. 1999, 35, 1687–1693. [Google Scholar] [CrossRef]
Materials 14 00864 g001 550
Figure 1. Schematic detailing the properties of pristine graphene.
Figure 1. Schematic detailing the properties of pristine graphene.
Materials 14 00864 g001
Materials 14 00864 g002 550
Figure 2. Schematic summarizing applications of graphene (Gr)-based materials.
Figure 2. Schematic summarizing applications of graphene (Gr)-based materials.
Materials 14 00864 g002
Materials 14 00864 g003 550
Figure 3. (a) SEM image showing clustered and poorly dispersed 2D Gr flakes (highlighted in red dashed oval) in a ceramic matrix; (b) schematic representing the discontinuous contact between poorly dispersed 2D Gr in the material matrix, impeding phonon and electron conduction; (c) SEM image of pristine GrF; and (d) schematic representing how the architecture of GrF offers a seamless pathway for phonon and electron conduction.
Figure 3. (a) SEM image showing clustered and poorly dispersed 2D Gr flakes (highlighted in red dashed oval) in a ceramic matrix; (b) schematic representing the discontinuous contact between poorly dispersed 2D Gr in the material matrix, impeding phonon and electron conduction; (c) SEM image of pristine GrF; and (d) schematic representing how the architecture of GrF offers a seamless pathway for phonon and electron conduction.
Materials 14 00864 g003
Materials 14 00864 g004 550
Figure 4. Data showing the number of scientific publications on GrF produced by chemical vapor disposition (CVD) since its advent in 2011.
Figure 4. Data showing the number of scientific publications on GrF produced by chemical vapor disposition (CVD) since its advent in 2011.
Materials 14 00864 g004
Materials 14 00864 g005 550
Figure 5. Schematic showing the production of 3D GrF by CVD technique [29,30,31].
Figure 5. Schematic showing the production of 3D GrF by CVD technique [29,30,31].
Materials 14 00864 g005
Materials 14 00864 g006a 550Materials 14 00864 g006b 550
Figure 6. (a,c) Schematic representing the influence of free energy of the slurry system on the length and diameter of the freezing ice crystal; (b) SEM image of 3D GNP foam with 8 wt.% GNP; and (d) SEM image of 3D GNP foam with 9.5 wt.% GNP (3D foam made in a circular Al mold of 20 mm length, 10 mm diameter, and 0.16 mm thick).
Figure 6. (a,c) Schematic representing the influence of free energy of the slurry system on the length and diameter of the freezing ice crystal; (b) SEM image of 3D GNP foam with 8 wt.% GNP; and (d) SEM image of 3D GNP foam with 9.5 wt.% GNP (3D foam made in a circular Al mold of 20 mm length, 10 mm diameter, and 0.16 mm thick).
Materials 14 00864 g006aMaterials 14 00864 g006b
Materials 14 00864 g007 550
Figure 7. (a) 3D GNP foam produced in a cylindrical mold (20 mm in length, 10 mm in diameter and mold thickness of 0.16 mm); and (b) 3D GNP foam made in a rectangular mold (50 mm in size, 20 mm in width, 10 mm in height diameter, and mold thickness of 0.16 mm). Both insets show the wall thickness (8 wt.% GNP slurry used in both the molds).
Figure 7. (a) 3D GNP foam produced in a cylindrical mold (20 mm in length, 10 mm in diameter and mold thickness of 0.16 mm); and (b) 3D GNP foam made in a rectangular mold (50 mm in size, 20 mm in width, 10 mm in height diameter, and mold thickness of 0.16 mm). Both insets show the wall thickness (8 wt.% GNP slurry used in both the molds).
Materials 14 00864 g007
Materials 14 00864 g008 550
Figure 8. SEM images showing: (a) FD 3D GNP foam with basal plane aligned; (b) reticulated 3D architecture of FD GNP foam with solid wall or strut; and (c) reticulated 3D architecture of CVD GrF with hollow wall or strut (inset showing hollow wall or strut).
Figure 8. SEM images showing: (a) FD 3D GNP foam with basal plane aligned; (b) reticulated 3D architecture of FD GNP foam with solid wall or strut; and (c) reticulated 3D architecture of CVD GrF with hollow wall or strut (inset showing hollow wall or strut).
Materials 14 00864 g008
Materials 14 00864 g009 550
Figure 9. (a) Compressive stress vs. strain plot comparison of FD GNP foam and CVD GrF; (b) thermal diffusivity of FD GNP foam vs. CVD GrF; and (c) Raman spectra of FD GNP foam and CVD GrF.
Figure 9. (a) Compressive stress vs. strain plot comparison of FD GNP foam and CVD GrF; (b) thermal diffusivity of FD GNP foam vs. CVD GrF; and (c) Raman spectra of FD GNP foam and CVD GrF.
Materials 14 00864 g009
Materials 14 00864 g010 550
Figure 10. 3D GNP foam produced in various sizes, shapes, and thicknesses by FD technique to highlight the scalability of the process with ease.
Figure 10. 3D GNP foam produced in various sizes, shapes, and thicknesses by FD technique to highlight the scalability of the process with ease.
Materials 14 00864 g010
Materials 14 00864 g011 550
Figure 11. Porous 3D Foam produced by FD: (a) GNP/CNT composite foam; (b) BNNP foam; and (c) BNNT foam. Insets show the wall features of the foams.
Figure 11. Porous 3D Foam produced by FD: (a) GNP/CNT composite foam; (b) BNNP foam; and (c) BNNT foam. Insets show the wall features of the foams.
Materials 14 00864 g011
Table
Table 1. The intensity of D, G, and 2D peak for FD GNP and CVD GrF foams.
Table 1. The intensity of D, G, and 2D peak for FD GNP and CVD GrF foams.
SampleIDIGI2DID/IGI2D/IG
CVD GrF873.55340853780.0811.453
FD GNP364.37451665650.261.57
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Thomas, T.; Agarwal, A. A Facile and Scalable Approach in the Fabrication of Tailored 3D Graphene Foam via Freeze Drying. Materials 2021, 14, 864. https://doi.org/10.3390/ma14040864

AMA Style

Thomas T, Agarwal A. A Facile and Scalable Approach in the Fabrication of Tailored 3D Graphene Foam via Freeze Drying. Materials. 2021; 14(4):864. https://doi.org/10.3390/ma14040864

Chicago/Turabian Style

Thomas, Tony, and Arvind Agarwal. 2021. "A Facile and Scalable Approach in the Fabrication of Tailored 3D Graphene Foam via Freeze Drying" Materials 14, no. 4: 864. https://doi.org/10.3390/ma14040864

APA Style

Thomas, T., & Agarwal, A. (2021). A Facile and Scalable Approach in the Fabrication of Tailored 3D Graphene Foam via Freeze Drying. Materials, 14(4), 864. https://doi.org/10.3390/ma14040864

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